A chromosome-level genome assembly of a free-living white-crowned sparrow (Zonotrichia leucophrys gambelii)

The white-crowned sparrow, Zonotrichia leucophrys, is a passerine bird with a wide distribution and it is extensively adapted to environmental changes. It has historically acted as a model species in studies on avian ecology, physiology and behaviour. Here, we present a high-quality chromosome-level genome of Zonotrichia leucophrys using PacBio and OmniC sequencing data. Gene models were constructed by combining RNA-seq and Iso-seq data from liver, hypothalamus, and ovary. In total a 1,123,996,003 bp genome was generated, including 31 chromosomes assembled in complete scaffolds along with other, unplaced scaffolds. This high-quality genome assembly offers an important genomic resource for the research community using the white-crowned sparrow as a model for understanding avian genome biology and development, and provides a genomic basis for future studies, both fundamental and applied.

Antibiotic-Induced Treatments Reveal Stress-Responsive Gene Expression in the Endangered Lichen Lobaria pulmonaria

Antibiotics are primarily found in the environment due to human activity, which has been reported to influence the structure of biotic communities and the ecological functions of soil and water ecosystems. Nonetheless, their effects in other terrestrial ecosystems have not been well studied. As a result of oxidative stress in organisms exposed to high levels of antibiotics, genotoxicity can lead to DNA damage and, potentially, cell death. In addition, in symbiotic organisms, removal of the associated microbiome by antibiotic treatment has been observed to have a big impact on the host, e.g., corals. The lung lichen Lobaria pulmonaria has more than 800 associated bacterial species, a microbiome which has been hypothesized to increase the lichen's fitness. We artificially exposed samples of L. pulmonaria to antibiotics and a stepwise temperature increase to determine the relative effects of antibiotic treatments vs. temperature on the mycobiont and photobiont gene expression and the viability and on the community structure of the lichen-associated bacteria. We found that the mycobiont and photobiont highly reacted to different antibiotics, independently of temperature exposure. We did not find major differences in bacterial community composition or alpha diversity between antibiotic treatments and controls. For these reasons, the upregulation of stress-related genes in antibiotic-treated samples could be caused by genotoxicity in L. pulmonaria and its photobiont caused by exposure to antibiotics, and the observed stress responses are reactions of the symbiotic partners to reduce damage to their cells. Our study is of great interest for the community of researchers studying symbiotic organisms as it represents one of the first steps to understanding gene expression in an endangered lichen in response to exposure to toxic environments, along with dynamics in its associated bacterial communities.

Comparative transcriptome reveals the thermal stress response differences between Heilongjiang population and Xinjiang population of Lota lota

Some cold-water fishes are particularly sensitive to the water temperature increasing caused by current global warming. However, the alterations in the physiology and behavior of infraspecific populations living in heterogeneous landscapes in response to water temperature increasing were significantly different. Consequently, understanding the impact of temperature increasing on different populations may be crucial for the conservation of cold-water fishes in the context of global warming. The burbot is the only freshwater specie in Gadiformes. To better understand the differences of different populations of burbot under similar thermal stress, Lota lota was selected as the research objects. Firstly, RNA-seq was applied to identify the transcriptomic responses of Heilongjiang population exposed to three temperature gradients (0 °C, 18 °C and 28 °C). Compared with 0 °C, 4216 and 12,657 genes were significantly differentially expressed at 18 °C and 28 °C, respectively. Meanwhile, 49 genes were significantly differentially expressed in three temperature pairs and these genes were presumed to involve in stress response process, immunologic process, reproductive process, development process, material metabolism process, signal transduction process, spermatogenesis process and cell apoptosis process. The response differences of two L. lota populations to similar thermal stress were compared and the results showed that they have different gene expression responses (the number of differentially expression genes and biological processes). The lower annual temperature of the Heilongjiang River might make it more sensitive to temperature increasing. Based on the comparative transcriptome analyses, 12 orthologous genes were considered as the potential regulators of L. lota preference for cold-water environment and these genes are potentially related to the immunologic process, reproductive process, development process, signal transduction process, and cell apoptosis process. Those results can provide basic information for the rational development of conservation strategies of different L. lota populations under the background of global warming.

Comparison between short-term stress and long-term adaptive responses reveal common paths to molecular adaptation

The phenotypic plasticity in responses to short-term stress can provide clues for understanding the adaptive fixation mechanism of genetic variation during long-term exposure to extreme environments. However, few studies have compared short-term stress responses with long-term evolutionary patterns; in particular, no interactions between the two processes have been evaluated in high-altitude environment. We performed RNA sequencing in embryo fibroblasts derived from great tits and mice to explore transcriptional responses after exposure to simulated high-altitude environmental stresses. Transcriptional changes of genes associated with metabolic pathways were identified in both bird and mice cells after short-term stress responses. Genomic comparisons among long-term highland tits and mammals and their lowland relatives revealed similar pathways (e.g., metabolic pathways) with that initiated under short-term stress transcriptional responses in vitro. These findings highlight the indicative roles of short-term stress in the long-term adaptation, and adopt common paths to molecular adaptation in mouse and bird cells.

•

Short-term stress and long-term adaptations share the common metabolic pathways

•

Phenotypic plasticity can promote adaptive evolution

•

Adopt common paths to molecular adaptation in mouse and bird cells

