De novo characterization of Larix gmelinii (Rupr.) Rupr. transcriptome and analysis of its gene expression induced by jasmonates

Analysis of the transcriptome of the needles and bark of Pinus radiata induced by bark stripping and methyl jasmonate

BACKGROUND

Plants are attacked by diverse insect and mammalian herbivores and respond with different physical and chemical defences. Transcriptional changes underlie these phenotypic changes. Simulated herbivory has been used to study the transcriptional and other early regulation events of these plant responses. In this study, constitutive and induced transcriptional responses to artificial bark stripping are compared in the needles and the bark of Pinus radiata to the responses from application of the plant stressor, methyl jasmonate. The time progression of the responses was assessed over a 4-week period.

RESULTS

Of the 6312 unique transcripts studied, 86.6% were differentially expressed between the needles and the bark prior to treatment. The most abundant constitutive transcripts were related to defence and photosynthesis and their expression did not differ between the needles and the bark. While no differential expression of transcripts were detected in the needles following bark stripping, in the bark this treatment caused an up-regulation and down-regulation of genes associated with primary and secondary metabolism. Methyl jasmonate treatment caused differential expression of transcripts in both the bark and the needles, with individual genes related to primary metabolism more responsive than those associated with secondary metabolism. The up-regulation of genes related to sugar break-down and the repression of genes related with photosynthesis, following both treatments was consistent with the strong down-regulation of sugars that has been observed in the same population. Relative to the control, the treatments caused a differential expression of genes involved in signalling, photosynthesis, carbohydrate and lipid metabolism as well as defence and water stress. However, non-overlapping transcripts were detected between the needles and the bark, between treatments and at different times of assessment. Methyl jasmonate induced more transcriptional responses in the bark than bark stripping, although the peak of expression following both treatments was detected 7 days post treatment application. The effects of bark stripping were localised, and no systemic changes were detected in the needles.

CONCLUSION

There are constitutive and induced differences in the needle and bark transcriptome of Pinus radiata. Some expression responses to bark stripping may differ from other biotic and abiotic stresses, which contributes to the understanding of plant molecular responses to diverse stresses. Whether the gene expression changes are heritable and how they differ between resistant and susceptible families identified in earlier studies needs further investigation.

Transcriptome Analysis Reveals Differentially Expressed Genes That Regulate Biosynthesis of the Active Compounds with Methyl Jasmonate in Rosemary Suspension Cells

To study the effects of Methyl jasmonates (MeJA) on rosemary suspension cells, the antioxidant enzymes’ change of activities under different concentrations of MeJA, including 0 (CK), 10 (M10), 50 (M50) and 100 μM MeJA (M100). The results demonstrated that MeJA treatments increased the activities of phenylalanine ammonla-lyase (PAL), superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) and polyphenol oxidase (PPO) and reduced the contents of hydrogen peroxide (H₂O₂) and malondialdehyde (MDA), thus accelerating the ROS scavenging. Comparative transcriptome analysis of different concentrations of MeJA showed that a total of 7836, 6797 and 8310 genes were differentially expressed in the comparisons of CKvsM10, CKvsM50, CKvsM100, respectively. The analysis of differentially expressed genes (DEGs) showed phenylpropanoid biosynthesis, vitamin B6, ascorbate and aldarate metabolism-related genes were significantly enriched. The transcripts of flavonoid and terpenoid metabolism pathways and plant hormone signal transduction, especially the jasmonic acid (JA) signal-related genes, were differentially expressed in CKvsM50 and CKvsM100 comparisons. In addition, the transcription factors (TFs), e.g., MYC2, DELLA, MYB111 played a key role in rosemary suspension cells under MeJA treatments. qRT-PCR of eleven DEGs showed a high correlation between the RNA-seq and the qRT-PCR result. Taken together, MeJA alleviated peroxidative damage of the rosemary suspension cells in a wide concentration range via concentration-dependent differential expression patterns. This study provided a transcriptome sequence resource responding to MeJA and a valuable resource for the genetic and genomic studies of the active compounds engineering in rosemary.

Concerted control of the LaRAV1-LaCDKB1;3 module by temperature during dormancy release and reactivation of larch

