Transcriptome analysis between invasive Pomacea canaliculata and indigenous Cipangopaludina cahayensis reveals genomic divergence and diagnostic microsatellite/SSR markers

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

Comparative transcriptome analysis on the mangrove Acanthus ilicifolius and its two terrestrial relatives provides insights into adaptation to intertidal habitats

Acanthus is a unique genus covering both mangroves and terrestrial species, and thus is an ideal system to comparatively analyze the mechanisms of mangrove adaptation to intertidal habitats. We performed RNA sequencing of the mangrove plant Acanthus ilicifolius and its two terrestrial relatives, Acanthus leucostachyus and Acanthus mollis. A total of 91,125, 118,290, and 141,640 unigenes were obtained. Simple sequence repeats (SSR) analysis showed that A. ilicifolius had more SSRs, the highest frequency of distribution, and higher in polymorphism potential compared to the two terrestrial relatives. Phylogenetic analyses suggested a relatively recent split between A. ilicifolius and A. leucostachyus, i.e., about 16.76 million years ago (Mya), after their ancestor divergence with A. mollis (32.11 Mya), indicating that speciation of three Acanthus species occurred in the Early to Middle Miocene. Gene Ontology (GO) enrichment revealed that the unique unigenes in A. ilicifolius are predominantly related to rhythmic process, reproductive process and response to stimuli. The accelerated evolution and positive selection analyses indicated that the genus Acanthus migrated from terrestrial to intertidal habitats, where 311 pairs may be under positive selection. Functional enrichment analysis revealed that these genes associated with essential metabolism and biosynthetic pathways such as oxidative phosphorylation, plant hormone signal transduction, photosynthetic carbon fixation and arginine and proline metabolism, are related to the adaptation of A. ilicifolius to intertidal habitats, which are characterized by high salinity and hypoxia. Our results indicate the evolutionary processes and the mechanisms underlying the adaptability of Acanthus to various harsh environments from the arid terrestrial to intertidal habitats.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

Transcriptomic responses of Artemia salina exposed to an environmentally relevant dose of Alexandrium minutum cells or Gonyautoxin2/3

Toxicities of the marine algae Alexandrium minutum and its excreted gonyautoxins (GTXs) to the marine crustacean Artemia salina were investigated. Mortality was observed for neither larvae nor adult A. salina exposed to A. minutum at a density of 5000 cells/mL or 0.5 μM GTX2/3. After exposure, the full transcriptome of adult A. salina was assembled and functionally annotated. A total of 599,286 transcripts were obtained, which were clustered into 515,196 unigenes. Results of the transcriptional effect level index revealed that direct exposure to the toxic algae A. minutum caused greater alterations in the transcriptome than did exposure to the extracellular product GTX2/3. Mechanisms of effects were different between exposure of A. salina to A. minutum cells or GTX2/3. Exposure to A. minutum modulated formation of the ribonucleoprotein complex and metabolism of amino acids and lipids in A. salina. Exposure to GTX2/3 exposure inhibited expression of genes related to metabolism of chitin, which might result in disruption of molting process or disturbed sheath morphogenesis. Overall, effects on transcription observed in this study represent the first report based on application of next generation sequencing techniques to investigate the transcriptomic response of A. salina exposed to an environmentally realistic level of A. minutum or GTX2/3.

Evolutionary adaptation analysis of immune defense and hypoxia tolerance in two closely related Marsupenaeus species based on comparative transcriptomics

Kuruma shrimp, a major farmed shrimp species in the world, includes two cryptic or sibling species, Form I (Marsupenaeus japonicus) and Form II (Marsupenaeus pulchricaudatus). Due to the lack of genomic resources, little is known about the molecular mechanisms associated with immune defense and hypoxia tolerance. Here, we sequenced the transcriptomes of two closely related Marsupenaeus species and compared genomic divergence. This study obtained 77049 and 84561 unigenes with N50 values of 1281bp and 1244bp for M. japonicus and M. pulchricaudatus, respectively, and 5036 pairs of putative orthologs were identified between two Marsupenaeus species. Estimation of Ka/Ks ratios indicated that 165 orthologous genes may be under positive selection (Ka/Ks > 0.5), including 49 pairs with a Ka/Ks ratio >1. According to the peak of synonymous rates, the divergence time between M. japonicus and M. pulchricaudatus was about 0.26-0.69 Mya. These positively selected orthologous genes related to the immune process mainly comprised single VWC domain protein, legumain, ras-related C3 botulinum, caspase, C-type lectin and were enriched in functions related to immune (Toll-like receptor and PI3K-Akt signaling) and hypoxia signaling (HIF-1 signaling and VEGF signaling). In this study, dozens of caspase-like unigenes were screened from two Marsupenaeus transcriptomes. Among these, the PjCaspase orthologous gene was subjected to positive selection (Ka/Ks = 1.22), which had different secondary and three-dimensional structure prediction. Based on the single copy caspase gene, eight populations of Marsupenaeus species were divided into two phylogeographic lineages from the East and South China. We characterized the transcriptomes of the two Marsupenaeus species and obtained several key orthologs associated with immune defense and hypoxia tolerance, which provides new insights into the immunity and genetic divergence of the two varieties. Moreover, this study will facilitate further comparative genomic studies of the two varieties.

