Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

Transcriptomic and Histological Analysis of the Greentail Prawn (Metapenaeus bennettae) Following Light Crude Oil Exposure

Oil spills pose a significant threat to marine biodiversity. Crude oil can partition into sediments where it may be persistent, placing benthic species such as decapods at particular risk of exposure. Transcriptomic and histological tools are often used to investigate the effects of hydrocarbon exposure on marine organisms following oil spill events, allowing for the identification of metabolic pathways impacted by oil exposure. However, there is limited information available for decapod crustaceans, many of which carry significant economic value. In the present study, we assess the sublethal impacts of crude oil exposure in the commercially important Australian greentail prawn (Metapenaeus bennettae) using transcriptomic and histological analyses. Prawns exposed to light, unweathered crude oil "spiked" sediments for 90 h were transferred to clean sediments for a further 72 h to assess recovery. Chemical analyses indicated that polycyclic aromatic hydrocarbons increased by approximately 65% and 91% in prawn muscle following 24 and 90 h of exposure, respectively, and significantly decreased during 24- and 72-h recovery periods. Transcriptomic responses followed an exposure and recovery pattern with innate immunity and nutrient metabolism transcripts significantly lowered in abundance after 24 h of exposure and were higher in abundance after 72 h of recovery. In addition, transcription/translation, cellular responses, and DNA repair pathways were significantly impacted after 24 h of exposure and recovered after 72 h of recovery. However, histological alterations such as tubule atrophy indicated an increase in severity after 24 and 72 h of recovery. The present study provides new insights into the sublethal impacts of crude oil exposure in greentail prawns and identifies molecular pathways altered by exposure. We expect these findings to inform future management associated with oil extraction activity and spills. Environ Toxicol Chem 2022;41:2162-2180. © 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Salivary and Intestinal Transcriptomes Reveal Differential Gene Expression in Starving, Fed and Trypanosoma cruzi-Infected Rhodnius neglectus

Rhodnius neglectus is a potential vector of Trypanosoma cruzi (Tc), the causative agent of Chagas disease. The salivary glands (SGs) and intestine (INT) are actively required during blood feeding. The saliva from SGs is injected into the vertebrate host, modulating immune responses and favoring feeding for INT digestion. Tc infection significantly alters the physiology of these tissues; however, studies that assess this are still scarce. This study aimed to gain a better understanding of the global transcriptional expression of genes in SGs and INT during fasting (FA), fed (FE), and fed in the presence of Tc (FE + Tc) conditions. In FA, the expression of transcripts related to homeostasis maintenance proteins during periods of stress was predominant. Therefore, the transcript levels of Tret1-like and Hsp70Ba proteins were increased. Blood appeared to be responsible for alterations found in the FE group, as most of the expressed transcripts, such as proteases and cathepsin D, were related to digestion. In FE + Tc group, there was a decreased expression of blood processing genes for insect metabolism (e.g., Antigen-5 precursor, Pr13a, and Obp), detoxification (Sult1) in INT and acid phosphatases in SG. We also found decreased transcriptional expression of lipocalins and nitrophorins in SG and two new proteins, pacifastin and diptericin, in INT. Several transcripts of unknown proteins with investigative potential were found in both tissues. Our results also show that the presence of Tc can change the expression in both tissues for a long or short period of time. While SG homeostasis seems to be re-established on day 9, changes in INT are still evident. The findings of this study may be used for future research on parasite-vector interactions and contribute to the understanding of food physiology and post-meal/infection in triatomines.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Transcriptional Sequencing and Gene Expression Analysis of Various Genes in Fruit Development of Three Different Black Pepper (Piper nigrum L.) Varieties

Black pepper (Piper nigrum) is a vital spice crop with uses ranging from culinary to pharmacological applications. However, limited genetic information has constrained the understanding of the molecular regulation of flower and fruit development in black pepper. In this study, a comparison among three different black pepper varieties, Semengok Aman (SA), Kuching (KC), and Semengok 1 (S1), with varying fruit characteristics was used to provide insight on the genetic regulation of flower and fruit development. Next-generation sequencing (NGS) technology was used to determine the flower and fruit transcriptomes by sequencing on an Illumina HiSeq 2500 platform followed by de novo assembly using SOAPdenovo-Trans. The high-quality assembly of 66,906 of unigenes included 64.4% of gene sequences (43,115) with similarity to one or more protein sequences from the GenBank database. Annotation with Blast2Go assigned 37,377 genes to one or more Gene Ontology terms. Of these genes, 5,874 genes were further associated with the biological pathways recorded in the KEGG database. Comparison of flower and fruit transcriptome data from the three different black pepper varieties revealed a large number of DEGs between flower and fruit of the SA variety. Gene Ontology (GO) enrichment analysis further supports functions of DEGs between flower and fruit in the categories of carbohydrate metabolic processes, embryo development, and DNA metabolic processes while the DEGs in fruit relate to biosynthetic process, secondary metabolic process, and catabolic process. The enrichment of DEGs in KEGG pathways was also investigated, and a large number of genes were found to belong to the nucleotide metabolism and carbohydrate metabolism categories. Gene expression profiling of flower formation-related genes reveals that other than regulating the flowering in black pepper, the flowering genes might also be implicated in the fruit development process. Transcriptional analysis of sugar transporter and carbohydrate metabolism genes in different fruit varieties suggested that the carbohydrate metabolism in black pepper fruit is developmentally regulated, and some genes might serve as potential genes for future crop quality improvement. Study on the piperine-related gene expression analysis suggested that lysine-derived products might present in all stages of fruit development, but the transportation was only active at the early stage of fruit development. These results indicate several candidate genes related to the development of flower and fruit in black pepper and provide a resource for future functional analysis and potentially for future crop improvement.

Transcriptomic responses to predator kairomones in embryos of the aquatic snail Radix balthica

The ability of organisms to respond to predation threat by exhibiting induced defenses is well documented, but studies on the potential mechanistic basis for such responses are scarce. Here, we examine the transcriptomic response to predator kairomones of two functionally distinct developmental stages in embryos of the aquatic snail Radix balthica: E8-the stage at which a range-finding trial indicated that kairomone-induced accelerated growth and development first occurred; and E9-the stage at which embryos switched from ciliary- to crawling-driven locomotion. We tested whether expression profiles were influenced by kairomones and whether this influence varied between stages. We also identified potential candidate genes for investigating mechanisms underpinning induced responses. There were 6,741 differentially expressed transcripts between developmental stages, compared to just five in response to predator kairomones. However, on examination of functional enrichment in the transcripts responding to predator kairomones and adopting a less stringent significance threshold, 206 transcripts were identified relating to muscle function, growth, and development, with this response being greater at the later E9 stage. Furthermore, these transcripts included putative annotations for genes identified as responding to predator kairomones in other taxa, including C1q, lectin, and actin domains. Globally, transcript expression appeared reduced in response to predator kairomones and we hypothesize that this might be a result of metabolic suppression, as has been reported in other taxa in response to predation threat.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BACKGROUND

Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.

