The role of gene expression in ecological speciation

The Role of Hybridization in Species Formation and Persistence

Hybridization, or interbreeding between different taxa, was traditionally considered to be rare and to have a largely detrimental impact on biodiversity, sometimes leading to the breakdown of reproductive isolation and even to the reversal of speciation. However, modern genomic and analytical methods have shown that hybridization is common in some of the most diverse clades across the tree of life, sometimes leading to rapid increase of phenotypic variability, to introgression of adaptive alleles, to the formation of hybrid species, and even to entire species radiations. In this review, we identify consensus among diverse research programs to show how the field has progressed. Hybridization is a multifaceted evolutionary process that can strongly influence species formation and facilitate adaptation and persistence of species in a rapidly changing world. Progress on testing this hypothesis will require cooperation among different subdisciplines.

Gene Expression Shifts in Emperor Penguin Adaptation to the Extreme Antarctic Environment

Gene expression can accelerate ecological divergence by rapidly tweaking the response of an organism to novel environments, with more divergent environments exerting stronger selection and supposedly, requiring faster adaptive responses. Organisms adapted to extreme environments provide ideal systems to test this hypothesis, particularly when compared to related species with milder ecological niches. The Emperor penguin (Aptenodytes forsteri) is the only endothermic vertebrate breeding in the harsh Antarctic winter, in stark contrast with the less cold-adapted sister species, the King penguin (A. patagonicus). Assembling the first de novo transcriptomes and analysing multi-tissue (brain, kidney, liver, muscle, skin) RNA-Seq data from natural populations of both species, we quantified the shifts in tissue-enhanced genes, co-expression gene networks, and differentially expressed genes characterising Emperor penguin adaptation to the extreme Antarctic. Our analyses revealed the crucial role played by muscle and liver in temperature homeostasis, fasting, and whole-body energy metabolism (glucose/insulin regulation, lipid metabolism, fatty acid beta-oxidation, and blood coagulation). Repatterning at the regulatory level appears as more important in the brain of the Emperor penguin, showing the lowest signature of differential gene expression, but the largest co-expression gene network shift. Nevertheless, over-expressed genes related to mTOR signalling in the brain and the liver support their central role in cold and fasting responses. Besides contributing to understanding the genetics underlying complex traits, like body energy reservoir management, our results provide a first insight into the role of gene expression in adaptation to one of the most extreme environmental conditions endured by an endotherm.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Comparative transcriptomics and gene expression divergence associated with homoploid hybrid speciation in Argyranthemum

Ecological isolation is increasingly thought to play an important role in speciation, especially for the origin and reproductive isolation of homoploid hybrid species. However, the extent to which divergent and/or transgressive gene expression changes are involved in speciation is not well studied. In this study, we employ comparative transcriptomics to investigate gene expression changes associated with the origin and evolution of two homoploid hybrid plant species, Argyranthemum sundingii and A. lemsii (Asteraceae). As there is no standard methodology for comparative transcriptomics, we examined five different pipelines for data assembly and analysing gene expression across the four species (two hybrid and two parental). We note biases and problems with all pipelines, and the approach used affected the biological interpretation of the data. Using the approach that we found to be optimal, we identify transcripts showing DE between the parental taxa and between the homoploid hybrid species and their parents; in several cases, putative functions of these DE transcripts have a plausible role in ecological adaptation and could be the cause or consequence of ecological speciation. Although independently derived, the homoploid hybrid species have converged on similar expression phenotypes, likely due to adaptation to similar habitats.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

A first annotated genome sequence for Haliotis midae with genomic insights into abalone evolution and traits of economic importance

Haliotis midae or "perlemoen" is one of five abalone species endemic to South Africa, and being palatable, the only commercially important abalone species with a high international demand. The higher demand for this abalone species has resulted in the decrease of natural stocks due to overexploitation by capture fisheries and poaching. Facilitating aquaculture production of H. midae should assist in minimising the pressure on the wild populations. Here, the draft genome of H. midae has been sequenced, assembled, and annotated. The draft assembly resulted in a total length of 1.5 Gb, contig N50 of 0.238 Mb, scaffold N50 of 0. 238 Mb and GC level of 40%. Gene annotation, combining ab initio and evidence-based pipelines identified 52,280 genes with protein coding potential. The genes identified were used to predict orthologous genes shared among the four other abalone species (H. laevigata, H. rubra, H. discus hannai and H. rufescens) and 4702 orthologous genes were shared across the five species. Among the orthologous genes in abalones, single copy genes were further analysed for signatures of selection and several molecular regulatory proteins involved in developmental functions were found to be under positive selection in specific abalone lineages. Furthermore, whole genome SNP-based phylogenomic assessment was performed to confirm the evolutionary relationship among the considered abalone species with draft genomes, reaffirming that H. midae is closely related to the Australian Greenlip (H. laevigata) and Blacklip (H. rubra). The study assists in the understanding of genes related to various biological systems underscoring the evolution and development of abalones, with potential applications for genetic improvement of commercial stocks.