Dormancy release and reactivation of temperate-zone trees involve the temperature-modulated expression of cell-cycle genes. However, information on the detailed regulatory mechanism is limited. Here, we compared the transcriptomes of the stems of active and dormant larch trees, emphasizing the expression patterns of cell-cycle genes and transcription factors and assessed their relationships and responses to temperatures. Twelve cell-cycle genes and 31 transcription factors were strongly expressed in the active stage. Promoter analysis suggested that these 12 genes might be regulated by transcription factors from 10 families. Altogether, 73 cases of regulation between 16 transcription factors and 12 cell-cycle genes were predicted, while the regulatory interactions between LaMYB20 and LaCYCB1;1, and LaRAV1 and LaCDKB1;3 were confirmed by yeast one-hybrid and dual-luciferase assays. Last, we found that LaRAV1 and LaCDKB1;3 had almost the same expression patterns during dormancy release and reactivation induced naturally or artificially by temperature, indicating that the LaRAV1-LaCDKB1;3 module functions in the temperature-modulated dormancy release and reactivation of larch trees. These results provide new insights into the link between temperature and cell-cycle gene expression, helping to understand the temperature control of tree growth and development in the context of climate change.

Enhanced Biosynthesis of Chlorogenic Acid and Its Derivatives in Methyl-Jasmonate-Treated Gardenia jasminoides Cells: A Study on Metabolic and Transcriptional Responses of Cells

Chlorogenic acid and its derivatives (CQAs) are considered as important bioactive secondary metabolites in Gardenia jasminoides Ellis (G. jasminoides). However, few studies have investigated the biosynthesis and regulation of CQAs in G. jasminoides. In this study, methyl jasmonate (MeJA) was used to enhance CQAs accumulation in cultured G. jasminoides cells. Moreover, the possible molecular mechanism of MeJA-mediated accumulation of CQAs is also explored. To this end, time-course transcriptional profiles of G. jasminoides cells responding to MeJA were used to investigate the mechanism from different aspects, including jasmonate (JAs) biosynthesis, signal transduction, biosynthesis of precursor, CQAs biosynthesis, transporters, and transcription factors (TFs). A total of 57,069 unigenes were assembled from the clean reads, in which 80.7% unigenes were successfully annotated. Furthermore, comparative transcriptomic results indicated that differentially expressed genes (DEGs) were mainly involved in JAs biosynthesis and signal transduction (25 DEGs), biosynthesis of precursor for CQAs (18 DEGs), CQAs biosynthesis (19 DEGs), and transporters (29 DEGs). Most of these DEGs showed continuously upregulated expressions over time, which might activate the jasmonic acid (JA) signal transduction network, boost precursor supply, and ultimately stimulate CQAs biosynthesis. Additionally, various TFs from different TF families also responded to MeJA elicitation. Interestingly, 38 DEGs from different subgroups of the MYB family might display positive or negative regulations on phenylpropanoids, especially on CQAs biosynthesis. Conclusively, our results provide insight into the possible molecular mechanism of regulation on CQAs biosynthesis, which led to a high CQAs yield in the G. jasminoides cells under MeJA treatment.

Methyl Jasmonate Affects Photosynthesis Efficiency, Expression of HvTIP Genes and Nitrogen Homeostasis in Barley

Jasmonates modulate many growth and developmental processes and act as stress hormones that play an important role in plant tolerance to biotic and abiotic stresses. Therefore, there is a need to identify the genes that are regulated through the jasmonate signalling pathway. Aquaporins, and among them the Tonoplast Intrinsic Proteins (TIPs), form the channels in cell membranes that are responsible for the precise regulation of the movement of water and other substrates between cell compartments. We identified the cis-regulatory motifs for the methyl jasmonate (MeJA)-induced genes in the promoter regions of all the HvTIP genes, which are active in barley seedlings, and thus we hypothesised that the HvTIP expression could be a response to jasmonate signalling. In the presented study, we determined the effect of methyl jasmonate on the growth parameters and photosynthesis efficiency of barley seedlings that had been exposed to different doses of MeJA (15–1000 µM × 120 h) in a hydroponic solution. All of the applied MeJA concentrations caused a significant reduction of barley seedling growth, which was most evident in the length of the first leaf sheath and dry leaf weight. The observed decrease of the PSII parameters after the exposure to high doses of MeJA (500 µM or higher) was associated with the downregulation of HvPsbR gene encoding one of the extrinsic proteins of the Oxygen Evolving Complex. The reduced expression of HvPsbR might lead to the impairment of the OEC action, manifested by the occurrence of the K-band in an analysis of fluorescence kinetics after MeJA treatment as well as reduced photosynthesis efficiency. Furthermore, methyl jasmonate treatment caused a decrease in the nitrogen content in barley leaves, which was associated with an increased expression the four tonoplast aquaporin genes (HvTIP1;2, HvTIP2;2, HvTIP4;1 and HvTIP4;2) predicted to transport the nitrogen compounds from the vacuole to the cytosol. The upregulation of the nitrogen-transporting HvTIPs might suggest their involvement in the vacuolar unloading of ammonia and urea, which both could be remobilised when the nitrogen content in the leaves decreases. Our research provides tips on physiological role of the individual TIP subfamily members of aquaporins under methyl jasmonate action.