RESULTS

Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.

CONCLUSIONS

Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Identification of Putative Precursor Genes for the Biosynthesis of Cannabinoid-Like Compound in Radula marginata

The liverwort Radula marginata belongs to the bryophyte division of land plants and is a prospective alternate source of cannabinoid-like compounds. However, mechanistic insights into the molecular pathways directing the synthesis of these cannabinoid-like compounds have been hindered due to the lack of genetic information. This prompted us to do deep sequencing, de novo assembly and annotation of R. marginata transcriptome, which resulted in the identification and validation of the genes for cannabinoid biosynthetic pathway. In total, we have identified 11,421 putative genes encoding 1,554 enzymes from 145 biosynthetic pathways. Interestingly, we have identified all the upstream genes of the central precursor of cannabinoid biosynthesis, cannabigerolic acid (CBGA), including its two first intermediates, stilbene acid (SA) and geranyl diphosphate (GPP). Expression of all these genes was validated using quantitative real-time PCR. We have characterized the protein structure of stilbene synthase (STS), which is considered as a homolog of olivetolic acid in R. marginata. Moreover, the metabolomics approach enabled us to identify CBGA-analogous compounds using electrospray ionization mass spectrometry (ESI-MS/MS) and gas chromatography mass spectrometry (GC-MS). Transcriptomic analysis revealed 1085 transcription factors (TF) from 39 families. Comparative analysis showed that six TF families have been uniquely predicted in R. marginata. In addition, the bioinformatics analysis predicted a large number of simple sequence repeats (SSRs) and non-coding RNAs (ncRNAs). Our results collectively provide mechanistic insights into the putative precursor genes for the biosynthesis of cannabinoid-like compounds and a novel transcriptomic resource for R. marginata. The large-scale transcriptomic resource generated in this study would further serve as a reference transcriptome to explore the Radulaceae family.

Characterization of a de novo assembled transcriptome of the Common Blackbird (Turdus merula)

BACKGROUND

In recent years, next generation high throughput sequencing technologies have proven to be useful tools for investigations concerning the genomics or transcriptomics also of non-model species. Consequently, ornithologists have adopted these technologies and the respective bioinformatics tools to survey the genomes and transcriptomes of a few avian non-model species. The Common Blackbird is one of the most common bird species living in European cities, which has successfully colonized urban areas and for which no reference genome or transcriptome is publicly available. However, to target questions like genome wide gene expression analysis, a reference genome or transcriptome is needed.

METHODS

Therefore, in this study two Common Blackbirds were sacrificed, their mRNA was isolated and analyzed by RNA-Seq to de novo assemble a transcriptome and characterize it. Illumina reads (125 bp paired-end) and a Velvet/Oases pipeline led to 162,158 transcripts. For the annotation (using Blast+), an unfiltered protein database was used. SNPs were identified using SAMtools and BCFtools. Furthermore, mRNA from three single tissues (brain, heart and liver) of the same two Common Blackbirds were sequenced by Illumina (75 bp single-end reads). The draft transcriptome and the three single tissues were compared by their BLAST hits with the package VennDiagram in R.

RESULTS

Following the annotation against protein databases, we found evidence for 15,580 genes in the transcriptome (all well characterized hits after annotation). On 18% of the assembled transcripts, 144,742 SNPs were identified which are, consequently, 0.09% of all nucleotides in the assembled transcriptome. In the transcriptome and in the single tissues (brain, heart and liver), 10,182 shared genes were found.

DISCUSSION

Using a next-generation technology and bioinformatics tools, we made a first step towards the genomic investigation of the Common Blackbird. The de novo assembled transcriptome is usable for downstream analyses such as differential gene expression analysis and SNP identification. This study shows the importance of the approach to sequence single tissues to understand functions of tissues, proteins and the phenotype.

Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions

The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube-dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non-coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

The transcriptome of a "sleeping" invader: de novo assembly and annotation of the transcriptome of aestivating Cornu aspersum

Background

Cornu aspersum is a quite intriguing species from the point of view of ecology and evolution and its potential use in medical and environmental applications. It is a species of economic importance since it is farmed and used for culinary purposes. However, the genomic tools that would allow a thorough insight into the ecology, evolution, nutritional and medical properties of this highly adaptable organism, are missing. In this work, using next-generation sequencing (NGS) techniques we assessed a significant portion of the transcriptome of this non-model organism.

Results

Out of the 9445 de novo assembled contigs, 2886 (30.6%) returned significant hits and for 2261 (24%) of them Gene Ontology (GO) terms associated to the hits were retrieved. A high percentage of the contigs (69.4%) produced no BLASTx hits. The GO terms were grouped to reflect biological processes, molecular functions and cellular components. Certain GO terms were dominant in all groups. After scanning the assembled transcriptome for microsatellites (simple sequence repeats, SSRs), a total of 563 SSRs were recovered. Among the identified SSRs, trinucleotide repeats were the predominant followed by tetranucleotide and dinucleotide repeats.

Conclusion

The annotation success of the transcriptome of C. aspersum was relatively low. This is probably due to the very limited number of annotated reference genomes existing for mollusc species, especially terrestrial ones. Several biological processes being active in the aestivating species were revealed through the association of the transcripts to enzymes relating to the pathways. The genomic tools provided herein will eventually aid in the study of the global genomic diversity of the species and the investigation of aspects of the ecology, evolution, behavior, nutritional and medical properties of this highly adaptable organism.

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.

An annotated draft genome for Radix auricularia (Gastropoda, Mollusca)

Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia . Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms

Background

Next-generation sequencing (NGS) technologies are arguably the most revolutionary technical development to join the list of tools available to molecular biologists since PCR. For researchers working with nonconventional model organisms one major problem with the currently dominant NGS platform (Illumina) stems from the obligatory fragmentation of nucleic acid material that occurs prior to sequencing during library preparation. This step creates a significant bioinformatic challenge for accurate de novo assembly of novel transcriptome data. This challenge becomes apparent when a variety of modern assembly tools (of which there is no shortage) are applied to the same raw NGS dataset. With the same assembly parameters these tools can generate markedly different assembly outputs.

Results

In this study we present an approach that generates an optimized consensus de novo assembly of eukaryotic coding transcriptomes. This approach does not represent a new assembler, rather it combines the outputs of a variety of established assembly packages, and removes redundancy via a series of clustering steps. We test and validate our approach using Illumina datasets from six phylogenetically diverse eukaryotes (three metazoans, two plants and a yeast) and two simulated datasets derived from metazoan reference genome annotations. All of these datasets were assembled using three currently popular assembly packages (CLC, Trinity and IDBA-tran). In addition, we experimentally demonstrate that transcripts unique to one particular assembly package are likely to be bioinformatic artefacts. For all eight datasets our pipeline generates more concise transcriptomes that in fact possess more unique annotatable protein domains than any of the three individual assemblers we employed. Another measure of assembly completeness (using the purpose built BUSCO databases) also confirmed that our approach yields more information.

