RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

Genetic and molecular dissection of ginseng (Panax ginseng Mey.) germplasm using high-density genic SNP markers, secondary metabolites, and gene expressions

Genetic and molecular knowledge of a species is crucial to its gene discovery and enhanced breeding. Here, we report the genetic and molecular dissection of ginseng, an important herb for healthy food and medicine. A mini-core collection consisting of 344 cultivars and landraces was developed for ginseng that represents the genetic variation of ginseng existing in its origin and diversity center. We sequenced the transcriptomes of all 344 cultivars and landraces; identified over 1.5 million genic SNPs, thereby revealing the genic diversity of ginseng; and analyzed them with 26,600 high-quality genic SNPs or a selection of them. Ginseng had a wide molecular diversity and was clustered into three subpopulations. Analysis of 16 ginsenosides, the major bioactive components for healthy food and medicine, showed that ginseng had a wide variation in the contents of all 16 ginsenosides and an extensive correlation of their contents, suggesting that they are synthesized through a single or multiple correlated pathways. Furthermore, we pair-wisely examined the relationships between the cultivars and landraces, revealing their relationships in gene expression, gene variation, and ginsenoside biosynthesis. These results provide new knowledge and new genetic and genic resources for advanced research and breeding of ginseng and related species.

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Acute thermal stress elicits interactions between gene expression and alternative splicing in a fish of conservation concern

Transcriptomic research provides a mechanistic understanding of an organism's response to environmental challenges such as increasing temperatures, which can provide key insights into the threats posed by thermal challenges associated with urbanization and climate change. Differential gene expression and alternative splicing are two elements of the transcriptomic stress response that may work in tandem, but relatively few studies have investigated these interactions in fishes of conservation concern. We studied the imperilled redside dace (Clinostomus elongatus) as thermal stress is hypothesised to be an important cause of population declines. We tested the hypothesis that gene expression-splicing interactions contribute to the thermal stress response. Wild fish exposed to acute thermal stress were compared with both handling controls and fish sampled directly from a river. Liver tissue was sampled to study the transcriptomic stress response. With a gene set enrichment analysis, we found that thermally stressed fish showed a transcriptional response related to transcription regulation and responses to unfolded proteins, and alternatively spliced genes related to gene expression regulation and metabolism. One splicing factor, prpf38b, was upregulated in the thermally stressed group compared to the other treatments. This splicing factor may have a role in the Jun/AP-1 cellular stress response, a pathway with wide-ranging and context-dependent effects. Given large gene interaction networks and the context-dependent nature of transcriptional responses, our results highlight the importance of understanding interactions between gene expression and splicing for understanding transcriptomic responses to thermal stress. Our results also reveal transcriptional pathways that can inform conservation breeding, translocation, and reintroduction programs for redside dace and other imperilled species by identifying appropriate source populations.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Little is known about how alternative splicing and other post-transcriptional molecular mechanisms impact plasticity. Steward et al. use transcriptomic and genomic data from the butterfly Bicyclus anynana, finding that splicing plasticity is likely to experience genetic constraints.