A transcriptome analysis uncovers Panax notoginseng resistance to Fusarium solani induced by methyl jasmonate

BACKGROUND

Panax notoginseng is a famous Chinese herbal medicine, but the root rot disease mainly caused by Fusarium solani severely reduces the yield and quality of its medicinal materials.

OBJECTIVE

The defense priming in P. notoginseng through exogenous application of signaling molecule will supply theoretical support for the exogenous regulation of disease resistance in P. notoginseng.

METHODS

In this study, the exogenous application of methyl jasmonate (MeJA) increased P. notoginseng's resistance to F. solani. Furthermore, the P. notoginseng transcriptome during F. solani infection was investigated through next-generation sequencing to uncover the resistance mechanism of P. notogingseng induced by MeJA.

RESULTS

The de novo assembly of transcriptome sequences produced 80,551 unigenes, and 36,771 of these unigenes were annotated by at least one database. A differentially expressed gene analysis revealed that a large number of genes related to terpenoid backbone biosynthesis, phenylalanine metabolism, and plant-pathogen interactions were predominantly up-regulated by MeJA. Moreover, jasmonic acid (JA) biosynthesis-related genes and the JA signaling pathway genes, such as linoleate 13S-lipoxygenase, allene oxide cyclase, allene oxide synthase, TIFY, defensin, and pathogenesis-related proteins, showed increased transcriptional levels after inoculation with F. solani. Notably, according to the gene expression analysis, JA and ethylene signaling pathways may act synergistically to positively regulate the defense responses of P. notoginseng to F. solani.

CONCLUSION

JA signaling appears to play a vital role in P. notoginseng responses to F. solani infection, which will be helpful in improving the disease resistance of P. notoginseng cultivars as well as in developing an environmentally friendly biological control method for root rot disease.

Transcriptome analysis of Hevea brasiliensis in response to exogenous methyl jasmonate provides novel insights into regulation of jasmonate-elicited rubber biosynthesis

The phytohomorne methyl jasmonate (MeJA) is known to trigger extensive reprogramming of gene expression leading to transcriptional activation of many secondary metabolic pathways. However, natural rubber is a commercially important secondary metabolite and little is known about the genetic and genomic basis of jasmonate-elicited rubber biosynthesis in rubber tree (Hevea brasiliensis). RNA sequencing (RNA-seq) of H. brasiliensis bark treated with 1 g lanolin paste containing 0.02% w/w MeJA for 24 h (M2) and 0.04% w/w MeJA for 24 h (M4) was performed. A total of 2950 and 2850 differentially expressed genes in M2 and M4 compared with control (C) were respectively detected. Key genes involved in 2-C-methyl-D-erythritol 4-phosphate, rubber biosynthesis, glycolysis and carbon fixation (Calvin cycle) pathway were found to be up-regulated by MeJA treatment. Particularly, the expression of 3-hydroxy-3-metylglutaryl coenzyme A reductase in MVA pathway was down-regulated by MeJA treatment, but the expression of farnesyl diphosphate synthase (FPS) and cis-prenyltransferase (CPT, or rubber transferase) in rubber biosynthesis pathway were up-regulated by MeJA treatment. Up-regulation of critical genes in JA biosynthesis in response to MeJA treatment exhibited the self-activation of JA biosynthesis. In addition, up-regulated genes of great regulatory importance in cross-talk between JA and other hormone signaling, and of transcriptional regulation were identified. The increased expression levels of FPS and CPT in rubber biosynthesis pathway possibly resulted in an increased latex production in rubber tree treated with MeJA. The present results provide insights into the mechanism by which MeJA activates the rubber biosynthesis and the transcriptome data can also serve as the foundation for future research into the molecular basis for MeJA regulation of other cellular processes.