The genome of the golden apple snail Pomacea canaliculata provides insight into stress tolerance and invasive adaptation

Background

The golden apple snail (Pomacea canaliculata) is a freshwater snail listed among the top 100 worst invasive species worldwide and a noted agricultural and quarantine pest that causes great economic losses. It is characterized by fast growth, strong stress tolerance, a high reproduction rate, and adaptation to a broad range of environments.

Results

Here, we used long-read sequencing to produce a 440-Mb high-quality, chromosome-level assembly of the P. canaliculata genome. In total, 50 Mb (11.4%) repeat sequences and 21,533 gene models were identified in the genome. The major findings of this study include the recent explosion of DNA/hAT-Charlie transposable elements, the expansion of the P450 gene family, and the constitution of the cellular homeostasis system, which contributes to ecological plasticity in stress adaptation. In addition, the high transcriptional levels of perivitelline genes in the ovary and albumen gland promote the function of nutrient supply and defense ability in eggs. Furthermore, the gut metagenome also contains diverse genes for food digestion and xenobiotic degradation.

Conclusions

These findings collectively provide novel insights into the molecular mechanisms of the ecological plasticity and high invasiveness.

AmpuBase: a transcriptome database for eight species of apple snails (Gastropoda: Ampullariidae)

BACKGROUND

Gastropoda, with approximately 80,000 living species, is the largest class of Mollusca. Among gastropods, apple snails (family Ampullariidae) are globally distributed in tropical and subtropical freshwater ecosystems and many species are ecologically and economically important. Ampullariids exhibit various morphological and physiological adaptations to their respective habitats, which make them ideal candidates for studying adaptation, population divergence, speciation, and larger-scale patterns of diversity, including the biogeography of native and invasive populations. The limited availability of genomic data, however, hinders in-depth ecological and evolutionary studies of these non-model organisms.

RESULTS

Using Illumina Hiseq platforms, we sequenced 1220 million reads for seven species of apple snails. Together with the previously published RNA-Seq data of two apple snails, we conducted de novo transcriptome assembly of eight species that belong to five genera of Ampullariidae, two of which represent Old World lineages and the other three New World lineages. There were 20,730 to 35,828 unigenes with predicted open reading frames for the eight species, with N50 (shortest sequence length at 50% of the unigenes) ranging from 1320 to 1803 bp. 69.7% to 80.2% of these unigenes were functionally annotated by searching against NCBI's non-redundant, Gene Ontology database and the Kyoto Encyclopaedia of Genes and Genomes. With these data we developed AmpuBase, a relational database that features online BLAST functionality for DNA/protein sequences, keyword searching for unigenes/functional terms, and download functions for sequences and whole transcriptomes.

CONCLUSIONS

In summary, we have generated comprehensive transcriptome data for multiple ampullariid genera and species, and created a publicly accessible database with a user-friendly interface to facilitate future basic and applied studies on ampullariids, and comparative molecular studies with other invertebrates.

RNA seq- and DEG-based comparison of developmental toxicity in fish embryos of two species exposed to Iranian heavy crude oil