Gene expression responses to thermal shifts in the endangered lichen Lobaria pulmonaria

Anthropogenic climate change has led to unprecedented shifts in temperature across many ecosystems. In a context of rapid environmental changes, acclimation is an important process as it may influence the capacity of organisms to survive under novel thermal conditions. Mechanisms of acclimation could involve upregulation of stress response genes involved in protein folding, DNA damage repair and the regulation of signal transduction genes, along with a simultaneous downregulation of genes involved in growth or the cell cycle, in order to maintain cellular functions and equilibria. We transplanted Lobaria pulmonaria lichens originating from different forests to determine the relative effects of long-term acclimation and genetic factors on the variability in expression of mycobiont and photobiont genes. We found a strong response of the mycobiont and photobiont to high temperatures, regardless of sample origin. The green-algal photobiont had an overall lower response than the mycobiont. Gene expression of both symbionts was also influenced by acclimation to transplantation sites and by genetic factors. L. pulmonaria seems to have evolved powerful molecular pathways to deal with environmental fluctuations and stress and can acclimate to new habitats by transcriptomic convergence. Although L. pulmonaria has the molecular machinery to counteract short-term thermal stress, survival of lichens such as L. pulmonaria depends mostly on their long-term positive carbon balance, which can be compromised by higher temperatures and reduced precipitation, and both these outcomes have been predicted for Central Europe in connection with global climate change.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kākā (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kākā in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

Natural history collections are critical resources for contemporary and future studies of urban evolution

Urban environments are among the fastest changing habitats on the planet, and this change has evolutionary implications for the organisms inhabiting them. Herein, we demonstrate that natural history collections are critical resources for urban evolution studies. The specimens housed in these collections provide great potential for diverse types of urban evolution research, and strategic deposition of specimens and other materials from contemporary studies will determine the resources and research questions available to future urban evolutionary biologists. As natural history collections are windows into the past, they provide a crucial historical timescale for urban evolution research. While the importance of museum collections for research is generally appreciated, their utility in the study of urban evolution has not been explicitly evaluated. Here, we: (a) demonstrate that museum collections can greatly enhance urban evolution studies, (b) review patterns of specimen use and deposition in the urban evolution literature, (c) analyze how urban versus rural and native versus nonnative vertebrate species are being deposited in museum collections, and (d) make recommendations to researchers, museum professionals, scientific journal editors, funding agencies, permitting agencies, and professional societies to improve archiving policies. Our analyses of recent urban evolution studies reveal that museum specimens can be used for diverse research questions, but they are used infrequently. Further, although nearly all studies we analyzed generated resources that could be deposited in natural history collections (e.g., collected specimens), a minority (12%) of studies actually did so. Depositing such resources in collections is crucial to allow the scientific community to verify, replicate, and/or re-visit prior research. Therefore, to ensure that adequate museum resources are available for future urban evolutionary biology research, the research community-from practicing biologists to funding agencies and professional societies-must make adjustments that prioritize the collection and deposition of urban specimens.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H₂ S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H₂ S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H₂ S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Experimental study of hypoxia-induced changes in gene expression in an Asian pika, Ochotona dauurica

Acclimation to environmental changes driven by alterations in gene expression will serve as an important response for some species facing rapid Anthropogenic climate change. Pikas, genus Ochotona, are particularly vulnerable to climate change and current trends suggest that only the highest, coldest elevations within their ranges may remain suitable habitat for these species. In this study we aimed to assess the role of changes in gene expression in potentially facilitating elevational movements in pikas by measuring gene expression in the only known captive pika population, Ochotona dauurica, in response to hypoxic conditions. Using a controlled experiment, we exposed four male pikas to oxygen concentrations characteristic of sea-level, 2,000 m, and 4,000 m for 5 days each. Using blood samples collected after each treatment, we used RNAseq to determine if candidate pathways were undergoing significant changes in gene expression at different levels of oxygen (~100%, ~77%, and ~61% of sea-level oxygen concentrations). Gene set enrichment analyses showed that gene sets associated with the oxidative phosphorylation pathway and electron transport chain were significantly enriched for up-regulated genes in the 4,000 m samples compared to samples from the same individuals at lower-elevation conditions. Up-regulation of these pathways is consistent with known mechanisms of oxygen compensation. Our results suggest that these pikas have the acclimation capacity to tolerate oxygen concentrations characteristic of any elevation within their species range and that gene expression can be changed in a matter of days to accommodate drastically different oxygen concentrations. Thus, rapid and radical elevational movements that may be required of some pika species to avoid warmer temperatures in the Anthropocene will likely not be limited by hypoxic stress.

Transcriptome analysis indicates a broad range of toxic effects of Deepwater Horizon oil on Seaside Sparrows

In marine species, the transcriptomic response to Deepwater Horizon (DWH) oil implicated many biochemical pathways, with corresponding adverse outcomes on organ development and physiological performance. Terrestrial organisms differ in their mechanisms of exposure to polycyclic aromatic hydrocarbons (PAHs) and their physiological challenges, and may reveal either distinct effects of oil on biochemical pathways or the generality of the responses to oil shown in marine species. Using a cross-species hybridization microarray approach, we investigated the transcriptomic response in the liver of Seaside Sparrows (Ammospiza maritima) exposed to DWH oil compared with birds from a control site. Our analysis identified 295 genes differentially expressed between birds exposed to oil and controls. Gene ontology (GO) and canonical pathway analysis suggested that the identified genes were involved in a coordinated response that promoted hepatocellular proliferation and liver regeneration while inhibiting apoptosis, necrosis, and liver steatosis. Exposure to oil also altered the expression of genes regulating energy homeostasis, including carbohydrate metabolism and gluconeogenesis, and the biosynthesis, transport and metabolism of lipids. These results provide a molecular mechanism for the long-standing observation of hepatic hypertrophy and altered lipid biosynthesis and transport in birds exposed to crude oil. Several of the activated pathways and pathological outcomes shown here overlap with the ones altered in fish species upon exposure to oil. Overall, our study shows that the path of oil contamination from the marine system into salt marshes can lead to similar responses in terrestrial birds to those described in marine organisms, suggesting similar adverse outcomes and shared machinery for detoxification.

Phenotypic plasticity as a long-term memory easing readaptations to ancestral environments

Phenotypic plasticity refers to environment-induced phenotypic changes without mutation and is present in all organisms. The role of phenotypic plasticity in organismal adaptations to novel environments has attracted much attention, but its role in readaptations to ancestral environments is understudied. To address this question, we use the reciprocal transplant approach to investigate the multitissue transcriptomes of chickens adapted to the Tibetan Plateau and adjacent lowland. While many genetic transcriptomic changes had occurred in the forward adaptation to the highland, plastic changes largely transform the transcriptomes to the preferred state when Tibetan chickens are brought back to the lowland. The same trend holds for egg hatchability, a key component of the chicken fitness. These findings, along with highly similar patterns in comparable experiments of guppies and Escherichia coli, demonstrate that organisms generally "remember" their ancestral environments via phenotypic plasticity and reveal a mechanism by which past experience affects future evolution.

Transcriptome analysis reveals plasticity in gene regulation due to environmental cues in Primula sikkimensis, a high altitude plant species

BACKGROUND

Studying plasticity in gene expression in natural systems is crucial, for predicting and managing the effects of climate change on plant species. To understand the contribution of gene expression level variations to abiotic stress compensation in a Himalaya plant (Primula sikkimensis), we carried out a transplant experiment within (Ambient), and beyond (Below Ambient and Above Ambient) the altitudinal range limit of species. We sequenced nine transcriptomes (three each from each altitudinal range condition) using Illumina sequencing technology. We compared the fitness variation of transplants among three transplant conditions.

RESULTS

A large number of significantly differentially expressed genes (DEGs) between below ambient versus ambient (109) and above ambient versus ambient (85) were identified. Transcripts involved in plant growth and development were mostly up-regulated in below ambient conditions. Transcripts involved in signalling, defence, and membrane transport were mostly up-regulated in above ambient condition. Pathway analysis revealed that most of the genes involved in metabolic processes, secondary metabolism, and flavonoid biosynthesis were differentially expressed in below ambient conditions, whereas most of the genes involved in photosynthesis and plant hormone signalling were differentially expressed in above ambient conditions. In addition, we observed higher reproductive fitness in transplant individuals at below ambient condition compared to above ambient conditions; contrary to what we expect from the cold adaptive P. sikkimensis plants.

CONCLUSIONS

We reveal P. sikkimensis's capacity for rapid adaptation to climate change through transcriptome variation, which may facilitate the phenotypic plasticity observed in morphological and life history traits. The genes and pathways identified provide a genetic resource for understanding the temperature stress (both the hot and cold stress) tolerance mechanism of P. sikkimensis in their natural environment.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

Fine-tuning of seasonal timing of breeding is regulated downstream in the underlying neuro-endocrine system in a small songbird

The timing of breeding is under selection in wild populations as a result of climate change, and understanding the underlying physiological processes mediating this timing provides insight into the potential rate of adaptation. Current knowledge on this variation in physiology is, however, mostly limited to males. We assessed whether individual differences in the timing of breeding in females are reflected in differences in candidate gene expression and, if so, whether these differences occur in the upstream (hypothalamus) or downstream (ovary and liver) parts of the neuroendocrine system. We used 72 female great tits from two generations of lines artificially selected for early and late egg laying, which were housed in climate-controlled aviaries and went through two breeding cycles within 1 year. In the first breeding season we obtained individual egg-laying dates, while in the second breeding season, using the same individuals, we sampled several tissues at three time points based on the timing of the first breeding attempt. For each tissue, mRNA expression levels were measured using qPCR for a set of candidate genes associated with the timing of reproduction and subsequently analysed for differences between generations, time points and individual timing of breeding. We found differences in gene expression between generations in all tissues, with the most pronounced differences in the hypothalamus. Differences between time points, and early- and late-laying females, were found exclusively in the ovary and liver. Altogether, we show that fine-tuning of the seasonal timing of breeding, and thereby the opportunity for adaptation in the neuroendocrine system, is regulated mostly downstream in the neuro-endocrine system.

Transcriptomic Responses of Two Ecologically Divergent Populations of Japanese Mantis Shrimp (Oratosquilla oratoria) under Thermal Stress

Simple Summary

Rising ocean temperature would change the seawater chemistry and affect the external and internal physiology of crustaceans due to their lack of certain efficient temperature regulators. In addition, the infraspecific populations of crustaceans might also have different response strategies to the rising of temperature. Therefore, we identified the transcriptomic variations to the same thermal stress between ecologically divergent populations of Oratosquilla oratoria. The aim of this study was to investigate the population-specific function genes and relevant pathways in response to thermal stress in O. oratoria. The results showed that gene-expressed variation was in a population-specific pattern, which indicated that the local environment could lead to the evolvement of changes in gene regulation, ultimately leading to adaptive divergences. Additionally, we found several genes with large pleiotropic effects in the Zhoushan population, which might indicate that the regulation mechanisms of the Zhoushan population were more efficient than those of the Qingdao population under same thermal stress. The results provided some novel insights into the local adaptive differences of the infraspecific populations of O. oratoria and other crustaceans.

Abstract

Crustaceans are generally considered more sensitive to ocean warming due to their lack of certain efficient regulators. However, the alterations in the physiology and behavior of crustaceans in response to thermal stress differ vastly even among the infraspecific populations of heterogeneous landscapes. Consequently, understanding the impact of temperature fluctuation on crustacean infraspecific populations might be essential for maintaining a sustainable persistence of populations at existing locations. In the present study, we chose the Japanese mantis shrimp (Oratosquilla oratoria) as the representative crustacean population, and conducted transcriptome analyses in two divergent O. oratoria populations (the Zhoushan and Qingdao populations) under same thermal stress (20–28 °C) to identify the population-specific expression response to thermal stress. The results showed significant differences in gene expressions, GO terms and metabolic pathways between the two populations. We hypothesized that intraspecific mutations in the same or different genes might lead to thermal adaptive divergences. Temperature increases from 20–28 °C produced significant enrichment in GO terms and altered the metabolic pathways in the Zhoushan population despite the lack of differentially expressed unigenes. Therefore, several functional genes with large pleiotropic effects may underlie the response to thermal stress in the Zhoushan population. Furthermore, the most significantly enriched biological processes of the Qingdao population were associated with the state or activity of cells and its significant enriched pathways with genetic information processing as well as immune and environmental information processing. In contrast, the differentially regulated unigenes of the Zhoushan population were primarily involved in the regulatory cellular and transcription processes and the most significant pathways found were metabolic and digestive. Consequently, the regulatory mechanisms of the Zhoushan population are probably more efficient than those of the Qingdao population under the same thermal stress.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

New Zealand Tree and Giant Wētā (Orthoptera) Transcriptomics Reveal Divergent Selection Patterns in Metabolic Loci

Exposure to low temperatures requires an organism to overcome physiological challenges. New Zealand wētā belonging to the genera Hemideina and Deinacrida are found across a wide range of thermal environments and therefore subject to varying selective pressures. Here we assess the selection pressures across the wētā phylogeny, with a particular emphasis on identifying genes under positive or diversifying selection. We used RNA-seq to generate transcriptomes for all 18 Deinacrida and Hemideina species. A total of 755 orthologous genes were identified using a bidirectional best-hit approach, with the resulting gene set encompassing a diverse range of functional classes. Analysis of ortholog ratios of synonymous to nonsynonymous amino acid changes found 83 genes that are under positive selection for at least one codon. A wide variety of Gene Ontology terms, enzymes, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways are represented among these genes. In particular, enzymes involved in oxidative phosphorylation, melanin synthesis, and free-radical scavenging are represented, consistent with physiological and metabolic changes that are associated with adaptation to alpine environments. Structural alignment of the transcripts with the most codons under positive selection revealed that the majority of sites are surface residues, and therefore have the potential to influence the thermostability of the enzyme, with the exception of prophenoloxidase where two residues near the active site are under selection. These proteins provide interesting candidates for further analysis of protein evolution.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus)

