Phenotypic plasticity and epigenetic marking: an assessment of evidence for genetic accommodation

The relationship between genotype (which is inherited) and phenotype (the target of selection) is mediated by environmental inputs on gene expression, trait development, and phenotypic integration. Phenotypic plasticity or epigenetic modification might influence evolution in two general ways: (1) by stimulating evolutionary responses to environmental change via population persistence or by revealing cryptic genetic variation to selection, and (2) through the process of genetic accommodation, whereby natural selection acts to improve the form, regulation, and phenotypic integration of novel phenotypic variants. We provide an overview of models and mechanisms for how such evolutionary influences may be manifested both for plasticity and epigenetic marking. We point to promising avenues of research, identifying systems that can best be used to address the role of plasticity in evolution, as well as the need to apply our expanding knowledge of genetic and epigenetic mechanisms to our understanding of how genetic accommodation occurs in nature. Our review of a wide variety of studies finds widespread evidence for evolution by genetic accommodation.

Developmental thermal plasticity among Drosophila melanogaster populations

Many biotic and abiotic variables influence the dispersal and distribution of organisms. Temperature has a major role in determining these patterns because it changes daily, seasonally and spatially, and these fluctuations have a significant impact on an organism's behaviour and fitness. Most ecologically relevant phenotypes that are adaptive are also complex and thus they are influenced by many underlying loci that interact with the environment. In this study, we quantified the degree of thermal phenotypic plasticity within and among populations by measuring chill-coma recovery times of lines reared from egg to adult at two different environmental temperatures. We used sixty genotypes from six natural populations of Drosophila melanogaster sampled along a latitudinal gradient in South America. We found significant variation in thermal plasticity both within and among populations. All populations exhibit a cold acclimation response, with flies reared at lower temperatures having increased resistance to cold. We tested a series of environmental parameters against the variation in population mean thermal plasticity and discovered the mean thermal plasticity was significantly correlated with altitude of origin of the population. Pairing our data with previous experiments on viability fitness assays in the same populations in fixed and variable environments suggests an adaptive role of this thermal plasticity in variable laboratory environments. Altogether, these data demonstrate abundant variation in adaptive thermal plasticity within and among populations.

Examining sources of variation in HPG axis function among individuals and populations of the dark-eyed junco

Gonadal steroids are important mediators of traits relevant to fitness, and thus may be targets of selection. However, more knowledge is needed about sources of variation along the endocrine axes that may contribute to functional variation in steroid levels. In a controlled captive environment, we studied males of two closely related subspecies of the dark-eyed junco (Junco hyemalis) that differ in testosterone-related phenotype, asking whether they also differ in testosterone (T), and assessing the contribution of the sequential links of the hypothalamic-pituitary-gonadal axis. When males of both subspecies were challenged with gonadotropin-releasing hormone (GnRH), they were similar in circulating luteinizing hormone (LH) and T responses. When challenged with exogenous LH, they again produced levels of T similar to one another, and to the levels produced in response to GnRH. However, the smaller, less ornamented, and less aggressive subspecies had greater abundance of mRNA for LH receptor in the testes and for androgen receptor in the rostral hypothalamus, suggesting potential differences in regulatory feedback. We suggest that circulating hormone levels may be less prone to evolutionary change than the responsiveness of individual hormone targets. Among individuals, T titers were highly repeatable whether males were challenged with GnRH or with LH, but LH produced in response to GnRH did not covary with T produced in response to LH. Testis mass, but not LH receptor transcript abundance, predicted individual variation in T responses. These data implicate the gonad, but not the pituitary, as an important source of individual variation in T production.

Genetic mapping of the female mimic morph locus in the ruff

BACKGROUND

Ruffs (Aves: Philomachus pugnax) possess a genetic polymorphism for male mating behaviour resulting in three permanent alternative male reproductive morphs: (i) territorial 'Independents', (ii) non-territorial 'Satellites', and (iii) female-mimicking 'Faeders'. Development into independent or satellite morphs has previously been shown to be due to a single-locus, two-allele autosomal Mendelian mode of inheritance at the Satellite locus. Here, we use linkage analysis to map the chromosomal location of the Faeder locus, which controls development into the Faeder morph, and draw further conclusions about candidate genes, assuming shared synteny with other birds.

RESULTS

