The influence of population structure on gene expression and flowering time variation in the ubiquitous weed Capsella bursa‐pastoris (Brassicaceae)

Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa‐pastoris. We used RNA‐Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping‐by‐sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.     

Phenotypic plasticity in gene expression contributes to divergence of locally adapted populations of Fundulus heteroclitus

We examine the interaction between phenotypic plasticity and evolutionary adaptation using muscle gene expression levels among populations of the fish Fundulus heteroclitus acclimated to three temperatures. Our analysis reveals shared patterns of phenotypic plasticity due to thermal acclimation as well as non‐neutral patterns of variation among populations adapted to different thermal environments. For the majority of significant differences in gene expression levels, phenotypic plasticity and adaptation operate on different suites of genes. The subset of genes that demonstrate both adaptive differences and phenotypic plasticity, however, exhibit countergradient variation of expression. Thus, expression differences among populations counteract environmental effects, reducing the phenotypic differentiation between populations. Finally, gene‐by‐environment interactions among genes with non‐neutral patterns of expression suggest that the penetrance of adaptive variation depends on the environmental conditions experienced by the individual.     

Selection of reference genes for qRT‐PCR and expression analysis of high‐altitude‐related genes in grassland caterpillars (Lepidoptera: Erebidae: Gynaephora) along an altitude gradient

Changes in gene expression patterns can reflect the adaptation of organisms to divergent environments. Quantitative real‐time PCR (qRT‐PCR) is an important tool for ecological adaptation studies at the gene expression level. The quality of the results of qRT‐PCR analysis largely depends on the availability of reliable reference genes (RGs). To date, reliable RGs have not been determined for adaptive evolution studies in insects using a standard approach. Here, we evaluated the reliability of 17 candidate RGs for five Gynaephora populations inhabiting various altitudes of the Tibetan Plateau (TP) using four independent (geNorm, NormFinder, BestKeeper, and the deltaCt method) and one comprehensive (RefFinder) algorithms. Our results showed that EF1‐α, RPS15, and RPS13 were the top three most suitable RGs, and a combination of these three RGs was the most optimal for normalization. Conversely, RPS2, ACT, and RPL27 were the most unstable RGs. The expression profiles of two target genes (HSP70 and HSP90) were used to confirm the reliability of the chosen RGs. Additionally, the expression patterns of four other genes (GPI, HIF1A, HSP20, and USP) associated with adaptation to extreme environments were assessed to explore the adaptive mechanisms of TP Gynaephora species to divergent environments. Each of these six target genes showed discrepant expression patterns among the five populations, suggesting that the observed expression differences may be associated with the local adaptation of Gynaephora to divergent altitudinal environments. This study is a useful resource for studying the adaptive evolution of TP Gynaephora to divergent environments using qRT‐PCR, and it also acts as a guide for selecting suitable RGs for ecological and evolutionary studies in insects.     

COMPARATIVE FUNCTIONAL ANALYSES OF ULTRABITHORAX REVEAL MULTIPLE STEPS AND PATHS TO DIVERSIFICATION OF LEGS IN THE ADAPTIVE RADIATION OF SEMI‐AQUATIC INSECTS

Invasion of new ecological habitats is often associated with lineage diversification, yet the genetic changes underlying invasions and radiations are poorly understood. Over 200 million years ago, the semi‐aquatic insects invaded water surface from a common terrestrial ancestor and diversified to exploit a wide array of niches. Here, we uncover the changes in regulation and function of the gene Ultrabithorax associated with both the invasion of water surface and the subsequent diversification of the group. In the common ancestor of the semi‐aquatic insects, a novel deployment of Ubx protein in the mid‐legs increased their length, thereby enhancing their role in water surface walking. In derived lineages that specialize in rowing on the open water, additional changes in the timing of Ubx expression further elongated the mid‐legs thereby facilitating their function as oars. In addition, Ubx protein function was selectively reversed to shorten specific rear‐leg segments, thereby enabling their function as rudders. These changes in Ubx have generated distinct niche‐specialized morphologies that account for the remarkable diversification of the semi‐aquatic insects. Therefore, changes in the regulation and function of a key developmental gene may facilitate both the morphological change necessary to transition to novel habitats and fuel subsequent morphological diversification.     

EVALUATING PATTERNS OF CONVERGENT EVOLUTION AND TRANS‐SPECIES POLYMORPHISM AT MHC IMMUNOGENES IN TWO SYMPATRIC STICKLEBACK SPECIES