Chromosome-level reference genome of Tetrastigma hemsleyanum (Vitaceae) provides insights into genomic evolution and the biosynthesis of phenylpropanoids and flavonoids

Here, we present a high-quality chromosome-scale genome assembly (2.19 Gb) and annotation of Tetrastigma hemsleyanum, a perennial herbaceous liana native to subtropical China with diverse medicinal applications. Approximately 73% of the genome was comprised of transposable elements (TEs), of which long terminal repeat retrotransposons (LTR-RTs) were a predominant group (69% of the genome). The genome size increase of T. hemsleyanum (relative to Vitis species) was mostly due to the proliferation of LTR-RTs. Of the different modes of gene duplication identified, transposed duplication (TRD) and dispersed duplication (DSD) were the predominant ones. Genes, particularly those involved in the phenylpropanoid-flavonoid (PF) pathway and those associated with therapeutic properties and environmental stress resistance, were significantly amplified through recent tandem duplications. We dated the divergence of two intraspecific lineages in Southwest (SW) versus Central-South-East (CSE) China to the late Miocene (approximately 5.2 million years ago). Of those, the former showed more upregulated genes and metabolites. Based on resequencing data of 38 individuals representing both lineages, we identified various candidate genes related to 'response to stimulus' and 'biosynthetic process', including ThFLS11, which is putatively involved in flavonoid accumulation. Overall, this study provides abundant genomic resources for future evolutionary, ecological, and functional genomics studies in T. hemsleyanum and related species.

The lack of genetic variation underlying thermal transcriptomic plasticity suggests limited adaptability of the Northern shrimp, Pandalus borealis

Introduction

Genetic variation underlies the populations’ potential to adapt to and persist in a changing environment, while phenotypic plasticity can play a key role in buffering the negative impacts of such change at the individual level.

Methods

We investigated the role of genetic variation in the thermal response of the northern shrimp Pandalus borealis, an ectotherm species distributed in the Arctic and North Atlantic Oceans. More specifically, we estimated the proportion transcriptomic responses explained by genetic variance of female shrimp from three origins after 30 days of exposure to three temperature treatments.

Results

We characterized the P. borealis transcriptome (170,377 transcripts, of which 27.48% were functionally annotated) and then detected a total of 1,607 and 907 differentially expressed transcripts between temperatures and origins, respectively. Shrimp from different origins displayed high but similar level of transcriptomic plasticity in response to elevated temperatures. Differences in transcript expression among origins were not correlated to population genetic differentiation or diversity but to environmental conditions at origin during sampling.

Discussion

The lack of genetic variation explaining thermal plasticity suggests limited adaptability in this species’ response to future environmental changes. These results together with higher mortality observed at the highest temperature indicate that the thermal niche of P. borealis will likely be restricted to higher latitudes in the future. This prediction concurs with current decreases in abundance observed at the southern edge of this species geographical distribution, as it is for other cold-adapted crustaceans.

Environment and genotype predict the genomic nature of domestication of salmonids as revealed by gene expression

Billions of salmonids are produced annually by artificial reproduction for harvest and conservation. Morphologically, behaviourally and physiologically these fish differ from wild-born fish, including in ways consistent with domestication. Unlike most studied domesticates, which diverged from wild ancestors millennia ago, salmonids offer a tractable model for early-stage domestication. Here, we review a fundamental mechanism for domestication-driven differences in early-stage domestication, differentially expressed genes (DEGs), in salmonids. We found 34 publications examining DEGs under domestication driven by environment and genotype, covering six species, over a range of life-history stages and tissues. Three trends emerged. First, domesticated genotypes have increased expression of growth hormone and related metabolic genes, with differences magnified under artificial environments with increased food. Regulatory consequences of these DEGs potentially drive overall DEG patterns. Second, immune genes are often DEGs under domestication and not simply owing to release from growth-immune trade-offs under increased food. Third, domesticated genotypes exhibit reduced gene expression plasticity, with plasticity further reduced in low-complexity environments typical of production systems. Recommendations for experimental design improvements, coupled with tissue-specific expression and emerging analytical approaches for DEGs present tractable avenues to understand the evolution of domestication in salmonids and other species.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Extensive Adaptive Variation in Gene Expression within and between Closely Related Horseshoe Bats (Chiroptera, Rhinolophus) Revealed by Three Organs