Segregation data on the Faeder locus were obtained from captive-bred pedigrees comprising 64 multi-generation families (N = 381). There was no evidence that the Faeder locus was linked to the Satellite locus, but it was linked with microsatellite marker Ppu020. Comparative mapping of ruff microsatellite markers against the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes places the Ppu020 and Faeder loci on a region of chromosome 11 that includes the Melanocortin-1 receptor (MC1R) gene, which regulates colour polymorphisms in numerous birds and other vertebrates. Melanin-based colouration varies with life-history strategies in ruffs and other species, thus the MC1R gene is a strong candidate to play a role in alternative male morph determination.

CONCLUSION

Two unlinked loci appear to control behavioural development in ruffs. The Faeder locus is linked to Ppu020, which, assuming synteny, is located on avian chromosome 11. MC1R is a candidate gene involved in alternative male morph determination in ruffs.

Nature's inordinate fondness for metabolic enzymes: why metabolic enzyme loci are so frequently targets of selection

Metabolic enzyme loci were some of the first genes accessible for molecular evolution and ecology research. New technologies now make the whole genome, transcriptome or proteome readily accessible, allowing unbiased scans for loci exhibiting significant differences in allele frequency or expression level and associated with phenotypes and/or responses to natural selection. With surprising frequency and in many cases in proportions greater than chance relative to other genes, glycolysis and TCA cycle enzyme loci appear among the genes with significant associations in these studies. Hence, there is an ongoing need to understand the basis for fitness effects of metabolic enzyme polymorphisms. Allele-specific effects on the binding affinity and catalytic rate of individual enzymes are well known, but often of uncertain significance because metabolic control theory and in vivo studies indicate that many individual metabolic enzymes do not affect pathway flux rate. I review research, so far little used in evolutionary biology, showing that metabolic enzyme substrates affect signalling pathways that regulate cell and organismal biology, and that these enzymes have moonlighting functions. To date there is little knowledge of how alleles in natural populations affect these phenotypes. I discuss an example in which alleles of a TCA enzyme locus associate with differences in a signalling pathway and development, organismal performance, and ecological dynamics. Ultimately, understanding how metabolic enzyme polymorphisms map to phenotypes and fitness remains a compelling and ongoing need for gaining robust knowledge of ecological and evolutionary processes.

Functional genomics of adaptation to hypoxic cold-stress in high-altitude deer mice: transcriptomic plasticity and thermogenic performance

In species that are distributed across steep environmental gradients, adaptive variation in physiological performance may be attributable to transcriptional plasticity in underlying regulatory networks. Here we report the results of common-garden experiments that were designed to elucidate the role of regulatory plasticity in evolutionary adaptation to hypoxic cold-stress in deer mice (Peromyscus maniculatus). We integrated genomic transcriptional profiles with measures of metabolic enzyme activities and whole-animal thermogenic performance under hypoxia in highland (4350 m) and lowland (430 m) mice from three experimental groups: (1) wild-caught mice that were sampled at their native elevations; (2) wild-caught/lab-reared mice that were deacclimated to low-elevation conditions in a common-garden lab environment; and (3) the F(1) progeny of deacclimated mice that were maintained under the same low-elevation common-garden conditions. In each experimental group, highland mice exhibited greater thermogenic capacities than lowland mice, and this enhanced performance was associated with upregulation of transcriptional modules that influence several hierarchical steps in the O(2) cascade, including tissue O(2) diffusion (angiogenesis) and tissue O(2) utilization (metabolic fuel use and cellular oxidative capacity). Most of these performance-related transcriptomic changes occurred over physiological and developmental timescales, suggesting that regulatory plasticity makes important contributions to fitness-related physiological performance in highland deer mice.

Stress hormone receptors change as range expansion progresses in house sparrows

As ranges expand, individuals encounter different environments at the periphery than at the centre of the range. Previously, we have shown that glucocorticoids (GCs) vary with range expansion: individuals at the range edge release more GCs in response to restraint. Here, we measured hippocampal mRNA expression of GC receptors (mineralocorticoid, MR and glucocorticoid, GR) in eight house sparrow (Passer domesticus) populations varying in age. We found that individuals closest to the range edge had the lowest expression of MR relative to GR; in all likelihood, this relationship was driven by a marginal reduction of MR mRNA at the range edge. Reduced MR (relative to GR) might allow enhanced GC binding to GR, the lower affinity receptor that would enhance a rapid physiological and behavioural response to stressors. The insights gained from this study are not only enlightening to introduced species, but may also predict how certain species will react as their ranges shift owing to anthropogenic changes.

RNA-seq analysis reveals extensive transcriptional plasticity to temperature stress in a freshwater fish species

Background

