Parallel evolution of opsin visual pigments in hawkmoths by tuning of spectral sensitivities during transition from a nocturnal to a diurnal ecology

Light environments differ dramatically between day and night. The transition between diurnal and nocturnal visual ecology has happened repeatedly throughout evolution in many species. However, the molecular mechanism underlying the evolution of vision in recent diurnal-nocturnal transition is poorly understood. Here, we focus on hawkmoths (Lepidoptera: Sphingidae) to address this question by investigating five nocturnal and five diurnal species. We performed RNA-sequencing analysis and identified opsin genes corresponding to the ultraviolet (UV), short-wavelength (SW) and long-wavelength (LW)-absorbing visual pigments. We found no significant differences in the expression patterns of opsin genes between the nocturnal and diurnal species. We then constructed the phylogenetic trees of hawkmoth species and opsins. The diurnal lineages had emerged at least three times from the nocturnal ancestors. The evolutionary rates of amino acid substitutions in the three opsins differed between the nocturnal and diurnal species. We found an excess number of parallel amino acid substitutions in the opsins in three independent diurnal lineages. The numbers were significantly more than those inferred from neutral evolution, suggesting that positive selection acted on these parallel substitutions. Moreover, we predicted the visual pigment absorption spectra based on electrophysiologically determined spectral sensitivity in two nocturnal and two diurnal species belonging to different clades. In the diurnal species, the LW pigments shift 10 nm towards shorter wavelengths, and the SW pigments shift 10 nm in the opposite direction. Taken together, our results suggest that parallel evolution of opsins may have enhanced the colour discrimination properties of diurnal hawkmoths in ambient light.

Molecular Phylogenetic Analysis of the AIG Family in Vertebrates

Androgen-inducible genes (AIGs), which can be regulated by androgen level, constitute a group of genes characterized by the presence of the AIG/FAR-17a domain in its protein sequence. Previous studies on AIGs demonstrated that one member of the gene family, AIG1, is involved in many biological processes in cancer cell lines and that ADTRP is associated with cardiovascular diseases. It has been shown that the numbers of AIG paralogs in humans, mice, and zebrafish are 2, 2, and 3, respectively, indicating possible gene duplication events during vertebrate evolution. Therefore, classifying subgroups of AIGs and identifying the homologs of each AIG member are important to characterize this novel gene family further. In this study, vertebrate AIGs were phylogenetically grouped into three major clades, ADTRP, AIG1, and AIG-L, with AIG-L also evident in an outgroup consisting of invertebrsate species. In this case, AIG-L, as the ancestral AIG, gave rise to ADTRP and AIG1 after two rounds of whole-genome duplications during vertebrate evolution. Then, the AIG family, which was exposed to purifying forces during evolution, lost or gained some of its members in some species. For example, in eutherians, Neognathae, and Percomorphaceae, AIG-L was lost; in contrast, Salmonidae and Cyprinidae acquired additional AIG copies. In conclusion, this study provides a comprehensive molecular phylogenetic analysis of vertebrate AIGs, which can be employed for future functional characterization of AIGs.

Exploring the Origin and Relatedness of Maternal Lineages Through Analysis of Mitochondrial DNA in the Holstein Horse

Maternal lineages are important for the breeding decision in the Holstein horse breed. To investigate the genetic diversity of the maternal lineages and the relationships between founder mares, the maternal inherited mitochondrial genome (except the repetitive part of the non-coding region) of 271 mares representing 75 lineages was sequenced. The sequencing predominantly revealed complete homology in the nucleotide sequences between mares from one lineage with exceptions in 13 lineages, where differences in one to three positions are probably caused by de novo mutations or alternate fixation of heteroplasmy. We found 78 distinct haplotypes that have not yet been described in other breeds. Six of these occurred in two or three different lineages indicating a common ancestry. Haplotypes can be divided into eight clusters with all mares from one lineage belonging to the same cluster. Within a cluster, the average number of pairwise differences ranged from zero to 16.49 suggesting close maternal relationships between these mares. The results showed that the current breeding population originated from at least eight ancestral founder mares.

Environment-Driven Adaptations of Leaf Cuticular Waxes Are Inheritable for Medicago ruthenica

Cuticular waxes covering the plant surface play pivotal roles in helping plants adapt to changing environments. However, it is still not clear whether the responses of plant cuticular waxes to their growing environments are inheritable. We collected seeds of Medicago ruthenica (a perennial legume) populations from 30 growing sites in northern China and examined the variations of leaf cuticular waxes in a common garden experiment. Four wax genes, MrFAR3-1, MrFAR3-2, MrCER1, and MrKCS1, involved in biosynthesis of predominant wax classes (primary alcohol and alkane) and wax precursors, were isolated to test the contributions of genetic variations of the coding sequences (CDS) and the promoter sequences and epigenetic modifications. The plasticity responses of the cuticular waxes were further validated by two stress-modeling experiments (drought and enhancing ultraviolet B). Great variations in total wax coverage and abundance of wax classes or wax compounds were observed among M. ruthenica populations in a common garden experiment. Stress-modeling experiments further validated that M. ruthenica would alter leaf wax depositions under changed growing conditions. The transcriptional levels of the wax genes were positively or negatively correlated with amounts of cuticular waxes. However, the analysis of promoter methylation showed that the methylation level of the promoter region was not associated with their expressions. Although both promoter sequences and CDS showed a number of polymorphic sites, the promoters were not naturally selected and insignificant difference could be observed in the numbers and types of acting elements of the four wax genes among populations. In contrast, the CDS of the wax genes were naturally selected, with a number of missense mutations resulting in alterations of the amino acid as well as their isoelectric points and polarities, which could impact on enzyme function/activity. We conclude that long-term adaptation under certain environments would induce genetic mutation of wax biosynthesis genes, resulting in inheritable alterations of cuticular wax depositions.

Self-Synthesizing Nanorods from Dynamic Combinatorial Libraries against Drug Resistant Cancer