In the process of species differentiation and adaption, the relative influence of natural selection on gene expression variation often remains unclear (especially its impact on phenotypic divergence). In this study, we used differentially expressed genes from brain, cochlea, and liver samples collected from two species of bats to determine the gene expression variation forced by natural selection when comparing at the interspecific (Rhinolophus siamensis and R. episcopus episcopus) and the intraspecific (R. e. episcopus and R. episcopus spp.) levels. In both cases, gene expression variation was extensively adaptive (>66.0%) and mainly governed by directional selection, followed by stabilizing selection, and finally balancing selection. The expression variation related to acoustic signals (resting frequency, RF) and body size (forearm length, FA) was also widely governed by natural selection (>69.1%). Different functional patterns of RF- or FA-related adaptive expression variation were found between the two comparisons, which manifested as abundant immune-related regulations between subspecies (indicating a relationship between immune response and phenotypic adaption). Our study verifies the extensive adaptive expression variation between both species and subspecies and provides insight into the effects of natural selection on species differentiation and adaptation as well as phenotypic divergence at the expression level.

Ecological significance of G protein-coupled receptors in the Pacific oyster (Crassostrea gigas): Pervasive gene duplication and distinct transcriptional response to marine environmental stresses

Marine ecosystems with ocean warming and industry pollution threaten the survival and adaptation of organisms. G protein-coupled receptors (GPCRs) play critical roles in various physiological and toxicological processes in vertebrates and invertebrates. The Pacific oyster (Crassostrea gigas) was widely used to study the adaptation of marine molluscs to coastal environments. In this work, we identified a total of 586 GPCRs in C. gigas genome. The C. gigas GPCRs were divided into five classes (including class A, B, C, E and F) with different degrees of expansion. Meta-analysis of multiple RNA-seq datasets revealed that transcriptional expression patterns of GPCRs in C. gigas were distinct in response to high temperature, salinity, air exposure, heavy metal, ostreid herpes virus 1 (OsHV-1) and Vibrio challenge. This work for the first time characterized the GPCR gene family and provided insights into the potential roles of GPCRs in adaptation of marine molluscs to stressful coastal environment.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

How to Accurately Delineate Morphologically Conserved Taxa and Diagnose Their Phenotypic Disparities: Species Delimitation in Cryptic Rhinolophidae (Chiroptera)

Systematics and taxonomy are the backbone of all components of biology and ecology, yet cryptic species present a major challenge for accurate species identification. This is especially problematic as they represent a substantial portion of undiscovered biodiversity, and have implications for not only species conservation, but even assaying potential risk of zoonotic spillover. Here, we use integrative approaches to delineate potential cryptic species in horseshoe bats (Rhinolophidae), evaluate the phenotypic disparities between cryptic species, and identify key traits for their identification. We tested the use of multispecies coalescent models (MSC) using Bayesian Phylogenetic and Phylogeography (BPP) and found that BPP was useful in delineating potential cryptic species, and consistent with acoustic traits. Our results show that around 40% of Asian rhinolophid species are potentially cryptic and have not been formally described. In order to avoid potential misidentification and allow species to be accurately identified, we identified quantitative noseleaf sella and acoustic characters as the most informative traits in delineating between potential cryptic species in Rhinolophidae. This highlights the physical differences between cryptic species that are apparent in noseleaf traits which often only qualitatively described but rarely measured. Each part of the noseleaf including the sella, lateral lappets, and lancet furrows, play roles in focusing acoustic beams and thus, provide useful characteristics to identify cryptic Rhinolophus species. Finally, species delimitation for cryptic species cannot rely on genetic data alone, but such data should be complemented by other evidence, including phenotypic, acoustic data, and geographic distributions to ensure accurate species identification and delineation.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

Mitochondrial volume density and evidence for its role in adaptive divergence in response to thermal tolerance in threespine stickleback

Phenotypic plasticity is predicted to permit persistence in new environments, and may subsequently evolve to enhance fitness. Colonizing environments with lower winter temperatures can lead to the evolution of lower critical thermal minima; the corresponding physiological traits associated with temperature tolerance are predicted to involve mitochondrial function. Threespine stickleback (Gasterosteus aculeatus) have colonized freshwater lakes along the Pacific Northwest. These freshwater populations are known to exhibit cold-induced increases in mitochondrial volume density in pectoral muscle, but whether such plasticity evolved before or after colonization is uncertain. Here, we measure critical thermal minima (CTmin) in one marine and one freshwater population of threespine stickleback, and mitochondrial volume density in pectoral and cardiac tissue of both populations acclimated to different temperature treatments (6.2, 14.5 and 20.6 ℃). Mitochondrial volume density increased with cold acclimation in pectoral muscle; cardiac muscle was non-plastic but had elevated mitochondrial volume densities compared to pectoral muscle across all temperature treatments. There were no differences in the levels of plasticity between marine and freshwater stickleback, but neither were there differences in CTmin. Importantly, marine stickleback exhibited plasticity under low-salinity conditions, suggesting that marine stickleback had at least one necessary phenotype for persistence in freshwater environments before colonization occurred.

