Transcriptome analyses and differential gene expression in a non-model fish species with alternative mating tactics

Sensitivity of transcriptomics: Different samples and methodology alter conclusions in Gulf pipefish (Syngnathus scovelli)

Transcriptome analysis has become a central tool in evolutionary and functional genomics. However, variation among biological samples and analysis techniques can greatly influence results, potentially compromising insights into the phenomenon under study. Here, we evaluate differences in the brain transcriptome between female and male Gulf pipefish (Syngnathus scovelli). We perform comparisons between results from entire pipelines for brain transcriptome assembly, quantification, and analysis. We also offer a unique biological comparison between two sampling instances (Redfish Bay: n = 15, Port Lavaca: n = 7). Our results demonstrate crucial shortcomings with current experimental approaches. We found high variation within our results that was driven by both technical differences between pipelines and biological differences between pipefish samples. In our analysis of highly expressed genes, we found that the choice of methods influenced the degree of contamination or noise included in the identified genes. Notably, genes identified within the same pipeline were more similar than any other comparison. Our differential expression analysis revealed that both methodology and sampling location influenced the quantity and consistency of statistically significant transcripts. In the context of these results, we offer modifications to current practices that may increase the robustness of transcriptome-based conclusions. In particular, the use of a reference-guided assembly and an increase in sample sizes are likely to improve resistance to noise or error.

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment-social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Genomics and transcriptomics of the Chinese mitten crabs (Eriocheir sinensis)

To gain a deeper understanding of the genetic factors influencing the growth and development of Eriocheir sinensis, a well-known species of hairy crab found in Yangcheng Lake, this study focused on the de novo genome and full-length transcriptome information of the selected subjects. Specifically, Yangcheng Lake hairy crabs were chosen as the experimental samples. Initially, a genome analysis was performed, resulting in the identification of gene fragments with a combined length of 1266,092,319 bp. Subsequently, a transcriptome analysis was conducted on a mixture of tissues from four different sites, namely muscle, brain, eye, and heart, to further investigate the genetic characteristics at the transcriptome level. The Pacific Biosciences (Pacio) single-molecule real-time sequencing system generated a total of 36.93 G sub-fragments and 175,90041 effective inserts. This research contributes to the indirect comprehension of genetic variations underlying individual traits. Furthermore, a comparison of the obtained data with relevant literature emphasizes the advantages of this study and establishes a basis for further investigations on the Chinese mitten crab.

Evaluating restriction enzyme selection for reduced representation sequencing in conservation genomics

Conservation genomic studies in non-model organisms generally rely on reduced representation sequencing techniques based on restriction enzymes to identify population structure as well as candidate loci for local adaptation. While the expectation is that the reduced representation of the genome is randomly distributed, the proportion of the genome sampled might depend on the GC content of the recognition site of the restriction enzyme used. Here, we evaluated the distribution and functional composition of loci obtained after a reduced representation approach using Genotyping-by-Sequencing (GBS). To do so, we compared experimental data from two endemic fish species (Symphodus ocellatus and Symphodus tinca, EcoT22I enzyme) and two ecosystem engineer sea urchins (Paracentrotus lividus and Arbacia lixula, ApeKI enzyme). In brief, we mapped the sequenced loci to the phylogenetically closest reference genome available (Labrus bergylta in the fish and Strongylocentrotus purpuratus in the sea urchin datasets), classified them as exonic, intronic and intergenic, and studied their function by using Gene Ontology (GO) terms. We also simulated the effect of using both enzymes in the two reference genomes. In both simulated and experimental data, we detected an enrichment towards exonic or intergenic regions depending on the restriction enzyme used and failed to detect differences between total loci and candidate loci for adaptation in the empirical dataset. Most of the functions assigned to the mapped loci were shared between the four species and involved a myriad of general functions. Our results highlight the importance of restriction enzyme selection and the need for high-quality annotated genomes in conservation genomic studies.

Immunoecology of Species with Alternative Reproductive Tactics and Strategies

Alternative reproductive tactics and strategies (ARTS) refer to polymorphic reproductive behaviours in which in addition to the usual two sexes, there are one or more alternative morphs, usually male, that have evolved the ability to circumvent direct intra-sexual competition. Each morph has its own morphological, ecological, developmental, behavioural, life-history, and physiological profile that shifts the balance between reproduction and self-maintenance, one aspect being immunity. Immunoecological work on species with ARTS, which is the topic of this review, is particularly interesting because the alternative morphs make it possible to separate the effects of sex per se from other factors that in other species are inextricably linked with sex. We first summarize the evolution, development, and maintenance of ARTS. We then review immunoecological hypotheses relevant to species with ARTS, dividing them into physiological, life-history, and ecological hypotheses. In context of these hypotheses, we critically review in detail all immunoecological studies we could find on species with ARTS. Several interesting patterns emerge. Oddly, there is a paucity of studies on insects, despite the many benefits that arise from working with insects: larger sample sizes, simple immune systems, and countless forms of alternative reproductive strategies and tactics. Of all the hypotheses considered, the immunocompetence handicap hypothesis has generated the greatest amount of work, but not necessarily the greatest level of understanding. Unfortunately, it is often used as a general guiding principle rather than a source of explicitly articulated predictions. Other hypotheses are usually considered a posteriori, but perhaps they should take centre stage. Whereas blanket concepts such as “immunocompetence” and “androgens” might be useful to develop a rationale, predictions need to be far more explicitly articulated. Integration so far has been a one-way street, with ecologists delving deeper into physiology, sometimes at the cost of ignoring their organisms’ evolutionary history and ecology. One possible useful framework is to divide ecological and evolutionary factors affecting immunity into those that stimulate the immune system, and those that depress it. Finally, the contributions of genomics to ecology are being increasingly recognized and sometimes applied to species with ARTS, but we must ensure that evolutionary and ecological hypotheses drive the effort, as there is no grandeur in the strict reductionist view of life.