Conclusions

Our approach yields coding transcriptome assemblies that are more likely to be closer to biological reality than any of the three individual assembly packages we investigated. This approach (freely available as a simple perl script) will be of use to researchers working with species for which there is little or no reference data against which the assembly of a transcriptome can be performed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1406-x) contains supplementary material, which is available to authorized users.

Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery

Satsuma myomphala is critically endangered through loss of natural habitats, predation by natural enemies, and indiscriminate collection. It is a protected species in Korea but lacks genomic resources for an understanding of varied functional processes attributable to evolutionary success under natural habitats. For assessing the genetic information of S. myomphala, we performed for the first time, de novo transcriptome sequencing and functional annotation of expressed sequences using Illumina Next-Generation Sequencing (NGS) platform and bioinformatics analysis. We identified 103,774 unigenes of which 37,959, 12,890, and 17,699 were annotated in the PANM (Protostome DB), Unigene, and COG (Clusters of Orthologous Groups) databases, respectively. In addition, 14,451 unigenes were predicted under Gene Ontology functional categories, with 4581 assigned to a single category. Furthermore, 3369 sequences with 646 having Enzyme Commission (EC) numbers were mapped to 122 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. The prominent protein domains included the Zinc finger (C2H2-like), Reverse Transcriptase, Thioredoxin-like fold, and RNA recognition motif domain. Many unigenes with homology to immunity, defense, and reproduction-related genes were screened in the transcriptome. We also detected 3120 putative simple sequence repeats (SSRs) encompassing dinucleotide to hexanucleotide repeat motifs from >1kb unigene sequences. A list of PCR primers of SSR loci have been identified to study the genetic polymorphisms. The transcriptome data represents a valuable resource for further investigations on the species genome structure and biology. The unigenes information and microsatellites would provide an indispensable tool for conservation of the species in natural and adaptive environments.

Deep, multi-stage transcriptome of the schistosomiasis vector Biomphalaria glabrata provides platform for understanding molluscan disease-related pathways

BACKGROUND

The gastropod mollusc Biomphalaria glabrata is well known as a vector for the tropical disease schistosomiasis, which affects nearly 200 million people worldwide. Despite intensive study, our understanding of the genetic basis of B. glabrata development, growth and disease resistance is constrained by limited genetic resources, constraints for which next-generation sequencing methods provide a ready solution.

METHODS

Illumina sequencing and de novo assembly using the Trinity program was used to generate a high-quality transcriptomic dataset spanning the entirety of in ovo development in schistosomiasis-free B. glabrata. This was subjected to automated (KEGG, BLAST2GO) and manual annotation efforts, allowing insight into the gene complements of this species in a number of contexts.

RESULTS

Excellent dataset recovery was observed, with 133,084 contigs produced of mean size 2219.48 bp. 80,952 (60.8 %) returned a BLASTx hit with an E value of less than 10^-3, and 74,492 (55.97 %) were either mapped or assigned a GO identity using the BLAST2GO program. The CEGMA set of core eukaryotic genes was found to be 99.6 % present, indicating exceptional transcriptome completeness. We were able to identify a wealth of disease-pathway related genes within our dataset, including the Wnt, apoptosis and Notch pathways. This provides an invaluable reference point for further work into molluscan development and evolution, for studying the impact of schistosomiasis in this species, and perhaps providing targets for the treatment of this widespread disease.

CONCLUSIONS

Here we present a deep transcriptome of an embryonic sample of schistosomiasis-free B. glabrata, presenting a comprehensive dataset for comparison to disease-affected specimens and from which conclusions can be drawn about the genetics of this widespread medical model. Furthermore, the dataset provided by this sequencing provides a useful reference point for comparison to other mollusc species, which can be used to better understand the evolution of this commercially, ecologically and medically important phylum.

De novo transcriptome assembly of the marine gastropod Reishia clavigera for supporting toxic mechanism studies

The intertidal whelk Reishia clavigera is commonly used as a biomonitor of chemical contamination in the marine environment along Western Pacific region, and as a model for mechanistic studies of organotin-mediated imposex development. However, limited genomic resources of R. clavigera have restricted its role for the investigation of molecular mechanisms of such endocrine disruptions. This study, therefore, aimed to establish tissue-specific transcriptomes of the digestive gland, gonad, head ganglia, penis and the remaining body part of the male and female R. clavigera. By combining the results, a global transcriptome was obtained. A total of 578,134,720 high-quality filtered reads were obtained using Illumina sequencing. The R. clavigera transcriptome comprised of 38,466 transcripts and 32,798 unigenes with predicted open reading frames. The average length of transcripts was 1,709bp with N50 of 2,236bp. Based on sequence similarity searches against public databases, 28,657 transcripts and 24,403 unigenes had at least one BLAST hit. There were 17,530 transcripts and 14,897 unigenes annotated with at least one Gene Ontology (GO) term. Moreover, 5,776 transcripts and 5,137 unigenes were associated with 333 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways. The numbers of unigenes were similar among the five target tissues and between sexes, but tissue-specific expression profiles were revealed by multivariate analyses. Based on the functional annotation, putative steroid hormone-associated unigenes were identified. In particular, we highlighted the presence of steroid hormone receptor homologues that could be the targets for mechanistic studies of the organotin-mediated imposex development in marine gastropods. This newly generated transcriptome assembly of R. clavigera provides a valuable molecular resource for ecotoxicological and environmental genomic studies.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

Transcriptomic and Proteomic Analysis of Arion vulgaris--Proteins for Probably Successful Survival Strategies?

The Spanish slug, Arion vulgaris, is considered one of the hundred most invasive species in Central Europe. The immense and very successful adaptation and spreading of A. vulgaris suggest that it developed highly effective mechanisms to deal with infections and natural predators. Current transcriptomic and proteomic studies on gastropods have been restricted mainly to marine and freshwater gastropods. No transcriptomic or proteomic study on A. vulgaris has been carried out so far, and in the current study, the first transcriptomic database from adult specimen of A. vulgaris is reported. To facilitate and enable proteomics in this non-model organism, a mRNA-derived protein database was constructed for protein identification. A gel-based proteomic approach was used to obtain the first generation of a comprehensive slug mantle proteome. A total of 2128 proteins were unambiguously identified; 48 proteins represent novel proteins with no significant homology in NCBI non-redundant database. Combined transcriptomic and proteomic analysis revealed an extensive repertoire of novel proteins with a role in innate immunity including many associated pattern recognition, effector proteins and cytokine-like proteins. The number and diversity in gene families encoding lectins point to a complex defense system, probably as a result of adaptation to a pathogen-rich environment. These results are providing a fundamental and important resource for subsequent studies on molluscs as well as for putative antimicrobial compounds for drug discovery and biomedical applications.