Comparative analysis of the transcriptomes of two rice subspecies during domestication

Two subspecies of rice, Oryza sativa ssp. indica and O. sativa ssp. japonica, with reproductive isolation and differences in morphology and phenotypic differences, were established during the process of rice domestication. To understand how domestication has changed the transcriptomes of the two rice subspecies and given rise to the phenotypic differences, we obtained approximately 700 Gb RNA-Seq data from 26 indica and 25 japonica accessions, and identified 97,005 transcribed fragments and 4579 novel transcriptionally active regions. The two rice subspecies had significantly different gene expression profiles, we identified 1,357 (3.3% in all genes) differentially expressed genes (DEGs) between indica and japonica rice. Combining existing gene function studies, it is found that some of these differential genes are related to the differentiation of the two subspecies, such as grain shape and cold tolerance, etc. Functional annotation of these DEGs indicates that they are involved in cell wall biosynthesis and reproductive processes. Furthermore, compared with the non-DEGs, the DEGs from both subspecies had more 5'flanking regions with low polymorphism to divergence ratios, indicating a stronger positive selection pressure on the regulation of the DEGs. This study improves our understanding of the rice genome by comparatively analyzing the transcriptomes of indica and japonica rice and identifies DEGs those may be responsible for the reproductive isolation and phenotypic differences between the two rice subspecies.

Using Reciprocal Transplants to Assess Local Adaptation, Genetic Rescue, and Sexual Selection in Newly Established Populations

Small populations establishing on colonization fronts have to adapt to novel environments with limited genetic variation. The pace at which they can adapt, and the influence of genetic variation on their success, are key questions for understanding intraspecific diversity. To investigate these topics, we performed a reciprocal transplant experiment between two recently founded populations of brown trout in the sub-Antarctic Kerguelen Islands. Using individual tagging and genetic assignment methods, we tracked the fitness of local and foreign individuals, as well as the fitness of their offspring over two generations. In both populations, although not to the same extent, gene flow occurred between local and foreign gene pools. In both cases, however, we failed to detect obvious footprints of local adaptation (which should limit gene flow) and only weak support for genetic rescue (which should enhance gene flow). In the population where gene flow from foreign individuals was low, no clear differences were observed between the fitness of local, foreign, and F1 hybrid individuals. In the population where gene flow was high, foreign individuals were successful due to high mating success rather than high survival, and F1 hybrids had the same fitness as pure local offspring. These results suggest the importance of considering sexual selection, rather than just local adaptation and genetic rescue, when evaluating the determinants of success in small and recently founded populations.

Mating-regulated atrial proteases control reinsemination rates in Anopheles gambiae females

Anopheles gambiae mosquitoes are the most important vectors of human malaria. The reproductive success of these mosquitoes relies on a single copulation event after which the majority of females become permanently refractory to further mating. This refractory behavior is at least partially mediated by the male-synthetized steroid hormone 20-hydroxyecdysone (20E), which is packaged together with other seminal secretions into a gelatinous mating plug and transferred to the female atrium during mating. In this study, we show that two 20E-regulated chymotrypsin-like serine proteases specifically expressed in the reproductive tract of An. gambiae females play an important role in modulating the female susceptibility to mating. Silencing these proteases by RNA interference impairs correct plug processing and slows down the release of the steroid hormone 20E from the mating plug. In turn, depleting one of these proteases, the Mating Regulated Atrial Protease 1 (MatRAP1), reduces female refractoriness to further copulation, so that a significant proportion of females mate again. Microscopy analysis reveals that MatRAP1 is localized on a previously undetected peritrophic matrix-like structure surrounding the mating plug. These data provide novel insight into the molecular mechanisms shaping the post-mating biology of these important malaria vectors.

Differential gene expression reveals mechanisms related to habitat divergence between hybridizing orchids from the Neotropical coastal plains

BACKGROUND

Closely related hybridizing species are ideal systems for identifying genomic regions underlying adaptive divergence. Although gene expression plays a central role in determining ecologically-based phenotypic differences, few studies have inferred the role of gene expression for adaptive divergence in Neotropical systems. In this study, we conduct genome-wide expression analysis alongside soil elemental analysis in sympatric and allopatric populations of Epidendrum fulgens and E. puniceoluteum (Orchidaceae), which occur in contrasting adjacent habitats in the Neotropical coastal plains.

RESULTS

These species were highly differentiated by their gene expression profiles, as determined by 18-21% of transcripts. Gene ontology (GO) terms associated with reproductive processes were enriched according to comparisons between species in both allopatric and sympatric populations. Species showed differential expression in genes linked to salt and waterlogging tolerance according to comparisons between species in sympatry, and biological processes related to environmental stimulus appeared as representative among those transcripts associated with edaphic characteristics in each sympatric zone. Hybrids, in their turn, were well differentiated from E. fulgens, but exhibited a similar gene expression profile to flooding-tolerant E. puniceolutem. When compared with parental species, hybrids showed no transcripts with additive pattern of expression and increased expression for almost all transgressive transcripts.