Alternative splicing plays key roles in response to stress across different stages of fighting in the fish Betta splendens

Background

Aggression is an evolutionarily conserved behavior critical for animal survival. In the fish Betta splendens, across different stages of fighting interactions, fighting opponents suffer from various stressors, especially from the great demand for oxygen. Using RNA sequencing, we profiled differential alternative splicing (DAS) events in the brains of fish collected before fighting, during fighting, and after fighting to study the involvement of alternative splicing (AS) in the response to stress during the fight.

Results

We found that fighting interactions induced the greatest increase in AS in the ‘during-fighting’ fish, followed by that of the ‘after-fighting’ fish. Intron retention (IR) was the most enriched type among all the basic AS events. DAS genes were mainly associated with synapse assembly, ion transport, and regulation of protein secretion. We further observed that IR events significantly differentiated between winners and losers for 19 genes, which were associated with messenger RNA biogenesis, DNA repair, and transcription machinery. These genes share many common features, including shorter intron length and higher GC content.

Conclusions

This study is the first comprehensive view of AS induced by fighting interactions in a fish species across different stages of those interactions, especially with respect to IR events in winners and losers. Together, these findings facilitate future investigations into transcriptome complexity and AS regulation in response to stress under the context of aggression in vertebrates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08609-2.

Neurogenomic divergence during speciation by reinforcement of mating behaviors in chorus frogs (Pseudacris)

BACKGROUND

Species interactions can promote mating behavior divergence, particularly when these interactions are costly due to maladaptive hybridization. Selection against hybridization can indirectly cause evolution of reproductive isolation within species, a process termed cascade reinforcement. This process can drive incipient speciation by generating divergent selection pressures among populations that interact with different species assemblages. Theoretical and empirical studies indicate that divergent selection on gene expression networks has the potential to increase reproductive isolation among populations. After identifying candidate synaptic transmission genes derived from neurophysiological studies in anurans, we test for divergence of gene expression in a system undergoing cascade reinforcement, the Upland Chorus Frog (Pseudacris feriarum).

RESULTS

Our analyses identified seven candidate synaptic transmission genes that have diverged between ancestral and reinforced populations of P. feriarum, including five that encode synaptic vesicle proteins. Our gene correlation network analyses revealed four genetic modules that have diverged between these populations, two possessing a significant concentration of neurotransmission enrichment terms: one for synaptic membrane components and the other for metabolism of the neurotransmitter nitric oxide. We also ascertained that a greater number of genes have diverged in expression by geography than by sex. Moreover, we found that more genes have diverged within females as compared to males between populations. Conversely, we observed no difference in the number of differentially-expressed genes within the ancestral compared to the reinforced population between the sexes.

CONCLUSIONS

This work is consistent with the idea that divergent selection on mating behaviors via cascade reinforcement contributed to evolution of gene expression in P. feriarum. Although our study design does not allow us to fully rule out the influence of environment and demography, the fact that more genes diverged in females than males points to a role for cascade reinforcement. Our discoveries of divergent candidate genes and gene networks related to neurotransmission support the idea that neural mechanisms of acoustic mating behaviors have diverged between populations, and agree with previous neurophysiological studies in frogs. Increasing support for this hypothesis, however, will require additional experiments under common garden conditions. Our work points to the importance of future replicated and tissue-specific studies to elucidate the relative contribution of gene expression divergence to the evolution of reproductive isolation during incipient speciation.

The tyranny of phylogeny-A plea for a less dogmatic stance on two-species comparisons: Funding bodies, journals and referees discourage two- or few-species comparisons, but such studies provide essential insights complementary to phylogenetic comparative studies

Phylogenetically controlled studies across multiple species correct for taxonomic confounds in physiological performance traits. Therefore, they are preferred over comparisons of two or few closely-related species. Funding bodies, referees and journal editors nowadays often even reject to consider detailed comparisons of two or few closely related species. Here, we plea for a less dogmatic stance on such comparisons, because phylogenetic studies come with their own limitations similar in magnitude as those of two-species comparisons. Two-species comparisons are particularly relevant and instructive for understanding physiological pathways and de novo mutations in three contexts: in a purely mechanistic context, when differences in the regulation of a trait are the focus of investigation, when a physiological trait lacks a direct connection to fitness, and when physiological measures cannot easily be standardized among laboratories. In conclusion, phylogenetic comparative and two-species studies have different strengths and weaknesses and combining these complementary approaches will help integrating biology.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Epigenetic regulation of gonadal and brain aromatase expression in a cichlid fish with environmental sex determination

