The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing

Multiple Tissues Transcriptome of Zig-Zag Eel (Mastacembelus armatus) with Different Growth Rates

In order to explore the main regulatory genes and related pathways of growth traits, transcriptome sequencing was first performed on the brain, liver, and muscle tissues of 3-month-old M. armatus with different growth rates. By comparative transcriptome analysis of fast-growing and slow-growing groups of M. armatus, a total of 2887 DEGs were screened, of which 59 up-regulated genes and 105 down-regulated genes were detected in the brain, 146 up-regulated genes and 202 down-regulated genes were detected in the liver, and 529 up-regulated genes and 1846 down-regulated genes were detected in muscle, including insulin-like growth factor binding protein 1a (IGFBP1A), insulin-like growth factor binding protein 1b (IGFBP1B), myosin, light chain 1 (MYL1), and myoglobin (MB). Through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, we identified a total of 288 significantly enriched GO entries and 68 significantly enriched KEGG pathways related to growth, such as skeletal muscle tissue development, insulin-like growth factor binding, and the mitotic cell cycle. These key genes and signaling pathways may play a key role in regulating the growth of M. armatus. Digging into the regulatory mechanisms of these key genes will provide a theoretical basis for further exploration of the molecular mechanisms related to the growth and development of M. armatus, and help to breed new varieties of M. armatus with rapid growth.

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.

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

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.

Transcriptome-Wide Patterns of the Genetic and Expression Variations in Two Sympatric Schizothoracine Fishes in a Tibetan Plateau Glacier Lake

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of F_ST = 0.25–0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations.

Molecular Response to High Hydrostatic Pressure: Time-Series Transcriptomic Analysis of Shallow-Water Sea Cucumber Apostichopus japonicus

Hydrostatic pressure is a key environmental factor constraining the benthic migration of shallow-water invertebrates. Although many studies have examined the physiological effects of high hydrostatic pressure on shallow-water invertebrates, the molecular response to high pressure is not fully understood. This question has received increasing attention because ocean warming is forcing the bathymetric migrations of shallow-water invertebrates. Here, we applied time-series transcriptomic analysis to high-pressure incubated and atmospheric pressure-recovered shallow-water sea cucumber (Apostichopus japonicus) to address this question. A total of 44 samples from 15 experimental groups were sequenced. Our results showed that most genes responded to pressure stress at the beginning when pressure was changed, but significant differences of gene expression appeared after 4 to 6 h. Transcription was the most sensitive biological process responding to high-pressure exposure, which was enriched among up-regulated genes after 2 h, followed by ubiquitination (4 h), endocytosis (6 h), stress response (6 h), methylation regulation (24 h), and transmembrane transportation (24 h). After high-pressure incubation, all these biological processes remained up-regulated within 4–6 h at atmospheric pressure. Overall, our results revealed the dynamic transcriptional response of A. japonicus to high-pressure exposure. Additionally, few quantitative or functional responses related to A. japonicus on transcriptional level were introduced by hydrostatic pressure changes after 1 h, and main biological responses were introduced after 4 h, suggesting that, when hydrostatic pressure is the mainly changed environmental factor, it will be better to fix sea cucumber samples for transcriptomic analysis within 1 h, but 4 h will be also acceptable.

Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations

Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes.

The transcriptional correlates of divergent electric organ discharges in Paramormyrops electric fish

Background

Understanding the genomic basis of phenotypic diversity can be greatly facilitated by examining adaptive radiations with hypervariable traits. In this study, we focus on a rapidly diverged species group of mormyrid electric fish in the genus Paramormyrops, which are characterized by extensive phenotypic variation in electric organ discharges (EODs). The main components of EOD diversity are waveform duration, complexity and polarity. Using an RNA-sequencing based approach, we sought to identify gene expression correlates for each of these EOD waveform features by comparing 11 specimens of Paramormyrops that exhibit variation in these features.

Results

