Transcriptome Analysis Reveals Changes in Whole Gene Expression, Biological Process, and Molecular Functions Induced by Nickel in Jack Pine (Pinus banksiana)

Understanding the genetic response of plants to nickel stress is a necessary step to improving the utility of plants in environmental remediation and restoration. The main objective of this study was to generate whole genome expression profiles of P. banksiana exposed to nickel ion toxicity compared to reference genotypes. Pinus banksiana seedlings were screened in a growth chamber setting using a high concentration of 1600 mg of nickel per 1 kg of soil. RNA was extracted and sequenced using the Illumina platform, followed by de novo transcriptome assembly. Overall, 25,552 transcripts were assigned gene ontology. The biological processes in water-treated samples were analyzed, and 55% of transcripts were distributed among five categories: DNA metabolic process (19.3%), response to stress (13.3%), response to chemical stimuli (8.7%), signal transduction (7.7%) and response to biotic stimulus (6.0%). For molecular function, the highest percentages of genes were involved in nucleotide binding (27.6%), nuclease activity (27.3%) and kinase activity (10.3%). Sixty-two percent of genes were associated with cellular compartments. Of these genes, 21.7% were found in the plasma membrane, 16.1% in the cytosol, 12.4% with the chloroplast and 11.9% in the extracellular region. Nickel ions induced changes in gene expression, resulting in the emergence of differentially regulated categories. Overall, there were significant changes in gene expression with a total 4128 genes upregulated and 3754 downregulated genes detected in nickel-treated genotypes compared to water-treated control plants. For biological processes, the highest percentage of upregulated genes in plants exposed to nickel were associated with the response to stress (15%), the response to chemicals (11,1%), carbohydrate metabolic processes (7.4%) and catabolic processes (7.4%). The largest proportions of downregulated genes were associated with the biosynthetic process (21%), carbohydrate metabolic process (14.3%), response to biotic stimulus (10.7%) and response to stress (10.7%). For molecular function, genes encoding for enzyme regulatory and hydrolase activities represented the highest proportion (61%) of upregulated gene. The majority of downregulated genes were involved in the biosynthetic processes. Overall, 58% of upregulated genes were located in the extracellular region and the nucleus, while 42% of downregulated genes were localized to the plasma membrane and 33% to the extracellular region. This study represents the first report of a transcriptome from a conifer species treated with nickel.