A fit animal must develop testes or ovaries, with brain and physiology to match. In species with alternative male morphs this coordination of development across tissues operates within sexes as well as between. For Pelvicachromis pulcher, an African cichlid in which early pH exposure influences both sex and alternative male morph, we sequence both copies of aromatase (cyp19a1), a key gene for sex determination. We analyze gene expression and epigenetic state, comparing gonad and brain tissue from females, alternative male morphs, and fry. Relative to brain, we find elevated expression of the A-copy in the ovaries but not testes. Methylation analysis suggests strong epigenetic regulation, with one region specifying sex and another specifying tissue. We find elevated brain expression of the B-copy with no sex or male morph differences. B-copy methylation follows that of the A-copy rather than corresponding to B-copy expression. In 30-day old fry, we see elevated B-copy expression in the head, but we do not see the expected elevated A-copy expression in the trunk that would reflect ovarian development. Interestingly, the A-copy epialleles that distinguish ovaries from testes are among the most explanatory patterns for variation among fry, suggesting epigenetic marking of sex prior to differentiation and thus laying the groundwork for mechanistic studies of epigenetic regulation of sex and morph differentiation.

SEGN: Inferring real-time gene networks mediating phenotypic plasticity

The capacity of an organism to alter its phenotype in response to environmental perturbations changes over developmental time and is a process determined by multiple genes that are co-expressed in intricate but organized networks. Characterizing the spatiotemporal change of such gene networks can offer insight into the genomic signatures underlying organismic adaptation, but it represents a major methodological challenge. Here, we integrate the holistic view of systems biology and the interactive notion of evolutionary game theory to reconstruct so-called systems evolutionary game networks (SEGN) that can autonomously detect, track, and visualize environment-induced gene networks along the time axis. The SEGN overcomes the limitations of traditional approaches by inferring context-specific networks, encapsulating bidirectional, signed, and weighted gene-gene interactions into fully informative networks, and monitoring the process of how networks topologically alter across environmental and developmental cues. Based on the design principle of SEGN, we perform a transcriptional plasticity study by culturing Euphrates poplar, a tree that can grow in the saline desert, in saline-free and saline-stress conditions. SEGN characterize previously unknown gene co-regulation that modulates the time trajectories of the trees' response to salt stress. As a marriage of multiple disciplines, SEGN shows its potential to interpret gene interdependence, predict how transcriptional co-regulation responds to various regimes, and provides a hint for exploring the mass, energetic, or signal basis that drives various types of gene interactions.

Transcriptome sequencing reveals signatures of positive selection in the Spot-Tailed Earless Lizard

The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.

Brain transcriptomics of agonistic behaviour in the weakly electric fish Gymnotus omarorum, a wild teleost model of non-breeding aggression

Differences in social status are often mediated by agonistic encounters between competitors. Robust literature has examined social status-dependent brain gene expression profiles across vertebrates, yet social status and reproductive state are often confounded. It has therefore been challenging to identify the neuromolecular mechanisms underlying social status independent of reproductive state. Weakly electric fish, Gymnotus omarorum, display territorial aggression and social dominance independent of reproductive state. We use wild-derived G. omarorum males to conduct a transcriptomic analysis of non-breeding social dominance relationships. After allowing paired rivals to establish a dominance hierarchy, we profiled the transcriptomes of brain sections containing the preoptic area (region involved in regulating aggressive behaviour) in dominant and subordinate individuals. We identified 16 differentially expressed genes (FDR < 0.05) and numerous genes that co-varied with behavioural traits. We also compared our results with previous reports of differential gene expression in other teleost species. Overall, our study establishes G. omarorum as a powerful model system for understanding the neuromolecular bases of social status independent of reproductive state.

Gill Transcriptome Sequencing and De Novo Annotation of Acanthogobius ommaturus in Response to Salinity Stress

Acanthogobius ommaturus is a euryhaline fish widely distributed in coastal, bay and estuarine areas, showing a strong tolerance to salinity. In order to understand the mechanism of adaptation to salinity stress, RNA-seq was used to compare the transcriptome responses of Acanthogobius ommaturus to the changes of salinity. Four salinity gradients, 0 psu, 15 psu (control), 30 psu and 45 psu were set to conduct the experiment. In total, 131,225 unigenes were obtained from the gill tissue of A. ommaturus using the Illumina HiSeq 2000 platform (San Diego, USA). Compared with the gene expression profile of the control group, 572 differentially expressed genes (DEGs) were screened, with 150 at 0 psu, 170 at 30 psu, and 252 at 45 psu. Additionally, among these DEGs, Gene Ontology (GO) analysis indicated that binding, metabolic processes and cellular processes were significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis detected 3, 5 and 8 pathways related to signal transduction, metabolism, digestive and endocrine systems at 0 psu, 30 psu and 45 psu, respectively. Based on GO enrichment analysis and manual literature searches, the results of the present study indicated that A. ommaturus mainly responded to energy metabolism, ion transport and signal transduction to resist the damage caused by salinity stress. Eight DEGs were randomly selected for further validation by quantitative real-time PCR (qRT-PCR) and the results were consistent with the RNA-seq data.

Identification of sex differences in zebrafish (Danio rerio) brains during early sexual differentiation and masculinization using 17α-methyltestoterone