The immunologically important major histocompatibility complex (MHC) harbors some of the most polymorphic genes in vertebrates. These genes presumably evolve under parasite‐mediated selection and frequently show inconsistent allelic genealogies, where some alleles are more similar between species than within species. This phenomenon is thought to arise either from convergent evolution under parallel selection or from the preservation of ancient allelic lineages beyond speciation events (trans‐species polymorphism, TSP). Here, we examine natural populations of two sympatric stickleback species (Gasterosteus aculeatus and Pungitius pungitius) to investigate the contribution of these two mechanisms to the evolution of inconsistent allelic genealogies at the MHC. Overlapping parasite taxa between the two host species in three different habitats suggest contemporary parallel selection on the MHC genes. Accordingly, we detected a lack of species‐specific phylogenetic clustering in the immunologically relevant antigen‐binding residues of the MHC IIB genes which contrasted with the rest of the coding and noncoding sequence. However, clustering was not habitat‐specific and a codon‐usage analysis revealed patterns of similarity by descent. In this light, common descent via TSP, in combination with intraspecies gene conversion, rather than convergent evolution is the more strongly supported scenario for the inconsistent genealogy at the MHC.     

Long‐term acclimation to reciprocal light conditions suggests depth‐related selection in the marine foundation species Posidonia oceanica

Phenotypic differences among populations of the same species reflect selective responses to ecological gradients produced by variations in abiotic and biotic factors. Moreover, they can also originate from genetic differences among populations, due to a reduced gene flow. In this study, we examined the extent of differences in photo‐acclimative traits of Posidonia oceanica (L.) Delile clones collected above and below the summer thermocline (i.e., ?5 and ?25 m) in a continuous population extending along the water depth gradient. During a reciprocal light exposure and subsequent recovery in mesocosms, we assessed degree of phenotypic plasticity and local adaptation of plants collected at different depths, by measuring changes in several traits, such as gene expression of target genes, photo‐physiological features, and other fitness‐related traits (i.e., plant morphology, growth, and mortality rates). Samples were also genotyped, using microsatellite markers, in order to evaluate the genetic divergence among plants of the two depths. Measures collected during the study have shown a various degree of phenotypic changes among traits and experimental groups, the amount of phenotypic changes observed was also dependent on the type of light environments considered. Overall plants collected at different depths seem to be able to acclimate to reciprocal light conditions in the experimental time frame, through morphological changes and phenotypic buffering, supported by the plastic regulation of a reduced number of genes. Multivariate analyses indicated that plants cluster better on the base of their depth origin rather than the experimental light conditions applied. The two groups were genetically distinct, but the patterns of phenotypic divergence observed during the experiment support the hypothesis that ecological selection can play a role in the adaptive divergence of P. oceanica clones along the depth gradient.     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

Adaptive brain size divergence in nine‐spined sticklebacks (Pungitius pungitius)?

Most studies seeking to provide evolutionary explanations for brain size variability have relied on interspecific comparisons, while intraspecific studies utilizing ecologically divergent populations to this effect are rare. We investigated the brain size and structure of first‐generation laboratory‐bred nine‐spined sticklebacks (Pungitius pungitius) from four geographically and genetically isolated populations originating from markedly different habitats. We found that the relative size of bulbus olfactorius and telencephalon was significantly larger in marine than in pond populations. Significant, but habitat‐independent population differences were also found in relative brain and cerebellum sizes. The consistent, habitat‐specific differences in the relative size of bulbus olfactorius and telencephalon suggest their adaptive reduction in response to reduced (biotic and abiotic) habitat complexity in pond environments. In general, the results suggest that genetically based brain size and structure differences can evolve relatively rapidly and in repeatable fashion with respect to habitat structure.     

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species‐wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well‐watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA‐seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species‐wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population‐specific responses. Weighted gene co‐expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.     

ASYMMETRIC RESPONSE PATTERNS TO SUBSPECIES‐SPECIFIC SONG DIFFERENCES IN ALLOPATRY AND PARAPATRY IN THE GRAY‐BREASTED WOOD‐WREN

Song divergence between closely related taxa may play a critical role in the evolutionary processes of speciation and hybridization. We explored song variation between two Ecuadorian subspecies of the gray‐breasted wood‐wren (Henicorhina leucophrys) and tested the impact of song divergence on response behaviors. Songs were significantly different between the two subspecies, even between two parapatric populations 10 km apart. Playback experiments revealed an asymmetric response pattern to these divergent subspecies specific songs; one subspecies responded more to songs of its own subspecies than to the other subspecies’ songs, whereas the second responded equally strongly to songs of both subspecies. While song parameters revealed a mixed pattern of divergence between allopatric and parapatric populations, the majority of spectral characteristics showed increased divergence in parapatry, suggestive of character displacement. This increased song divergence in parapatry appeared to affect behavioral responses to playback as discriminating responses were most prominent in parapatry and against parapatric songs. The clear behavioral impact of subspecies‐specific song differences supports a potential role for song as an acoustic barrier to gene flow. The asymmetric nature of the responses suggests that song divergence could affect the direction of gene flow and the position of the subspecies‐specific transition.