De novo Transcriptome Generation and Annotation for Two Korean Endemic Land Snails, Aegista chejuensis and Aegista quelpartensis, Using Illumina Paired-End Sequencing Technology

Aegista chejuensis and Aegista quelpartensis (Family-Bradybaenidae) are endemic to Korea, and are considered vulnerable due to declines in their population. The limited genetic resources for these species restricts the ability to prioritize conservation efforts. We sequenced the transcriptomes of these species using Illumina paired-end technology. Approximately 257 and 240 million reads were obtained and assembled into 198,531 and 230,497 unigenes for A. chejuensis and A. quelpartensis, respectively. The average and N50 unigene lengths were 735.4 and 1073 bp, respectively, for A. chejuensis, and 705.6 and 1001 bp, respectively, for A. quelpartensis. In total, 68,484 (34.5%) and 77,745 (33.73%) unigenes for A. chejuensis and A. quelpartensis, respectively, were annotated to databases. Gene Ontology terms were assigned to 23,778 (11.98%) and 26,396 (11.45) unigenes, for A. chejuensis and A. quelpartensis, respectively, while 5050 and 5838 unigenes were mapped to 117 and 124 pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, we identified and annotated 9542 and 10,395 putative simple sequence repeats (SSRs) in unigenes from A. chejuensis and A. quelpartensis, respectively. We designed a list of PCR primers flanking the putative SSR regions. These microsatellites may be utilized for future phylogenetics and conservation initiatives.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.

Analysis of RNA-Seq Data Using TopHat and Cufflinks

The recent advances in high throughput RNA sequencing (RNA-Seq) have generated huge amounts of data in a very short span of time for a single sample. These data have required the parallel advancement of computing tools to organize and interpret them meaningfully in terms of biological implications, at the same time using minimum computing resources to reduce computation costs. Here we describe the method of analyzing RNA-seq data using the set of open source software programs of the Tuxedo suite: TopHat and Cufflinks. TopHat is designed to align RNA-seq reads to a reference genome, while Cufflinks assembles these mapped reads into possible transcripts and then generates a final transcriptome assembly. Cufflinks also includes Cuffdiff, which accepts the reads assembled from two or more biological conditions and analyzes their differential expression of genes and transcripts, thus aiding in the investigation of their transcriptional and post transcriptional regulation under different conditions. We also describe the use of an accessory tool called CummeRbund, which processes the output files of Cuffdiff and gives an output of publication quality plots and figures of the user's choice. We demonstrate the effectiveness of the Tuxedo suite by analyzing RNA-Seq datasets of Arabidopsis thaliana root subjected to two different conditions.

Transcriptome Characterization for Non-Model Endangered Lycaenids, Protantigius superans and Spindasis takanosis, Using Illumina HiSeq 2500 Sequencing

The Lycaenidae butterflies, Protantigius superans and Spindasis takanosis, are endangered insects in Korea known for their symbiotic association with ants. However, necessary genomic and transcriptomics data are lacking in these species, limiting conservation efforts. In this study, the P. superans and S. takanosis transcriptomes were deciphered using Illumina HiSeq 2500 sequencing. The P. superans and S. takanosis transcriptome data included a total of 254,340,693 and 245,110,582 clean reads assembled into 159,074 and 170,449 contigs and 107,950 and 121,140 unigenes, respectively. BLASTX hits (E-value of 1.0 × 10⁻⁵) against the known protein databases annotated a total of 46,754 and 51,908 transcripts for P. superans and S. takanosis. Approximately 41.25% and 38.68% of the unigenes for P. superans and S. takanosis found homologous sequences in Protostome DB (PANM-DB). BLAST2GO analysis confirmed 18,611 unigenes representing Gene Ontology (GO) terms and a total of 5259 unigenes assigned to 116 pathways for P. superans. For S. takanosis, a total of 6697 unigenes were assigned to 119 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 382,164 and 390,516 Simple Sequence Repeats (SSRs) were compiled from the unigenes of P. superans and S. takanosis, respectively. This is the first report to record new genes and their utilization for conservation of lycaenid species population and as a reference information for closely related species.

Data set for transcriptome analysis of the Chinese giant salamander (Andrias davidianus )

The Chinese giant salamander (Andrias davidianus) occupies a seat at the phylogenetic and species evolution process, which makes it an invaluable model for genetics; however, the genetic information and gene sequences about the Chinese giant salamander in public databases are scanty. Hence, we aimed to perform transcriptome analysis with the help of high-throughput sequencing. In this data, 61,317,940 raw reads were acquired from Chinese giant salamander mRNA using Illumina paired-end sequencing platform. After de novo assembly, a total of 72,072 unigenes were gained, in which 33,834 (46.95%) and 29,479 (40.91%) transcripts exhibited homology to sequences in the Nr database and Swiss-Prot database, (E-value <10(-5)), respectively. In the obtained unigenes, 18,019 (25%) transcripts were assigned with at least one Gene Ontology term, of which 1218 (6.8%) transcripts were assigned to immune system processes. In addition, a total of 17,572 assembled sequences were assigned into 241 predicted KEGG metabolic pathways. Among these, 2552 (14.5%) transcripts were assigned to the immune system relevant pathway and 5 transcripts were identified as potential antimicrobial peptides (AMPs).

High throughput de novo RNA sequencing elucidates novel responses in Penicillium chrysogenum under microgravity

In this study, the transcriptional alterations in Penicillium chrysogenum under simulated microgravity conditions were analyzed for the first time using an RNA-Seq method. The increasing plethora of eukaryotic microbial flora inside the spaceship demands the basic understanding of fungal biology in the absence of gravity vector. Penicillium species are second most dominant fungal contaminant in International Space Station. Penicillium chrysogenum an industrially important organism also has the potential to emerge as an opportunistic pathogen for the astronauts during the long-term space missions. But till date, the cellular mechanisms underlying the survival and adaptation of Penicillium chrysogenum to microgravity conditions are not clearly elucidated. A reference genome for Penicillium chrysogenum is not yet available in the NCBI database. Hence, we performed comparative de novo transcriptome analysis of Penicillium chrysogenum grown under microgravity versus normal gravity. In addition, the changes due to microgravity are documented at the molecular level. Increased response to the environmental stimulus, changes in the cell wall component ABC transporter/MFS transporters are noteworthy. Interestingly, sustained increase in the expression of Acyl-coenzyme A: isopenicillin N acyltransferase (Acyltransferase) under microgravity revealed the significance of gravity in the penicillin production which could be exploited industrially.

Transcriptome sequencing, annotation and polymorphism detection in the hop bush, Dodonaea viscosa

Background

