Transcriptome Profiling of Louisiana iris Root and Identification of Genes Involved in Lead-Stress Response

Transcriptomics Combined with Physiology and Metabolomics Reveals the Mechanism of Tolerance to Lead Toxicity in Maize Seedling

Lead (Pb) exposure can induce molecular changes in plants, disrupt metabolites, and impact plant growth. Therefore, it is essential to comprehend the molecular mechanisms involved in Pb tolerance in plants to evaluate the long-term environmental consequences of Pb exposure. This research focused on maize as the test subject to study variations in biomass, root traits, genes, and metabolites under hydroponic conditions under Pb conditions. The findings indicate that high Pb stress significantly disrupts plant growth and development, leading to a reduction in catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activities by 17.12, 5.78, and 19.38%, respectively. Conversely, Pb stress led to increase malondialdehyde (MDA) contents, ultimately impacting the growth of maize. The non-targeted metabolomics analysis identified 393 metabolites categorized into 12 groups, primarily consisting of organic acids and derivatives, organ heterocyclic compounds, lipids and lipid-like molecules and benzenoids. Further analysis indicated that Pb stress induced an accumulation of 174 metabolites mainly enriched in seven metabolic pathways, for example phenylpropanoid biosynthesis and flavonoid biosynthesis. Transcriptome analysis revealed 1933 shared differentially expressed genes (DEGs), with 1356 upregulated and 577 downregulated genes across all Pb treatments. Additionally, an integrated analysis identified several DEGs and differentially accumulated metabolites (DAMs), including peroxidase, alpha-trehalose, and D-glucose 6-phosphate, which were linked to cell wall biosynthesis. These findings imply the significance of this pathway in Pb detoxification. This comprehensive investigation, employing multiple methodologies, provides a detailed molecular-level insight into maize's response to Pb stress.

Pb uptake, accumulation, and translocation in plants: Plant physiological, biochemical, and molecular response: A review

Lead (Pb) is a highly toxic contaminant that is ubiquitously present in the ecosystem and poses severe environmental issues, including hazards to soil-plant systems. This review focuses on the uptake, accumulation, and translocation of Pb metallic ions and their toxicological effects on plant morpho-physiological and biochemical attributes. We highlight that the uptake of Pb metal is controlled by cation exchange capacity, pH, size of soil particles, root nature, and other physio-chemical limitations. Pb toxicity obstructs seed germination, root/shoot length, plant growth, and final crop-yield. Pb disrupts the nutrient uptake through roots, alters plasma membrane permeability, and disturbs chloroplast ultrastructure that triggers changes in respiration as well as transpiration activities, creates the reactive oxygen species (ROS), and activates some enzymatic and non-enzymatic antioxidants. Pb also impairs photosynthesis, disrupts water balance and mineral nutrients, changes hormonal status, and alters membrane structure and permeability. This review provides consolidated information concentrating on the current studies associated with Pb-induced oxidative stress and toxic conditions in various plants, highlighting the roles of different antioxidants in plants mitigating Pb-stress. Additionally, we discussed detoxification and tolerance responses in plants by regulating different gene expressions, protein, and glutathione metabolisms to resist Pb-induced phytotoxicity. Overall, various approaches to tackle Pb toxicity have been addressed; the phytoremediation techniques and biochar amendments are economical and eco-friendly remedies for improving Pb-contaminated soils.

Hybrid RNA Sequencing Strategy for the Dynamic Transcriptomes of Winter Dormancy in an Evergreen Herbaceous Perennial, Iris japonica

Japanese iris (Iris japonica) is a popular perennial ornamental that originated in China; it has a long display period and remains green outdoors throughout the year. winter dormancy characteristics contribute greatly to the evergreenness of herbaceous perennials. Thus, it is crucial to explore the mechanism of winter dormancy in this evergreen herbaceous perennial. Here, we used the hybrid RNA-seq strategy including single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies to generate large-scale Full-length transcripts to examine the shoot apical meristems of Japanese iris. A total of 10.57 Gb clean data for SMRT and over 142 Gb clean data for NGS were generated. Using hybrid error correction, 58,654 full-length transcripts were acquired and comprehensively analysed, and their expression levels were validated by real-time qPCR. This is the first full-length RNA-seq study in the Iris genus; our results provide a valuable resource and improve understanding of RNA processing in this genus, for which little genomic information is available as yet. In addition, our data will facilitate in-depth analyses of winter dormancy mechanisms in herbaceous perennials, especially evergreen monocotyledons.