Molecular self-assembly has been widely used to develop nanocarriers for drug delivery. However, most of them have unsatisfactory drug loading capacity (DLC) and the dilemma between stimuli-responsiveness and stability, stagnating their translational process. Herein, we overcame these drawbacks using dynamic combinatorial chemistry. A carrier molecule was spontaneously and quantitatively synthesized, aided by co-self-assembly with a template molecule and an anti-cancer drug doxorubicin (DOX) from a dynamic combinatorial library that was operated by disulfide exchange under thermodynamic control. The highly selective synthesis guaranteed a stable yet pH- and redox- responsive nanocarrier with a maximized DLC of 40.1 % and an enhanced drug potency to fight DOX resistance in vitro and in vivo. Our findings suggested that harnessing the interplay between synthesis and self-assembly in complex chemical systems could yield functional nanomaterials for advanced applications.

Evolution of the Highly Repetitive PEVK Region of Titin Across Mammals

The protein titin plays a key role in vertebrate muscle where it acts like a giant molecular spring. Despite its importance and conservation over vertebrate evolution, a lack of high quality annotations in non-model species makes comparative evolutionary studies of titin challenging. The PEVK region of titin—named for its high proportion of Pro-Glu-Val-Lys amino acids—is particularly difficult to annotate due to its abundance of alternatively spliced isoforms and short, highly repetitive exons. To understand PEVK evolution across mammals, we developed a bioinformatics tool, PEVK_Finder, to annotate PEVK exons from genomic sequences of titin and applied it to a diverse set of mammals. PEVK_Finder consistently outperforms standard annotation tools across a broad range of conditions and improves annotations of the PEVK region in non-model mammalian species. We find that the PEVK region can be divided into two subregions (PEVK-N, PEVK-C) with distinct patterns of evolutionary constraint and divergence. The bipartite nature of the PEVK region has implications for titin diversification. In the PEVK-N region, certain exons are conserved and may be essential, but natural selection also acts on particular codons. In the PEVK-C, exons are more homogenous and length variation of the PEVK region may provide the raw material for evolutionary adaptation in titin function. The PEVK-C region can be further divided into a highly repetitive region (PEVK-CA) and one that is more variable (PEVK-CB). Taken together, we find that the very complexity that makes titin a challenge for annotation tools may also promote evolutionary adaptation.

ALLOZYME VARIATION IN A SNAIL (LITTORINA SAXATILIS)-DECONFOUNDING THE EFFECTS OF MICROHABITAT AND GENE FLOW

It is commonly observed that a restricted gene flow among populations of a species generates genetic differentiation in, for example, allozyme markers. However, recent studies suggest that microhabitat-specific variation may contribute to the total differentiation. To appreciate the relative contributions of geographic variation and habitat-specific variation, we sampled 42 subpopulations of the intertidal snail Littorina saxatilis from three different microhabitats (boulders, low and high rocky intertidal) on five small islands within a distance of 15 km. We used a modified orthogonal version of Nei's gene diversity analysis with a modified analysis of variance (ANOVA) that estimated the significance of habitat and geographic separation and the interaction between them. Between subpopulation differentiation (G_ST ) was usually in the range of 5% to 10% but was exceptionally high in one locus (Aat; 53%). Genetic differentiation attributable to different habitats accounted for 10% to 81% (mean 35%) of the between subpopulation variation and was significant (P < 0.05) in six loci. Differentiation due to geographic separation accounted for 11% to 61% (mean 36%) and was significant (P < 0.05) in seven loci. Furthermore, three loci showed interactions between habitat and island, suggesting varying effects of habitat in different islands. Microhabitat-specific variation, probably through spatially varying fitness, seems particularly likely in Aat and Pgm-2. Moderate levels of habitat associated variation added to the observed differentiation due to gene flow in Pgi, Pnp, and Pgm-1, whereas in the remaining three loci either the habitat effect was confused by strong habitat-island interaction (Ark) or was virtually absent (Pep and Mpi).

HISTORICAL DEMOGRAPHY AND PRESENT DAY POPULATION STRUCTURE OF THE GREENFINCH, CARDUEUS CHLORIS-AN ANALYSIS OF mtDNA CONTROL-REGION SEQUENCES

Genetic variability within and among 10 geographically distinct populations of Greenfinches (Carduelis chloris) was assayed by directly sequencing a 637 BP part of the mtDNA control region from 194 individuals. Thirteen variable positions defined 18 haplotypes with a maximum sequence divergence of 0.8%. Haplotype (h = 0.28-0.77) and nucleotide (π = 0.058-0.17%) diversities within populations were low, and decreased with increasing latitude (h:r_s = -0.81; π: r_s = -0.89). The distribution of pairwise nucleotide differences fit better with expectations of a "sudden expansion" than of an "equilibrium" model, and the estimates of long term effective population sizes were considerably lower than current census estimates, especially in northern European samples. Selection is an unlikely cause of observed patterns because the distribution of variability conformed to expectations of neutral infinite alleles model and haplotype diversity across populations was positively correlated with heterozygosity (H_E ) in nuclear genes (r_s = 0.74, P < 0.05). Hence, a recent bottleneck, followed by serial bottlenecking during the process of post-Pleistocene recolonization of northern Europe, together with recent population expansion provide a plausible explanation for the low genetic diversity in the north. Genetic distances among populations showed a clear pattern of isolation-by-distance, and 14% of the haplotypic variation was among populations, the rest being distributed among individuals within populations. In accordance with allozyme and morphological data, a hierarchical analysis of nucleotide diversity recognized southern European populations as distinct from northern European ones. However, the magnitude of divergence in mtDNA, allozymes and morphology were highly dissimilar (morphology > mtDNA > allozymes).

Diversifying selection of the anthocyanin biosynthetic downstream gene UFGT accelerates floral diversity of island Scutellaria species

BACKGROUND

Adaptive divergence, which usually explains rapid diversification within island species, might involve the positive selection of genes. Anthocyanin biosynthetic pathway (ABP) genes are important for floral diversity, and are related to stress resistance and pollination, which could be responsible for species diversification. Previous studies have shown that upstream genes of ABP are subject to selective constraints and have a slow evolutionary rate, while the constraints on downstream genes are lower.

RESULTS