Identifying genes of adaptive significance in a changing environment is a major focus of ecological genomics. Such efforts were restricted, until recently, to researchers studying a small group of model organisms or closely related taxa. With the advent of next generation sequencing (NGS), genomes and transcriptomes of virtually any species are now available for studies of adaptive evolution. We experimentally manipulated temperature conditions for two groups of crimson spotted rainbowfish (Melanotaenia duboulayi) and measured differences in RNA transcription between them. This non-migratory species is found across a latitudinal thermal gradient in eastern Australia and is predicted to be negatively impacted by ongoing environmental and climatic change.

Results

Using next generation RNA-seq technologies on an Illumina HiSeq2000 platform, we assembled a de novo transcriptome and tested for differential expression across the treatment groups. Quality of the assembly was high with a N50 length of 1856 bases. Of the 107,749 assembled contigs, we identified 4251 that were differentially expressed according to a consensus of four different mapping and significance testing approaches. Once duplicate isoforms were removed, we were able to annotate 614 up-regulated transfrags and 349 that showed reduced expression in the higher temperature group.

Conclusions

Annotated blast matches reveal that differentially expressed genes correspond to critical metabolic pathways previously shown to be important for temperature tolerance in other fish species. Our results indicate that rainbowfish exhibit predictable plastic regulatory responses to temperature stress and the genes we identified provide excellent candidates for further investigations of population adaptation to increasing temperatures.

Continental phylogeography of an ecologically and morphologically diverse Neotropical songbird, Zonotrichia capensis

BACKGROUND

The Neotropics are exceptionally diverse, containing roughly one third of all extant bird species on Earth. This remarkable species richness is thought to be a consequence of processes associated with both Andean orogenesis throughout the Tertiary, and climatic fluctuations during the Quaternary. Phylogeographic studies allow insights into how such events might have influenced evolutionary trajectories of species and ultimately contribute to a better understanding of speciation. Studies on continentally distributed species are of particular interest because different populations of such taxa may show genetic signatures of events that impacted the continent-wide biota. Here we evaluate the genealogical history of one of the world's most broadly-distributed and polytypic passerines, the rufous-collared sparrow (Zonotrichia capensis).

RESULTS

We obtained control region DNA sequences from 92 Zonotrichia capensis individuals sampled across the species' range (Central and South America). Six additional molecular markers, both nuclear and mitochondrial, were sequenced for a subset of individuals with divergent control region haplotypes. Median-joining network analysis, and Bayesian and maximum parsimony phylogenetic analyses all recovered three lineages: one spanning Middle America, the Dominican Republic, and north-western South America; one encompassing the Dominican Republic, Roraima (Venezuela) and La Paz (Bolivia) south to Tierra del Fuego, Argentina; and a third, including eastern Argentina and Brazil. Phylogenetic analyses suggest that the Middle American/north-western South American clade is sister to the remaining two. Bayesian and maximum likelihood coalescent simulations used to study lineage demographic history, diversification times, migration rates and population expansion together suggested that diversification of the three lineages occurred rapidly during the Pleistocene, with negligible gene flow, leaving genetic signatures of population expansions.

CONCLUSIONS

The Pleistocene history of the rufous-collared sparrow involved extensive range expansion from a probable Central American origin. Its remarkable morphological and behavioral diversity probably represents recent responses to local conditions overlying deeper patterns of lineage diversity, which are themselves produced by isolation and the history of colonization of South America.

RNA-seq transcriptome analysis of male and female zebra finch cell lines

The derivation of stably cultured cell lines has been critical to the advance of molecular biology. We profiled gene expression in the first two generally available cell lines derived from the zebra finch. Using Illumina RNA-seq, we generated ~93 million reads and mapped the majority to the recently assembled zebra finch genome. Expression of most Ensembl-annotated genes was detected, but over half of the mapped reads aligned outside annotated genes. The male-derived G266 line expressed Z-linked genes at a higher level than did the female-derived ZFTMA line, indicating persistence in culture of the distinctive lack of avian sex chromosome dosage compensation. Although these cell lines were not derived from neural tissue, many neurobiologically relevant genes were expressed, although typically at lower levels than in a reference sample from auditory forebrain. These cell lines recapitulate fundamental songbird biology and will be useful for future studies of songbird gene regulation and function.

Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and nonparallel physiological divergence