Identification of differentially expressed genes and pathways in isonuclear kenaf genotypes under salt stress

Salinity is a potential abiotic stress and globally threatens crop productivity. However, the molecular mechanisms underlying salt stress tolerance with respect to cytoplasmic effect, gene expression, and metabolism pathway under salt stress have not yet been reported in isonuclear kenaf genotypes. To fill this knowledge gap, growth, physiological, biochemical, transcriptome, and cytoplasm changes in kenaf cytoplasmic male sterile (CMS) line (P3A) and its iso-nuclear maintainer line (P3B) under 200 mM sodium chloride (NaCl) stress and control conditions were analyzed. Salt stress significantly reduced leaf soluble protein, soluble sugars, proline, chlorophyll content, antioxidant enzymatic activity, and induced oxidative damage in terms of higher MDA contents in both genotypes. The reduction of these parameters resulted in a reduced plant growth compared with control. However, P3A was relatively more tolerant to salt stress than P3B. This tolerance of P3A was further confirmed by improved physio-biochemical traits under salt stress conditions. Transcriptome analysis showed that 4256 differentially expressed genes (DEGs) between the two genotypes under salt stress were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that photosynthesis, photosynthesis antenna-protein, and plant hormone signal transduction pathways might be associated with the improved NaCl stress tolerance in P3A. Conclusively, P3A cytoplasmic male sterile could be a potential salt-tolerant material for future breeding program of kenaf and can be used for phytoremediation of salt-affected soils. These data provide a valuable resource on the cytoplasmic effect of salt-responsive genes in kenaf and salt stress tolerance in kenaf.

Evaluation of drought resistance and transcriptome analysis for the identification of drought-responsive genes in Iris germanica

Iris germanica, a species with very high ornamental value, exhibits the strongest drought resistance among the species in the genus Iris, but the molecular mechanism underlying its drought resistance has not been evaluated. To investigate the gene expression profile changes exhibited by high-drought-resistant I. germanica under drought stress, 10 cultivars with excellent characteristics were included in pot experiments under drought stress conditions, and the changes in the chlorophyll (Chl) content, plasma membrane relative permeability (RP), and superoxide dismutase (SOD), malondialdehyde (MDA), free proline (Pro), and soluble protein (SP) levels in leaves were compared among these cultivars. Based on their drought-resistance performance, the 10 cultivars were ordered as follows: 'Little Dream' > 'Music Box' > 'X'Brassie' > 'Blood Stone' > 'Cherry Garden' > 'Memory of Harvest' > 'Immortality' > 'White and Gold' > 'Tantara' > 'Clarence'. Using the high-drought-resistant cultivar 'Little Dream' as the experimental material, cDNA libraries from leaves and rhizomes treated for 0, 6, 12, 24, and 48 h with 20% polyethylene glycol (PEG)-6000 to simulate a drought environment were sequenced using the Illumina sequencing platform. We obtained 1, 976, 033 transcripts and 743, 982 unigenes (mean length of 716 bp) through a hierarchical clustering analysis of the resulting transcriptome data. The unigenes were compared against the Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG, and gene ontology (GO) databases for functional annotation, and the gene expression levels in leaves and rhizomes were compared between the 20% PEG-6000 stress treated (6, 12, 24, and 48 h) and control (0 h) groups using DESeq2. 7849 and 24,127 differentially expressed genes (DEGs) were obtained from leaves and rhizomes, respectively. GO and KEGG enrichment analyses of the DEGs revealed significantly enriched KEGG pathways, including ribosome, photosynthesis, hormone signal transduction, starch and sucrose metabolism, synthesis of secondary metabolites, and related genes, such as heat shock proteins (HSPs), transcription factors (TFs), and active oxygen scavengers. In conclusion, we conducted the first transcriptome sequencing analysis of the I. germanica cultivar 'Little Dream' under drought stress and generated a large amount of genetic information. This study lays the foundation for further exploration of the molecular mechanisms underlying the responses of I. germanica to drought stress and provides valuable genetic resources for the breeding of drought-resistant plants.

De novo transcriptome characterization of Iris atropurpurea (the Royal Iris) and phylogenetic analysis of MADS-box and R2R3-MYB gene families