In this study, we confirmed these earlier observations of heterogeneous evolutionary rate in upstream gene CHS and the downstream gene UFGT, both of which were expressed in Scutellaria sp. inflorescence buds. We found a higher evolutionary rate and positive selection for UFGT. The codons under positive selection corresponded to the diversified regions, and the presence or absence of an α-helix might produce conformation changes in the Rossmann-like fold structure. The significantly high evolutionary rates for UFGT genes in Taiwanese lineages suggested rapid accumulation of amino acid mutations in island species. The results showed positive selection in closely related species and explained the high diversification of floral patterns in these recently diverged species. In contrast, non-synonymous mutation rate decreases in long-term divergent species could reduce mutational load and prevent the accumulation of deleterious mutations.

CONCLUSIONS

Together with the positive selection of transcription factors for ABP genes described in a previous study, these results confirmed that positive selection takes place at a translational level and suggested that the high divergence of ABP genes is related to the floral diversity in island Scutellaria species.

A Mutation Model from First Principles of the Genetic Code

The paper presents a neutral Codons Probability Mutations (CPM) model of molecular evolution and genetic decay of an organism. The CPM model uses a Markov process with a 20-dimensional state space of probability distributions over amino acids. The transition matrix of the Markov process includes the mutation rate and those single point mutations compatible with the genetic code. This is an alternative to the standard Point Accepted Mutation (PAM) and BLOcks of amino acid SUbstitution Matrix (BLOSUM). Genetic decay is quantified as a similarity between the amino acid distribution of proteins from a (group of) species on one hand, and the equilibrium distribution of the Markov chain on the other. Amino acid data for the eukaryote, bacterium, and archaea families are used to illustrate how both the CPM and PAM models predict their genetic decay towards the equilibrium value of 1. A family of bacteria is studied in more detail. It is found that warm environment organisms on average have a higher degree of genetic decay compared to those species that live in cold environments. The paper addresses a new codon-based approach to quantify genetic decay due to single point mutations compatible with the genetic code. The present work may be seen as a first approach to use codon-based Markov models to study how genetic entropy increases with time in an effectively neutral biological regime. Various extensions of the model are also discussed.

Comprehensive annotation of Glossina pallidipes salivary gland hypertrophy virus from Ethiopian tsetse flies: a proteogenomics approach

Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; family Hytrosaviridae) can establish asymptomatic and symptomatic infection in its tsetse fly host. Here, we present a comprehensive annotation of the genome of an Ethiopian GpSGHV isolate (GpSGHV-Eth) compared with the reference Ugandan GpSGHV isolate (GpSGHV-Uga; GenBank accession number EF568108). GpSGHV-Eth has higher salivary gland hypertrophy syndrome prevalence than GpSGHV-Uga. We show that the GpSGHV-Eth genome has 190 291 nt, a low G+C content (27.9 %) and encodes 174 putative ORFs. Using proteogenomic and transcriptome mapping, 141 and 86 ORFs were mapped by transcripts and peptides, respectively. Furthermore, of the 174 ORFs, 132 had putative transcriptional signals [TATA-like box and poly(A) signals]. Sixty ORFs had both TATA-like box promoter and poly(A) signals, and mapped by both transcripts and peptides, implying that these ORFs encode functional proteins. Of the 60 ORFs, 10 ORFs are homologues to baculovirus and nudivirus core genes, including three per os infectivity factors and four RNA polymerase subunits (LEF4, 5, 8 and 9). Whereas GpSGHV-Eth and GpSGHV-Uga are 98.1 % similar at the nucleotide level, 37 ORFs in the GpSGHV-Eth genome had nucleotide insertions (n = 17) and deletions (n = 20) compared with their homologues in GpSGHV-Uga. Furthermore, compared with the GpSGHV-Uga genome, 11 and 24 GpSGHV ORFs were deleted and novel, respectively. Further, 13 GpSGHV-Eth ORFs were non-canonical; they had either CTG or TTG start codons instead of ATG. Taken together, these data suggest that GpSGHV-Eth and GpSGHV-Uga represent two different lineages of the same virus. Genetic differences combined with host and environmental factors possibly explain the differential GpSGHV pathogenesis observed in different G. pallidipes colonies.

Identification of dopamine receptors across the extant avian family tree and analysis with other clades uncovers a polyploid expansion among vertebrates

Dopamine is an important central nervous system transmitter that functions through two classes of receptors (D1 and D2) to influence a diverse range of biological processes in vertebrates. With roles in regulating neural activity, behavior, and gene expression, there has been great interest in understanding the function and evolution dopamine and its receptors. In this study, we use a combination of sequence analyses, microsynteny analyses, and phylogenetic relationships to identify and characterize both the D1 (DRD1A, DRD1B, DRD1C, and DRD1E) and D2 (DRD2, DRD3, and DRD4) dopamine receptor gene families in 43 recently sequenced bird genomes representing the major ordinal lineages across the avian family tree. We show that the common ancestor of all birds possessed at least seven D1 and D2 receptors, followed by subsequent independent losses in some lineages of modern birds. Through comparisons with other vertebrate and invertebrate species we show that two of the D1 receptors, DRD1A and DRD1B, and two of the D2 receptors, DRD2 and DRD3, originated from a whole genome duplication event early in the vertebrate lineage, providing the first conclusive evidence of the origin of these highly conserved receptors. Our findings provide insight into the evolutionary development of an important modulatory component of the central nervous system in vertebrates, and will help further unravel the complex evolutionary and functional relationships among dopamine receptors.

Evolution and Structural Analyses of Glossina morsitans (Diptera; Glossinidae) Tetraspanins