Sexual behavior in teleost fish is highly plastic. It can be attributed to the relatively few sex differences found in adult brain transcriptomes. Environmental and hormonal factors can influence sex-specific behavior. Androgen treatment stimulates behavioral masculinization. Sex dimorphic gene expression in developing teleost brains and the molecular basis for androgen-induced behavioral masculinization are poorly understood. In this study, juvenile zebrafish (Danio rerio) were treated with 100 ng/L of 17 alpha-methyltestosterone (MT) during sexual development from 20 days post fertilization to 40 days and 60 days post fertilization. We compared brain gene expression patterns in MT-treated zebrafish with control males and females using RNA-Seq to shed light on the dynamic changes in brain gene expression during sexual development and how androgens affect brain gene expression leading to behavior masculinization. We found modest differences in gene expression between juvenile male and female zebrafish brains. Brain aromatase (cyp19a1b), prostaglandin 3a synthase (ptges3a), and prostaglandin reductase 1 (ptgr1) were among the genes with sexually dimorphic expression patterns. MT treatment significantly altered gene expression relative to both male and female brains. Fewer differences were found among MT-treated brains and male brains compared to female brains, particularly at 60 dpf. MT treatment upregulated the expression of hydroxysteroid 11-beta dehydrogenase 2 (hsd11b2), deiodinase, iodothyronine, type II (dio2), and gonadotrophin releasing hormones (GnRH) 2 and 3 (gnrh2 and gnrh3) suggesting local synthesis of 11-ketotestosterone, triiodothyronine, and GnRHs in zebrafish brains which are influenced by androgens. Androgen, estrogen, prostaglandin, thyroid hormone, and GnRH signaling pathways likely interact to modulate teleost sexual behavior.

A novel integrative approach elucidates fine-scale dispersal patchiness in marine populations

Dispersal is one of the main determining factors of population structure. In the marine habitat, well-connected populations with large numbers of reproducing individuals are common but even so population structure can exist on a small-scale. Variation in dispersal patterns between populations or over time is often associated to geographic distance or changing oceanographic barriers. Consequently, detecting structure and variation in dispersal on a fine-scale within marine populations still remains a challenge. Here we propose and use a novel approach of combining a clustering model, early-life history trait information from fish otoliths, spatial coordinates and genetic markers to detect very fine-scale dispersal patterns. We collected 1573 individuals (946 adults and 627 juveniles) of the black-faced blenny across a small-scale (2 km) coastline as well as at a larger-scale area (<50 kms). A total of 178 single nucleotide polymorphism markers were used to evaluate relatedness patterns within this well-connected population. In our clustering models we categorized SHORT-range dispersers to be potential local recruits based on their high relatedness within and low relatedness towards other spatial clusters. Local retention and/or dispersal of this potential local recruitment varied across the 2 km coastline with higher frequency of SHORT-range dispersers towards the southwest of the area for adults. An inverse pattern was found for juveniles, showing an increase of SHORT-range dispersers towards the northeast. As we rule out selective movement and mortality from one year to the next, this pattern reveals a complex but not full genetic mixing, and variability in coastal circulation is most likely the main driver of this fine-scale chaotic genetic patchiness within this otherwise homogeneous population. When focusing on the patterns within one recruitment season, we found large differences in temperatures (from approx. 17 °C to 25 °C) as well as pelagic larval duration (PLD) for juveniles from the beginning of the season and the end of the season. We were able to detect fine-scale differences in LONG-range juvenile dispersers, representing distant migrants, depending on whether they were born at the beginning of the season with a longer PLD, or at the end of the reproductive season. The ability to detect such fine-scale dispersal patchiness will aid in our understanding of the underlying mechanisms of population structuring and chaotic patchiness in a wide range of species even with high potential dispersal abilities.

CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features

With advances in next-generation sequencing technologies, numerous novel transcripts in a large number of organisms have been identified. With the goal of fast, accurate assessment of the coding ability of RNA transcripts, we upgraded the coding potential calculator CPC1 to CPC2. CPC2 runs ∼1000 times faster than CPC1 and exhibits superior accuracy compared with CPC1, especially for long non-coding transcripts. Moreover, the model of CPC2 is species-neutral, making it feasible for ever-growing non-model organism transcriptomes. A mobile-friendly web server, as well as a downloadable standalone package, is freely available at http://cpc2.cbi.pku.edu.cn.

DNA Methylation Patterns in the Round Goby Hypothalamus Support an On-The-Spot Decision Scenario for Territorial Behavior

The question as to how early life experiences are stored on a molecular level and affect traits later in life is highly topical in ecology, medicine, and epigenetics. In this study, we use a fish model to investigate whether DNA methylation mediates early life experiences and predetermines a territorial male reproductive phenotype. In fish, adult reproductive phenotypes frequently depend on previous life experiences and are often associated with distinct morphological traits. DNA methylation is an epigenetic mechanism which is both sensitive to environmental conditions and stably inherited across cell divisions. We therefore investigate early life predisposition in the round goby Neogobius melanostomus by growth back-calculations and then study DNA methylation by MBD-Seq in the brain region controlling vertebrate reproductive behavior, the hypothalamus. We find a link between the territorial reproductive phenotype and high growth rates in the first year of life. However, hypothalamic DNA methylation patterns reflect the current behavioral status independently of early life experiences. Together, our data suggest a non-predetermination scenario in the round goby, in which indeterminate males progress to a non-territorial status in the spawning season, and in which some males then assume a specialized territorial phenotype if current conditions are favorable.

Environmental and Evolutionary Drivers of the Modular Gene Regulatory Network Underlying Phenotypic Plasticity for Stress Resistance in the Nematode Caenorhabditis remanei