CONCLUSIONS

This study sheds light on general mechanisms promoting ecological differentiation and assortative mating, and suggests candidate genes, such as those encoding catalase and calcium-dependent protein kinase, underling adaptation to harsh edaphic conditions in the Neotropical coastal plains. Moreover, it demonstrates that differential gene expression plays a central role in determining ecologically-based phenotypic differences among co-occurring species and their hybrids.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus)

BACKGROUND

Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression.

RESULTS

Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring.

CONCLUSIONS

This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

Transcriptomic divergence between upland and lowland ecotypes contributes to rice adaptation to a drought-prone agroecosystem

INTRODUCTION

Transcriptomic divergence drives plant ecological adaptation. Upland rice is differentiated in drought tolerance from lowland rice during its adaptation to the drought-prone environment. They provide a good system to learn the evolution of drought tolerance in rice.

METHODS AND RESULTS

We estimate morphological differences between the two rice ecotypes under well-watered and drought conditions, as well as their genetic and transcriptomic divergences by the high-throughput sequencing. Upland rice possesses higher expression diversity than lowland rice does. Thousands of genes exhibit expression divergences between the two rice ecotypes, which contributes to their morphological differences in drought tolerance. These transcriptomic divergences contribute to drought adaptation of upland rice during its domestication. Mutations in transcriptional regulatory regions, which cause presence and absence of cis-elements, are the cause of expression divergence. About 15.3% transcriptionally selected genes also receive sequence-based selection in upland or lowland ecotype. Some highly differentiated genes promote the transcriptomic divergence between rice ecotypes via gene co-expression network. In addition, we also detected transcriptomic trade-offs between drought tolerance and productivity.

DISCUSSION

Many key genes, which promote transcriptomic adaptation to drought in upland rice, have great prospective in breeding water-saving and drought-resistant rice. Meanwhile, appropriate strategies are required in breeding to overcome the potential transcriptomic trade-off.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.

Restoring ancestral phenotypes is a general pattern in gene expression evolution during adaptation to new environments in Tribolium castaneum

Plasticity and evolution are two processes allowing populations to respond to environmental changes, but how both are related and impact each other remains controversial. We studied plastic and evolutionary responses in gene expression of Tribolium castaneum after exposure of the beetles to new environments that differed from ancestral conditions in temperature, humidity or both. Using experimental evolution with 10 replicated lines per condition, we were able to demonstrate adaptation after 20 generations. We measured whole-transcriptome gene expression with RNA-sequencing to infer evolutionary and plastic changes. We found more evidence for changes in mean expression (shift in the intercept of reaction norms) in adapted lines than for changes in plasticity (shifts in slopes). Plasticity was mainly preserved in selected lines and was responsible for a large part of the phenotypic divergence in expression between ancestral and new conditions. However, we found that genes with the largest evolutionary changes in expression also evolved reduced plasticity and often showed expression levels closer to the ancestral stage. Results obtained in the three different conditions were similar, suggesting that restoration of ancestral expression levels during adaptation is a general evolutionary pattern. With a larger sample in the most stressful condition, we were able to detect a positive correlation between the proportion of genes with reversion of the ancestral plastic response and mean fitness per selection line.

Transcriptome profiling revealed potentially important roles of defensive gene expression in the divergence of insect biotypes: a case study with the cereal aphid Sitobion avenae

Background

Many insects can develop differential biotypes on variable host plants, but the underlying molecular factors and mechanisms are not well understood. To address this issue, transcriptome profiling analyses were conducted for two biotypes of the cereal aphid, Sitobion avenae (Fabricius), on both original and alternative plants.

Results

Comparisons between both biotypes generated 4174 differentially expressed unigenes (DEGs). In their response to host plant shift, 39 DEGs were shared by both biotypes, whereas 126 and 861 DEGs occurred only in biotypes 1 and 3, respectively. MMC (modulated modularity clustering) analyses showed that specific DEGs of biotypes 1 and 3 clustered into five and nine transcriptional modules, respectively. Among these DEGs, defense-related genes underwent intensive expression restructuring in both biotypes. However, biotype 3 was found to have relatively lower gene transcriptional plasticity than biotype 1. Gene enrichment analyses of the abovementioned modules showed functional divergence in defensive DEGs for the two biotypes in response to host transfer. The expression plasticity for some defense related genes was showed to be directly related to fecundity of S. avenae biotypes on both original and alternative plants, suggesting that expression plasticity of key defensive genes could have significant impacts on the adaptive potential and differentiation of S. avenae biotypes on different plants.