Across multiple lakes in North America, lake whitefish (Coregonus clupeaformis) have independently evolved 'dwarf' and 'normal' sympatric species pairs that exhibit pronounced phenotypic and genetic divergence. In particular, traits associated with metabolism have been shown to be highly differentiated between whitefish species. Here, we examine the transcription of genes associated with the five mitochondrial and nuclear genome-encoded oxidative phosphorylation (OXPHOS) complexes, the primary physiological mechanism responsible for the production of ATP, in whitefish species pairs from Cliff Lake and Webster Lake in Maine, USA. We observed OXPHOS gene transcription divergence between dwarf and normal whitefish in each of the two lakes, with the former exhibiting transcription upregulation for genes associated with each of the OXPHOS complexes. We also observed a significant influence of lake on transcription levels for some of the genes, indicating that inter-lake ecological or genetic differences are contributing to variation in OXPHOS gene transcription levels. Together, our results support the hypothesis that metabolic divergence is a critical adaptation involved in whitefish speciation and implicate OXPHOS gene upregulation as a factor involved in meeting the enhanced energetic demands of dwarf whitefish. Further studies are now needed to evaluate the contribution of genetically vs. plasticity driven variation in transcription associated with this critical physiological pathway.

Genomic insights into adaptation to high-altitude environments

Elucidating the molecular genetic basis of adaptive traits is a central goal of evolutionary genetics. The cold, hypoxic conditions of high-altitude habitats impose severe metabolic demands on endothermic vertebrates, and understanding how high-altitude endotherms cope with the combined effects of hypoxia and cold can provide important insights into the process of adaptive evolution. The physiological responses to high-altitude stress have been the subject of over a century of research, and recent advances in genomic technologies have opened up exciting opportunities to explore the molecular genetic basis of adaptive physiological traits. Here, we review recent literature on the use of genomic approaches to study adaptation to high-altitude hypoxia in terrestrial vertebrates, and explore opportunities provided by newly developed technologies to address unanswered questions in high-altitude adaptation at a genomic scale.

Comparative genomics in ecological physiology: toward a more nuanced understanding of acclimation and adaptation

Organisms that live in variable environments must adjust their physiology to compensate for environmental change. Modern functional genomics technologies offer global top-down discovery-based tools for identifying and exploring the mechanistic basis by which organisms respond physiologically to a detected change in the environment. Given that populations and species from different niches may exhibit different acclimation abilities, comparative genomic approaches may offer more nuanced understanding of acclimation responses, and provide insight into the mechanistic and genomic basis of variable acclimation. The physiological genomics literature is large and growing, as is the comparative evolutionary genomics literature. Yet, expansion of physiological genomics experiments to exploit taxonomic variation remains relatively undeveloped. Here, recent advances in the emerging field of comparative physiological genomics are considered, including examples of plants, bees and fish, and opportunities for further development are outlined particularly in the context of climate change research. Elements of robust experimental design are discussed with emphasis on the phylogenetic comparative approach. Understanding how acclimation ability is partitioned among populations and species in nature, and knowledge of the relevant genes and mechanisms, will be important for characterizing and predicting the ecological and evolutionary consequences of human-accelerated environmental change.

The effect of short-term hypoxic exposure on metabolic gene expression

The long-term effect of hypoxia is to decrease both the production and use of ATP and thus decrease the reliance on mitochondrial oxidative energy production. Yet, recent studies include more immediate affects of hypoxia on gene expression and these data suggest the maintenance of mitochondrial function. To better understand the short-term physiological response to hypoxia, we quantified metabolic mRNA expression in the heart ventricles and livers of the teleost fish Fundulus grandis exposed to partial oxygen pressure of 2.8 kPa (-13.5% air saturation).Twenty-eight individuals from a single population were exposed to hypoxia for 0, 4, 8, 12, 24, 48, and 96  hr. Liver and cardiac tissues were sampled from the same individuals at 0-48  hr. At 96  hr, only cardiac tissue was assayed. Gene expression was significantly different (ANOVA, P < 0.05) for 17 of 226 metabolic genes (7.5%) in cardiac tissue and for 20 of 256 (7.8%) metabolic genes in hepatic tissue. For the two tissues examined in this study, the maximum response occurred at different times. For cardiac tissue, using Dunnett's post hoc test, most of these significant differences occurred at 96  hr of exposure. For liver, all but one significant difference occurred at 4  hr. Surprisingly, too many (relative to random expectations) of the genes with significant increase in mRNA are involved in the oxidative phosphorylation pathway: 44% of the significant genes at 96  hr in the heart and 33% of the significant genes at 4  hr in the liver are involved in the oxidative phosphorylation pathway. These data indicate that there are tissue-specific differences in the timing of the response to hypoxia, yet both cardiac and hepatic tissues have increases in mRNA that code for enzyme in the oxidative phosphorylation pathway. If these changes in mRNA produce a similar change in protein, then these data suggest that the initial response to hypoxia involves an increase in the oxidative pathway potentially as a mechanism to maintain ATP production.