The Royal Irises (section Oncocyclus) are a Middle-Eastern group of irises, characterized by extremely large flowers with a huge range of flower colors and a unique pollination system. The Royal Irises are considered to be in the course of speciation and serve as a model for evolutionary processes of speciation and pollination ecology. However, no transcriptomic and genomic data are available for these plants. Transcriptome sequencing is a valuable resource for determining the genetic basis of ecological-meaningful traits, especially in non-model organisms. Here we describe the de novo transcriptome assembly of Iris atropurpurea, an endangered species endemic to Israel's coastal plain. We sequenced and analyzed the transcriptomes of roots, leaves, and three stages of developing flower buds. To identify genes involved in developmental processes we generated phylogenetic gene trees for two major gene families, the MADS-box and MYB transcription factors, which play an important role in plant development. In addition, we identified 1503 short sequence repeats that can be developed for molecular markers for population genetics in irises. This first reported transcriptome for the Royal Irises, and the data generated, provide a valuable resource for this non-model plant that will facilitate gene discovery, functional genomic studies, and development of molecular markers in irises, to complete the intensive eco-evolutionary studies of this group.

Lead exposure-induced defense responses result in low lead translocation from the roots to aerial tissues of two contrasting poplar species

To explore whether lead (Pb)-induced defense responses are responsible for the low root-to-shoot Pb translocation, we exposed saplings of the two contrasting poplar species, Populus × canescens with relatively high root-to-shoot Pb translocation and P. nigra with low Pb translocation, to 0 or 8 mM PbCl₂. Pb translocation from the roots to aboveground tissues was lower by 57% in P. nigra than that in P. × canescens. Lower Pb concentrations in the roots and aerial tissues, greater root biomass, and lower ROS overproduction in the roots were found in P. nigra than those in P. × canescens treated with Pb. P. nigra roots had higher proportions of cell walls (CWs)-bound Pb and water insoluble Pb compounds, and higher transcript levels of some pivotal genes related to Pb vacuolar sequestration, such as phytochelatin synthetase 1.1 (PCS1.1), ATP-binding cassette transporter C1.1 (ABCC1.1) and ABCC3.1 than P. × canescens roots. Pb exposure induced defense responses including increases in the contents of pectin and hemicellulose, and elevated oxalic acid accumulation, and the transcriptional upregulation of PCS1.1, ABCC1.1 and ABCC3.1 in the roots of P. nigra and P. × canescens. These results suggest that the stronger defense barriers in P. nigra roots are probably associated with the lower Pb translocation from the roots to aerial tissues, and that Pb exposure-induced defense responses can enhance the barriers against Pb translocation in poplar roots.

Transcriptome -wide modulation combined with morpho-physiological analyses of Typha orientalis roots in response to lead challenge

Lead (Pb) is a common pollutant in many environments, including in the soil, water, and/or air. Typha orientalis Presl, a large emergent aquatic plant, has been reported to function as a Pb-tolerant and Pb-accumulating plant; however, very little molecular information regarding the tolerance of T. orientalis towards Pb is known. In this study, Pb accumulation and key factors involved in the Pb stress response at different Pb concentrations were investigated. Pb was primarily accumulated in the roots and was mainly located in the cell wall and membrane systems. Differentially expressed genes (DEGs) were identified in T. orientalis roots after Pb exposure via RNA-seq analyses. In the 0.10 mM and 0.25 mM Pb²⁺-treated groups, a total of 3275 DEGs were detected relative to the control. Many of these genes were associated with oxidation-reduction processes, metal transport, protein kinase/phosphorylation, and DNA binding transcription factors, which were shown to be Pb-responsive DEGs. Mapping Kyoto Encyclopedia of Genes and Genomes (KEGG) database, "phenylpropanoid biosynthesis" was analyzed as the major pathway of the important modules of overlapping DEGs of 0.10 mM and 0.25 mM Pb²⁺ treatments. Furthermore, a lead response gene named ToLR1 with unknown function was of particular interest. The full-length of ToLR1 sequence was cloned using rapid amplification of cDNA ends (RACE) and overexpressed in Arabidopsis thaliana, which resulted in enhanced resistance to Pb stress. This is the first report providing genomic information detailing Pb responsive genes in T. orientalis. Moreover, this study provides novel insights into the molecular mechanisms underlying the response of T. orientalis and other accumulators towards Pb stress. The key genes identified in this study may serve as potential targets for genetic engineering targeting phytoremediation.