Profiling methyl jasmonate-responsive transcriptome for understanding induced systemic resistance in whitebark pine (Pinus albicaulis)

RNA-seq analysis on whitebark pine needles demonstrated that methyl jasmonate (MeJA)-triggered transcriptome re-programming substantially overlapped with defense responses against insects and fungal pathogens in Pinus species, increasing current knowledge regarding induced systemic resistance (ISR) to pathogens and pests in whitebark pine. Many whitebark pine populations are in steep decline due to high susceptibility to mountain pine beetle and the non-native white pine blister rust (WPBR). Resistance, including induced systemic resistance (ISR), is not well characterized in whitebark pine, narrowing the current options for increasing the success of restoration and breeding programs. Exogenous jasmonates are known to trigger ISR by activating the plant's immune system through regulation of gene expression to produce chemical defense compounds. This study reports profiles of whitebark pine needle transcriptomes, following methyl jasmonate (MeJA) treatment using RNA-seq. A MeJA-responsive transcriptome was de novo assembled and transcriptome profiling identified a set of differentially expressed genes (DEGs), revealing 1422 up- and 999 down-regulated transcripts with at least twofold change (FDR corrected p < 0.05) in needle tissues in response to MeJA application. GO analysis revealed that these DEGs have putative functions in plant defense signalling, transcription regulation, biosyntheses of secondary metabolites, and other biological processes. Lineage-specific expression of defense-related genes was characterized through comparison with MeJA signalling in model plants. In particular, MeJA-triggered transcriptome re-programming substantially overlapped with defense responses against WPBR and insects in related Pinus species, suggesting that MeJA may be used to improve whitebark pine resistance to pathogens/pests. Our study provides new insights into molecular mechanisms and metabolic pathways involved in whitebark pine ISR. DEGs identified in this study can be used as candidates to facilitate identification of genomic variation contributing to host resistance and aid in breeding selection of elite genotypes with better adaptive fitness to environmental stressors in this endangered tree species.

Studying Secondary Growth and Bast Fiber Development: The Hemp Hypocotyl Peeks behind the Wall

Cannabis sativa L. is an annual herbaceous crop grown for the production of long extraxylary fibers, the bast fibers, rich in cellulose and used both in the textile and biocomposite sectors. Despite being herbaceous, hemp undergoes secondary growth and this is well exemplified by the hypocotyl. The hypocotyl was already shown to be a suitable model to study secondary growth in other herbaceous species, namely Arabidopsis thaliana and it shows an important practical advantage, i.e., elongation and radial thickening are temporally separated. This study focuses on the mechanisms marking the transition from primary to secondary growth in the hemp hypocotyl by analysing the suite of events accompanying vascular tissue and bast fiber development. Transcriptomics, imaging and quantification of phytohormones were carried out on four representative developmental stages (i.e., 6-9-15-20 days after sowing) to provide a comprehensive overview of the events associated with primary and secondary growth in hemp. This multidisciplinary approach provides cell wall-related snapshots of the growing hemp hypocotyl and identifies marker genes associated with the young (expansins, β-galactosidases, and transcription factors involved in light-related processes) and the older hypocotyl (secondary cell wall biosynthetic genes and transcription factors).

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.

Identification of genes regulated by histone acetylation during root development in Populus trichocarpa

Background

Histone deacetylases (HDACs) are key enzymes catalyzing the removal of acetyl groups from histones. HDACs act in concert with histone acetyltransferases (HATs) to regulate histone acetylation status, which modifies chromatin structure, affecting gene transcription and thus regulating multiple biological processes such as plant growth and development. Over a decade, certain HDACs in herbaceous plants have been deeply studied. However, functions of HDACs in woody plants are not well understood.

Results