Transcriptomic resilience to global warming in the seagrass Zostera marina, a marine foundation species

Large-scale transcription profiling via direct cDNA sequencing provides important insights as to how foundation species cope with increasing climatic extremes predicted under global warming. Species distributed along a thermal cline, such as the ecologically important seagrass Zostera marina, provide an opportunity to assess temperature effects on gene expression as a function of their long-term adaptation to heat stress. We exposed a southern and northern European population of Zostera marina from contrasting thermal environments to a realistic heat wave in a common-stress garden. In a fully crossed experiment, eight cDNA libraries, each comprising ~125 000 reads, were obtained during and after a simulated heat wave, along with nonstressed control treatments. Although gene-expression patterns during stress were similar in both populations and were dominated by classical heat-shock proteins, transcription profiles diverged after the heat wave. Gene-expression patterns in southern genotypes returned to control values immediately, but genotypes from the northern site failed to recover and revealed the induction of genes involved in protein degradation, indicating failed metabolic compensation to high sea-surface temperature. We conclude that the return of gene-expression patterns during recovery provides critical information on thermal adaptation in aquatic habitats under climatic stress. As a unifying concept for ecological genomics, we propose transcriptomic resilience, analogous to ecological resilience, as an important measure to predict the tolerance of individuals and hence the fate of local populations in the face of global warming.

Elevated performance: the unique physiology of birds that fly at high altitudes

Birds that fly at high altitudes must support vigorous exercise in oxygen-thin environments. Here I discuss the characteristics that help high fliers sustain the high rates of metabolism needed for flight at elevation. Many traits in the O(2) transport pathway distinguish birds in general from other vertebrates. These include enhanced gas-exchange efficiency in the lungs, maintenance of O(2) delivery and oxygenation in the brain during hypoxia, augmented O(2) diffusion capacity in peripheral tissues and a high aerobic capacity. These traits are not high-altitude adaptations, because they are also characteristic of lowland birds, but are nonetheless important for hypoxia tolerance and exercise capacity. However, unique specializations also appear to have arisen, presumably by high-altitude adaptation, at every step in the O(2) pathway of highland species. The distinctive features of high fliers include an enhanced hypoxic ventilatory response, an effective breathing pattern, larger lungs, haemoglobin with a higher O(2) affinity, further augmentation of O(2) diffusion capacity in the periphery and multiple alterations in the metabolic properties of cardiac and skeletal muscle. These unique specializations improve the uptake, circulation and efficient utilization of O(2) during high-altitude hypoxia. High-altitude birds also have larger wings than their lowland relatives to reduce the metabolic costs of staying aloft in low-density air. High fliers are therefore unique in many ways, but the relative roles of adaptation and plasticity (acclimatization) in high-altitude flight are still unclear. Disentangling these roles will be instrumental if we are to understand the physiological basis of altitudinal range limits and how they might shift in response to climate change.

Rapid evolution of disease resistance is accompanied by functional changes in gene expression in a wild bird

Wild organisms are under increasing pressure to adapt rapidly to environmental changes. Predicting the impact of these changes on natural populations requires an understanding of the speed with which adaptive phenotypes can arise and spread, as well as of the underlying mechanisms. However, our understanding of these parameters is poor in natural populations. Here we use experimental and molecular approaches to investigate the recent emergence of resistance in eastern populations of North American house finches (Carpodacus mexicanus) to Mycoplasma galliseptum (MG), a severe conjunctivitis-causing bacterium. Two weeks following an experimental infection that took place in 2007, finches from eastern US populations with a 12-y history of exposure to MG harbored 33% lower MG loads in their conjunctivae than finches from western US populations with no prior exposure to MG. Using a cDNA microarray, we show that this phenotypic difference in resistance was associated with differences in splenic gene expression, with finches from the exposed populations up-regulating immune genes postinfection and those from the unexposed populations generally down-regulating them. The expression response of western US birds to experimental infection in 2007 was more similar to that of the eastern US birds studied in 2000, 7 y earlier in the epizootic, than to that of eastern birds in 2007. These results support the hypothesis that resistance has evolved by natural selection in the exposed populations over the 12 y of the epizootic. We hypothesize that host resistance arose and spread from standing genetic variation in the eastern US and highlight that natural selection can lead to rapid phenotypic evolution in populations when acting on such variation.

Role of DNA methylation of AHR1 and AHR2 promoters in differential sensitivity to PCBs in Atlantic Killifish, Fundulus heteroclitus