Transcriptome profiling of Fagopyrum tataricum leaves in response to lead stress

BACKGROUND

Lead (Pb) pollution is a widespread environmental problem that is harmful to living organisms. Tartary buckwheat (Fagopyrum tataricum), a member of the family Polygonaceae, exhibits short growth cycles and abundant biomass production, could be an ideal plant for phytoremediation due to its high Pb tolerance. Here, we aimed to explore the molecular basis underlying the responses of this plant to Pb stress.

RESULTS

In our study, ultrastructural localization assays revealed that Pb ions primarily accumulate in leaf vacuoles. RNA deep sequencing (RNA-Seq) of tartary buckwheat leaves was performed on two Pb-treated samples, named Pb1 (2000 mg/kg Pb (NO₃)₂) and Pb2 (10,000 mg/kg Pb (NO₃)₂), and a control (CK). A total of 88,977 assembled unigenes with 125,203,555 bases were obtained. In total, 2400 up-regulated and 3413 down-regulated differentially expressed genes (DEGs) were identified between CK and Pb1, and 2948 up-regulated DEGs and 3834 down-regulated DEGs were generated between CK and Pb2, respectively. Gene Ontology (GO) and pathway enrichment analyses showed that these DEGs were primarily associated with 'cell wall', 'binding', 'transport', and 'lipid and energy' metabolism. The results of quantitative real-time PCR (qRT-PCR) analyses of 15 randomly selected candidate DEGs and 6 regulated genes were consistent with the results of the transcriptome analysis. Heterologous expression assays in the yeast strain Δycf1 indicated that overexpressing CCCH-type zinc finger protein 14 (ZFP14) enhanced sensitivity to Pb²⁺, while 5 other genes, namely, metal transporter protein C2 (MTPC2), phytochelatin synthetase-like family protein (PCSL), vacuolar cation/proton exchanger 1a (VCE1a), natural resistance-associated macrophage protein 3 (Nramp3), and phytochelatin synthetase (PCS), enhanced the Pb tolerance of the mutant strain.

CONCLUSION

Combining our findings with those of previous studies, we generated a schematic model that shows the metabolic processes of tartary buckwheat under Pb stress. This study provides important data for further genomic analyses of the biological and molecular mechanisms of Pb tolerance and accumulation in tartary buckwheat.

De novo sequencing, assembly, and analysis of Iris lactea var. chinensis roots' transcriptome in response to salt stress

As a halophyte, Iris lactea var. chinensis (I. lactea var. chinensis) is widely distributed and has good drought and heavy metal resistance. Moreover, it is an excellent ornamental plant. I. lactea var. chinensis has extensive application prospects owing to the global impacts of salinization. To better understand its molecular mechanism involved in salt resistance, the de novo sequencing, assembly, and analysis of I. lactea var. chinensis roots' transcriptome in response to salt-stress conditions was performed. On average, 74.17% of the clean reads were mapped to unigenes. A total of 121,093 unigenes were constructed and 56,398 (46.57%) were annotated. Among these, 13,522 differentially expressed genes (DEGs) were identified between salt-treated and control samples Compared to the transcriptional level of control, 7037 DEGs were up-regulated and 6539 down-regulated. In addition, 129 up-regulated and 1609 down-regulated genes were simultaneously detected in all three pairwise comparisons between control and salt-stressed libraries. At least 247 and 250 DEGs encoding transcription factors and transporter proteins were identified. Meanwhile, 130 DEGs regarding reactive oxygen species (ROS) scavenging system were also summarized. Based on real-time quantitative RT-PCR, we verified the changes in the expression patterns of 10 unigenes. Our study identified potential salt-responsive candidate genes and increased the understanding of halophyte responses to salinity stress.

Transcriptome assembly and expression profiling of the molecular responses to cadmium toxicity in cerebral ganglia of wolf spider Pardosa pseudoannulata (Araneae: Lycosidae)