Tetraspanins are important conserved integral membrane proteins expressed in many organisms. Although there is limited knowledge about the full repertoire, evolution and structural characteristics of individual members in various organisms, data obtained so far show that tetraspanins play major roles in membrane biology, visual processing, memory, olfactory signal processing, and mechanosensory antennal inputs. Thus, these proteins are potential targets for control of insect pests. Here, we report that the genome of the tsetse fly, Glossina morsitans (Diptera: Glossinidae) encodes at least seventeen tetraspanins (GmTsps), all containing the signature features found in the tetraspanin superfamily members. Whereas six of the GmTsps have been previously reported, eleven could be classified as novel because their amino acid sequences do not map to characterized tetraspanins in the available protein data bases. We present a model of the GmTsps by using GmTsp42Ed, whose presence and expression has been recently detected by transcriptomics and proteomics analyses of G. morsitans. Phylogenetically, the identified GmTsps segregate into three major clusters. Structurally, the GmTsps are largely similar to vertebrate tetraspanins. In view of the exploitation of tetraspanins by organisms for survival, these proteins could be targeted using specific antibodies, recombinant large extracellular loop (LEL) domains, small-molecule mimetics and siRNAs as potential novel and efficacious putative targets to combat African trypanosomiasis by killing the tsetse fly vector.

Reconstructing a B-Cell Clonal Lineage. II. Mutation, Selection, and Affinity Maturation

Affinity maturation of the antibody response is a fundamental process in adaptive immunity during which B-cells activated by infection or vaccination undergo rapid proliferation accompanied by the acquisition of point mutations in their rearranged immunoglobulin (Ig) genes and selection for increased affinity for the eliciting antigen. The rate of somatic hypermutation at any position within an Ig gene is known to depend strongly on the local DNA sequence, and Ig genes have region-specific codon biases that influence the local mutation rate within the gene resulting in increased differential mutability in the regions that encode the antigen-binding domains. We have isolated a set of clonally related natural Ig heavy chain-light chain pairs from an experimentally infected influenza patient, inferred the unmutated ancestral rearrangements and the maturation intermediates, and synthesized all the antibodies using recombinant methods. The lineage exhibits a remarkably uniform rate of improvement of the effective affinity to influenza hemagglutinin (HA) over evolutionary time, increasing 1000-fold overall from the unmutated ancestor to the best of the observed antibodies. Furthermore, analysis of selection reveals that selection and mutation bias were concordant even at the level of maturation to a single antigen. Substantial improvement in affinity to HA occurred along mutationally preferred paths in sequence space and was thus strongly facilitated by the underlying local codon biases.

Detecting natural selection in genomic data

The past fifty years have seen the development and application of numerous statistical methods to identify genomic regions that appear to be shaped by natural selection. These methods have been used to investigate the macro- and microevolution of a broad range of organisms, including humans. Here, we provide a comprehensive outline of these methods, explaining their conceptual motivations and statistical interpretations. We highlight areas of recent and future development in evolutionary genomics methods and discuss ongoing challenges for researchers employing such tests. In particular, we emphasize the importance of functional follow-up studies to characterize putative selected alleles and the use of selection scans as hypothesis-generating tools for investigating evolutionary histories.

Genome-scale analysis of demographic history and adaptive selection

One of the main topics in population genetics is identification of adaptive selection among populations. For this purpose, population history should be correctly inferred to evaluate the effect of random drift and exclude it in selection identification. With the rapid progress in genomics in the past decade, vast genome-scale variations are available for population genetic analysis, which however requires more sophisticated models to infer species' demographic history and robust methods to detect local adaptation. Here we aim to review what have been achieved in the fields of demographic modeling and selection detection. We summarize their rationales, implementations, and some classical applications. We also propose that some widely-used methods can be improved in both theoretical and practical aspects in near future.

Adaptive evolution of M3 lysin--a candidate gamete recognition protein in the Mytilus edulis species complex

Marine invertebrate gamete recognition proteins (GRPs) are classic examples of rapid adaptive evolution of reproductive proteins, and hybridizing Mytilus blue mussels allow us to study the evolution of GRPs during speciation following secondary contact. Even with frequent hybridization, positive selection drives divergence of M7 lysin, one of the three Mytilus egg vitelline envelope (VE) lysins. Mytilus trossulus and M. edulis form a broad hybrid zone in the Canadian Maritimes and eastern Maine, isolated by strong (but partial) gamete incompatibility. M7 lysin, however, is an unlikely GRP controlling this gametic incompativility, as earlier studies showed either weak or no positive selection and extensive introgression between the two species. We used reverse transcriptase-polymerase chain reaction and cloned several alleles of M3 lysin, a potent VE lysin encoded by a nonhomologous gene whose evolution has not been studied. McDonald-Kreitman and HKA tests reveal strong positive selection, which PAML branch-site models detect in 19.7% of the codons. Protein structure predictions show that replacements map exclusively to one face of the carbohydrate recognition domain (CRD) of this C-type lectin, with codons under positive selection localizing to CRD regions known to control ligand specificity. Polymorphism/divergence analyses show that selective sweep has purged M. edulis but not M. trossulus of polymorphism, and unique to M3 is an absence of fixed substitutions and broad haplotype sharing between M. edulis and Mediterranean M. galloprovincialis. Taken together, these results suggest that different lysins serve as GRPs in different Mytilus hybrid zones, with M3 likely co-opted to play this role in the western Atlantic.

Perspectives of genomic diversification and molecular recombination towards R-gene evolution in plants

Plants are under strong evolutionary pressure in developing new and noble R genes to recognize pathogen avirulence (avr) determinants and bring about stable defense for generation after generations. Duplication, sequence variation by mutation, disparity in the length and structure of leucine rich repeats etc., causes tremendous variations within and among R genes in a plant thereby developing diverse recognitional specificity suitable enough for defense against new pathogens. Recent studies on genome sequencing, diversity and population genetics in different plants have thrown new insights on the molecular evolution of these genes. Tandem and segmental duplication are important factors in R gene abundance as inferred from the distribution of major nucleotide binding site-leucine rich repeats (NBS-LRRs) type R-genes in plant genomes. Likewise, R-gene evolution is also thought to be facilitated by cluster formation thereby causing recombination and sequence exchange and resulting in haplotypic diversity. Population studies have further proven that balancing selection is responsible for the maintenance of allelic diversity in R genes. In this review, we emphasize and discuss on improved perspectives towards the molecular mechanisms and selection pressure responsible for the evolution of NBS-LRR class resistance genes in plants.

Statistical tests of neutrality in the age of weak selection