Histone deacetylase specific inhibitor trichostatin A (TSA) was used to investigate the role of HDACs in organogenesis of roots and root development in Populus trochocarpa. The adventitious roots were regenerated and grown on medium supplemented with 0, 1, and 2.5 μM TSA. TSA treatment delayed root regeneration and inhibited primary root growth. To examine the genes modified by TSA in the regenerated roots, tag-based digital gene expression (DGE) analysis was performed using Illumina HiSeqTM 2000. Approximately 4.5 million total clean tags were mapped per library. The distinct clean tags for the three libraries corresponding to 0, 1 and 2.5 μM TSA treatment were 166167, 143103 and 153507, from which 38.45 %, 31.84 % and 38.88 % were mapped unambiguously to the unigene database, respectively. Most of the tags were expressed at similar levels, showing a < 5-fold difference after 1 μM and 2.5 μM TSA treatments and the maximum fold-change of the tag copy number was around 20. The expression levels of many genes in roots were significantly altered by TSA. A total of 36 genes were up-regulated and 1368 genes were down-regulated after 1 μM TSA treatment, while 166 genes were up-regulated and 397 genes were down-regulated after 2.5 μM TSA treatment. Gene ontology (GO) and pathway analyses indicated that the differentially expressed genes were related to many kinds of molecular functions and biological processes. The genes encoding key enzymes catalyzing gibberellin biosynthesis were significantly down-regulated in the roots exposed to 2.5 μM TSA and their expression changes were validated by using real-time PCR.

Conclusions

HDACs were required for de novo organogenesis and normal growth of populus roots. DGE data provides the gene profiles in roots probably regulated by histone acetylation during root growth and development, which will lead to a better understanding of the mechanism controlling root development.

Electronic supplementary material

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

Conifer genomics and adaptation: at the crossroads of genetic diversity and genome function

Conifers have been understudied at the genomic level despite their worldwide ecological and economic importance but the situation is rapidly changing with the development of next generation sequencing (NGS) technologies. With NGS, genomics research has simultaneously gained in speed, magnitude and scope. In just a few years, genomes of 20-24 gigabases have been sequenced for several conifers, with several others expected in the near future. Biological insights have resulted from recent sequencing initiatives as well as genetic mapping, gene expression profiling and gene discovery research over nearly two decades. We review the knowledge arising from conifer genomics research emphasizing genome evolution and the genomic basis of adaptation, and outline emerging questions and knowledge gaps. We discuss future directions in three areas with potential inputs from NGS technologies: the evolutionary impacts of adaptation in conifers based on the adaptation-by-speciation model; the contributions of genetic variability of gene expression in adaptation; and the development of a broader understanding of genetic diversity and its impacts on genome function. These research directions promise to sustain research aimed at addressing the emerging challenges of adaptation that face conifer trees.

Transcriptomic Response to Nitric Oxide Treatment in Larix olgensis Henry

Larix olgensis Henry is an important coniferous species found in plantation forests in northeastern China, but it is vulnerable to pathogens. Nitric oxide (NO) is an important molecule involved in plant resistance to pathogens. To study the regulatory role of NO at the transcriptional level, we characterized the transcriptomic response of L. olgensis seedlings to sodium nitroprusside (SNP, NO donor) using Illumina sequencing and de novo transcriptome assembly. A significant number of putative metabolic pathways and functions associated with the unique sequences were identified. Genes related to plant pathogen infection (FLS2, WRKY33, MAPKKK, and PR1) were upregulated with SNP treatment. This report describes the potential contribution of NO to disease resistance in L. olgensis as induced by biotic stress. Our results provide a substantial contribution to the genomic and transcriptomic resources for L. olgensis, as well as expanding our understanding of the involvement of NO in defense responses at the transcriptional level.

Combination of transcriptomic and metabolomic analyses reveals a JAZ repressor in the jasmonate signaling pathway of Salvia miltiorrhiza