Conclusions

The divergence patterns of transcriptional plasticity in defense related genes may play important roles in the phenotypic evolution and differentiation of S. avenae biotypes. Our results can provide insights into the role of gene expression plasticity in the divergence of insect biotypes and adaptive evolution of insect populations.

Speciation and changes in male gene expression in Drosophila

It has long been acknowledged that changes in the regulation of gene expression may account for major organismal differences. However, we still do not fully understand how changes in gene expression evolve and how do such changes influence organisms' differences. We are even less aware of the impact such changes might have in restricting gene flow between species. Here, we focus on studies of gene expression and speciation in the Drosophila model. We review studies that have identified gene interactions in post-mating reproductive isolation and speciation, particularly those that modulate male gene expression. We also address studies that have experimentally manipulated changes in gene expression to test their effect in post-mating reproductive isolation. We highlight the need for a more in-depth analysis of the role of selection causing disrupted gene expression of such candidate genes in sterile/inviable hybrids. Moreover, we discuss the relevance to incorporate more routinely assays that simultaneously evaluate the potential effects of environmental factors and genetic background in modulating plastic responses in male genes and their potential role in speciation.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Comparative Transcriptomics Reveals Patterns of Adaptive Evolution Associated with Depth and Age Within Marine Rockfishes (Sebastes)

The genetic underpinnings that contribute to ecological adaptation and speciation are not completely understood, especially within marine ecosystems. These evolutionary processes can be elucidated by studying adaptive radiations, because they provide replicates of divergence within a given environment or time-frame. Marine rockfishes (genus Sebastes) are an adaptive radiation and unique model system for studying adaptive evolution in the marine realm. We investigated molecular evolution associated with ecological (depth) and life history (lifespan) divergence in 2 closely related clades of Sebastes. Brain transcriptomes were sequenced via RNA-Seq from 3 species within the subgenus Pteropodus and a pair of related congeners from the subgenus Sebastosomus in order to identify patterns of adaptive evolution. De novo assemblies from these transcriptomes were used to identify 3867 orthologous clusters, and genes subject to positive selection were identified based on all 5 species, depth, and lifespan. Within all our analyses, we identified hemoglobin subunit α to be under strong positive selection and is associated with the depth of occurrence. In our lifespan analysis we identified immune function genes under positive selection in association with maximum lifespan. This study provides insight on the molecular evolution of rockfishes and these candidate genes may provide a better understanding of how these subgenera radiated within the Northeast Pacific.

Codon Pairs are Phylogenetically Conserved: A comprehensive analysis of codon pairing conservation across the Tree of Life

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques.

Additive and mostly adaptive plastic responses of gene expression to multiple stress in Tribolium castaneum

Gene expression is known to be highly responsive to the environment and important for adjustment of metabolism but there is also growing evidence that differences in gene regulation contribute to species divergence and differences among locally adapted populations. However, most studies so far investigated populations when divergence had already occurred. Selection acting on expression levels at the onset of adaptation to an environmental change has not been characterized. Understanding the mechanisms is further complicated by the fact that environmental change is often multivariate, meaning that organisms are exposed to multiple stressors simultaneously with potentially interactive effects. Here we use a novel approach by combining fitness and whole-transcriptome data in a large-scale experiment to investigate responses to drought, heat and their combination in Tribolium castaneum. We found that fitness was reduced by both stressors and their combined effect was almost additive. Expression data showed that stressor responses were acting independently and did not interfere physiologically. Since we measured expression and fitness within the same individuals, we were able to estimate selection on gene expression levels. We found that variation in fitness can be attributed to gene expression variation and that selection pressures were environment dependent and opposite between control and stress conditions. We could further show that plastic responses of expression were largely adaptive, i.e. in the direction that should increase fitness.

Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes

Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Local adaptation fuels cryptic speciation in terrestrial annelids

Uncovering the genetic and evolutionary basis of cryptic speciation is a major focus of evolutionary biology. Next Generation Sequencing (NGS) allows the identification of genome-wide local adaptation signatures, but has rarely been applied to cryptic complexes - particularly in the soil milieu - as it is the case with integrative taxonomy. The earthworm genus Carpetania, comprising six previously suggested putative cryptic lineages, is a promising model to study the evolutionary phenomena shaping cryptic speciation in soil-dwelling lineages. Genotyping-By-Sequencing (GBS) was used to provide genome-wide information about genetic variability between 17 populations, and geometric morphometrics analyses of genital chaetae were performed to investigate unexplored cryptic morphological evolution. Genomic analyses revealed the existence of three cryptic species, with half of the previously-identified potential cryptic lineages clustering within them. Local adaptation was detected in more than 800 genes putatively involved in a plethora of biological functions (most notably reproduction, metabolism, immunological response and morphogenesis). Several genes with selection signatures showed shared mutations for each of the cryptic species, and genes under selection were enriched in functions related to regulation of transcription, including SNPs located in UTR regions. Finally, geometric morphometrics approaches partially confirmed the phylogenetic signal of relevant morphological characters such as genital chaetae. Our study therefore unveils that local adaptation and regulatory divergence are key evolutionary forces orchestrating genome evolution in soil fauna.