Organisms can cope with stressful environments via a combination of phenotypic plasticity at the individual level and adaptation at the population level. Changes in gene expression can play an important role in both. Significant advances in our understanding of gene regulatory plasticity and evolution have come from comparative studies in the field and laboratory. Experimental evolution provides another powerful path by which to learn about how differential regulation of genes and pathways contributes to both acclimation and adaptation. Here we present results from one such study using the nematode Caenorhabditis remanei We selected one set of lines to withstand heat stress and another oxidative stress. We then compared transcriptional responses to acute heat stress of both and an unselected control to the ancestral population using a weighted gene coexpression network analysis, finding that the transcriptional response is primarily dominated by a plastic response that is shared in the ancestor and all evolved populations. In addition, we identified several modules that respond to artificial selection by (1) changing the baseline level of expression, (2) altering the magnitude of the plastic response, or (3) a combination of the two. Our findings therefore reveal that while patterns of transcriptional response can be perturbed with short bouts of intense selection, the overall ancestral structure of transcriptional plasticity is largely maintained over time.

Female Mimicry by Sneaker Males Has a Transcriptomic Signature in Both the Brain and the Gonad in a Sex-Changing Fish

Phenotypic plasticity represents an elegant adaptive response of individuals to a change in their environment. Bluehead wrasses (Thalassoma bifasciatum) exhibit astonishing sexual plasticity, including female-to-male sex change and discrete male morphs that differ strikingly in behavior, morphology, and gonadal investment. Using RNA-seq transcriptome profiling, we examined the genes and physiological pathways underlying flexible behavioral and gonadal differences among female, dominant (bourgeois) male, and female-mimic (sneaker) male blueheads. For the first time in any organism, we find that female mimicry by sneaker males has a transcriptional signature in both the brain and the gonad. Sneaker males shared striking similarity in neural gene expression with females, supporting the idea that males with alternative reproductive phenotypes have "female-like brains." Sneaker males also overexpressed neuroplasticity genes, suggesting that their opportunistic reproductive strategy requires a heightened capacity for neuroplasticity. Bourgeois males overexpressed genes associated with socio-sexual behaviors (e.g., isotocin), but also neuroprotective genes and biomarkers of oxidative stress and aging, indicating a hitherto unexplored cost to these males of attaining the reproductively privileged position at the top of the social hierarchy. Our novel comparison of testicular transcriptomes in a fish with male sexual polymorphism associates greater gonadal investment by sneaker males with overexpression of genes involved in cell proliferation and sperm quality control. We propose that morphological female-mimicry by sneaker male teleosts entails pervasive downregulation of androgenesis genes, consistent with low androgen production in males lacking well-developed secondary sexual characters.

Behavioural tactic predicts preoptic-hypothalamic gene expression more strongly than developmental morph in fish with alternative reproductive tactics

Reproductive success relies on the coordination of social behaviours, such as territory defence, courtship and mating. Species with extreme variation in reproductive tactics are useful models for identifying the neural mechanisms underlying social behaviour plasticity. The plainfin midshipman (Porichthys notatus) is a teleost fish with two male reproductive morphs that follow widely divergent developmental trajectories and display alternative reproductive tactics (ARTs). Type I males defend territories, court females and provide paternal care, but will resort to cuckoldry if they cannot maintain a territory. Type II males reproduce only through cuckoldry. We sought to disentangle gene expression patterns underlying behavioural tactic, in this case ARTs, from those solely reflective of developmental morph. Using RNA-sequencing, we investigated differential transcript expression in the preoptic area-anterior hypothalamus (POA-AH) of courting type I males, cuckolding type I males and cuckolding type II males. Unexpectedly, POA-AH differential expression was more strongly coupled to behavioural tactic than morph. This included a suite of transcripts implicated in hormonal regulation of vertebrate social behaviour. Our results reveal that divergent expression patterns in a conserved neuroendocrine centre known to regulate social-reproductive behaviours across vertebrate lineages may be uncoupled from developmental history to enable plasticity in the performance of reproductive tactics.

Full-length transcriptome of Misgurnus anguillicaudatus provides insights into evolution of genus Misgurnus

Reconstruction and annotation of transcripts, particularly for a species without reference genome, plays a critical role in gene discovery, investigation of genomic signatures, and genome annotation in the pre-genomic era. This study generated 33,330 full-length transcripts of diploid M. anguillicaudatus using PacBio SMRT Sequencing. A total of 6,918 gene families were identified with two or more isoforms, and 26,683 complete ORFs with an average length of 1,497 bp were detected. Totally, 1,208 high-confidence lncRNAs were identified, and most of these appeared to be precursor transcripts of miRNAs or snoRNAs. Phylogenetic tree of the Misgurnus species was inferred based on the 1,905 single copy orthologous genes. The tetraploid and diploid M. anguillicaudatus grouped into a clade, and M. bipartitus showed a closer relationship with the M. anguillicaudatus. The overall evolutionary rates of tetraploid M. anguillicaudatus were significantly higher than those of other Misgurnus species. Meanwhile, 28 positively selected genes were identified in M. anguillicaudatus clade. These positively selected genes may play critical roles in the adaptation to various habitat environments for M. anguillicaudatus. This study could facilitate further exploration of the genomic signatures of M. anguillicaudatus and provide potential insights into unveiling the evolutionary history of tetraploid loach.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Transcriptome response of Atlantic salmon (Salmo salar) to competition with ecologically similar non-native species