Why review statistical tests of neutrality at a time when pan-selectionists and pan-neutralists alike seem to have been replaced by weak selectionists? First, we still don't actually know how variation is maintained at the molecular level; and second, tests of neutrality have a utility for evolutionary biologists beyond the neutralist/selectionist debate. New tests and variations on the existing tests are arising practically every month. From the complementary viewpoints of an empiricist and a theoretician, we sample the recent literature on tests of statistical neutrality and discuss the motivations, applications, assumptions, interpretations and future directions of these tests.

Weighing the evidence for adaptation at the molecular level

The abundance of genome polymorphism and divergence data has provided unprecedented insight into how mutation, drift and natural selection shape genome evolution. Application of the McDonald-Kreitman (MK) test to such data indicates a pervasive influence of positive selection, particularly in Drosophila species. However, evidence for positive selection in other species ranging from yeast to humans is often weak or absent. Although evidence for positive selection could be obscured in some species, there is also reason to believe that the frequency of adaptive substitutions could be overestimated as a result of epistatic fitness effects or hitchhiking of deleterious mutations. Based on these considerations it is argued that the common assumption of independence among sites must be relaxed before abandoning the neutral theory of molecular evolution.

Fusion-fission experiments in Aphidius: evolutionary split without isolation in response to environmental bimodality

Studying host-based divergence naturally maintained by a balance between selection and gene flow can provide valuable insights into genetic underpinnings of host adaptation and ecological speciation in parasites. Selection-gene flow balance is often postulated in sympatric host races, but direct experimental evidence is scarce. In this study, we present such evidence obtained in host races of Aphidius ervi, an important hymenopteran agent of biological control of aphids in agriculture, using a novel fusion-fission method of gene flow perturbation. In our study, between-race genetic divergence was obliterated by means of advanced hybridisation, followed by a multi-generation exposure of the resulting genetically uniform hybrid swarm to a two-host environment. This fusion-fission procedure was implemented under two contrasting regimes of between-host gene flow in two replicated experiments involving different racial pairs. Host-based genetic fission in response to environmental bimodality occurred in both experiments in as little as six generations of divergent adaptation despite continuous gene flow. We demonstrate that fission recovery of host-based divergence evolved faster and hybridisation-induced linkage disequilibrium decayed slower under restricted (6.7%) compared with unrestricted gene flow, directly pointing at a balance between gene flow and divergent selection. We also show, in four separate tests, that random drift had no or little role in the observed genetic split. Rates and patterns of fission divergence differed between racial pairs. Comparative linkage analysis of these differences is currently under way to test for the role of genomic architecture of adaptation in ecology-driven divergent evolution.

Pervasive natural selection in the Drosophila genome?

Over the past four decades, the predominant view of molecular evolution saw little connection between natural selection and genome evolution, assuming that the functionally constrained fraction of the genome is relatively small and that adaptation is sufficiently infrequent to play little role in shaping patterns of variation within and even between species. Recent evidence from Drosophila, reviewed here, suggests that this view may be invalid. Analyses of genetic variation within and between species reveal that much of the Drosophila genome is under purifying selection, and thus of functional importance, and that a large fraction of coding and noncoding differences between species are adaptive. The findings further indicate that, in Drosophila, adaptations may be both common and strong enough that the fate of neutral mutations depends on their chance linkage to adaptive mutations as much as on the vagaries of genetic drift. The emerging evidence has implications for a wide variety of fields, from conservation genetics to bioinformatics, and presents challenges to modelers and experimentalists alike.

Direct evidence that genetic variation in glycerol-3-phosphate and malate dehydrogenase genes (Gpdh and Mdh1) affects adult ethanol tolerance in Drosophila melanogaster

Many studies of alcohol adaptation in Drosophila melanogaster have focused on the Adh polymorphism, yet the metabolic elimination of alcohol should involve many enzymes and pathways. Here we evaluate the effects of glycerol-3-phosphate dehydrogenase (Gpdh) and cytosolic malate dehydrogenase (Mdh1) genotype activity on adult tolerance to ethanol. We have created a set of P-element-excision-derived Gpdh, Mdh1, and Adh alleles that generate a range of activity phenotypes from full to zero activity. Comparisons of paired Gpdh genotypes possessing 10 and 60% normal activity and 66 and 100% normal activity show significant effects where higher activity increases tolerance. Mdh1 null allele homozygotes show reductions in tolerance. We use piggyBac FLP-FRT site-specific recombination to create deletions and duplications of Gpdh. Duplications show an increase of 50% in activity and an increase of adult tolerance to ethanol exposure. These studies show that the molecular polymorphism associated with GPDH activity could be maintained in natural populations by selection related to adaptation to alcohols. Finally, we examine the interactions between activity genotypes for Gpdh, Mdh1, and Adh. We find no significant interlocus interactions. Observations on Mdh1 in both Gpdh and Adh backgrounds demonstrate significant increases in ethanol tolerance with partial reductions (50%) in cytosolic MDH activity. This observation strongly suggests the operation of pyruvate-malate and, in particular, pyruvate-citrate cycling in adaptation to alcohol exposure. We propose that an understanding of the evolution of tolerance to alcohols will require a system-level approach, rather than a focus on single enzymes.

Turning a hobby into a job: how duplicated genes find new functions

Gene duplication provides raw material for functional innovation. Recent advances have shed light on two fundamental questions regarding gene duplication: which genes tend to undergo duplication? And how does natural selection subsequently act on them? Genomic data suggest that different gene classes tend to be retained after single-gene and whole-genome duplications. We also know that functional differences between duplicate genes can originate in several different ways, including mutations that directly impart new functions, subdivision of ancestral functions and selection for changes in gene dosage. Interestingly, in many cases the 'new' function of one copy is a secondary property that was always present, but that has been co-opted to a primary role after the duplication.

Genetic hitchhiking in a subdivided population of Mytilus edulis