Cadmium (Cd) is a heavy metal that can cause irreversible toxicity to animals, and is an environmental pollutant in farmlands. Spiders are considered to be an excellent model for investigating the impacts of heavy metals on the environment. To date, the changes at the molecular level in the cerebral ganglia of spiders are poorly understood. Cd exposure leads to strong damage in the nervous system, such as apoptosis and necrosis of nerve cells, therefore we conducted a transcriptomic analysis of Pardosa pseudoannulata cerebral ganglia under Cd stress to profile differential gene expression (DGE). We obtained a total of 123,328 assembled unigenes, and 1441 Cd stress-associated DEGs between the Cd-treated and control groups. Expression profile analysis demonstrated that many genes involved in calcium signaling, cGMP-PKG signaling, tyrosine metabolism, phototransduction-fly, melanogenesis and isoquinoline alkaloid biosynthesis were up-regulated under Cd stress, whereas oxidative phosphorylation-related, nervous disease-associated, non-alcoholic fatty liver disease-associated, and ribosomal-associated genes were down-regulated. Here, we provide a comprehensive set of DEGs influenced by Cd stress, and heavy metal stress, and provide new information for elucidating the neurotoxic mechanisms of Cd stress in spiders.

Transcriptome-Wide Identification of Differentially Expressed Genes and Long Non-coding RNAs in Aluminum-Treated Rat Hippocampus

Aluminum (Al) is an environmental neurotoxicant with a wide exposure, but the molecular mechanism underlying its toxicity remains unclear. We used RNA sequencing (RNA-seq) in the hippocampus of Al-treated rats to identify 96 upregulated and 652 downregulated mRNAs, and 37 dysregulated long non-coding (lnc)RNAs. Gene ontology analysis showed that dysregulated genes were involved in glial cell differentiation, neural transmission, and vesicle trafficking. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed clustering of differentially expressed mRNAs and lncRNA target genes in several pathways, including the "adenosine monophosphate-activated protein kinase signaling pathway," "extracellular matrix receptor interaction," "the phosphatidylinositol 3 kinase-protein kinase B signaling pathway," and "focal adhesion" signaling pathway. RNA-seq results were validated by reverse transcription (RT)-PCR. Additionally, Al induced changes to the number and morphology of glial cells in the hippocampus of rats, as shown by glial fibrillary acidic protein (GFAP) and ionized calcium binding adaptor molecule 1 (Iba1) immunochemistry. RT-PCR and western blotting validated the significant increase in expression of glial cell-related genes GFAP and SOX10 following Al exposure compared with control rats, consistent with RNA-seq results. Collectively, these results suggest that aberrant mRNAs and lncRNAs respond to Al neurotoxicity, and that glial cell-related genes play important roles in the Al neurotoxicity mechanism. These findings provide the basis for designing targeted approaches for the treatment or prevention of Al-induced neurotoxicity.

De novo characterization of the Iris lactea var. chinensis transcriptome and an analysis of genes under cadmium or lead exposure

Iris lactea var. chinensis (I. lactea var. chinensis) is tolerant to accumulations of cadmium (Cd) and lead (Pb). In this study, the transcriptome of I. lactea var. chinensis was investigated under Cd or Pb stresses. Using the gene ontology database, 31,974 unigenes were classified into biological process, cellular component and molecular function. In total, 13,132 unigenes were involved in enriched Encyclopedia of Genes and Genomes (KEGG) metabolic pathways, and the expression levels of 5904 unigenes were significantly changed after exposure to Cd or Pb stresses. Of these, 974 were co-up-regulated and 1281 were co-down-regulated under the two stresses. The transcriptome expression profiles of I. lactea var. chinensis under Cd or Pb stresses obtained in this study provided a resource for identifying common mechanisms in the detoxification of different heavy metals. Furthermore, the identified unigenes may be used for the genetic breeding of heavy-metal tolerant plants.

Transcriptomic and physiological analyses of Medicago sativa L. roots in response to lead stress

Lead (Pb) is one of the nonessential and toxic metals that threaten the environment and human health. Medicago sativa L. is a legume with high salt tolerance and high biomass production. It is not only a globally important forage crop but is also an ideal plant for phytoremediation. However, the biological and molecular mechanisms that respond to heavy metals are still not well defined in M. sativa. In this study, de novo and strand-specific RNA-sequencing was performed to identify genes involved in the Pb stress response in M. sativa roots. A total of 415,350 unigenes were obtained from the assembled cDNA libraries, among which 5,416 were identified as significantly differentially expressed genes (DEGs) (false discovery rate < 0.005) between cDNA libraries from control and Pb-treated plants. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs showed they mainly clustered with terms associated with binding, transport, membranes, and the pathways related to signal and energy metabolism. Moreover, a number of candidate genes included antioxidant enzymes, metal transporters, and transcription factors involved in heavy metal response were upregulated under Pb stress. Quantitative real-time PCR(qRT-PCR) validation of the expression patterns of 10 randomly selected candidate DEGs were consistent with the transcriptome analysis results. Thus, this study offers new information towards the investigation of biological changes and molecular mechanisms related to Pb stress response in plants.