Patterns of gene expression among Paramormyrops are highly correlated, and 3274 genes (16%) were differentially expressed. Using our most restrictive criteria, we detected 145–183 differentially expressed genes correlated with each EOD feature, with little overlap between them. The predicted functions of several of these genes are related to extracellular matrix, cation homeostasis, lipid metabolism, and cytoskeletal and sarcomeric proteins. These genes are of significant interest given the known morphological differences between electric organs that underlie differences in the EOD waveform features studied.

Conclusions

In this study, we identified plausible candidate genes that may contribute to phenotypic differences in EOD waveforms among a rapidly diverged group of mormyrid electric fish. These genes may be important targets of selection in the evolution of species-specific differences in mate-recognition signals.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

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.

Differences in Effective Ploidy Drive Genome-Wide Endosperm Expression Polarization and Seed Failure in Wild Tomato Hybrids

Parental imbalances in the endosperm leading to impaired development and eventual hybrid seed failure are common causes of postzygotic isolation in flowering plants. Endosperm sensitivity to parental dosage is reflected by canonical phenotypes of "parental excess" in reciprocal interploid crosses. Moreover, parental-excess traits are also evident in many homoploid interspecific crosses, potentially reflecting among-lineage variation in "effective ploidy" driven by endosperm properties. However, the genetic basis of effective ploidy is unknown and genome-wide expression perturbations in parental-excess endosperms from homoploid crosses have yet to be reported. The tomato clade (Solanum section Lycopersicon), encompassing closely related diploids with partial-to-complete hybrid seed failure, provides outstanding opportunities to study these issues. Here, we compared replicated endosperm transcriptomes from six crosses within and among three wild tomato lineages. Strikingly, strongly inviable hybrid crosses displayed conspicuous, asymmetric expression perturbations that mirror previously characterized parental-excess phenotypes. Solanum peruvianum, the species inferred to have evolved higher effective ploidy than the other two, drove expression landscape polarization between maternal and paternal roles. This global expression divergence was mirrored in functionally important gene families such as MADS-box transcription factors and E3 ubiquitin ligases, and revealed differences in cell cycle tuning that match phenotypic differences in developing endosperm and mature seed size between reciprocal crosses. Our work starts to uncover the complex interactions between expression divergence, parental conflict, and hybrid seed failure that likely contributed to plant diversity.

Good Citizenship Made Easy: A Step-by-Step Guide to Submitting RNA-Seq Data to NCBI

The analysis of transcriptome data from non-model organisms contributes to our understanding of diverse aspects of evolutionary biology, including developmental processes, speciation, adaptation, and extinction. Underlying this diversity is one shared feature, the generation of enormous amounts of sequence data. Data availability requirements in most journals oblige researchers to make their raw transcriptome data publicly available, and the databases housed at the National Center for Biotechnology Information (NCBI) are a popular choice for data deposition. Unfortunately, the successful submission of raw sequences to the Sequence Read Archive (SRA) and transcriptome assemblies to the Transcriptome Shotgun Assembly (TSA) can be challenging for novice users, significantly delaying data availability and publication. Here we present two comprehensive protocols for submitting RNA-Seq data to NCBI databases, accompanied by an easy-to-use website that facilitates the timely submission of data by researchers of any experience level. © 2018 by John Wiley & Sons, Inc.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

Holobionts and ecological speciation: the intestinal microbiota of lake whitefish species pairs

BACKGROUND

It is well established that symbionts have considerable impact on their host, yet the investigation of the possible role of the holobiont in the host's speciation process is still in its infancy. In this study, we compared the intestinal microbiota among five sympatric pairs of dwarf (limnetic) and normal (benthic) lake whitefish Coregonus clupeaformis representing a continuum in the early stage of ecological speciation. We sequenced the 16s rRNA gene V3-V4 regions of the intestinal microbiota present in a total of 108 wild sympatric dwarf and normal whitefish as well as the water bacterial community from five lakes to (i) test for differences between the whitefish intestinal microbiota and the water bacterial community and (ii) test for parallelism in the intestinal microbiota of dwarf and normal whitefish.

RESULTS