Jasmonates (JAs) are plant-specific key signaling molecules that respond to various stimuli and are involved in the synthesis of secondary metabolites. However, little is known about the JA signal pathway, especially in economically significant medicinal plants. To determine the functions of novel genes that participate in the JA-mediated accumulation of secondary metabolites, we examined the metabolomic and transcriptomic signatures from Salvia miltiorrhiza. For the metabolome, 35 representative metabolites showing significant changes in rates of accumulation were extracted and identified. We also screened out 2131 differentially expressed unigenes, of which 30 were involeved in the phenolic secondary metabolic pathway, while 25 were in the JA biosynthesis and signal pathways. Among several MeJA-induced novel genes, SmJAZ8 was selected for detailed functional analysis. Transgenic plants over-expressing SmJAZ8 exhibited a JA-insensitive phenotype, suggesting that the gene is a transcriptional regulator in the JA signal pathway of S. miltiorrhiza. Furthermore, this transgenic tool revealed that JAZ genes have novel function in the constitutive accumulation of secondary metabolites. Based on these findings, we propose that the combined strategy of transcriptomic and metabolomic analyses is valuable for efficient discovery of novel genes in plants.

Development and Characterization of Novel EST-SSRs from Larix gmelinii and Their Cross-Species Transferability

A set of 899 L. gmelinii expression sequence tags (ESTs), available at the National Center of Biotechnology Information (NCBI), was employed to address the feasibility on development of simple sequence repeat (SSR) markers for Larch species. Totally, 634 non-redundant unigenes including 145 contigs and 489 singletons were finally identified and mainly involved in biosynthetic, metabolic processes and response to stress according to BLASTX results, gene ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) maps. Approximately 11.7% (74) unigenes contained 90 candidate SSRs, which were mainly trinucleotides (29, 32.2%) and dinucleotides (26, 28.9%). A relatively high frequency of SSRs was respectively found in the Open Reading Frame (ORF, about 54.4%) and 5'-untranslated region (5'-UTR, 31.2%), while a low frequency was observed in the 3'-untranslated region (3'-UTR, about 14.4%). Of the 45 novel EST-SSRs markers, nine were found to be polymorphic at two L. gmelinii populations. The number of alleles per locus (Na) ranged from two to four, and the observed (Ho) and expected (He) heterozygosity values were 0.200-0.733 and 0.408-0.604, respectively. The inbreeding coefficients (FIS) for all loci were more than zero except Lg41. Most of these 9EST-SSR markers were transferable to its related species L. kaempferi, L. principis-rupprechtii and L. olgensis. These novel EST-SSRs will be useful for further research on comparative genomics, genetic resources conservation and molecular breeding in larch trees.

De novo transcriptome analysis of Perna viridis highlights tissue-specific patterns for environmental studies

BACKGROUND

The tropical green-lipped mussel Perna viridis is a common biomonitor throughout the Indo-Pacific region that is used for environmental monitoring and ecotoxicological investigations. However, there is limited molecular data available regarding this species. We sought to establish a global transcriptome database from the tissues of adductor muscle, gills and the hepatopancreas of P. viridis in an effort to advance our understanding of the molecular aspects involved during specific toxicity responses in this sentinel species.

RESULTS

Illumina sequencing results yielded 544,272,542 high-quality filtered reads. After de novo assembly using Trinity, 233,257 contigs were generated with an average length of 1,264 bp and an N50 length of 2,868 bp; 192,879 assembled transcripts and 150,111 assembled unigenes were obtained after clustering. A total of 93,668 assembled transcripts (66,692 assembled genes) with putative functions for protein domains were predicted based on InterProScan analysis. Based on similarity searches, 44,713 assembled transcripts and 25,319 assembled unigenes were annotated with at least one BLAST hit. A total of 21,262 assembled transcripts (11,947 assembled genes) were annotated with at least one well-defined Gene Ontology (GO) and 5,131 assembled transcripts (3,181 assembled unigenes) were assigned to 329 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The quantity of assembled unigenes and transcripts obtained from male and female mussels were similar but varied among the three studied tissues, with the highest numbers recorded in the gills, followed by the hepatopancreas, and then the adductor muscle. Multivariate analyses revealed strong tissue-specific patterns among the three different tissues, but not between sexes in terms of expression profiles for annotated genes in various GO terms, and genes associated with stress responses and degradation of xenobiotics. The expression profiles of certain selected genes in each tissue type were further validated using real-time quantitative polymerase chain reaction assays and a similar tissue-specific trend was seen.

CONCLUSIONS

The extensive sequence data generated from this study will provide a valuable molecular resource for facilitating environmental studies with P. viridis, and highlight the importance of tissue-specific approaches in the future.