De Novo Assembly, Annotation, and Characterization of Root Transcriptomes of Three Caladium Cultivars with a Focus on Necrotrophic Pathogen Resistance/Defense-Related Genes

Roots are vital to plant survival and crop yield, yet few efforts have been made to characterize the expressed genes in the roots of non-model plants (root transcriptomes). This study was conducted to sequence, assemble, annotate, and characterize the root transcriptomes of three caladium cultivars (Caladium × hortulanum) using RNA-Seq. The caladium cultivars used in this study have different levels of resistance to Pythiummyriotylum, the most damaging necrotrophic pathogen to caladium roots. Forty-six to 61 million clean reads were obtained for each caladium root transcriptome. De novo assembly of the reads resulted in approximately 130,000 unigenes. Based on bioinformatic analysis, 71,825 (52.3%) caladium unigenes were annotated for putative functions, 48,417 (67.4%) and 31,417 (72.7%) were assigned to Gene Ontology (GO) and Clusters of Orthologous Groups (COG), respectively, and 46,406 (64.6%) unigenes were assigned to 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. A total of 4518 distinct unigenes were observed only in Pythium-resistant "Candidum" roots, of which 98 seemed to be involved in disease resistance and defense responses. In addition, 28,837 simple sequence repeat sites and 44,628 single nucleotide polymorphism sites were identified among the three caladium cultivars. These root transcriptome data will be valuable for further genetic improvement of caladium and related aroids.

Unveiling gene trait relationship by cross-platform meta-analysis on Chinese hamster ovary cell transcriptome

In the past few years, transcriptome analysis has been increasingly employed to better understand the physiology of Chinese hamster ovary (CHO) cells at a global level. As more transcriptome data accumulated, meta-analysis on data sets collected from various sources can potentially provide better insights on common properties of those cells. Here, we performed meta-analysis on transcriptome data of different CHO cell lines obtained using NimbleGen or Affymetrix microarray platforms. Hierarchical clustering, non-negative matrix factorization (NMF) analysis, and principal component analysis (PCA) accordantly showed the samples were clustered into two groups: one consists of adherent cells in serum-containing medium, and the other suspension cells in serum-free medium. Genes that were differentially expressed between the two clusters were enriched in a few functional classes by Database for Annotation, Visualization, and Integrated Discovery (DAVID) of which many were common with the enriched gene sets identified by Gene Set Enrichment Analysis (GSEA), including extracellular matrix (ECM) receptor interaction, cell adhesion molecules (CAMs), and lipid related metabolism pathways. Despite the heterogeneous sources of the cell samples, the adherent and suspension growth characteristics and serum-supplementation appear to be a dominant feature in the transcriptome. The results demonstrated that meta-analysis of transcriptome could uncover features in combined data sets that individual data set might not reveal. As transcriptome data sets accumulate over time, meta-analysis will become even more revealing. Biotechnol. Bioeng. 2017;114: 1583-1592. © 2017 Wiley Periodicals, Inc.

Transcriptome Profiling Analysis of Wolf Spider Pardosa pseudoannulata (Araneae: Lycosidae) after Cadmium Exposure

Pardosa pseudoannulata is one of the most common wandering spiders in agricultural fields and a potentially good bioindicator for heavy metal contamination. However, little is known about the mechanisms by which spiders respond to heavy metals at the molecular level. In the present study, high-throughput transcriptome sequencing was employed to characterize the de novo transcriptome of the spiders and to identify differentially expressed genes (DEGs) after cadmium exposure. We obtained 60,489 assembled unigenes, 18,773 of which were annotated in the public databases. A total of 2939 and 2491 DEGs were detected between the libraries of two Cd-treated groups and the control. Functional enrichment analysis revealed that metabolism processes and digestive system function were predominately enriched in response to Cd stress. At the cellular and molecular levels, significantly enriched pathways in lysosomes and phagosomes as well as replication, recombination and repair demonstrated that oxidative damage resulted from Cd exposure. Based on the selected DEGs, certain critical genes involved in defence and detoxification were analysed. These results may elucidate the molecular mechanisms underlying spiders' responses to heavy metal stress.