The hop bush, Dodonaea viscosa, is a trans-oceanic species distributed oversix continents. It evolved in Australia where it is found over a wide range of habitat types and is an ecologically important species. Limited genomic resources are currently available for this species, thus our understanding of its evolutionary history and ecological adaptation is restricted. Here, we present a comprehensive transcriptome dataset for future genomic studies into this species.

Methods

We performed Illumina sequencing of cDNA prepared from leaf tissue collected from seven populations of D. viscosa ssp. angustissima and spatulata distributed along an environmental gradient in South Australia. Sequenced reads were assembled to provide a transcriptome resource. Contiguous sequences (contigs) were annotated using BLAST searches against the NCBI non-redundant database and gene ontology definitions were assigned. Single nucleotide polymorphisms were detected for the establishment of a genetic marker set. A comparison between the two subspecies was also carried out.

Results

Illumina sequencing returned 268,672,818 sequence reads, which were de novoassembled into 105,125 contigs. Contigs with significant BLAST alignments (E value < 1e^-5)numbered at 44,191, with 38,311 of these having their most significant hits to sequences from land plant species. Gene Ontology terms were assigned to 28,440 contigs and KEGG analysis identified 146 pathways that the gene products from 5,070 contigs are potentially involved in. The subspecies comparison identified 8,494 fixed SNP differences across 3,979 contiguous sequences, indicating a level of genetic differentiation between them. Across all samples, 248,235 SNPs were detected.

Conclusions

We have established a significant genomic data resource for D. viscosa,providing a comprehensive transcriptomic reference. Genetic differences among morphologically distinct subspecies were found. A wide range of putative gene regions were identified along with a large set of variable SNP markers, providing a basis for studies into the evolution and ecological adaptation of D. viscosa.

Electronic supplementary material

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

A Transcriptomic Analysis of Cave, Surface, and Hybrid Isopod Crustaceans of the Species Asellus aquaticus

Cave animals, compared to surface-dwelling relatives, tend to have reduced eyes and pigment, longer appendages, and enhanced mechanosensory structures. Pressing questions include how certain cave-related traits are gained and lost, and if they originate through the same or different genetic programs in independent lineages. An excellent system for exploring these questions is the isopod, Asellus aquaticus. This species includes multiple cave and surface populations that have numerous morphological differences between them. A key feature is that hybrids between cave and surface individuals are viable, which enables genetic crosses and linkage analyses. Here, we advance this system by analyzing single animal transcriptomes of Asellus aquaticus. We use high throughput sequencing of non-normalized cDNA derived from the head of a surface-dwelling male, the head of a cave-dwelling male, the head of a hybrid male (produced by crossing a surface individual with a cave individual), and a pooled sample of surface embryos and hatchlings. Assembling reads from surface and cave head RNA pools yielded an integrated transcriptome comprised of 23,984 contigs. Using this integrated assembly as a reference transcriptome, we aligned reads from surface-, cave- and hybrid- head tissue and pooled surface embryos and hatchlings. Our approach identified 742 SNPs and placed four new candidate genes to an existing linkage map for A. aquaticus. In addition, we examined SNPs for allele-specific expression differences in the hybrid individual. All of these resources will facilitate identification of genes and associated changes responsible for cave adaptation in A. aquaticus and, in concert with analyses of other species, will inform our understanding of the evolutionary processes accompanying adaptation to the subterranean environment.

Transcriptome analysis of the pearl oyster (Pinctada fucata) hemocytes in response to Vibrio alginolyticus infection