Atlantic killifish (Fundulus heteroclitus) inhabiting the PCB-contaminated Superfund site in New Bedford Harbor (MA, USA) have evolved genetic resistance to the toxic effects of these compounds. They also lack induction of cytochrome P4501A (CYP1A) and other aryl hydrocarbon receptor (AHR)-dependent responses after exposure to AHR agonists, suggesting an overall down-regulation of the AHR signaling pathway. In this study, we hypothesized that the genetic resistance is due to altered AHR expression resulting from hypermethylation of DNA in the promoter region of AHR genes in fish inhabiting New Bedford Harbor. To test this hypothesis, we cloned and sequenced AHR1 and AHR2 promoter regions and employed bisulfite conversion-polymerase chain reaction (BS-PCR) followed by clonal analysis to compare the methylation status of CpG islands of AHR1 and AHR2 in livers of adult killifish collected from New Bedford Harbor and a reference site (Scorton Creek, MA). No significant differences in methylation profiles were observed in either AHR1 or AHR2 promoter regions between NBH and SC fish. However, hypermethylation of the AHR1 promoter correlated with low expression of transcripts in the liver in both populations. In comparison to AHR1, hepatic mRNA expression of AHR2 is high and its promoter is hypomethylated. Taken together, our results suggest that genetic resistance to contaminants in NBH fish is not due to altered methylation of AHR promoter regions, but that promoter methylation may control tissue-specific expression of AHR genes in killifish.

Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates

High-altitude environments provide ideal testing grounds for investigations of mechanism and process in physiological adaptation. In vertebrates, much of our understanding of the acclimatization response to high-altitude hypoxia derives from studies of animal species that are native to lowland environments. Such studies can indicate whether phenotypic plasticity will generally facilitate or impede adaptation to high altitude. Here, we review general mechanisms of physiological acclimatization and genetic adaptation to high-altitude hypoxia in birds and mammals. We evaluate whether the acclimatization response to environmental hypoxia can be regarded generally as a mechanism of adaptive phenotypic plasticity, or whether it might sometimes represent a misdirected response that acts as a hindrance to genetic adaptation. In cases in which the acclimatization response to hypoxia is maladaptive, selection will favor an attenuation of the induced phenotypic change. This can result in a form of cryptic adaptive evolution in which phenotypic similarity between high- and low-altitude populations is attributable to directional selection on genetically based trait variation that offsets environmentally induced changes. The blunted erythropoietic and pulmonary vasoconstriction responses to hypoxia in Tibetan humans and numerous high-altitude birds and mammals provide possible examples of this phenomenon. When lowland animals colonize high-altitude environments, adaptive phenotypic plasticity can mitigate the costs of selection, thereby enhancing prospects for population establishment and persistence. By contrast, maladaptive plasticity has the opposite effect. Thus, insights into the acclimatization response of lowland animals to high-altitude hypoxia can provide a basis for predicting how altitudinal range limits might shift in response to climate change.

Ecological genomics--changing perspectives on Darwin's basic concerns

Ecological Genomics is an interdisciplinary field that seeks to understand the genetic and physiological basis of species interactions for evolutionary inferences. At the 7th annual Ecological Genomics Symposium, November 13-15, 2009, members of the Ecological Genomics program at Kansas State University invited 13 speakers and 56 poster presentations.

Molecular evolution of cytochrome C oxidase underlies high-altitude adaptation in the bar-headed goose

Bar-headed geese (Anser indicus) fly at up to 9,000 m elevation during their migration over the Himalayas, sustaining high metabolic rates in the severe hypoxia at these altitudes. We investigated the evolution of cardiac energy metabolism and O(2) transport in this species to better understand the molecular and physiological mechanisms of high-altitude adaptation. Compared with low-altitude geese (pink-footed geese and barnacle geese), bar-headed geese had larger lungs and higher capillary densities in the left ventricle of the heart, both of which should improve O(2) diffusion during hypoxia. Although myoglobin abundance and the activities of many metabolic enzymes (carnitine palmitoyltransferase, citrate synthase, 3-hydroxyacyl-coA dehydrogenase, lactate dehydrogenase, and pyruvate kinase) showed only minor variation between species, bar-headed geese had a striking alteration in the kinetics of cytochrome c oxidase (COX), the heteromeric enzyme that catalyzes O(2) reduction in oxidative phosphorylation. This was reflected by a lower maximum catalytic activity and a higher affinity for reduced cytochrome c. There were small differences between species in messenger RNA and protein expression of COX subunits 3 and 4, but these were inconsistent with the divergence in enzyme kinetics. However, the COX3 gene of bar-headed geese contained a nonsynonymous substitution at a site that is otherwise conserved across vertebrates and resulted in a major functional change of amino acid class (Trp-116 → Arg). This mutation was predicted by structural modeling to alter the interaction between COX3 and COX1. Adaptations in mitochondrial enzyme kinetics and O(2) transport capacity may therefore contribute to the exceptional ability of bar-headed geese to fly high.