Simple Summary

Micromys minutus is a small rodent species that has a wide range of vertical distribution in Taiwan. By comparing the gene expression profile of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site, the Tnfrsf12a gene was demonstrated to have a differential expression pattern. Although this finding may be correlated with the altitude acclimation, the observation of only one gene transcript with significant alteration leads us to suggest that genetic response may not play a major role in altitude acclimation in M. minutus. Future comparative functional genomics studies involving other organ systems (in addition to skeletal muscles) and alarger sample size are warranted for better insight into the altitude acclimation of this small rodent species.

Abstract

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < −4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Genomic Insights into the Adaptive Convergent Evolution

Adaptive convergent evolution, which refers to the same or similar phenotypes produced by species from independent lineages under similar selective pressures, has been widely examined for a long time. Accumulating studies on the adaptive convergent evolution have been reported from many different perspectives (cellular, anatomical, morphological, physiological, biochemical, and behavioral). Recent advances in the genomic technologies have demonstrated that adaptive convergence can arise from specific genetic mechanisms in different hierarchies, ranging from the same nucleotide or amino acid substitutions to the biological functions or pathways. Among these genetic mechanisms, the same amino acid changes in protein-coding genes play an important role in adaptive phenotypic convergence. Methods for detecting adaptive convergence at the protein sequence level have been constantly debated and developed. Here, we review recent progress on using genomic approaches to evaluate the genetic mechanisms of adaptive convergent evolution, summarize the research methods for identifying adaptive amino acid convergence, and discuss the future perspectives for researching adaptive convergent evolu-tion.