The pearl oyster Pinctada fucata is cultured widely for production of marine pearls in China, while mass mortalities, likely related to pathogenic infections, have occurred frequently in juvenile, mother and operated oysters. To address this issue, understanding host defense mechanisms of P. fucata against pathogenic challenge is extremely important. In the present study, a comparative analysis of hemocyte transcriptomes of P. fucata before and after Vibrio alginolyticus infection was conducted using the Illumina/Hiseq-2000 RNA-Seq technology. A total of 56,345,139 clean reads were generated and then assembled into 74,007 unigenes with an average length of 680 bp and an N50 of 1197 bp. Unigenes were annotated by comparing against non-redundant protein sequence (nr), non-redundant nucleotide (nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 29,615 unigenes (40.01%) were annotated in at least one database. There were 636 genes (518 up-regulated and 118 down-regulated) that were significantly differentially expressed after bacterial challenge, and among which 369 were associated with 122 pathways, including classical immune-related pathways, such as 'MAPK signaling pathway', 'Chemokine signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. These findings provide information on the pearl oyster innate immunity and may contribute to developing strategies for management of diseases and long-term sustainability of P. fucata culture.

Rapid evolution of chemosensory receptor genes in a pair of sibling species of orchid bees (Apidae: Euglossini)

Background

Insects rely more on chemical signals (semiochemicals) than on any other sensory modality to find, identify, and choose mates. In most insects, pheromone production is typically regulated through biosynthetic pathways, whereas pheromone sensory detection is controlled by the olfactory system. Orchid bees are exceptional in that their semiochemicals are not produced metabolically, but instead male bees collect odoriferous compounds (perfumes) from the environment and store them in specialized hind-leg pockets to subsequently expose during courtship display. Thus, the olfactory sensory system of orchid bees simultaneously controls male perfume traits (sender components) and female preferences (receiver components). This functional linkage increases the opportunities for parallel evolution of male traits and female preferences, particularly in response to genetic changes of chemosensory detection (e.g. Odorant Receptor genes). To identify whether shifts in pheromone composition among related lineages of orchid bees are associated with divergence in chemosensory genes of the olfactory periphery, we searched for patterns of divergent selection across the antennal transcriptomes of two recently diverged sibling species Euglossa dilemma and E. viridissima.

Results

We identified 3185 orthologous genes including 94 chemosensory loci from five different gene families (Odorant Receptors, Ionotropic Receptors, Gustatory Receptors, Odorant Binding Proteins, and Chemosensory Proteins). Our results revealed that orthologs with signatures of divergent selection between E. dilemma and E. viridissima were significantly enriched for chemosensory genes. Notably, elevated signals of divergent selection were almost exclusively observed among chemosensory receptors (i.e. Odorant Receptors).

Conclusions

Our results suggest that rapid changes in the chemosensory gene family occurred among closely related species of orchid bees. These findings are consistent with the hypothesis that strong divergent selection acting on chemosensory receptor genes plays an important role in the evolution and diversification of insect pheromone systems.

Electronic supplementary material

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

Characterizing a hybrid zone between a cryptic species pair of freshwater snails

Characterizing hybrid zones and their dynamics is a central goal in evolutionary biology, but this is particularly challenging for morphologically cryptic species. The lack of conspicuous divergence between parental types means intermediate hybrid forms often go undetected. We aimed to detect and characterize a suspected hybrid zone between a pair of morphologically cryptic lineages of the freshwater snail, Radix. We sampled Radix from across a contact zone between two mitochondrial lineages (Radix balthica and an undescribed lineage termed 'MOTU3') and detected admixture between two nuclear genotype clusters, which were significantly but not categorically associated with the mitochondrial lineages. Using a model selection approach, we show that the admixture cline is best explained by an interaction between precipitation and temperature gradients over the area, rather than geographic distance. We thus hypothesize that the correlation with climatic gradients suggests environmental selection has played a role in maintaining the hybrid zone. In a 2050 climate change scenario, we furthermore predict an expansion of one of the nuclear clusters and a widening of the hybrid zone as the climate warms and dries.

Positive selection in development and growth rate regulation genes involved in species divergence of the genus Radix

Background

Life history traits like developmental time, age and size at maturity are directly related to fitness in all organisms and play a major role in adaptive evolution and speciation processes. Comparative genomic or transcriptomic approaches to identify positively selected genes involved in species divergence can help to generate hypotheses on the driving forces behind speciation. Here we use a bottom-up approach to investigate this hypothesis by comparative analysis of orthologous transcripts of four closely related European Radix species.

Results

Snails of the genus Radix occupy species specific distribution ranges with distinct climatic niches, indicating a potential for natural selection driven speciation based on ecological niche differentiation. We then inferred phylogenetic relationships among the four Radix species based on whole mt-genomes plus 23 nuclear loci. Three different tests to infer selection and changes in amino acid properties yielded a total of 134 genes with signatures of positive selection. The majority of these genes belonged to the functional gene ontology categories “reproduction” and “genitalia” with an overrepresentation of the functions “development” and “growth rate”.

Conclusions

We show here that Radix species divergence may be primarily enforced by selection on life history traits such as (larval-) development and growth rate. We thus hypothesise that life history differences may confer advantages under the according climate regimes, e.g., species occupying warmer and dryer habitats might have a fitness advantage with fast developing susceptible life stages, which are more tolerant to habitat desiccation.

Electronic supplementary material

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

An embryonic transcriptome of the pulmonate snail Radix balthica

The pond snail, Radix balthica (Linnaeus 1758), is an emerging model species within ecological developmental biology. While its development has been characterised in detail, genomic resources for embryonic stages are lacking. We applied Illumina MiSeq RNA-seq to RNA isolated from pools of embryos at two points during development. Embryos were cultured in either the presence or absence of predator kariomones to increase the diversity of the transcripts assembled. Sequencing produced 47.2M paired-end reads, assembled into 54,360 contigs of which 73% were successfully annotated. This transcriptome provides an invaluable resource to build a mechanistic understanding of developmental plasticity.

Transcriptomics Analysis of Crassostrea hongkongensis for the Discovery of Reproduction-Related Genes

BACKGROUND

The reproductive mechanisms of mollusk species have been interesting targets in biological research because of the diverse reproductive strategies observed in this phylum. These species have also been studied for the development of fishery technologies in molluscan aquaculture. Although the molecular mechanisms underlying the reproductive process have been well studied in animal models, the relevant information from mollusks remains limited, particularly in species of great commercial interest. Crassostrea hongkongensis is the dominant oyster species that is distributed along the coast of the South China Sea and little genomic information on this species is available. Currently, high-throughput sequencing techniques have been widely used for investigating the basis of physiological processes and facilitating the establishment of adequate genetic selection programs.

RESULTS

The C.hongkongensis transcriptome included a total of 1,595,855 reads, which were generated by 454 sequencing and were assembled into 41,472 contigs using de novo methods. Contigs were clustered into 33,920 isotigs and further grouped into 22,829 isogroups. Approximately 77.6% of the isogroups were successfully annotated by the Nr database. More than 1,910 genes were identified as being related to reproduction. Some key genes involved in germline development, sex determination and differentiation were identified for the first time in C.hongkongensis (nanos, piwi, ATRX, FoxL2, β-catenin, etc.). Gene expression analysis indicated that vasa, nanos, piwi, ATRX, FoxL2, β-catenin and SRD5A1 were highly or specifically expressed in C.hongkongensis gonads. Additionally, 94,056 single nucleotide polymorphisms (SNPs) and 1,699 simple sequence repeats (SSRs) were compiled.

CONCLUSIONS

Our study significantly increased C.hongkongensis genomic information based on transcriptomics analysis. The group of reproduction-related genes identified in the present study constitutes a new tool for research on bivalve reproduction processes. The large group of molecular markers discovered in this study will be useful for population screening and marker assisted selection programs in C.hongkongensis aquaculture.

De novo sequencing and characterization of the Bradysia odoriphaga (Diptera: Sciaridae) larval transcriptome

The most serious pestilent threat to the Chinese chive, Allium tuberosum Rottle ex Spreng (Liliaceae) is the Bradysia odoriphaga Yang and Zhang. There is limited genetic research focused on B. odoriphaga, partially due to the lack of genomic resources. The advent of high-throughput sequencing technologies has enabled generation of genomic resources in a short time frame and at minimal costs. In this study, we performed, for the first time, de novo transcriptome sequencing of the B. odoriphaga. Here, 16,829 unigenes were assembled from the total reads, 12,024 of these unigenes were annotated in the NCBI NR protein database, and 9784 were annotated in the Swiss-Prot database. Of these annotated unigenes, 7903 and 5060 unigenes have been assigned to gene ontology categories and clusters of orthologous groups, respectively. Furthermore, 8647 unigenes were mapped to 257 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. We found that 408 unigenes were related to insecticide resistance and metabolism. In addition, 23,122 simple sequence repeats (SSRs) were identified in 11,009 unigenes, and 100 PCR primers of SSR loci were used to validate the assembly quality and polymorphisms. These results provide a good platform for further investigations into the insecticide resistance of B. odoriphaga. Finally, the SSRs identified in B. odoriphaga may be a useful genomic resource.

Optimized deep-targeted proteotranscriptomic profiling reveals unexplored Conus toxin diversity and novel cysteine frameworks

Cone snails are predatory marine gastropods characterized by a sophisticated venom apparatus responsible for the biosynthesis and delivery of complex mixtures of cysteine-rich toxin peptides. These conotoxins fold into small highly structured frameworks, allowing them to potently and selectively interact with heterologous ion channels and receptors. Approximately 2,000 toxins from an estimated number of >70,000 bioactive peptides have been identified in the genus Conus to date. Here, we describe a high-resolution interrogation of the transcriptomes (available at www.ddbj.nig.ac.jp) and proteomes of the diverse compartments of the Conus episcopatus venom apparatus. Using biochemical and bioinformatic tools, we found the highest number of conopeptides yet discovered in a single Conus specimen, with 3,305 novel precursor toxin sequences classified into 9 known superfamilies (A, I1, I2, M, O1, O2, S, T, Z), and identified 16 new superfamilies showing unique signal peptide signatures. We were also able to depict the largest population of venom peptides containing the pharmacologically active C-C-CC-C-C inhibitor cystine knot and CC-C-C motifs (168 and 44 toxins, respectively), as well as 208 new conotoxins displaying odd numbers of cysteine residues derived from known conotoxin motifs. Importantly, six novel cysteine-rich frameworks were revealed which may have novel pharmacology. Finally, analyses of codon usage bias and RNA-editing processes of the conotoxin transcripts demonstrate a specific conservation of the cysteine skeleton at the nucleic acid level and provide new insights about the origin of sequence hypervariablity in mature toxin regions.

The venomous cocktail of the vampire snail Colubraria reticulata (Mollusca, Gastropoda)

Background

Hematophagy arose independently multiple times during metazoan evolution, with several lineages of vampire animals particularly diversified in invertebrates. However, the biochemistry of hematophagy has been studied in a few species of direct medical interest and is still underdeveloped in most invertebrates, as in general is the study of venom toxins. In cone snails, leeches, arthropods and snakes, the strong target specificity of venom toxins uniquely aligns them to industrial and academic pursuits (pharmacological applications, pest control etc.) and provides a biochemical tool for studying biological activities including cell signalling and immunological response. Neogastropod snails (cones, oyster drills etc.) are carnivorous and include active predators, scavengers, grazers on sessile invertebrates and hematophagous parasites; most of them use venoms to efficiently feed. It has been hypothesized that trophic innovations were the main drivers of rapid radiation of Neogastropoda in the late Cretaceous.

We present here the first molecular characterization of the alimentary secretion of a non-conoidean neogastropod, Colubraria reticulata. Colubrariids successfully feed on the blood of fishes, throughout the secretion into the host of a complex mixture of anaesthetics and anticoagulants. We used a NGS RNA-Seq approach, integrated with differential expression analyses and custom searches for putative secreted feeding-related proteins, to describe in detail the salivary and mid-oesophageal transcriptomes of this Mediterranean vampire snail, with functional and evolutionary insights on major families of bioactive molecules.

Results

A remarkably low level of overlap was observed between the gene expression in the two target tissues, which also contained a high percentage of putatively secreted proteins when compared to the whole body. At least 12 families of feeding-related proteins were identified, including: 1) anaesthetics, such as ShK Toxin-containing proteins and turripeptides (ion-channel blockers), Cysteine-rich secretory proteins (CRISPs), Adenosine Deaminase (ADA); 2) inhibitors of primary haemostasis, such as novel vWFA domain-containing proteins, the Ectonucleotide pyrophosphatase/phosphodiesterase family member 5 (ENPP5) and the wasp Antigen-5; 3) anticoagulants, such as TFPI-like multiple Kunitz-type protease inhibitors, Peptidases S1 (PS1), CAP/ShKT domain-containing proteins, Astacin metalloproteases and Astacin/ShKT domain-containing proteins; 4) additional proteins, such the Angiotensin-Converting Enzyme (ACE: vasopressive) and the cytolytic Porins.