BACKGROUND

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average FST~0.007), a strong genetic differentiation is observed at a single locus (FST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.

RESULTS

We first showed that introgression of M. galloprovincialis alleles was very weak in both populations and did not significantly affect their differentiation. Secondly, we observed the genetic signature of a selective sweep within both M. edulis populations in the form of a star-shaped clade of alleles. This clade was nearly fixed in the North Sea and was segregating at a moderate frequency in the Bay of Biscay, explaining their genetic differentiation. Incomplete fixation reveals that selection was not direct on the locus but that the studied sequence recombined with a positively selected allele at a linked locus while it was on its way to fixation. Finally, using a deterministic model we showed that the wave of advance of a favourable allele at a linked locus, when crossing a strong enough barrier to gene flow, generates a step in neutral allele frequencies comparable to the step observed between the two M. edulis populations at the outlier locus. In our case, the position of the barrier is now materialised by a large patch of heterospecific M. galloprovincialis populations.

CONCLUSION

High FST outlier loci are usually interpreted as being the consequence of ongoing divergent local adaptation. Combining models and data we show that among-population differentiation can also dramatically increase following a selective sweep in a structured population. Our study illustrates how a striking geographical pattern of neutral diversity can emerge from past indirect hitchhiking selection in a structured population. NOTE: Nucleotide sequences reported in this paper are available in the GenBanktrade mark database under the accession numbers EU684165 - EU684228.

Phylogeography of stable fly (Diptera: Muscidae) estimated by diversity at ribosomal 16S and cytochrome oxidase I mitochondrial genes

The blood-feeding cosmopolitan stable fly, Stomoxys calcitrans L. (Diptera: Muscidae), is thought to disperse rapidly and widely, and earlier studies of allozyme variation were consistent with high vagility in this species. The geographic origins of New World populations are unknown. Diversity at mitochondrial loci r16S and cytochrome oxidase I was examined in 277 stable flies from 11 countries, including five zoogeographical regions. Of 809 nucleotides, 174 were polymorphic and 133 were parsimony informative. Seventy-six haplotypes were found in frequencies consistent with the Wright-Fisher infinite allele model. None were shared among four or more zoogeographical regions. The null hypothesis of mutation neutrality was not rejected, thereby validating the observed distribution. Fifty-nine haplotypes were singular, eight were private and confined to the Old World, and three of 76 haplotypes were shared between the Old and New World. Only 19 haplotypes were found in the New World, 14 of which were singletons. Haplotype and nucleotide diversities were heterogeneous among countries and regions. The most diversity was observed in sub-Saharan Africa. Regional differentiation indices were C(RT) = 0.26 and N(RT) = 0.31, indicating populations were highly structured macrogeographically. Palearctic and New World flies were the least differentiated from each other. There were strong genetic similarities among populations in the Nearctic, Neotropical, and Palearctic regions, and it is most likely that New World populations were derived from the Palearctic after 1492 CE, in the colonial era.

Adaptive evolution of color vision as seen through the eyes of butterflies

Butterflies and primates are interesting for comparative color vision studies, because both have evolved middle- (M) and long-wavelength- (L) sensitive photopigments with overlapping absorbance spectrum maxima (lambda(max) values). Although positive selection is important for the maintenance of spectral variation within the primate pigments, it remains an open question whether it contributes similarly to the diversification of butterfly pigments. To examine this issue, we performed epimicrospectrophotometry on the eyes of five Limenitis butterfly species and found a 31-nm range of variation in the lambda(max) values of the L-sensitive photopigments (514-545 nm). We cloned partial Limenitis L opsin gene sequences and found a significant excess of replacement substitutions relative to polymorphisms among species. Mapping of these L photopigment lambda(max) values onto a phylogeny revealed two instances within Lepidoptera of convergently evolved L photopigment lineages whose lambda(max) values were blue-shifted. A codon-based maximum-likelihood analysis indicated that, associated with the two blue spectral shifts, four amino acid sites (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) have evolved substitutions in parallel and exhibit significant d(N)/d(S) >1. Homology modeling of the full-length Limenitis arthemis astyanax L opsin placed all four substitutions within the chromophore-binding pocket. Strikingly, the Ser137Ala substitution is in the same position as a site that in primates is responsible for a 5- to 7-nm blue spectral shift. Our data show that some of the same amino acid sites are under positive selection in the photopigments of both butterflies and primates, spanning an evolutionary distance >500 million years.

Molecular paleoscience: systems biology from the past

Experimental paleomolecular biology, paleobiochemistry, and paleogenetics are closely related emerging fields that infer the sequences of ancient genes and proteins from now-extinct organisms, and then resurrect them for study in the laboratory. The goal of paleogenetics is to use information from natural history to solve the conundrum of modern genomics: How can we understand deeply the function of biomolecular structures uncovered and described by modern chemical biology? Reviewed here are the first 20 cases where biomolecular resurrections have been achieved. These show how paleogenetics can lead to an understanding of the function of biomolecules, analyze changing function, and put meaning to genomic sequences, all in ways that are not possible with traditional molecular biological studies.

Hypotheses for the origin and early evolution of triterpenoid cyclases

Hopanes and steranes are found almost universally in the sedimentary rock record where they often are used as proxies for aerobic organisms, metabolisms, and environments. In order to interpret ancient lipid signatures confidently we require a complementary understanding of how these modern biochemical pathways evolved since their conception. For example, generally it has been assumed that hopanoid biosynthesis was an evolutionary predecessor to steroid biosynthesis. Here we re-evaluate this assumption. Using a combined phylogenetic and biochemical perspective, we address the evolution of polycyclic triterpenoid biosynthesis and suggest several constraints on using these molecules as aerobic biomarkers. Amino acid sequence data show that the enzymes responsible for polycyclic triterpenoid biosynthesis (i.e. squalene and 2,3-oxidosqualene cyclases) are homologous. Numerous conserved domains correspond to active sites in the enzymes that are required to complete the complex cyclization reaction. From these sites we develop an evolutionary analysis of three independent characters to explain the evolution of the major classes of polycyclic triterpenoids. These characters are: (i) the number of unfavourable anti-Markovnikov ring closures, (ii) all-chair (CCC) or chair-boat-chair (CBC) substrate conformation, and (iii) the choice between squalene and 2,3-oxidosqualene as the substrate. We use these characters to construct four competing phylogenies to describe the evolution of polycyclic triterpenoid biosynthesis. The analysis suggests that malabaricanoids would be the most ancient polycyclic triterpenoids. The two most parsimonious evolutionary trees are the ones in which hopanoid and steroid cyclases diverged from a common ancestor. The transition from a CCC- to CBC-fold marks the major divergence in the evolution of these pathways, and it is diagnosable in the geological record. However, this transition does not require the simultaneous adoption of the aerobic substrate, 2,3-oxidosqualene, because these characters are controlled by independent parts of the enzyme.

Evolution of a tumorigenic property conferred by glycophosphatidyl-inositol membrane anchors of carcinoembryonic antigen gene family members during the primate radiation

GPI membrane anchors of cell surface glycoproteins have been shown to confer functional properties that are different from their transmembrane (TM)-anchored counterparts. For the human carcinoembryonic antigen (CEA) family, a subfamily of the immunoglobulin superfamily, conversion of the mode of membrane linkage from TM to GPI confers radical changes in function: from tumor suppression or neutrality toward inhibition of differentiation and anoikis and distortion of tissue architecture, thereby contributing to tumorigenesis. We show here that GPI anchorage in the CEA family evolved twice independently in primates, very likely from more primitive TM anchors, by different packages of mutations. Both mutational packages, one package found in many primates, including humans, and a second, novel package found only in the Cebidae radiation of New World monkeys, give rise to efficiently processed GPI-linked proteins. Both types of GPI anchors mediate inhibition of cell differentiation. The estimated rate of nonsynonymous mutations (Ka) in the anchor-determining domain for conversion from TM to GPI anchorage in the CEA family that were fixed during evolution in these primates is 7 times higher than the average Ka in primates, indicating positive selection. These results suggest therefore that the functional changes mediated by CEA GPI anchors, including the inhibition of differentiation and anoikis, could be adaptive and advantageous.

Patterns of genetic diversity and differentiation in the tsetse fly Glossina morsitans morsitans Westwood populations in East and southern Africa

Genetic diversity and differentiation within and among nine G. morsitans morsitans populations from East and southern Africa was assessed by examining variation at seven microsatellite loci and a mitochondrial locus, cytochrome oxidase (COI). Mean COI diversity within populations was 0.63+/-0.33 and 0.81 taken over all populations. Diversities averaged over microsatellite loci were high (mean number of alleles/locus>or=7.4; mean HE>or=65%) in all populations. Diversities averaged across populations were greater in East Africa (mean number of alleles=22+/-2.6; mean he=0.773+/-0.033) than in southern Africa (mean number of alleles=18.7+/-4.0; mean he=0.713+/-0.072). Differentiation among all populations was highly significant (RST=0.25, FST=0.132). Nei's Gij statistics were 0.09 and 0.19 within regions for microsatellites and mitochondria, respectively; between regions, Gij was 0.14 for microsatellites and 0.23 for mitochondria. GST among populations was 0.23 for microsatellite loci and 0.40 for mitochondria. The F, G and R statistics indicate highly restricted gene flow among G. m. morsitans populations separated over geographic scales of 12-917 km.

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.

Integrating biogeographic and genetic approaches to investigate the history of bioluminescent colour alleles in the Jamaican click beetle, Pyrophorus plagiophthalamus

Bioluminescent colour in the Jamaican click beetle, Pyrophorus plagiophthalamus, is an ideal system for studies moving from gene to landscape to gain a holistic understanding of the molecular, ecological, and historical bases for adaptation. Previous studies have established the genetics of bioluminescent colour variation in the beetle to the level of the nucleotide base pair in the target gene luciferase. Three different luciferase colour alleles affecting ventral light organ colour [yellow-green (vYG), yellow (vYE), and orange (vOR)] were found segregating in P. plagiophthalamus populations. These alleles differ from each other in a number of replacement mutations (14 total), the majority of which (11) have a measurable effect on colour. Phylogenetic analysis revealed a long-term adaptive trend on Jamaica towards longer wavelength bioluminescence, culminating in the most recently derived vOR allele. Here, we further investigate the historical and geographic context of adaptive colour evolution by testing a vicariance model for the origins of the extant ventral light organ polymorphism: that the vOR allele arose and differentiated in an isolated deme on the east side of Jamaica before spreading westward. Comparisons of colour phenotypes, luciferase coding sequences, the third intron of the gene, mtDNA, and microsatellite data provided evidence for past population subdivision on Jamaica and ongoing gene flow, as has been found for other island endemics. However, the pattern of differentiation supported the allopatric divergence of vYG and vYE alleles. The vOR gene appears to have arisen relatively recently from a vYE precursor and postdates the period of major biogeographic isolation. We discuss the implications of the results for discerning ecological causation in the adaptive sequence from nucleotide to landscape to population change for bioluminescent colour.

Genetic and morphological diversity in populations of Nucella lapillus (L.; neogastropoda) in response to tributyltin contamination

During the past three decades, North Atlantic populations of dogwhelks have been severely reduced in numbers, due to imposex and female sterility caused by TBT. We examined the relationship between the known history of female sterility and the present genetic diversity (measured by five allozyme-loci) and phenotypic variance (measured by nine quantitative characters) in six populations of Nucella lapillus along the coast of Norway. The environment of one of the populations was severely polluted with a vas deferens sequence index (VDSI) of 4.31, resulting in 24% sterile females. The environment of four of the other five populations was moderately polluted (VDSI 3.31-3.97, 0-1.6% sterile females), while the last population was unpolluted (VDSI 0.02, 0% sterile females). The six populations did not differ significantly with respect to genetic diversity. However, the population from the most polluted area had greater phenotypic variation in shell size, compared to all other populations. The results indicate that although TBT has had severe local negative effects on N. lapillus populations, it has not caused a reduction in long-term adaptive potential.

Estimating the genomewide rate of adaptive protein evolution in Drosophila

When polymorphism and divergence data are available for multiple loci, extended forms of the McDonald-Kreitman test can be used to estimate the average proportion of the amino acid divergence due to adaptive evolution--a statistic denoted alpha. But such tests are subject to many biases. Most serious is the possibility that high estimates of alpha reflect demographic changes rather than adaptive substitution. Testing for between-locus variation in alpha is one possible way of distinguishing between demography and selection. However, such tests have yielded contradictory results, and their efficacy is unclear. Estimates of alpha from the same model organisms have also varied widely. This study clarifies the reasons for these discrepancies, identifying several method-specific biases in widely used estimators and assessing the power of the methods. As part of this process, a new maximum-likelihood estimator is introduced. This estimator is applied to a newly compiled data set of 115 genes from Drosophila simulans, each with each orthologs from D. melanogaster and D. yakuba. In this way, it is estimated that alpha approximately 0.4+/-0.1, a value that does not vary substantially between different loci or over different periods of divergence. The implications of these results are discussed.

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.

Triglyceride pools, flight and activity variation at the Gpdh locus in Drosophila melanogaster

We have created a set of P-element excision-derived Gpdh alleles that generate a range of GPDH activity phenotypes ranging from zero to full activity. By placing these synthetic alleles in isogenic backgrounds, we characterize the effects of minor and major activity variation on two different aspects of Gpdh function: the standing triglyceride pool and glycerol-3-phosphate shuttle-assisted flight. We observe small but statistically significant reductions in triglyceride content for adult Gpdh genotypes possessing 33-80% reductions from normal activity. These small differences scale to a notable proportion of the observed genetic variation in triglyceride content in natural populations. Using a tethered fly assay to assess flight metabolism, we observed that genotypes with 100 and 66% activity exhibited no significant difference in wing-beat frequency (WBF), while activity reductions from 60 to 10% showed statistically significant reductions of approximately 7% in WBF. These studies show that the molecular polymorphism associated with GPDH activity could be maintained in natural populations by selection in the triglyceride pool.

A comparison of variation between a MHC pseudogene and microsatellite loci of the little greenbul (Andropadus virens)

Background

We investigated genetic variation of a major histcompatibility complex (MHC) pseudogene (Anvi-DAB1) in the little greenbul (Andropadus virens) from four localities in Cameroon and one in Ivory Coast, West Africa. Previous microsatellite and mitochondrial DNA analyses had revealed little or no genetic differentiation among Cameroon localities but significant differentiation between localities in Cameroon and Ivory Coast.

Results

Levels of genetic variation, heterozygosity, and allelic diversity were high for the MHC pseudogene in Cameroon. Nucleotide diversity of the MHC pseudogene in Cameroon and Ivory Coast was comparable to levels observed in other avian species that have been studied for variation in nuclear genes. An excess of rare variants for the MHC pseudogene was found in the Cameroon population, but this excess was not statistically significant. Pairwise measures of population differentiation revealed high divergence between Cameroon and Ivory Coast for microsatellites and the MHC locus, although for the latter distance measures were much higher than the comparable microsatellite distances.

Conclusion

We provide the first ever comparison of variation in a putative MHC pseudogene to variation in neutral loci in a passerine bird. Our results are consistence with the action of neutral processes on the pseudogene and suggest they can provide an independent perspective on demographic history and population substructure.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

The rank ordering of genotypic fitness values predicts genetic constraint on natural selection on landscapes lacking sign epistasis

Sewall Wright's genotypic fitness landscape makes explicit one mechanism by which epistasis for fitness can constrain evolution by natural selection. Wright distinguished between landscapes possessing multiple fitness peaks and those with only a single peak and emphasized that the former class imposes substantially greater constraint on natural selection. Here I present novel formalism that more finely partitions the universe of possible fitness landscapes on the basis of the rank ordering of their genotypic fitness values. In this report I focus on fitness landscapes lacking sign epistasis (i.e., landscapes that lack mutations the sign of whose fitness effect varies epistatically), which constitute a subset of Wright's single peaked landscapes. More than one fitness rank ordering lacking sign epistasis exists for L > 2 (where L is the number of interacting loci), and I find that a highly statistically significant effect exists between landscape membership in fitness rank-ordering partition and two different proxies for genetic constraint, even within this subset of landscapes. This statistical association is robust to population size, permitting general inferences about some of the characteristics of fitness rank orderings responsible for genetic constraint on natural selection.

Translational effects of differential codon usage among intragenic domains of new genes in Drosophila

Evolved codon usages often pose a technical challenge over the expressing of eukaryotic genes in microbial systems because of changed translational machinery. In the present study, we investigated the translational effects of intragenic differential codon usage on the expression of the new Drosophila gene, jingwei (jgw), a chimera derived from two unrelated parental genes: Ymp and Adh. We found that jgw possesses a strong intragenic differential usage of synonymous codons, i.e. the Adh-derived C-domain has a significantly higher codon bias than that of the Ymp-derived N-domain (P=0.0023 by t-test). Additional evolutionary analysis revealed the heterogeneous distribution of rare codons, implicating its role in gene regulation and protein translation. The in vitro expression of jgw further demonstrated that the heterogeneous distribution of rare codons has played a role in regulating gene expression, particularly, affecting the quality of protein translation.

Planetary systems biology

Combining paleogenetics, protein engineering, synthetic biology, and metabolic modeling, a planetary biology perspective is brought to bear on adaptive evolutionary events in ancient bacteria.

Robustness, evolvability, and neutrality

Biological systems, from macromolecules to whole organisms, are robust if they continue to function, survive, or reproduce when faced with mutations, environmental change, and internal noise. I focus here on biological systems that are robust to mutations and ask whether such systems are more or less evolvable, in the sense that they can acquire novel properties. The more robust a system is, the more mutations in it are neutral, that is, without phenotypic effect. I argue here that such neutral change--and thus robustness--can be a key to future evolutionary innovation, if one accepts that neutrality is not an essential feature of a mutation. That is, a once neutral mutation may cause phenotypic effects in a changed environment or genetic background. I argue that most, if not all, neutral mutations are of this sort, and that the essentialist notion of neutrality should be abandoned. This perspective reconciles two opposing views on the forces dominating organismal evolution, natural selection and random drift: neutral mutations occur and are especially abundant in robust systems, but they do not remain neutral indefinitely, and eventually become visible to natural selection, where some of them lead to evolutionary innovations.