Gene expression is implicated in the ability of pikas to occupy Himalayan elevational gradient

Species are shifting their ranges due to climate change, many moving to cooler and higher locations. However, with elevation increase comes oxygen decline, potentially limiting a species’ ability to track its environment depending on what mechanisms it has available to compensate for hypoxic stress. Pikas (Family Ochotonidae), cold-specialist small mammal species, are already undergoing elevational range shifts. We collected RNA samples from one population of Ochotona roylei in the western Himalaya at three sites– 3,600, 4,000, and 5,000 meters–and found no evidence of significant population genetic structure nor positive selection among sites. However, out of over 10,000 expressed transcripts, 26 were significantly upregulated at the 5,000 m site and were significantly enriched for pathways consistent with physiological compensation for limited oxygen. These results suggest that differences in gene expression may play a key role in enabling hypoxia tolerance on this local scale, indicating elevational flexibility that may facilitate successful range shifts in response to climate change.

Chromosomal polymorphism and comparative chromosome painting in the rufous-collared sparrow (Zonotrichia capensis)

Zonotrichia capensis is widely distributed in the Neotropics. Previous cytogenetic studies demonstrated the presence of polymorphisms in two chromosome pairs (ZCA2 and ZCA4). Here, we report results based on comparative chromosome painting, using probes derived from Gallus gallus and Leucopternis albicollis, focused on characterizing the chromosome organization of Z. capensis. Our results demonstrate the conservation of ancestral syntenies as observed previously in other species of passerine. Syntenies were rearranged by a series of inversions in the second chromosome as described in other Passeriformes, but in this species, by using probes derived from L. albicollis we observed an extra inversion in the second chromosome that had not previously been reported. We also report a paracentric inversion in pair 3; this chromosome corresponds to the second chromosome in Zonotrichia albicollis and may indicate the presence of ancestral inversions in the genus. The chromosomal inversions we found might be important for understanding the phenotypic variation that exists throughout the distribution of Z. capensis.