Conclusions

Colubraria feeding physiology seems to involve inhibitors of both primary and secondary haemostasis, anaesthetics, a vasoconstrictive enzyme to reduce feeding time and tissue-degrading proteins such as Porins and Astacins. The complexity of Colubraria venomous cocktail and the divergence from the arsenal of the few neogastropods studied to date (mostly conoideans) suggest that biochemical diversification of neogastropods might be largely underestimated and worth of extensive investigation.

Electronic supplementary material

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

RNA-Seq analysis and gene discovery of Andrias davidianus using Illumina short read sequencing

The Chinese giant salamander, Andrias davidianus, is an important species in the course of evolution; however, there is insufficient genomic data in public databases for understanding its immunologic mechanisms. High-throughput transcriptome sequencing is necessary to generate an enormous number of transcript sequences from A. davidianus for gene discovery. In this study, we generated more than 40 million reads from samples of spleen and skin tissue using the Illumina paired-end sequencing technology. De novo assembly yielded 87,297 transcripts with a mean length of 734 base pairs (bp). Based on the sequence similarities, searching with known proteins, 38,916 genes were identified. Gene enrichment analysis determined that 981 transcripts were assigned to the immune system. Tissue-specific expression analysis indicated that 443 of transcripts were specifically expressed in the spleen and skin. Among these transcripts, 147 transcripts were found to be involved in immune responses and inflammatory reactions, such as fucolectin, β-defensins and lymphotoxin beta. Eight tissue-specific genes were selected for validation using real time reverse transcription quantitative PCR (qRT-PCR). The results showed that these genes were significantly more expressed in spleen and skin than in other tissues, suggesting that these genes have vital roles in the immune response. This work provides a comprehensive genomic sequence resource for A. davidianus and lays the foundation for future research on the immunologic and disease resistance mechanisms of A. davidianus and other amphibians.

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

RNA-seq analysis of Macrobrachium rosenbergii hepatopancreas in response to Vibrio parahaemolyticus infection

Background

The Malaysian giant freshwater prawn, Macrobrachium rosenbergii, is an economically important crustacean worldwide. However, production of this prawn is facing a serious threat from Vibriosis disease caused by Vibrio species such as Vibrio parahaemolyticus. Unfortunately, the mechanisms involved in the immune response of this species to bacterial infection are not fully understood. We therefore used a high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the hepatopancreas from this freshwater prawn infected with V. parahaemolyticus to gain an increased understanding of the molecular mechanisms underlying the species’ immune response to this pathogenic bacteria.

Result

A total of 59,122,940 raw reads were obtained from the control group, and 58,385,094 reads from the Vibrio-infected group. Via de novo assembly by Trinity assembler, 59,050 control unigenes and 73,946 Vibrio-infected group unigenes were obtained. By clustering unigenes from both libraries, a total of 64,411 standard unigenes were produced. The standard unigenes were annotated against the NCBI non-redundant, Swiss-Prot, Kyoto Encyclopaedia of Genes and Genome pathway (KEGG) and Orthologous Groups of Proteins (COG) databases, with 19,799 (30.73%), 16,832 (26.13%), 14,706 (22.83%) and 7,856 (12.19%) hits respectively, giving a final total of 22,455 significant hits (34.86% of all unigenes). A Gene Ontology (GO) analysis search using the Blast2GO program resulted in 6,007 unigenes (9.32%) being categorized into 55 functional groups. A differential gene expression analysis produced a total of 14,569 unigenes aberrantly expressed, with 11,446 unigenes significantly up-regulated and 3,103 unigenes significantly down-regulated. The differentially expressed immune genes fall under various processes of the animal immune system.

Conclusion

This study provided an insight into the antibacterial mechanism in M. rosenbergii and the role of differentially expressed immune genes in response to V. parahaemolyticus infection. Furthermore, this study has generated an abundant list of transcript from M.rosenbergii which will provide a fundamental basis for future genomics research in this field.

Electronic supplementary material

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

Pathways associated with lignin biosynthesis in lignomaniac jute fibres