Molecular responses to freshwater limitation in the mangrove tree Avicennia germinans (Acanthaceae)

Environmental variation along the geographical space can shape populations by natural selection. In the context of global warming and changing precipitation regimes, it is crucial to understand the role of environmental heterogeneity in tropical trees adaptation, given their disproportional contribution to water and carbon biogeochemical cycles. Here, we investigated how heterogeneity in freshwater availability along tropical wetlands has influenced molecular variations of the black mangrove (Avicennia germinans). A total of 57 trees were sampled at seven sites differing markedly in precipitation regime and riverine freshwater inputs. Using 2,297 genome-wide single nucleotide polymorphic markers, we found signatures of natural selection by the association between variations in allele frequencies and environmental variables, including the precipitation of the warmest quarter and the annual precipitation. Additionally, we found candidate loci for selection based on statistical deviations from neutral expectations of interpopulation differentiation. Most candidate loci within transcribed sequences were functionally associated with central aspects of drought tolerance or plant response to drought. Moreover, our results suggest the occurrence of the rapid evolution of a population, probably in response to sudden and persistent limitations in plant access to soil water, following a road construction in 1974. Observations supporting rapid evolution included the reduction in tree size and changes in allele frequencies and in transcript expression associated with increased drought tolerance through the accumulation of osmoprotectants and antioxidants, biosynthesis of cuticles, protection against protein degradation, stomatal closure, photorespiration and photosynthesis. We describe a major role of spatial heterogeneity in freshwater availability in the specialization of this typically tropical tree.

Evolutionary genomics of host plant adaptation: insights from Drosophila

Variation in gene expression in response to the use of alternate host plants can reveal genetic and physiological mechanisms explaining why insect-host relationships vary from host specialism to generalism. Interpreting transcriptome variation relies on well-annotated genomes, making drosophilids valuable model systems, particularly those species with tractable ecological associations. Patterns of whole genome expression and alternate gene splicing in response to growth on different hosts have revealed expression of gene networks of known detoxification genes as well as novel functionally enriched genes of diverse metabolic and structural functions. Integrating trancriptomic responses with fitness differences and levels of phenotypic plasticity in response to alternate hosts will help to reveal the general nature of genotype-phenotype relationships.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

MicroRNA Gene Regulation in Extremely Young and Parallel Adaptive Radiations of Crater Lake Cichlid Fish

Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Selection of reference genes for flowering pathway analysis in the masting plants, Celmisia lyallii and Chionochloa pallens, under variable environmental conditions

Mast flowering is characterised by mass synchronised flowering at irregular intervals over a wide geographical area. An understanding of the molecular drivers of mast flowering requires expression analysis of key developmentally regulated gene(s). Reverse transcription-quantitative PCR is the gold standard technique used to assess expression of target gene(s) and to validate high-throughput sequencing data. Selection and validation of appropriate reference gene(s), used as normalisation factors in transcript abundance analysis, is an essential step to avoid ambiguous expression results. Eight candidate reference genes were assessed to select the best internal normalisation factors in naturally growing masting plants Chionochloa pallens and Celmisia lyallii. Statistical packages geNorm, Normfinder, BestKeeper, ΔCt and RefFinder were used to determine the expression stability in plants translocated to different altitudes and sampled across the season. GAPDH and PP2a in Celmisia and ExP and THP in Chionochloa were found to be the best pairs of reference genes for normalisation of the gene expression data. Our study revealed environmentally-induced changes in reference gene expression, information that will be utilised as we investigate flowering phenology of masting plants under global climatic change.

Mate choice and gene expression signatures associated with nutritional adaptation in the medfly (Ceratitis capitata)

Evolutionary responses to nutrition are key to understanding host shifts and the resulting potential for reproductive isolation. Experimental evolution has previously been used to describe the responses of the medfly (Ceratitis capitata) to larval diets with different nutritional properties. Within 30 generations this led to divergence in larval development time, egg to adult survival and adaptation in adult body size. Here we used mRNA-seq to identify differences in gene expression patterns in these same populations, using males from the 60^th generation of nutritional selection. We validated differential expression by using qRT-PCR and found that genes linked to metabolism, oxidative phosphorylation and proteolysis were significantly over-represented among the differentially expressed genes. The results provide the first genome-wide survey of the putative mechanisms underpinning evolved responses to nutritional adaptation. In addition, we tested the hypothesis that nutritional adaptation can alter mating patterns. We found evidence for assortative mating by diet at generation 60, but not 90. Hence, the pattern was variable across generations and there was no evidence overall for any isolating mating divergence between the lines. Overall, the results provide insight into the mechanisms underpinning dietary adaptation and extend our knowledge of which traits represent core responses to nutritional selection.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