Analysis of mitochondrial DNA sequence and copy number variation across five high-altitude species and their low-altitude relatives

High-altitude inhospitable environments impose a formidable life challenge for the local animals. Training and exposure to high-altitude environments produce both distinct physiological and phenotypic characteristics. The mitochondrion, an organelle crucial for the energy production, plays an important role in hypoxia adaptation. In this study, we investigated the mitochondrial DNA (mtDNA) polymorphism and copy number variation between the population pairs from distinct altitudes across the multi-species. Higher mitochondrial DNA control region's genetic diversity is conspicuous in high-altitude animals versus low-altitude relatives. We also found an accordant decrease of mtDNA copy number in most of the tissues from high-altitude animals. Compared to mammals, chickens have significantly distinct mitogenomic characteristics, and more significant changes in the skeletal muscle mtDNA copy number between high- and low-altitude individuals. Our study catches a snapshot of the biological similarities and differences in the mitochondrial high-altitude acclimation across the species.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Gene expression variations in high-altitude adaptation: a case study of the Asiatic toad (Bufo gargarizans)

Background

Genome-wide investigation of molecular mechanisms for high-altitude adaptation has attracted great attention in the last few years. In order to understand the contribution of gene expression level variations to high-altitude adaptation in Asiatic toads (Bufo gargarizans), we implemented a reciprocal transplant experiment between low- and high-altitude sites and sequenced 12 transcriptomes from brain, heart, and liver tissues.

Results

A large number of genes with expression differences (DEGs) between high- and low-altitude individuals (193 fixed and 844 plastic) were identified, and the majority of them were tissue specific. Heart displayed the largest number of DEGs, both plastic and fixed. Fixed DEGs were particularly concentrated in functions associated with muscle contraction, and the majority of them were down-regulated in high-altitude individuals. Plastic DEGs were highly concentrated in several energy metabolism related functional categories, and the majority of them were also down-regulated at high-altitude environments. In liver samples, genes associated with nutrient metabolism experienced a broad-scale expression down-regulation in high-altitude toads.

Conclusions

These broadly suppressed expression patterns at high altitudes are in strong contrast to those of endothermic homeotherms, suggesting poikilothermic vertebrates may have adopted different strategies at high altitudes. Our results strongly support that both genotypic specialization and phenotypic plasticity play crucial role in adaptation to high altitude for Asiatic toads. Poikilothermic vertebrates are among the most hypoxia-tolerant animals known, and many molecular mechanisms remain elusive. We hope that our results will provide useful directions for future research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0529-z) contains supplementary material, which is available to authorized users.

Evolutionary Transcriptomics and Proteomics: Insight into Plant Adaptation

Comparative transcriptomics and proteomics (T&P) have brought biological insight into development, gene function, and physiological stress responses. However, RNA-seq and high-throughput proteomics remain underutilised in studies of plant adaptation. These methodologies have created discovery tools with the potential to significantly advance our understanding of adaptive diversification. We outline experimental recommendations for their application. We discuss analysis models and approaches that accelerate the identification of adaptive gene sets and integrate transcriptome, proteome, phenotypic, and environmental data. Finally, we encourage widespread uptake and future developments in T&P that will advance our understanding of evolution and adaptation.

Population transcriptomes reveal synergistic responses of DNA polymorphism and RNA expression to extreme environments on the Qinghai-Tibetan Plateau in a predatory bird