Non-native species may be introduced either intentionally or unintentionally, and their impact can range from benign to highly disruptive. Non-native salmonids were introduced into Lake Ontario, Canada, to provide recreational fishing opportunities; however, the establishment of those species has been proposed as a significant barrier to the reintroduction of native Atlantic salmon (Salmo salar) due to intense interspecific competition. In this study, we compared population differences of Atlantic salmon in transcriptome response to interspecific competition. We reared Atlantic salmon from two populations (LaHave River and Sebago Lake) with fish of each of three non-native salmonids (Chinook salmon Oncorhynchus tshawytscha, rainbow trout O. mykiss, and brown trout S. trutta) in artificial streams. We used RNA-seq to assess transcriptome differences between the Atlantic salmon populations and the responses of these populations to the interspecific competition treatments after 10 months of competition in the stream tanks. We found that population differences in gene expression were generally greater than the effects of interspecific competition. Interestingly, we found that the two Atlantic salmon populations exhibited similar responses to interspecific competition based on functional gene ontologies, but the specific genes within those ontologies were different. Our transcriptome analyses suggest that the most stressful competitor (as measured by the highest number of differentially expressed genes) differs between the two study populations. Our transcriptome characterization highlights the importance of source population selection for conservation applications, as organisms with different evolutionary histories can possess different transcriptional responses to the same biotic stressors. The results also indicate that generalized predictions of the response of native species to interactions with introduced species may not be appropriate without incorporating potential population-specific response to introduced species.

Temporal variation in brain transcriptome is associated with the expression of female mimicry as a sequential male alternative reproductive tactic in fish

Distinct patterns of gene expression often underlie intra- and intersexual differences, and the study of this set of coregulated genes is essential to understand the emergence of complex behavioural phenotypes. Here, we describe the development of a de novo transcriptome and brain gene expression profiles of wild-caught peacock blenny, Salaria pavo, an intertidal fish with sex-role reversal in courtship behaviour (i.e., females are the courting sex) and sequential alternative reproductive tactics in males (i.e., larger and older nest-holder males and smaller and younger sneaker males occur). Sneakers mimic both female's courtship behaviour and nuptial coloration to get access to nests and sneak fertilizations, and later in life transition into nest-holder males. Thus, this species offers the unique opportunity to study how the regulation of gene expression can contribute to intersex phenotypes and to the sequential expression of male and female behavioural phenotypes by the same individual. We found that at the whole brain level, expression of the sneaker tactic was paralleled by broader and divergent gene expression when compared to either females or nest-holder males, which were more similar between themselves. When looking at sex-biased transcripts, sneaker males are intersex rather than being either nest-holder or female-like, and their transcriptome is simultaneously demasculinized for nest-holder-biased transcripts and feminized for female-biased transcripts. These results indicate that evolutionary changes in reproductive plasticity can be achieved through regulation of gene expression, and in particular by varying the magnitude of expression of sex-biased genes, throughout the lifetime of the same individual.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Comparative transcriptomic analysis of the evolution and development of flower size in Saltugilia (Polemoniaceae)

BACKGROUND

Flower size varies dramatically across angiosperms, representing innovations over the course of >130 million years of evolution and contributing substantially to relationships with pollinators. However, the genetic underpinning of flower size is not well understood. Saltugilia (Polemoniaceae) provides an excellent non-model system for extending the genetic study of flower size to interspecific differences that coincide with variation in pollinators.

RESULTS

Using targeted gene capture methods, we infer phylogenetic relationships among all members of Saltugilia to provide a framework for investigating the genetic control of flower size differences via RNA-Seq de novo assembly. Nuclear concatenation and species tree inference methods provide congruent topologies. The inferred evolutionary trajectory of flower size is from small flowers to larger flowers. We identified 4 to 10,368 transcripts that are differentially expressed during flower development, with many unigenes associated with cell wall modification and components of the auxin and gibberellin pathways.

CONCLUSIONS

Saltugilia is an excellent model for investigating covarying floral and pollinator evolution. Four candidate genes from model systems (BIG BROTHER, BIG PETAL, GASA, and LONGIFOLIA) show differential expression during development of flowers in Saltugilia, and four other genes (FLOWERING-PROMOTING FACTOR 1, PECTINESTERASE, POLYGALACTURONASE, and SUCROSE SYNTHASE) fit into hypothesized organ size pathways. Together, these gene sets provide a strong foundation for future functional studies to determine their roles in specifying interspecific differences in flower size.

Differential gene expression in narrow- and broad-headed European glass eels (Anguilla anguilla) points to a transcriptomic link of head shape dimorphism with growth rate and chemotaxis

One of the major challenges in evolutionary biology is to understand the mechanisms underlying morphological dimorphism and plasticity, including the genomic basis of traits and links to ecology. At the yellow eel stage of the European eel (Anguilla anguilla), two morphotypes are found: broad- and narrow-heads. This dimorphism has been linked to dietary differences, with broad-heads feeding on harder, larger prey than narrow-heads. However, recent research showed that both morphotypes could be distinguished at the glass eel stage, the nonfeeding predecessor of the yellow eel stage, implying that nondietary factors play a role in the development of this head shape dimorphism. Here, we used transcriptome profiling (RNAseq) to identify differentially expressed genes between broad- and narrow-headed glass eels. We found 260 significantly differentially expressed genes between the morphotypes, of which most were related to defence and immune responses. Interestingly, two genes involved in growth (soma and igf2) were significantly upregulated in narrow-heads, while nine genes involved in chemotaxis showed significant differential expression. Thus, we found support for the observation that head shape is associated with somatic growth, with fast-growing eels developing a narrower head. Additionally, observations in the wild have shown that slow-growers prefer freshwater, while fast-growers prefer brackish water. The differential expression of genes involved in chemotaxis seems to indicate that glass eel growth rate and habitat choice are linked. We hypothesize that two levels of segregation could take place in the European eel: first according to habitat choice and second according to feeding preference.