Speciation genetics: current status and evolving approaches

The view of species as entities subjected to natural selection and amenable to change put forth by Charles Darwin and Alfred Wallace laid the conceptual foundation for understanding speciation. Initially marred by a rudimental understanding of hereditary principles, evolutionists gained appreciation of the mechanistic underpinnings of speciation following the merger of Mendelian genetic principles with Darwinian evolution. Only recently have we entered an era where deciphering the molecular basis of speciation is within reach. Much focus has been devoted to the genetic basis of intrinsic postzygotic isolation in model organisms and several hybrid incompatibility genes have been successfully identified. However, concomitant with the recent technological advancements in genome analysis and a newfound interest in the role of ecology in the differentiation process, speciation genetic research is becoming increasingly open to non-model organisms. This development will expand speciation research beyond the traditional boundaries and unveil the genetic basis of speciation from manifold perspectives and at various stages of the splitting process. This review aims at providing an extensive overview of speciation genetics. Starting from key historical developments and core concepts of speciation genetics, we focus much of our attention on evolving approaches and introduce promising methodological approaches for future research venues.

Nucleotide divergence vs. gene expression differentiation: comparative transcriptome sequencing in natural isolates from the carrion crow and its hybrid zone with the hooded crow

Recent advances in sequencing technology promise to provide new strategies for studying population differentiation and speciation phenomena in their earliest phases. We focus here on the black carrion crow (Corvus [corone] corone), which forms a zone of hybridization and overlap with the grey coated hooded crow (Corvus [corone] cornix). However, although these semispecies are taxonomically distinct, previous analyses based on several types of genetic markers did not reveal significant molecular differentiation between them. We here corroborate this result with sequence data obtained from a set of 25 nuclear intronic loci. Thus, the system represents a case of a very early phase of species divergence that requires new molecular approaches for its description. We have therefore generated RNAseq expression profiles using barcoded massively parallel pyrosequencing of brain mRNA from six individuals of the carrion crow and five individuals from a hybrid zone with the hooded crow. We obtained 856 675 reads from two runs, with average read length of 270 nt and coverage of 8.44. Reads were assembled de novo into 19 552 contigs, 70% of which could be assigned to annotated genes in chicken and zebra finch. This resulted in a total of 7637 orthologous genes and a core set of 1301 genes that could be compared across all individuals. We find a clear clustering of expression profiles for the pure carrion crow animals and disperse profiles for the animals from the hybrid zone. These results suggest that gene expression differences may indeed be a sensitive indicator of initial species divergence.

Evolving gene expression: from G to E to GxE

Analyses of gene expression data sets for multiple individuals and species promise to shed light on the mode of evolution of gene expression. However, complementary complexities challenge this goal. Characterization of the genetic variation underlying gene expression can easily be compromised by lack of environmental control. Conversely, the breadth of conclusions from studies of environmental effects has been limited by the use of single strains. Controlled studies have hinted at extensive genexenvironment interaction. Thus, both genetics and environment are key components in models of the evolution of gene expression. We review the literature on the evolution of gene expression in terms of genetics (G), environmental response (E) and GxE interactions to make this conceptual point.

Transcriptomic variation in a coral reveals pathways of clonal organisation

A microarray study was undertaken to examine the potential for clonal gene expression variation in a branching reef building coral, Acropora millepora. The role of small-scale gradients in light and water flow was examined by comparing gene expression levels between branch elevation (tip and base) and position (centre and edge) of replicate coral colonies (n=3). Analyses of variance revealed that almost 60% of variation in gene expression was present between colonies and 34 genes were considered differentially expressed between colonies (minimum P=6.5×10(-4)). These genes are associated with energy metabolism, protein biosynthesis and cell-cell recognition representing either genotypic variation in gene expression or the effects of specific environmental conditions that affect patterns of energy acquisition, growth and pathogen encounters. Less variation was present between central and peripheral branches (7%) and only a single gene was deemed differentially expressed (P=1.493×10(-3)). The function of this gene, a phosphatidylserine decarboxylase, suggests different growth patterns between branch positions within colonies and is consistent with the usual higher growth rates on the perimeter of corymbose-like branching coral colonies such as A. millepora. Four genes were differentially expressed between the tip and base of branches (P=3.239×10(-4)) and were associated with lysosome lipase activity and fluorescence, suggesting that branch tips may encounter higher pathogen loads or levels of mechanical stress and require greater levels of photo-protection associated with higher water flow and light levels. This study therefore confirms transcriptomic variation in response to small-scale environmental gradients consistent with differential resource allocation in clonal coral colonies.