Low oxygen and temperature pose key physiological challenges for endotherms living on the Qinghai-Tibetan Plateau (QTP). Molecular adaptations to high-altitude living have been detected in the genomes of Tibetans, their domesticated animals and a few wild species, but the contribution of transcriptional variation to altitudinal adaptation remains to be determined. Here we studied a top QTP predator, the saker falcon, and analysed how the transcriptome has become modified to cope with the stresses of hypoxia and hypothermia. Using a hierarchical design to study saker populations inhabiting grassland, steppe/desert and highland across Eurasia, we found that the QTP population is already distinct despite having colonized the Plateau <2000 years ago. Selection signals are limited at the cDNA level, but of only seventeen genes identified, three function in hypoxia and four in immune response. Our results show a significant role for RNA transcription: 50% of upregulated transcription factors were related to hypoxia responses, differentiated modules were significantly enriched for oxygen transport, and importantly, divergent EPAS1 functional variants with a refined co-expression network were identified. Conservative gene expression and relaxed immune gene variation may further reflect adaptation to hypothermia. Our results exemplify synergistic responses between DNA polymorphism and RNA expression diversity in coping with common stresses, underpinning the successful rapid colonization of a top predator onto the QTP. Importantly, molecular mechanisms underpinning highland adaptation involve relatively few genes, but are nonetheless more complex than previously thought and involve fine-tuned transcriptional responses and genomic adaptation.

Transcriptome analysis of a wild bird reveals physiological responses to the urban environment

Identifying the molecular basis of environmentally induced phenotypic variation presents exciting opportunities for furthering our understanding of how ecological processes and the environment can shape the phenotype. Urban and rural environments present free-living organisms with different challenges and opportunities, which have marked consequences for the phenotype, yet little is known about responses at the molecular level. We characterised transcriptomes from an urban and a rural population of great tits Parus major, demonstrating striking differences in gene expression profiles in both blood and liver tissues. Differentially expressed genes had functions related to immune and inflammatory responses, detoxification, protection against oxidative stress, lipid metabolism, and regulation of gene expression. Many genes linked to stress responses were expressed at higher levels in the urban birds, in accordance with our prediction that urban animals are exposed to greater environmental stress. This is one of the first studies to reveal transcriptional differences between urban- and rural-dwelling animals and suggests an important role for epigenetics in mediating environmentally induced physiological variation. The study provides valuable resources for developing further in-depth studies of the mechanisms driving phenotypic variation in the urban context at larger spatial and temporal scales.

Differential phenotypic and genetic expression of defence compounds in a plant-herbivore interaction along elevation

Elevation gradients impose large differences in abiotic and biotic conditions over short distances, in turn, likely driving differences in gene expression more than would genetic variation per se, as natural selection and drift are less likely to fix alleles at such a narrow spatial scale. As elevation increases, the pressure exerted on plants by herbivores and on arthropod herbivores by predators decreases, and organisms spanning the elevation gradient are thus expected to show lower levels of defence at high elevation. The alternative hypothesis, based on the optimal defence theory, is that defence allocation should be higher in low-resource habitats such as those at high elevation, due to higher costs associated with tissue replacement. In this study, we analyse variation with elevation in (i) defence compound content in the plant Lotus corniculatus and (ii) gene expression associated with defence against predators in the specific phytophagous moth, Zygaena filipendulae. Both species produce cyanogenic glycosides (CNglcs) such as lotaustralin and linamarin as defence mechanisms, with the moth, in addition, being able to sequester CNglcs from its host plant. Specifically, we tested the assumption that the defence-associated phenotype in plants and the gene expression in the insect herbivore should covary between low- and high-elevation environments. We found that L. corniculatus accumulated more CNglcs at high elevation, a result in agreement with the optimal defence theory. By contrast, we found that the levels of expression in the defence genes of Z. filipendulae larvae were not related to the CNglc content of their host plant. Overall, expression levels were not correlated with elevation either, with the exception of the UGT33A1 gene, which showed a marginally significant trend towards higher expression at high elevation when using a simple statistical framework. These results suggest that the defence phenotype of plants against herbivores, and subsequent herbivore sequestration machineries and de novo production, are based on a complex network of interactions.

Fungal and cyanobacterial gene expression in a lichen symbiosis: Effect of temperature and location

Organisms have evolved different cellular mechanisms to deal with environmental stress, primarily through complex molecular mechanisms including protein refolding and DNA repair. As mutualistic symbioses, lichens offer the possibility of analyzing molecular stress responses in a particularly tight interspecific relationship. We study the widespread cyanolichen Peltigera membranacea, a key player in carbon and nitrogen cycling in terrestrial ecosystems at northern latitudes. We ask whether increased temperature is reflected in mRNA levels of selected damage control genes, and do the response patterns show geographical associations? Using real-time PCR quantification of 38 transcripts, differential expression was demonstrated for nine cyanobacterial and nine fungal stress response genes (plus the fungal symbiosis-related lec2 gene) when the temperature was increased from 5 °C to 15 °C and 25 °C. Principle component analysis (PCA) revealed two gene groups with different response patterns. Whereas a set of cyanobacterial DNA repair genes and the fungal lec2 (PC1 group) showed an expression drop at 15 °C vs. 5 °C, most fungal candidates (PC2 group) showed increased expression at 25 °C vs. 5 °C. PC1 responses also correlated with elevation. The correlated downregulation of lec2 and cyanobacterial DNA repair genes suggests a possible interplay between the symbionts warranting further studies.