Increased Testosterone Decreases Medial Cortical Volume and Neurogenesis in Territorial Side-Blotched Lizards (Uta stansburiana)

Variation in an animal's spatial environment can induce variation in the hippocampus, an area of the brain involved in spatial cognitive processing. Specifically, increased spatial area use is correlated with increased hippocampal attributes, such as volume and neurogenesis. In the side-blotched lizard (Uta stansburiana), males demonstrate alternative reproductive tactics and are either territorial—defending large, clearly defined spatial boundaries—or non-territorial—traversing home ranges that are smaller than the territorial males' territories. Our previous work demonstrated cortical volume (reptilian hippocampal homolog) correlates with these spatial niches. We found that territorial holders have larger medial cortices than non-territory holders, yet these differences in the neural architecture demonstrated some degree of plasticity as well. Although we have demonstrated a link among territoriality, spatial use, and brain plasticity, the mechanisms that underlie this relationship are unclear. Previous studies found that higher testosterone levels can induce increased use of the spatial area and can cause an upregulation in hippocampal attributes. Thus, testosterone may be the mechanistic link between spatial area use and the brain. What remains unclear, however, is if testosterone can affect the cortices independent of spatial experiences and whether testosterone differentially interacts with territorial status to produce the resultant cortical phenotype. In this study, we compared neurogenesis as measured by the total number of doublecortin-positive cells and cortical volume between territorial and non-territorial males supplemented with testosterone. We found no significant differences in the number of doublecortin-positive cells or cortical volume among control territorial, control non-territorial, and testosterone-supplemented non-territorial males, while testosterone-supplemented territorial males had smaller medial cortices containing fewer doublecortin-positive cells. These results demonstrate that testosterone can modulate medial cortical attributes outside of differential spatial processing experiences but that territorial males appear to be more sensitive to alterations in testosterone levels compared with non-territorial males.

De novo transcriptome sequencing and analysis of freshwater snail (Radix balthica) to discover genes and pathways affected by exposure to oxazepam

Pharmaceuticals are increasingly found in aquatic ecosystems due to the non-efficiency of waste water treatment plants. Therefore, aquatic organisms are frequently exposed to a broad diversity of pharmaceuticals. Freshwater snail Radix balthica has been chosen as model to study the effects of oxazepam (psychotropic drug) on developmental stages ranging from trochophore to hatching. In order to provide a global insight of these effects, a transcriptome deep sequencing has been performed on exposed embryos. Eighteen libraries were sequenced, six libraries for three conditions: control, exposed to the lowest oxazepam concentration with a phenotypic effect (delayed hatching) (TA) and exposed to oxazepam concentration found in freshwater (TB). A total of 39,759,772 filtered raw reads were assembled into 56,435 contigs having a mean length of 1579.68 bp and mean depth of 378.96 reads. 44.91% of the contigs have at least one annotation. The differential expression analysis between the control condition and the two exposure conditions revealed 146 contigs differentially expressed of which 144 for TA and two for TB. 34.0% were annotated with biological function. There were four mainly impacted processes: two cellular signalling systems (Notch and JNK) and two biosynthesis pathways (Polyamine and Catecholamine pathways). This work reports a large-scale analysis of differentially transcribed genes of R. balthica exposed to oxazepam during egg development until hatching. In addition, these results enriched the de novo database of potential ecotoxicological models.

Brain Transcriptional Profiles of Male Alternative Reproductive Tactics and Females in Bluegill Sunfish

Bluegill sunfish (Lepomis macrochirus) are one of the classic systems for studying male alternative reproductive tactics (ARTs) in teleost fishes. In this species, there are two distinct life histories: parental and cuckolder, encompassing three reproductive tactics, parental, satellite, and sneaker. The parental life history is fixed, whereas individuals who enter the cuckolder life history transition from sneaker to satellite tactic as they grow. For this study, we used RNAseq to characterize the brain transcriptome of the three male tactics and females during spawning to identify gene ontology (GO) categories and potential candidate genes associated with each tactic. We found that sneaker males had higher levels of gene expression differentiation compared to the other two male tactics. Sneaker males also had higher expression in ionotropic glutamate receptor genes, specifically AMPA receptors, compared to other males, which may be important for increased spatial working memory while attempting to cuckold parental males at their nests. Larger differences in gene expression also occurred among male tactics than between males and females. We found significant expression differences in several candidate genes that were previously identified in other species with ARTs and suggest a previously undescribed role for cAMP-responsive element modulator (crem) in influencing parental male behaviors during spawning.

De novo transcriptome assembly for the lobster Homarus americanus and characterization of differential gene expression across nervous system tissues

BACKGROUND

The American lobster, Homarus americanus, is an important species as an economically valuable fishery, a key member in marine ecosystems, and a well-studied model for central pattern generation, the neural networks that control rhythmic motor patterns. Despite multi-faceted scientific interest in this species, currently our genetic resources for the lobster are limited. In this study, we de novo assemble a transcriptome for Homarus americanus using central nervous system (CNS), muscle, and hybrid neurosecretory tissues and compare gene expression across these tissue types. In particular, we focus our analysis on genes relevant to central pattern generation and the identity of the neurons in a neural network, which is defined by combinations of genes distinguishing the neuronal behavior and phenotype, including ion channels, neurotransmitters, neuromodulators, receptors, transcription factors, and other gene products.