Time matters! Developmental shift in gene expression between the head and the trunk region of the cichlid fish Astatotilapia burtoni

BACKGROUND

Differential gene expression can be translated into differing phenotypic traits. Especially during embryogenesis, specific gene expression networks regulate the development of different body structures. Cichlid fishes, with their impressive phenotypic diversity and propensity to radiate, are an emerging model system in the genomics era. Here we set out to investigate gene expression throughout development in the well-studied cichlid fish Astatotilapia burtoni, native to Lake Tanganyika and its affluent rivers.

RESULTS

Combining RNA-sequencing from different developmental time points as well as integrating adult gene expression data, we constructed a new genome annotation for A. burtoni comprising 103,253 transcripts (stemming from 52,584 genomic loci) as well as a new reference transcriptome set. We compared our transcriptome to the available reference genome, redefining transcripts and adding new annotations. We show that about half of these transcripts have coding potential. We also characterize transcripts that are not present in the genome assembly. Next, using our newly constructed comprehensive reference transcriptome, we characterized differential gene expression through time and showed that gene expression is shifted between different body parts. We constructed a gene expression network that identified connected genes responsible for particular phenotypes and made use of it to focus on genes under potential positive selection in A. burtoni, which were implicated in fin development and vision.

CONCLUSIONS

We provide new genomic resources for the cichlid fish Astatotilapia burtoni, which will contribute to its further establishment as a model system. Tracing gene expression through time, we identified gene networks underlying particular functions, which will help to understand the genetic basis of phenotypic diversity in cichlids.

Climatic niche evolution in the viviparous Sceloporus torquatus group (Squamata: Phrynosomatidae)

The cold-climate hypothesis maintains that viviparity arose as a means to prevent increased egg mortality in nests owing to low temperatures, and this hypothesis represents the primary and most strongly supported explanation for the evolution of viviparity in reptiles. In this regard, certain authors have stated that viviparous species will exhibit speciation via climatic niche conservatism, with similar climatic niches being observed in allopatric sister species. However, this prediction remains to be tested with bioclimatic variables relevant to each viviparous group. In the present study, we examined climatic niche evolution in a group of North American viviparous lizards to determine whether their diversification is linked to phylogenetic niche conservatism (PNC). We evaluated the phylogenetic signal and trait evolution of individual bioclimatic variables and principal component (PC) scores of a PC analysis, along with reconstructions of ancestral climate tolerances. The results suggest that diversification of the Sceloporus torquatus group species is associated with both niche differentiation and PNC. Furthermore, we did not observe PNC across nearly all bioclimatic variables and in PC2 and PC3. However, in Precipitation Seasonality (Bio15), in Precipitation of Coldest Quarter (Bio19) and in PC1 (weakly associated with variability of temperature), we did observe PNC. Additionally, variation of the scores along the phylogeny and Pagel's delta (δ) >1 of PC3 suggests a fast, recent evolution to dry conditions in the clade that sustains S. serrifer.

RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry

Adaptive divergence with gene flow often results in complex patterns of variation within taxa exhibiting substantial ecological differences among populations. One example where this may have occurred is the parallel evolution of freshwater-resident nonparasitic lampreys from anadromous-parasitic ancestors. Previous studies have focused on transitions between these two phenotypic extremes, but here, we considered more complex evolutionary scenarios where an intermediate freshwater form that remains parasitic is found sympatrically with the other two ecotypes. Using population genomic analysis (restriction-associated DNA sequencing), we found that a freshwater-parasitic ecotype was highly distinct from an anadromous-parasitic form (Q_lake-P  = 96.8%, F_st  = 0.154), but that a freshwater-nonparasitic form was almost completely admixed in Loch Lomond, Scotland. Demographic reconstructions indicated that both freshwater populations likely derived from a common freshwater ancestor. However, while the nonparasitic ecotype has experienced high levels of introgression from the anadromous-parasitic ecotype (Q_anad-P  = 37.7%), there is no evidence of introgression into the freshwater-parasitic ecotype. Paradoxically, mate choice experiments predicted high potential for gene flow: Males from all ecotypes were stimulated to spawn with freshwater-parasitic females, which released gametes in response to all ecotypes. Differentially fixed single nucleotide polymorphisms identified genes associated with growth and development, which could possibly influence the timing of metamorphosis, resulting in significant ecological differences between forms. This suggests that multiple lamprey ecotypes can persist in sympatry following shifts in adaptive peaks, due to environmental change during their repeated colonization of post-glacial regions, followed by periods of extensive gene flow among such diverging populations.