The water bacterial community was distinct from the intestinal microbiota, indicating that intestinal microbiota did not reflect the environment, but rather the intrinsic properties of the host microbiota. Our results revealed a strong influence of the host (dwarf or normal) on the intestinal microbiota with pronounced conservation of the core intestinal microbiota (mean ~ 44% of shared genera). However, no clear evidence for parallelism was observed, whereby non-parallel differences between dwarf and normal whitefish were observed in three of the lakes while similar taxonomic composition was observed for the two other species pairs.

CONCLUSIONS

This absence of parallelism across dwarf vs. normal whitefish microbiota highlighted the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota.

De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species

Background

Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date.

Results

We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species.

Conclusion

New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.

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

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

Heritable gene expression differences between lake and stream stickleback include both parallel and antiparallel components

The repeated phenotypic patterns that characterize populations undergoing parallel evolution provide support for a deterministic role of adaptation by natural selection. Determining the level of parallelism also at the genetic level is thus central to our understanding of how natural selection works. Many studies have looked for repeated genomic patterns in natural populations, but work on gene expression is less common. The studies that have examined gene expression have found some support for parallelism, but those studies almost always used samples collected from the wild that potentially confounds the effects of plasticity with heritable differences. Here we use two independent pairs of lake and stream threespine stickleback (Gasterosteus aculeatus) raised in common garden conditions to assess both parallel and antiparallel (that is, similar versus different directions of lake–stream expression divergence in the two watersheds) heritable gene expression differences as measured by total RNA sequencing. We find that more genes than expected by chance show either parallel (22 genes, 0.18% of expressed genes) or antiparallel (24 genes, 0.20% of expressed genes) lake–stream expression differences. These results correspond well with previous genomic studies in stickleback ecotype pairs that found similar levels of parallelism. We suggest that parallelism might be similarly constrained at the genomic and transcriptomic levels.

Large-scale SNP screenings identify markers linked with GCRV resistant traits through transcriptomes of individuals and cell lines in Ctenopharyngodon idella

Grass carp (Ctenopharyngodon idella) is an important economic species in freshwater aquaculture and its industry has been confined due to variety degeneration and frequent diseases. Marker-assisted selection is a feasible method for selective breeding of new varieties. Transcriptome data have greatly facilitated high-throughput single nucleotide polymorphism (SNP) marker discovery and phenotype association study. In this study, we gained a total of 25,981 and 5,775 high quality SNPs in two transcriptomes from individuals and cell lines, respectively. Comparative transcriptome analysis identified 413 and 832 grass carp reovirus (GCRV)-resistant-association SNPs as well as 1,381 and 1,606 GCRV-susceptible-association SNPs in individuals and cell lines, respectively. Integrated analysis indicated 22 genes with single SNP share common resistant/susceptible traits in two transcriptomes. Furthermore, we infected grass carp with GCRV, genotyping and association analyses were performed, and 9 in 22 SNPs were confirmed by PCR-RFLP. Meanwhile, mRNA expression profiles of 6 genes containing confirmed SNPs were examined by qRT-PCR. The results demonstrated that mRNA expressions were significant differences in resistant/susceptible individuals and cell lines. The present study develops an important strategy for high throughput screening of phenotype association genetic markers and the results will serve in grass carp breeding for GCRV resistance.

Transcriptomic comparison of invasive bigheaded carps (Hypophthalmichthys nobilis and Hypophthalmichthys molitrix) and their hybrids

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix), collectively called bigheaded carps, are invasive species in the Mississippi River Basin (MRB). Interspecific hybridization between bigheaded carps has been considered rare within their native rivers in China; however, it is prevalent in the MRB. We conducted de novo transcriptome analysis of pure and hybrid bigheaded carps and obtained 40,759 to 51,706 transcripts for pure, F₁ hybrid, and backcross bigheaded carps. The search against protein databases resulted in 20,336–28,133 annotated transcripts (over 50% of the transcriptome) with over 13,000 transcripts mapped to 23 Gene Ontology biological processes and 127 KEGG metabolic pathways. More transcripts were detected in silver carp than in bighead carp; however, comparable numbers of transcripts were annotated. Transcriptomic variation detected between two F₁ hybrids may indicate a potential loss of fitness in hybrids. The neighbor‐joining distance tree constructed using over 2,500 one‐to‐one orthologous sequences suggests transcriptomes could be used to infer the history of introgression and hybridization. Moreover, we detected 24,792 candidate SNPs that can be used to identify different species. The transcriptomes, orthologous sequences, and candidate SNPs obtained in this study should provide further knowledge of interspecific hybridization and introgression.