Divergence in gene expression within and between two closely related flycatcher species

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock‐like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between‐species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression P_ST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20‐kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated.

Transcriptome profiling of immune tissues reveals habitat-specific gene expression between lake and river sticklebacks

The observation of habitat-specific phenotypes suggests the action of natural selection. The three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to diverse freshwater habitats across the northern hemisphere since the last glaciation, while giving rise to recurring phenotypes associated with specific habitats. Parapatric lake and river populations of sticklebacks harbour distinct parasite communities, a factor proposed to contribute to adaptive differentiation between these ecotypes. However, little is known about the transcriptional response to the distinct parasite pressure of those fish in a natural setting. Here, we sampled wild-caught sticklebacks across four geographical locations from lake and river habitats differing in their parasite load. We compared gene expression profiles between lake and river populations using 77 whole-transcriptome libraries from two immune-relevant tissues, the head kidney and the spleen. Differential expression analyses revealed 139 genes with habitat-specific expression patterns across the sampled population pairs. Among the 139 differentially expressed genes, eight are annotated with an immune function and 42 have been identified as differentially expressed in previous experimental studies in which fish have been immune challenged. Together, these findings reinforce the hypothesis that parasites contribute to adaptation of sticklebacks in lake and river habitats.

Genetic responses to seasonal variation in altitudinal stress: whole-genome resequencing of great tit in eastern Himalayas

Species that undertake altitudinal migrations are exposed to a considerable seasonal variation in oxygen levels and temperature. How they cope with this was studied in a population of great tit (Parus major) that breeds at high elevations and winters at lower elevations in the eastern Himalayas. Comparison of population genomics of high altitudinal great tits and those living in lowlands revealed an accelerated genetic selection for carbohydrate energy metabolism (amino sugar, nucleotide sugar metabolism and insulin signaling pathways) and hypoxia response (PI3K-akt, mTOR and MAPK signaling pathways) in the high altitudinal population. The PI3K-akt, mTOR and MAPK pathways modulate the hypoxia-inducible factors, HIF-1α and VEGF protein expression thus indirectly regulate hypoxia induced angiogenesis, erythropoiesis and vasodilatation. The strategies observed in high altitudinal great tits differ from those described in a closely related species on the Tibetan Plateau, the sedentary ground tit (Parus humilis). This species has enhanced selection in lipid-specific metabolic pathways and hypoxia-inducible factor pathway (HIF-1). Comparative population genomics also revealed selection for larger body size in high altitudinal great tits.

Comparative transcriptomic analysis revealed adaptation mechanism of Phrynocephalus erythrurus, the highest altitude Lizard living in the Qinghai-Tibet Plateau

BACKGROUND

Organisms living at high altitudes must overcome three major environmental challenges: hypoxia, cold, and intense UV radiation. The molecular mechanisms that enable these challenges to be overcome have mainly been studied in endothermic organisms; relatively little attention has been paid to poikilothermic species. Here, we present deep transcriptome sequencing in two closely related lizards, the high altitude-dwelling Phrynocephalus erythrurus and the lowland-dwelling P. putjatia, to identify candidate genes under positive selection and to explore the convergent evolutionary adaptation of poikilothermic animals to high altitude life.

RESULTS

More than 70 million sequence reads were generated for each species via Illumina sequencing. De novo assembly produced 56,845 and 63,140 transcripts for P. erythrurus and P. putjatia, respectively. P. erythrurus had higher Ka/Ks ratios than P. putjatia, implying an accelerated evolutionary rate in the high altitude lizard lineage. 206 gene ontology (GO) categories with accelerated evolutionary rates and 43 candidate positively selected genes were detected along the P. erythrurus lineage. Some of these GO categories have functions associated with responses to hypoxia, energy metabolism and responses to UV damage. We also found that the high-altitude ranid frog R. kukunoris had higher Ka/Ks ratios than the closely related low-altitude frog R. chensinensis, and that the functional categories with accelerated evolutionary rates in R. kukunoris overlapped extensively with those detected along the P. erythrurus lineage.

CONCLUSIONS

The mechanisms of high altitude adaptation in P. erythrurus were tentatively inferred. By comparing two pairs of low- and high-altitude poikilothermic species, we found that similar functional categories had undergone positive selection in high altitude-dwelling Phrynocephalus and Rana lineages, indicating that similar mechanisms of adaptation to high altitude might have evolved in both genera. Our findings provide important guidance for future functional studies on high altitude adaptation in poikilothermic animals.

High-altitude ancestry and hypoxia acclimation have distinct effects on exercise capacity and muscle phenotype in deer mice