We generated the bast transcriptomes of a deficient lignified phloem fibre mutant and its wild-type jute (Corchorus capsularis) using Illumina paired-end sequencing. A total of 34,163 wild-type and 29,463 mutant unigenes, with average lengths of 1442 and 1136 bp, respectively, were assembled de novo, ~77-79 % of which were functionally annotated. These annotated unigenes were assigned to COG (~37-40 %) and GO (~22-28 %) classifications and mapped to 189 KEGG pathways (~19-21 %). We discovered 38 and 43 isoforms of 16 and 10 genes of the upstream shikimate-aromatic amino acid and downstream monolignol biosynthetic pathways, respectively, rendered their sequence similarities, confirmed the identities of 22 of these candidate gene families by phylogenetic analyses and reconstructed the pathway leading to lignin biosynthesis in jute fibres. We also identified major genes and bast-related transcription factors involved in secondary cell wall (SCW) formation. The quantitative RT-PCRs revealed that phenylalanine ammonia-lyase 1 (CcPAL1) was co-down-regulated with several genes of the upstream shikimate pathway in mutant bast tissues at an early growth stage, although its expression relapsed to the normal level at the later growth stage. However, cinnamyl alcohol dehydrogenase 7 (CcCAD7) was strongly down-regulated in mutant bast tissues irrespective of growth stages. CcCAD7 disruption at an early growth stage was accompanied by co-up-regulation of SCW-specific genes cellulose synthase A7 (CcCesA7) and fasciclin-like arabinogalactan 6 (CcFLA6), which was predicted to be involved in coordinating the S-layers' deposition in the xylan-type jute fibres. Our results identified CAD as a promising target for developing low-lignin jute fibres using genomics-assisted molecular approaches.

Inferring bona fide transfrags in RNA-Seq derived-transcriptome assemblies of non-model organisms

BACKGROUND

De novo transcriptome assembly of short transcribed fragments (transfrags) produced from sequencing-by-synthesis technologies often results in redundant datasets with differing levels of unassembled, partially assembled or mis-assembled transcripts. Post-assembly processing intended to reduce redundancy typically involves reassembly or clustering of assembled sequences. However, these approaches are mostly based on common word heuristics and often create clusters of biologically unrelated sequences, resulting in loss of unique transfrags annotations and propagation of mis-assemblies.

RESULTS

Here, we propose a structured framework that consists of a few steps in pipeline architecture for Inferring Functionally Relevant Assembly-derived Transcripts (IFRAT). IFRAT combines 1) removal of identical subsequences, 2) error tolerant CDS prediction, 3) identification of coding potential, and 4) complements BLAST with a multiple domain architecture annotation that reduces non-specific domain annotation. We demonstrate that independent of the assembler, IFRAT selects bona fide transfrags (with CDS and coding potential) from the transcriptome assembly of a model organism without relying on post-assembly clustering or reassembly. The robustness of IFRAT is inferred on RNA-Seq data of Neurospora crassa assembled using de Bruijn graph-based assemblers, in single (Trinity and Oases-25) and multiple (Oases-Merge and additive or pooled) k-mer modes. Single k-mer assemblies contained fewer transfrags compared to the multiple k-mer assemblies. However, Trinity identified a comparable number of predicted coding sequence and gene loci to Oases pooled assembly. IFRAT selects bona fide transfrags representing over 94% of cumulative BLAST-derived functional annotations of the unfiltered assemblies. Between 4-6% are lost when orphan transfrags are excluded and this represents only a tiny fraction of annotation derived from functional transference by sequence similarity. The median length of bona fide transfrags ranged from 1.5kb (Trinity) to 2kb (Oases), which is consistent with the average coding sequence length in fungi. The fraction of transfrags that could be associated with gene ontology terms ranged from 33-50%, which is also high for domain based annotation. We showed that unselected transfrags were mostly truncated and represent sequences from intronic, untranslated (5' and 3') regions and non-coding gene loci.

CONCLUSIONS

IFRAT simplifies post-assembly processing providing a reference transcriptome enriched with functionally relevant assembly-derived transcripts for non-model organism.

Lectin-like molecules in transcriptome of Littorina littorea hemocytes

The common periwinkle Littorina littorea was introduced in the list of models for comparative immunobiology as a representative of phylogenetically important taxon Caenogastropoda. Using Illumina sequencing technology, we de novo assembled the transcriptome of Littorina littorea hemocytes from 182 million mRNA-Seq pair-end 100 bp reads into a total of 15,526 contigs clustered in 4472 unigenes. The transcriptome profile was analyzed for presence of carbohydrate-binding molecules in a variety of architectural contexts. Hemocytes' repertoire of lectin-like proteins bearing conserved carbohydrate-recognition domains (CRDs) is highly diversified, including 11 of 15 lectin families earlier described in animals, as well as the novel members of lectin family found for the first time in mollusc species. The new molluscan lineage-specific domain combinations were confirmed by cloning and sequencing, including the fuco-lectin related molecules (FLReMs) composed of N-terminal region with no sequence homology to any known protein, a middle Fucolectin Tachylectin-4 Pentaxrin (FTP) domain, and a C-terminal epidermal growth factor (EGF) repeat region. The repertoire of lectin-like molecules is discussed in terms of their potential participation in the receptor phase of immune response. In total, immune-associated functions may be attributed to 70 transcripts belonging to 6 lectin families. These lectin-like genes show low overlap between species of invertebrates, suggesting relatively rapid evolution of immune-associated genes in the group. The repertoire provides valuable candidates for further characterization of the gene functions in mollusc immunity.

Evaluation of de novo transcriptome assemblies from RNA-Seq data

De novo RNA-Seq assembly facilitates the study of transcriptomes for species without sequenced genomes, but it is challenging to select the most accurate assembly in this context. To address this challenge, we developed a model-based score, RSEM-EVAL, for evaluating assemblies when the ground truth is unknown. We show that RSEM-EVAL correctly reflects assembly accuracy, as measured by REF-EVAL, a refined set of ground-truth-based scores that we also developed. Guided by RSEM-EVAL, we assembled the transcriptome of the regenerating axolotl limb; this assembly compares favorably to a previous assembly. A software package implementing our methods, DETONATE, is freely available at http://deweylab.biostat.wisc.edu/detonate.

Evaluation of de novo transcriptome assemblies from RNA-Seq data

De novo RNA-Seq assembly facilitates the study of transcriptomes for species without sequenced genomes, but it is challenging to select the most accurate assembly in this context. To address this challenge, we developed a model-based score, RSEM-EVAL, for evaluating assemblies when the ground truth is unknown. We show that RSEM-EVAL correctly reflects assembly accuracy, as measured by REF-EVAL, a refined set of ground-truth-based scores that we also developed. Guided by RSEM-EVAL, we assembled the transcriptome of the regenerating axolotl limb; this assembly compares favorably to a previous assembly. A software package implementing our methods, DETONATE, is freely available at http://deweylab.biostat.wisc.edu/detonate.