Comparison of Ribosomal RNA Removal Methods for Transcriptome Sequencing Workflows in Teleost Fish

RNA sequencing (RNA-Seq) is becoming the standard for transcriptome analysis. Removal of contaminating ribosomal RNA (rRNA) is a priority in the preparation of libraries suitable for sequencing. These methods have been well documented in mammals but typically require some optimization for lower vertebrates. Three commercial kits, including Dynabeads mRNA Purification Kit, RiboMinus Eukaryote System v2, and Ribo-Zero Gold rRNA Removal Kit were examined for the ability to remove rRNAs from rainbow trout (Oncorhynchus mykiss) RNA isolations. Total RNA was isolated from liver and muscle tissue samples (n = 24) and rRNAs removed using one of the three kits. Samples were analyzed visually on the Agilent Bioanalyzer and by Illumina RNA-seq, screening for Oncorhynchus rRNAs. There were significant differences between the kits in regards to their ability to remove rRNA, ranging from 2.74% - 10.94% rRNA sequences left behind per kit on average. Using the Bioanalyzer to evaluate ribosomal contamination in rRNA-depleted samples for RNA-Seq was good for detecting samples with higher concentrations of rRNA (>5%), but not very accurate at lower levels. Although all three kits were able to remove a substantial portion of the rRNA from different fish tissues, the Ribo-Zero Gold rRNA Removal Kit eliminated significantly more contaminating ribosomal RNAs than the others.

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations - 'outlier loci' (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false-positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated FST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated FST . Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of FST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased FST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated FST in geographically distant locations. We discuss possible explanations for this unexpected result.

Embryonic gene expression of Coregonus palaea (whitefish) under pathogen stress as analyzed by high-throughput RNA-sequencing

Most fishes produce free-living embryos that are exposed to environmental stressors immediately following fertilization, including pathogenic microorganisms. Initial immune protection of embryos involves the chorion, as a protective barrier, and maternally-allocated antimicrobial compounds. At later developmental stages, host-genetic effects influence susceptibility and tolerance, suggesting a direct interaction between embryo genes and pathogens. So far, only a few host genes could be identified that correlate with embryonic survival under pathogen stress in salmonids. Here, we utilized high-throughput RNA-sequencing in order to describe the transcriptional response of a non-model fish, the Alpine whitefish Coregonus palaea, to infection, both in terms of host genes that are likely manipulated by the pathogen, and those involved in an early putative immune response. Embryos were produced in vitro, raised individually, and exposed at the late-eyed stage to a virulent strain of the opportunistic fish pathogen Pseudomonas fluorescens. The pseudomonad increased embryonic mortality and affected gene expression substantially. For example, essential, upregulated metabolic pathways in embryos under pathogen stress included ion binding pathways, aminoacyl-tRNA-biosynthesis, and the production of arginine and proline, most probably mediated by the pathogen for its proliferation. Most prominently downregulated transcripts comprised the biosynthesis of unsaturated fatty acids, the citrate cycle, and various isoforms of b-cell transcription factors. These factors have been shown to play a significant role in host blood cell differentiation and renewal. With regard to specific immune functions, differentially expressed transcripts mapped to the complement cascade, MHC class I and II, TNF-alpha, and T-cell differentiation proteins. The results of this study reveal insights into how P. fluorescens impairs the development of whitefish embryos and set a foundation for future studies investigating host pathogen interactions in fish embryos.

Widespread and Adaptive Alterations in Genome-Wide Gene Expression Associated with Ecological Divergence of Two Oryza Species

Ecological speciation is a common mechanism by which new species arise. Despite great efforts, the role of gene expression in ecological divergence and speciation is poorly understood. Here, we conducted a genome-wide gene expression investigation of two Oryza species that are evolutionarily young and distinct in ecology and morphology. Using digital gene expression technology and the paired-end RNA sequencing method, we obtained 21,415 expressed genes across three reproduction-related tissues. Of them, approximately 8% (1,717) differed significantly in expression levels between the two species and these differentially expressed genes are randomly distributed across the genome. Moreover, 62% (1,064) of the differentially expressed genes exhibited a signature of directional selection in at least one species. Importantly, the genes with differential expression between species evolved more rapidly at the 5' flanking sequences than the genes without differential expression relative to coding sequences, suggesting that cis-regulatory changes are likely adaptive and play an important role in the ecological divergence of the two species. Finally, we showed evidence of significant differentiation between species in phenotype traits and observed that genes with differential expression were overrepresented with functional terms involving phenotypic and ecological differentiation between the two species, including reproduction- and stress-related characteristics. Our findings demonstrate that ecological speciation is associated with widespread and adaptive alterations in genome-wide gene expression and provide new insights into the importance of regulatory evolution in ecological speciation in plants.

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size-matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming "dwarf" ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I-IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Divergence in expression of candidate genes for the smoltification process between juvenile resident rainbow and anadromous steelhead trout

Rainbow and steelhead trout (Oncorhynchus mykiss), among other salmonid fishes, exhibit tremendous life history diversity, foremost of which is variation in migratory propensity. While some individuals possess the ability to undertake an anadromous marine migration, others remain resident in freshwater throughout their life cycle. Those that will migrate undergo tremendous physiological, morphological, and behavioral transformations in a process called smoltification which transitions freshwater-adapted parr to marine-adapted smolts. While the behavior, ecology, and physiology of smoltification are well described, our understanding of the proximate genetic mechanisms that trigger the process are not well known. Quantitative genetic analyses have identified several genomic regions associated with smoltification and migration-related traits within this species. Here we investigate the divergence in gene expression of 18 functional and positional candidate genes for the smoltification process in the brain, gill, and liver tissues of migratory smolts, resident parr, and precocious mature male trout at the developmental stage of out-migration. Our analysis reveals several genes differentially expressed between life history classes and validates the candidate nature of several genes in the parr-smolt transformation including Clock1α, FSHβ, GR, GH2, GHR1, GHR2, NDK7, p53, SC6a7, Taldo1, THRα, THRβ, and Vdac2.

Differences in the metabolic response to temperature acclimation in nine-spined stickleback (Pungitius pungitius) populations from contrasting thermal environments

Metabolic responses to temperature changes are crucial for maintaining the energy balance of an individual under seasonal temperature fluctuations. To understand how such responses differ in recently isolated populations (<11,000 years), we studied four Baltic populations of the nine-spined stickleback (Pungitius pungitius L.) from coastal locations (seasonal temperature range, 0-29°C) and from colder, more thermally stable spring-fed ponds (1-19°C). Salinity and predation pressure also differed between these locations. We acclimatized wild-caught fish to 6, 11, and 19°C in common garden conditions for 4-6 months and determined their aerobic scope and hepatosomatic index (HSI). The freshwater fish from the colder (2-14°C), predator-free pond population exhibited complete temperature compensation for their aerobic scope, whereas the coastal populations underwent metabolic rate reduction during the cold treatment. Coastal populations had higher HSI than the colder pond population at all temperatures, with cold acclimation accentuating this effect. The metabolic rates and HSI for freshwater fish from the pond with higher predation pressure were more similar to those of the coastal ones. Our results suggest that ontogenic effects and/or genetic differentiation are responsible for differential energy storage and metabolic responses between these populations. This work demonstrates the metabolic versatility of the nine-spined stickleback and the pertinence of an energetic framework to better understand potential local adaptations. It also demonstrates that instead of using a single acclimation temperature thermal reaction norms should be compared when studying individuals originating from different thermal environments in a common garden setting.

Transcriptome analysis of a petal anthocyanin polymorphism in the arctic mustard, Parrya nudicaulis

Angiosperms are renown for their diversity of flower colors. Often considered adaptations to pollinators, the most common underlying pigments, anthocyanins, are also involved in plants’ stress response. Although the anthocyanin biosynthetic pathway is well characterized across many angiosperms and is composed of a few candidate genes, the consequences of blocking this pathway and producing white flowers has not been investigated at the transcriptome scale. We take a transcriptome-wide approach to compare expression differences between purple and white petal buds in the arctic mustard, Parrya nudicaulis, to determine which genes’ expression are consistently correlated with flower color. Using mRNA-Seq and de novo transcriptome assembly, we assembled an average of 722 bp per gene (49.81% coding sequence based on the A. thaliana homolog) for 12,795 genes from the petal buds of a pair of purple and white samples. Our results correlate strongly with qRT-PCR analysis of nine candidate genes in the anthocyanin biosynthetic pathway where chalcone synthase has the greatest difference in expression between color morphs (P/W = ∼7×). Among the most consistently differentially expressed genes between purple and white samples, we found 3× more genes with higher expression in white petals than in purple petals. These include four unknown genes, two drought-response genes (CDSP32, ERD5), a cold-response gene (GR-RBP2), and a pathogen defense gene (DND1). Gene ontology analysis of the top 2% of genes with greater expression in white relative to purple petals revealed enrichment in genes associated with stress responses including cold, drought and pathogen defense. Unlike the uniform downregulation of chalcone synthase that may be directly involved in the loss of petal anthocyanins, the variable expression of several genes with greater expression in white petals suggest that the physiological and ecological consequences of having white petals may be microenvironment-dependent.

Phenotype-environment association of the oxygen transport system in trimorphic European whitefish (Coregonus lavaretus) populations

Replicated adaptive radiation events, typified by phenotypic divergence across resource axes, provide important insight into the eco-evolutionary dynamics that lead to the formation of new species. Here, we show that in trimorphic adaptive radiations of European whitefish (Coregonus lavaretus), divergence of the oxygen transport system has occurred across the pelagic/littoral (shallow)-profundal (deep) resource axis, and at multiple biological scales. Profundal whitefish exhibited significantly larger red blood cells (RBCs), a greater proportion of cathodic hemoglobin protein components, and higher hemoglobin transcript abundance in kidney compared to littoral and pelagic morphs. Hemoglobin transcript abundance in brain and gill, but not kidney, and anodic hemoglobin protein component diversity in blood were also linked to variation at an intronic single nucleotide polymorphism (SNP). As the whitefish morphs differ in population genetic structure at this SNP, hemoglobin transcript and protein divergence between profundal and pelagic/littoral morphs is likely being driven by genetic divergence. Our findings, along with our previous work on lake whitefish, highlight the importance of the oxygen transport system to the postglacial colonization of novel lacustrine environments by whitefish throughout the northern hemisphere.

Characterization of the house sparrow (Passer domesticus) transcriptome: a resource for molecular ecology and immunogenetics

The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.

RNA-Seq technology and its application in fish transcriptomics

High-throughput sequencing technologies, also known as next-generation sequencing (NGS) technologies, have revolutionized the way that genomic research is advancing. In addition to the static genome, these state-of-art technologies have been recently exploited to analyze the dynamic transcriptome, and the resulting technology is termed RNA sequencing (RNA-seq). RNA-seq is free from many limitations of other transcriptomic approaches, such as microarray and tag-based sequencing method. Although RNA-seq has only been available for a short time, studies using this method have completely changed our perspective of the breadth and depth of eukaryotic transcriptomes. In terms of the transcriptomics of teleost fishes, both model and non-model species have benefited from the RNA-seq approach and have undergone tremendous advances in the past several years. RNA-seq has helped not only in mapping and annotating fish transcriptome but also in our understanding of many biological processes in fish, such as development, adaptive evolution, host immune response, and stress response. In this review, we first provide an overview of each step of RNA-seq from library construction to the bioinformatic analysis of the data. We then summarize and discuss the recent biological insights obtained from the RNA-seq studies in a variety of fish species.

Accuracy of allele frequency estimation using pooled RNA-Seq

For nonmodel organisms, genome-wide information that describes functionally relevant variation may be obtained by RNA-Seq following de novo transcriptome assembly. While sequencing has become relatively inexpensive, the preparation of a large number of sequencing libraries remains prohibitively expensive for population genetic analyses of nonmodel species. Pooling samples may be then an attractive alternative. To test whether pooled RNA-Seq accurately predicts true allele frequencies, we analysed the liver transcriptomes of 10 bank voles. Each sample was sequenced both as an individually barcoded library and as a part of a pool. Equal amounts of total RNA from each vole were pooled prior to mRNA selection and library construction. Reads were mapped onto the de novo assembled reference transcriptome. High-quality genotypes for individual voles, determined for 23,682 SNPs, provided information on 'true' allele frequencies; allele frequencies estimated from the pool were then compared with these values. 'True' frequencies and those estimated from the pool were highly correlated. Mean relative estimation error was 21% and did not depend on expression level. However, we also observed a minor effect of interindividual variation in gene expression and allele-specific gene expression influencing allele frequency estimation accuracy. Moreover, we observed strong negative relationship between minor allele frequency and relative estimation error. Our results indicate that pooled RNA-Seq exhibits accuracy comparable with pooled genome resequencing, but variation in expression level between individuals should be assessed and accounted for. This should help in taking account the difference in accuracy between conservatively expressed transcripts and these which are variable in expression level.

Probing functional polymorphisms in the dengue vector, Aedes aegypti

BACKGROUND

Dengue is the most prevalent arboviral disease world-wide and its primary vector is the mosquito Aedes aegypti. The current lack of commercially-available vaccines makes control of vector populations the only effective strategy to prevent dengue transmission. Aedes aegypti geographic populations exhibit great variability in insecticide resistance and susceptibility to dengue infection. The characterization of single nucleotide polymorphisms (SNPs) as molecular markers to study quantitatively this variation is needed greatly because this species has a low abundance of microsatellite markers and limited known restriction fragments length polymorphisms (RFLPs) and single-strand conformation polymorphism (SSCP) markers.

RESULTS

We used RNA-seq to characterize SNPs in three Ae. aegypti strains, including the Liverpool (LVP) strain, from which the current genome annotation is derived. We identified 131,764 unique genome locations with at least one alternative nucleotide to what is reported in the reference annotation. These comprised changes in both open-reading frames (ORFs) and untranslated regions (UTRs) of transcripts. An in depth-look at sequence variation in immunity genes revealed that those associated with autophagy, MD2-like receptors and Peptidoglycan Recognition Proteins had more sequence variation in their 3'UTRs than mutations associated with non-synonymous changes. This supports the conclusion that these genes had maintained their functional specificity while being adapted to different regulatory domains. In contrast, a number of peroxidases, serpins and Clip-domain serine proteases exhibited conservation of putative UTR regulatory sequences while displaying diversification of the ORFs. Transcriptome evidence also was found for ~2500 novel transcriptional units (NTUs) not annotated in the reference genome.

CONCLUSIONS

The transcriptome-wide assessment of within and inter-strain polymorphisms in Ae. aegypti adds considerably to the number of molecular markers available for genetic studies in this mosquito. Additionally, data supporting NTU discovery emphasizes the need for continuous amendments of the reference genome annotation.

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

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

Effect of sexual maturation on muscle gene expression of rainbow trout: RNA-Seq approach

Muscle degradation occurs as a response to various physiological states that are regulated by specific molecular mechanisms. Previously, we characterized the metabolic changes of muscle deterioration of the female rainbow trout at full sexual maturity and spawning (Salem et al., Physiol. Genomics 2006;28:33–45; J. Proteomics 2010;73:778–789). Muscle deterioration in this model represents nutrient mobilization as a response to the energetic overdemands of the egg/ovarian growth phase. Our recent studies showed that most of the changes in muscle growth and quality start 2–3 months before spawning. Gravid fish exhibited reduced intramuscular fat that is lower in saturated and monounsaturated fatty acids and higher in polyunsaturated fatty acids compared to sterile fish. In this study, RNA-Seq was used to explain the mechanisms underlying changes during this phase of sexual maturity. Furthermore, to minimize changes due to nutrient deficits, fish were fed on a high-plane of nutrition. The RNA-Seq technique identified a gene expression signature that is consistent with metabolic changes of gravid fish. Gravid fish exhibited increased abundance of transcripts in metabolic pathways of fatty acid degradation and up-regulated expression of genes involved in biosynthesis of unsaturated fatty acids. In addition, increased expression of genes involved in the citric acid cycle and oxidative phosphorylation was observed for gravid fish. This muscle transcriptomic signature of fish fed on a high nutritional plane is quite distinct from that previously described for fish at terminal stages of maturity and suggest that female rainbow trout approaching spawning, on high nutritional planes, likely mobilize intramuscular fat rather than protein to support gonadal maturation.

Functional and population genomic divergence within and between two species of killifish adapted to different osmotic niches

Adaptation to salinity affects species distributions, promotes speciation, and guides many evolutionary patterns in fishes. To uncover the basis of a complex trait like osmoregulation, genome-level analyses are sensible. We combine population genomic scans with genome expression profiling to discover candidate genes and pathways associated with divergence between osmotic environments. We compared transcriptome sequence divergence between multiple freshwater and saltwater populations of the rainwater killifish, Lucania parva. We also compared sequence divergence between L. parva and its sister species, Lucania goodei, a freshwater specialist. We found highly differentiated single nucleotide polymorphisms (SNPs) between freshwater and saltwater L. parva populations in cell junction and ion transport genes, including V-type H(+) ATPase. Between species, we found divergence in reproduction and osmotic stress genes. Genes that were differentially expressed between species during osmotic acclimation included genes involved in ion transport and cell volume regulation. Gene sets that were divergent in coding sequence and divergent in expression did not overlap, although they did converge in function. Like many studies using genomic scans, our approach may miss some loci that contribute to adaptation but have complicated patterns of allelic variation. Our study suggests that gene expression and coding sequence may evolve independently as populations adapt to a complex physiological challenge.

Targeted sequence capture and resequencing implies a predominant role of regulatory regions in the divergence of a sympatric lake whitefish species pair (Coregonus clupeaformis)

Latest technological developments in evolutionary biology bring new challenges in documenting the intricate genetic architecture of species in the process of divergence. Sympatric populations of lake whitefish represent one of the key systems to investigate this issue. Despite the value of random genotype-by-sequencing methods and decreasing cost of sequencing technologies, it remains challenging to investigate variation in coding regions, especially in the case of recently duplicated genomes as in salmonids, as this greatly complicates whole genome resequencing. We thus designed a sequence capture array targeting 2773 annotated genes to document the nature and the extent of genomic divergence between sympatric dwarf and normal whitefish. Among the 2728 genes successfully captured, a total of 2182 coding and 10,415 noncoding putative single-nucleotide polymorphisms (SNPs) were identified after applying a first set of basic filters. A genome scan with a quality-refined selection of 2203 SNPs identified 267 outlier SNPs in 210 candidate genes located in genomic regions potentially involved in whitefish divergence and reproductive isolation. We found highly heterogeneous FST estimates among SNP loci. There was an overall low level of coding polymorphism, with a predominance of noncoding mutations among outliers. The heterogeneous patterns of divergence among loci confirm the porous nature of genomes during speciation with gene flow. Considering that few protein-coding mutations were identified as highly divergent, our results, along with previous transcriptomic studies, imply that changes in regulatory regions most likely had a greater role in the process of whitefish population divergence than protein-coding mutations. This study is the first to demonstrate the efficiency of large-scale targeted resequencing for a nonmodel species with such a large and unsequenced genome.