The hypoxic and cold environment at high altitudes requires that small mammals sustain high rates of O2 transport for exercise and thermogenesis while facing a diminished O2 availability. We used laboratory-born and -raised deer mice (Peromyscus maniculatus) from highland and lowland populations to determine the interactive effects of ancestry and hypoxia acclimation on exercise performance. Maximal O₂consumption (V̇o(2max)) during exercise in hypoxia increased after hypoxia acclimation (equivalent to the hypoxia at ∼4,300 m elevation for 6-8 wk) and was consistently greater in highlanders than in lowlanders. V̇o(2max) during exercise in normoxia was not affected by ancestry or acclimation. Highlanders also had consistently greater capillarity, oxidative fiber density, and maximal activities of oxidative enzymes (cytochrome c oxidase and citrate synthase) in the gastrocnemius muscle, lower lactate dehydrogenase activity in the gastrocnemius, and greater cytochrome c oxidase activity in the diaphragm. Hypoxia acclimation did not affect any of these muscle traits. The unique gastrocnemius phenotype of highlanders was associated with higher mRNA and protein abundances of peroxisome proliferator-activated receptor γ (PPARγ). Vascular endothelial growth factor (VEGFA) transcript abundance was lower in highlanders, and hypoxia acclimation reduced the expression of numerous genes that regulate angiogenesis and energy metabolism, in contrast to the observed population differences in muscle phenotype. Lowlanders exhibited greater increases in blood hemoglobin content, hematocrit, and wet lung mass (but not dry lung mass) than highlanders after hypoxia acclimation. Genotypic adaptation to high altitude, therefore, improves exercise performance in hypoxia by mechanisms that are at least partially distinct from those underlying hypoxia acclimation.

Integrating evolutionary and functional tests of adaptive hypotheses: a case study of altitudinal differentiation in hemoglobin function in an Andean Sparrow, Zonotrichia capensis

In air-breathing vertebrates, the physiologically optimal blood-O2 affinity is jointly determined by the prevailing partial pressure of atmospheric O2, the efficacy of pulmonary O2 transfer, and internal metabolic demands. Consequently, genetic variation in the oxygenation properties of hemoglobin (Hb) may be subject to spatially varying selection in species with broad elevational distributions. Here we report the results of a combined functional and evolutionary analysis of Hb polymorphism in the rufous-collared sparrow (Zonotrichia capensis), a species that is continuously distributed across a steep elevational gradient on the Pacific slope of the Peruvian Andes. We integrated a population genomic analysis that included all postnatally expressed Hb genes with functional studies of naturally occurring Hb variants, as well as recombinant Hb (rHb) mutants that were engineered through site-directed mutagenesis. We identified three clinally varying amino acid polymorphisms: Two in the α(A)-globin gene, which encodes the α-chain subunits of the major HbA isoform, and one in the α(D)-globin gene, which encodes the α-chain subunits of the minor HbD isoform. We then constructed and experimentally tested single- and double-mutant rHbs representing each of the alternative α(A)-globin genotypes that predominate at different elevations. Although the locus-specific patterns of altitudinal differentiation suggested a history of spatially varying selection acting on Hb polymorphism, the experimental tests demonstrated that the observed amino acid mutations have no discernible effect on respiratory properties of the HbA or HbD isoforms. These results highlight the importance of experimentally validating the hypothesized effects of genetic changes in protein function to avoid the pitfalls of adaptive storytelling.

Control of breathing and the circulation in high-altitude mammals and birds

Hypoxia is an unremitting stressor at high altitudes that places a premium on oxygen transport by the respiratory and cardiovascular systems. Phenotypic plasticity and genotypic adaptation at various steps in the O2 cascade could help offset the effects of hypoxia on cellular O2 supply in high-altitude natives. In this review, we will discuss the unique mechanisms by which ventilation, cardiac output, and blood flow are controlled in high-altitude mammals and birds. Acclimatization to high altitudes leads to some changes in respiratory and cardiovascular control that increase O2 transport in hypoxia (e.g., ventilatory acclimatization to hypoxia). However, acclimatization or development in hypoxia can also modify cardiorespiratory control in ways that are maladaptive for O2 transport. Hypoxia responses that arose as short-term solutions to O2 deprivation (e.g., peripheral vasoconstriction) or regional variation in O2 levels in the lungs (i.e., hypoxic pulmonary vasoconstriction) are detrimental at in chronic high-altitude hypoxia. Evolved changes in cardiorespiratory control have arisen in many high-altitude taxa, including increases in effective ventilation, attenuation of hypoxic pulmonary vasoconstriction, and changes in catecholamine sensitivity of the heart and systemic vasculature. Parallel evolution of some of these changes in independent highland lineages supports their adaptive significance. Much less is known about the genomic bases and potential interactive effects of adaptation, acclimatization, developmental plasticity, and trans-generational epigenetic transfer on cardiorespiratory control. Future work to understand these various influences on breathing and circulation in high-altitude natives will help elucidate how complex physiological systems can be pushed to their limits to maintain cellular function in hypoxia.

Reproducibility and consistency of proteomic experiments on natural populations of a non-model aquatic insect

Population proteomics has a great potential to address evolutionary and ecological questions, but its use in wild populations of non-model organisms is hampered by uncontrolled sources of variation. Here we compare the response to temperature extremes of two geographically distant populations of a diving beetle species (Agabus ramblae) using 2-D DIGE. After one week of acclimation in the laboratory under standard conditions, a third of the specimens of each population were placed at either 4 or 27°C for 12 h, with another third left as a control. We then compared the protein expression level of three replicated samples of 2-3 specimens for each treatment. Within each population, variation between replicated samples of the same treatment was always lower than variation between treatments, except for some control samples that retained a wider range of expression levels. The two populations had a similar response, without significant differences in the number of protein spots over- or under-expressed in the pairwise comparisons between treatments. We identified exemplary proteins among those differently expressed between treatments, which proved to be proteins known to be related to thermal response or stress. Overall, our results indicate that specimens collected in the wild are suitable for proteomic analyses, as the additional sources of variation were not enough to mask the consistency and reproducibility of the response to the temperature treatments.