To determine and compare the toxic effects of Iranian heavy crude oil (IHCO) on the embryonic development of two fish species, we examined transcriptome profiles using RNA-seq. The assembled contigs were 66,070 unigenes in olive flounder embryos and 76,498 unigenes in spotted seabass embryos. In the differential gene expression (DEG) profiles, olive flounder embryos showed different up- and down-regulated patterns than spotted seabass embryos in response to fresh IHCO (FIHCO) and weathered IHCO (WIHCO). In this work, we categorized DEG profiles into six pathways: ribosome, oxidative phosphorylation, Parkinson's disease, Alzheimer's disease, Huntington's disease, and cardiac muscle contraction, validating the expression patterns of 13 DEGs using real-time quantitative RT-PCR. The expression of the CYP1A, CYP1B1, and CYP1C1 genes in spotted seabass embryos was higher than in olive flounder embryos, whereas genes related to cell processing, development, and the immune system showed the opposite trend. Orthologous gene cluster analysis showed that olive flounder embryos were sensitive (fold change of genes with cutoff P<0.05) to both FIHCO and WIHCO, but spotted seabass embryos exhibited higher sensitivity to WIHCO than FIHCO, indicating that species-specific differences are likely to be reflected in population levels after oil spills. Overall, our study provides new insight on the different embryonic susceptibilities of two marine fish species to FIHCO and WIHCO and a better understanding of the underlying molecular mechanisms via RNA-seq and DEGs.

Comparative profiling of hepatopancreas transcriptomes in satiated and starving Pomacea canaliculata

Background

Although Pomacea canaliculata is native to South and Central America, it has become one of the most abundant invasive species worldwide and causes extensive damage to agriculture and horticulture. Conventional physical and chemical techniques have been used to eliminate P. canaliculata, but the effects are not ideal. Therefore, it is important to devise a new method based on a greater understanding of the biology of P. canaliculata. However, the molecular mechanisms underlying digestion and absorption in P. canaliculata are not well understood due to the lack of available genomic information for this species, particularly for digestive enzyme genes.

Results

In the present study, hepatopancreas transcriptome sequencing produced over 223 million high-quality reads, and a global de novo assembly generated a total of 87,766 unique transcripts (unigenes), of which 19,942 (22.7%) had significant similarities to proteins in the UniProt database. In addition, 296,675 annotated sequences were associated with Gene Ontology (GO) terms. A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment was performed for the unique unigenes, and 262 pathways (p-value < 10⁻⁵) in P. canaliculata were found to be predominantly related to plant consumption and coarse fiber digestion and absorption. These transcripts were classified into four large categories: hydrolase, transferase, isomerase and cytochrome P450. The Reads Per Kilobase of transcript per Million mapped reads (RPKM) analysis showed that there were 523 down-regulated unigenes and 406 up-regulated unigenes in the starving apple snails compared with the satiated apple snails. Several important genes are associated with digestion and absorption in plants: endo-beta-1, 4-glucanase, xylanase, cellulase, cellulase EGX1, cellulase EGX3 and G-type lysozyme genes were identified. The qRT-PCR results confirmed that the expression patterns of these genes (except for the longipain gene) were consistent with the RNA-Seq results.

Conclusions

Our results provide a more comprehensive understanding of the molecular genes associated with hepatopancreas functioning. Differentially expressed genes corresponding to critical metabolic pathways were detected in the transcriptome of starving apple snails compared with satiated apple snails. In addition to the cellulase gene, other genes were identified that may be important factors in plant matter metabolism in P. canaliculata, and this information has the potential to expedite the study of digestive physiology in apple snails.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0485-7) contains supplementary material, which is available to authorized users.

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.

Comparative Transcriptome and Chloroplast Genome Analyses of Two Related Dipteronia Species

Dipteronia (order Sapindales) is an endangered genus endemic to China and has two living species, D.sinensis and D. dyeriana. The plants are closely related to the genus Acer, which is also classified in the order Sapindales. Evolutionary studies on Dipteronia have been hindered by the paucity of information on their genomes and plastids. Here, we used next generation sequencing to characterize the transcriptomes and complete chloroplast genomes of both Dipteronia species. A comparison of the transcriptomes of both species identified a total of 7814 orthologs. Estimation of selection pressures using Ka/Ks ratios showed that only 30 of 5435 orthologous pairs had a ratio significantly >1, i.e., showing positive selection. However, 4041 orthologs had a Ka/Ks < 0.5 (p < 0.05), suggesting that most genes had likely undergone purifying selection. Based on orthologous unigenes, 314 single copy nuclear genes (SCNGs) were identified. Through a combination of de novo and reference guided assembly, plastid genomes were obtained; that of D. sinensis was 157,080 bp and that of D. dyeriana was 157,071 bp. Both plastid genomes encoded 87 protein coding genes, 40 tRNAs, and 8 rRNAs; no significant differences were detected in the size, gene content, and organization of the two plastomes. We used the whole chloroplast genomes to determine the phylogeny of D. sinensis and D. dyeriana and confirmed that the two species were highly divergent. Overall, our study provides comprehensive transcriptomic and chloroplast genomic resources, which will be valuable for future evolutionary studies of Dipteronia.