RESULTS

Using samples from the central nervous system (brain, abdominal ganglia), abdominal muscle, and heart (cardiac ganglia, pericardial organs, muscle), we used RNA-Seq to characterize gene expression patterns across tissues types. We also compared control tissues with those challenged with the neuropeptide proctolin in vivo. Our transcriptome generated 34,813 transcripts with known protein annotations. Of these, 5,000-10,000 of annotated transcripts were significantly differentially expressed (DE) across tissue types. We found 421 transcripts for ion channels and identified receptors and/or proteins for over 20 different neurotransmitters and neuromodulators. Results indicated tissue-specific expression of select neuromodulator (allostatin, myomodulin, octopamine, nitric oxide) and neurotransmitter (glutamate, acetylcholine) pathways. We also identify differential expression of ion channel families, including kainite family glutamate receptors, inward-rectifying K(+) (IRK) channels, and transient receptor potential (TRP) A family channels, across central pattern generating tissues.

CONCLUSIONS

Our transcriptome-wide profiles of the rhythmic pattern generating abdominal and cardiac nervous systems in Homarus americanus reveal candidates for neuronal features that drive the production of motor output in these systems.

Natural variation in brain gene expression profiles of aggressive and nonaggressive individual sticklebacks

Within many species, some individuals are consistently more aggressive than others. We examine whether there are differences in brain gene expression between aggressive versus nonaggressive behavioural types of individuals within a natural population of male three-spined sticklebacks (Gasterosteus aculeatus). We compared gene expression profiles of aggressive male sticklebacks to nonaggressive males in four regions of the brain (brainstem, cerebellum, diencephalon and telencephalon). Relatively few genes were differentially expressed between behavioural types in telencephalon, cerebellum and diencephalon, but hundreds of genes were differentially expressed in brainstem, a brain area involved in detecting threats. Six genes that were differentially expressed in response to a territorial intrusion in a previous study were also differentially expressed between behavioural types in this study, implying primarily non-shared but some shared molecular mechanisms. Our findings offer new insights into the molecular causes and correlates of behavioural plasticity and individual variation in behaviour.

In Silico Comparative Transcriptome Analysis of Two Color Morphs of the Common Coral Trout (Plectropomus Leopardus)

The common coral trout is one species of major importance in commercial fisheries and aquaculture. Recently, two different color morphs of Plectropomus leopardus were discovered and the biological importance of the color difference is unknown. Since coral trout species are poorly characterized at the molecular level, we undertook the transcriptomic characterization of the two color morphs, one black and one red coral trout, using Illumina next generation sequencing technologies. The study produced 55162966 and 54588952 paired-end reads, for black and red trout, respectively. De novo transcriptome assembly generated 95367 and 99424 unique sequences in black and red trout, respectively, with 88813 sequences shared between them. Approximately 50% of both trancriptomes were functionally annotated by BLAST searches against protein databases. The two trancriptomes were enriched into 25 functional categories and showed similar profiles of Gene Ontology category compositions. 34110 unigenes were grouped into 259 KEGG pathways. Moreover, we identified 14649 simple sequence repeats (SSRs) and designed primers for potential application. We also discovered 130524 putative single nucleotide polymorphisms (SNPs) in the two transcriptomes, supplying potential genomic resources for the coral trout species. In addition, we identified 936 fast-evolving genes and 165 candidate genes under positive selection between the two color morphs. Finally, 38 candidate genes underlying the mechanism of color and pigmentation were also isolated. This study presents the first transcriptome resources for the common coral trout and provides basic information for the development of genomic tools for the identification, conservation, and understanding of the speciation and local adaptation of coral reef fish species.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Optimizing and benchmarking de novo transcriptome sequencing: from library preparation to assembly evaluation

Background

RNA-seq enables gene expression profiling in selected spatiotemporal windows and yields massive sequence information with relatively low cost and time investment, even for non-model species. However, there remains a large room for optimizing its workflow, in order to take full advantage of continuously developing sequencing capacity.

Method

Transcriptome sequencing for three embryonic stages of Madagascar ground gecko (Paroedura picta) was performed with the Illumina platform. The output reads were assembled de novo for reconstructing transcript sequences. In order to evaluate the completeness of transcriptome assemblies, we prepared a reference gene set consisting of vertebrate one-to-one orthologs.

Result

To take advantage of increased read length of >150 nt, we demonstrated shortened RNA fragmentation time, which resulted in a dramatic shift of insert size distribution. To evaluate products of multiple de novo assembly runs incorporating reads with different RNA sources, read lengths, and insert sizes, we introduce a new reference gene set, core vertebrate genes (CVG), consisting of 233 genes that are shared as one-to-one orthologs by all vertebrate genomes examined (29 species)., The completeness assessment performed by the computational pipelines CEGMA and BUSCO referring to CVG, demonstrated higher accuracy and resolution than with the gene set previously established for this purpose. As a result of the assessment with CVG, we have derived the most comprehensive transcript sequence set of the Madagascar ground gecko by means of assembling individual libraries followed by clustering the assembled sequences based on their overall similarities.

Conclusion

Our results provide several insights into optimizing de novo RNA-seq workflow, including the coordination between library insert size and read length, which manifested in improved connectivity of assemblies. The approach and assembly assessment with CVG demonstrated here would be applicable to transcriptome analysis of other species as well as whole genome analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2007-1) contains supplementary material, which is available to authorized users.