Microarray analysis reveals transcriptional plasticity in the reef building coral Acropora millepora

We investigated variation in transcript abundance in the scleractinian coral, Acropora millepora, within and between populations characteristically exposed to different turbidity regimes and hence different levels of light and suspended particulate matter. We examined phenotypic plasticity by comparing levels of gene expression between source populations and following 10 days of acclimatization to a laboratory environment. Analyses of variance revealed that 0.05% of genes were differentially expressed between source populations, 1.32% following translocation into a common laboratory and 0.07% in the interaction (source population-dependent responses to translocation). Functional analyses identified an over-representation of differentially expressed genes associated with metabolism and fluorescence categories (primarily downregulated), and environmental information processing (primarily upregulated) following translocation to a lower light and turbidity environment. Such metabolic downregulation may indicate nonoxidative stress, hibernation or caloric restriction associated with the changed environmental conditions. Green fluorescent protein-related genes were the most differentially expressed and were exclusively downregulated; however, green fluorescent protein levels remained unchanged following translocation. Photophysiological responses of corals from both locations were characterized by a decline when introduced to the common laboratory environment but remained healthy (F(v)/F(m) > 0.6). Declines in total lipid content following translocation were the greatest for inshore corals, suggesting that turbid water corals have a strong reliance on heterotrophic feeding.

A microarray's view of life in the desert: adding a powerful evolutionary genomics tool to the packrat's midden

Identifying the genetic architecture of adaptive traits is fundamental to understanding how organisms respond to their environment, over both ecological and evolutionary timeframes. Microarray technology that allows us to capture the simultaneous expression of thousands of genes provides unparalleled insight into how organisms cope with their environment at the transcriptional level. Recent studies in Molecular Ecology demonstrate how microarrays can rapidly identify which genes and pathways allow organisms to face some of the most fundamental physiological challenges posed by the environment, including compensation for the hypoxic and thermal stress of high-altitudes (Cheviron et al. 2008) and, in this issue, the biotransformation of toxic plant secondary compounds by mammals (Magnanou et al. 2009). Microarrays (Ekins et al. 1989; Fodor et al. 1991) are glass slides affixed with hundreds to thousands of oligonucleotide or cDNA sequences (probes). Messenger RNA transcripts (typically reverse transcribed to cDNA) are isolated from a tissue/sample of interest and hybridized to the array. Binding to specific probes indicates that a particular gene was transcriptionally active at or near the time of sampling and thus provides a potentially comprehensive measure of gene expression. Although a tremendously powerful tool, commercially produced oligonucleotide arrays are only available for a handful of model organisms. Nonetheless, evolutionary ecologists have exploited this resource by using a cross-species hybridization approach (e.g. Saetre et al. 2004), that is, hybridizing a model organism array with a nonmodel sample (Bar-Or et al. 2007). Magnanou et al. (2009) present a novel example of using a model muroid microarray (Agilent Technologies, Rattus) to study physiological response in a wild, nonmodel muroid, Neotoma.

Migration-selection balance and local adaptation of mitochondrial haplotypes in rufous-collared sparrows (Zonotrichia capensis) along an elevational gradient

Variable selection pressures across heterogeneous landscapes can lead to local adaptation of populations. The extent of local adaptation depends on the interplay between natural selection and gene flow, but the nature of this relationship is complex. Gene flow can constrain local adaptation by eroding differentiation driven by natural selection, or local adaptation can itself constrain gene flow through selection against maladapted immigrants. Here we test for evidence that natural selection constrains gene flow among populations of a widespread passerine bird (Zonotrichia capensis) that are distributed along an elevational gradient in the Peruvian Andes. Using multilocus sequences and microsatellites screened in 142 individuals collected along a series of replicate transects, we found that mitochondrial gene flow was significantly reduced along elevational transects relative to latitudinal control transects. Nuclear gene flow, however, was not similarly reduced. Clines in mitochondrial haplotype frequency were strongly associated with transitions in environmental variables along the elevational transects, but this association was not observed for the nuclear markers. These results suggest that natural selection constrains mitochondrial gene flow along elevational gradients and that the mitonuclear discrepancy may be due to local adaptation of mitochondrial haplotypes.