Next Article in Journal
Comparative Metagenomic and Metatranscriptomic Analyses Reveal the Response of Black Soldier Fly (Hermetia illucens) Larvae Intestinal Microbes and Reduction Mechanisms to High Concentrations of Tetracycline
Previous Article in Journal
Polychlorinated Biphenyls (PCBs) Impact Prostatic Collagen Density and Bladder Volume in Young Adult Mice Exposed during in Utero and Lactational Development
Previous Article in Special Issue
Change in Air Quality during 2014–2021 in Jinan City in China and Its Influencing Factors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 11
Issue 7
10.3390/toxics11070610
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole-Transcriptome Analysis on the Leaves of Rosa chinensis Jacq. under Exposure to Polycyclic Aromatic Hydrocarbons

by
Shili Tian
1,
Qingyang Liu
2,
Jingming Qu
1,
Ming Yang
1,
Qiaoyun Ma
1,
Jia Liu
1,
Peng Shao
1,* and
Yanju Liu
1,*
1
Beijing Center for Physical and Chemical Analysis, Institute of Analysis and Testing, Beijing Academy of Science and Technology, Beijing 100089, China
2
College of Biology and the Environment, Nanjing Forestry University, Nanjing 210037, China
*
Authors to whom correspondence should be addressed.
Toxics 2023, 11(7), 610; https://doi.org/10.3390/toxics11070610
Submission received: 25 June 2023 / Revised: 10 July 2023 / Accepted: 11 July 2023 / Published: 13 July 2023
(This article belongs to the Special Issue State-of-the-Art Techniques in Air Quality, Climate Change, and Pollutant Transport)
Browse Figures
Versions Notes

Abstract

:
The leaves of plants can be recommended as a cheap and sustainable environmental protection tool to mitigate PAHs with high toxicity in the ambient environment because they can serve as a reactor to remove ambient PAHs. Although previous studies have demonstrated that PAHs exhibit toxicological features, our knowledge about how ambient PAHs influence the leaves of plants is limited regarding the leaves of plants reducing ambient PAHs as a reactor. In this study, 1-year-old Rosa chinensis Jacq. with good growth potential was selected as a model plant. The leaves of Rosa chinensis Jacq. were exposed to 16 types of PAHs in the environmental concentration exposure group (0.1 μg L−1) and high-concentration exposure group (5 μg L−1) for seven days. In comparison, the leaves of Rosa chinensis Jacq. were exposed to de-ionized water and were chosen as the control group. During the exposure periods, the physiological parameters of leaves including, chlorophyll value, water content, temperature and nitrogen, were monitored using a chlorophyll meter. After 7 days of exposure, the leaves in the control and exposure groups were collected and used for whole-transcriptome analysis. Our results demonstrate that significant differentially expressed genes were observed in the leaves of Rosa chinensis Jacq. exposed to individual PAHs at 5 μg L−1 compared to the control group. These differentially expressed genes were involved in seven main pathways using bioinformatic analyses. In contrast, the levels of PAHs at environmentally relevant concentrations had negligible impacts on the physiological parameters and the gene transcription levels of the leaves of Rosa chinensis Jacq. Our results may provide direct evidence to remove ambient PAHs using terrestrial trees without considering the risk of PAHs at environmentally relevant concentrations on the leaves of terrestrial plants.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are recognized as teratogenic, carcinogenic and mutagenic traits by the World Health Organization (WHO) [1,2]. PAHs originate from the incomplete combustion of fossil fuel and hydrogen-containing substances, including wood, crop straw, tobacco, etc. Most PAHs are emitted from mobile sources, biomass burning and coal combustion [3,4,5]. The emissions of PAHs in China account for ~20% of the total PAH emissions globally [6]. In North China, mean levels of particulate BaP (a representative PAH) vary from 1.1 to 14.3 ng m−3 annually, which is higher than the recommended threshold value of 1.0 ng m−3 by the WHO [7,8]. The high BaP concentration in particle matter could pose a health risk to public health and exhibit ecological toxicological effects on the ecosystem [9,10]. It is estimated that the inhalation of particulate BaP ranging from 1.1 to 14.3 ng m−3 could lead to an incremental lifetime cancer risk (ILCR) as high as 3.1 × 10−5, greater than the recommended safe level (10−6) [8].
The ecological remediation strategy is recommended as a cheap and sustainable environmental protection tool for mitigating PAHs with high toxicity in the ambient environment [11,12,13]. Trees could serve as filters for PAHs, and the leaves of plants could be treated as a reactor to remove ambient PAHs [14,15,16]. A previous study indicated that trees can remove approximately 710,000 tons of air pollutants, including NO2, SO2 and air pollutants, in particulate matter in 55 cities in the United States [17,18,19]. It is estimated that plant leaves can absorb 1.71 ± 0.05 g m−2 of air pollutants in particulate matter per week, from 1960 to 2016 based on 150 field studies across 15 countries [20,21]. Klingberg et al. [12] studied the PAH accumulation in Quercus palustris and Pinus nigra in the urban landscape of Gothenburg, Sweden. They found a strong association between gaseous PAH concentrations in leaves and the air. The concentrations of PAH were observed to be higher in 3-year-old black pine needles relative to those in 1-year-old black pine needles. In oak leaves, a significant decrease in the concentration of low-molecular-mass PAHs was found between June and September. In contrast, there was a significant increase in high-molecular-mass PAHs from June to September. Once PAHs are released into the atmosphere, they can reach out to the surfaces of leaves in the form of gas and in a particulate state [22]. The waxes of leaves are considered to be essential substances for ambient PAH enrichment [23,24]. Yang et al. [22] demonstrated the relationship between foliage uptake and the inner-leaf translocation of PAHs by Cinnamomum camphora. They found a negative correlation between the wax contents and the total concentration of 16 PAHs. The transportation of PAHs from foliar dust to cuticular wax was the primary pathway of leaf accumulation. The values of the translocation factor for PAHs from foliar dust to cuticular wax were observed to be highly dependent on an increasing tendency of low-molecular-weight PAHs and a decreasing tendency of high-molecular-weight PAHs. Tian et al. [24] assessed the differences in the uptake and accumulation of PAHs by leaves across eight plants in Shanghai, China. They showed the differences in the ability to uptake and absorb PAHs across eight plants due to variations in the morphology and physiological characteristics of leaves. The retention of low-molecular-weight PAHs in eight plant leaves was associated with leaf morphology and physiological characteristics, including surface roughness, stomatal density, polar components, etc. For medium- and high-molecular-weight PAHs, wax content and adsorption were found to be the dominant factors in the accumulation ability of eight plants.
Most PAHs can enter plant leaves through the stomata [23]. Lipid components in the waxes of leaves can interact with ambient PAHs through adsorption [25]. Then, the absorbed PAHs can penetrate the surfaces of leaves through the stomata and migrate to the internal tissues of the leaves [23]. Prigioniero et al. [25] used optical microscopy and infrared spectroscopy to explore the relationships between the uptake rates of PAHs and leaf surface functional traits in four Mediterranean evergreen trees, including Chamaerops humilis, Citrus × aurantium, Magnolia grandiflora and Quercus ilex during adry month. They suggested that cutin in the dewaxed leaves of evergreen trees is the main contributor to the uptake of PAHs.
Prior studies have documented the mitigation effects of PAHs by trees and the processes of PAHs from the ambient environment to the internal tissues of leaves [26,27,28]. De Nicola et al. [27] compared the ability to capture PAHs from ambient air between evergreen (P. pinaster) and deciduous trees (Q. robur). The results of the findings indicated that the uptake of PAHs increases with the diameter at the breast height of trees. The evergreen tree (P. pinaster) enhanced the ability to capture higher amounts of low- and medium-molecular-weight PAHs compared to that of the deciduous tree (Q. robur). Anna et al. [28] linked the ability to remove PAHs and the canopy water storage capacity of coniferous trees. They measured 18 types of PAHs in the needles of three coniferous trees, including Pinus sylvestris L., Picea abies (L.) H. Karst and Abies alba mill. Samples were collected in three locations in Czarna Rózga Forest Reserve, Poland. They found that increases in the canopy water storage capacity were associated with the total content of hydrophobic PAHs in needles within the same species. This study demonstrated that the canopy water storage capacity is an essential parameter for capturing hydrophobic PAHs.
However, limited studies have illustrated the toxicological effects of PAHs on plant leaves during the migrating processes from the ambient environment to the internal tissues [21]. Thus, this study aims to investigate the responses of plant leaves under exposure to PAHs at an environmentally relevant level and higher level withwhole-transcriptome analysis. This study could accumulate evidence on the physiological and pathological responses of plants under exposure to ambient PAHs.

2. Experimental Section

2.1. Exposure Experiments

In this experiment, 1-year-old Rosa chinensis Jacq. with good growth potential and no diseases were selected. Each individual of Rosa chinensis Jacq. was planted in a 15 cm-diameter flowerpot with a height of 15 cm. The cultivation environment had a photoperiod of 12/12 h (light/darkness). The exposure experiments included three groups, which were the control group, the low-concentration exposure group and the high-concentration exposure group. The exposure chemical reagents were a mixture of 16 polycyclic aromatic hydrocarbons [13]. The 16 types of PAHs were naphthalene, anthracene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, pyrene, benzo (a) anthracene, benzo (b) fluoranthene, benzo (k) fluoranthene, picene, indeno (1, 2, 3-cd) pyrene, benzo (g, h, i)perylene, dibenzo (a, i)pyrene and dibenzo (a, h)pyrene. We used the mixtures of 16 PAHs as the exposure chemicals because these 16 types of PAHs are frequently detected in the ambient environment [2,14,15,16]. The choice of the mixtures of PAHs rather than individual PAHs with high toxicity in the exposure experiment was in accord with the real environmental conditions [11,12,13]. An amount of 200 μg mL−1 of PAH mixture dissolved in a methanol solution was selected as the standard solution. The levels of the low-concentration exposure group and high-concentration exposure group were set to be 0.1 and 5 μg L−1 for individual PAHs, respectively. The level of 5 μg L−1 for individual PAHs was prepared using 100 μL of PAHs at 200 μg mL−1 and 4000 mL of de-ionized water. Then, the level of 0.1 μg L−1 for individual PAHs was prepared using 1 mL of PAHs at 5 μg mL−1 and 50 mL of de-ionized water. Thus, due to methanol being used as the solvent in the standard solution, it was estimated that the levels of methanol were 25 μL L−1 in the low-concentration exposure group and 0.5 μL L−1 in the high-concentration exposure group. The production of methanol in leaves occurred with leaf development and expansion [29]. It was found that the methanol emission of leaves ranged from 10.0 to 26.8 μg g−1 [29]. Because methanol existed in the leaves, it was expected that trace amounts of methanol existed in the exposure solution, ranging from 0.5 to 25 μL L−1, which have negligible effects on the physiological conditions of the leaves. Each group included four replicates.
Before the exposure experiment, the leaves of Rosa chinensis Jacq. were washed with de-ionized water three times. Then, in the low-concentration exposure group, the levels of individual PAHs at 0.1 μg L−1 were sprayed on the leaves of Rosa chinensis Jacq. once a day, and in the high-concentration exposure group, 5 μg L−1 of individual PAHs was sprayed on the leaves once per day. The control group was carried out with the same procedure using de-ionized water on four replicated healthy leaves with similar ages and leaf areas [30]. During the experiment, the levels of chlorophyll, leaf water content and nitrogen in the leaves were measured using a chlorophyll meter (YT-YC, YunTang Technology, Shandong, China) simultaneously every day [31,32]. For each measurement, four leaf replicates in each group were measured. The mean levels of chlorophyll, leaf water content, and nitrogen in the leaves were treated as the daily mean concentrations.Prior studies have shown that the overall degradation half-lives of PAHs on leaves varies from approximately 0.6 to 7 days at room temperature [33]. It is expected that the trace levels of PAHs are decomposed due to the photo- and biodegradation effect within one week. To conduct the exposure experiment in real environmental conditions, the entire exposure experiment was preceded by 7 days at a temperature ranging from 18 to 20 °C and humidity ranging from 40 to 50%. After the exposure experiment, 5 g of leaves of each sample from the low-concentration exposure group and high-concentration exposure group was collected and stored at −80 °C for transcriptomic sequencing analysis.

2.2. Transcriptome Sequencing

RNA quantification and the qualification of RNA integrity in leaves were assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA) [34]. Total RNA was used as input material for the RNA sample preparations. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under an elevated temperature in First-Strand Synthesis Reaction Buffer (5X). First-strand cDNA was synthesized using a random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-). Second-strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After the adenylation of the 3′ends of DNA fragments, adaptors with hairpin loop structures were ligated to prepare for hybridization. To select cDNA fragments that were preferentially 370~420 bp in length, the library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, MA, USA). Then, PCR was performed with Phusion High-Fidelity DNA polymerase, universal PCR primers and index (X) primer. Last, PCR products were purified (AMPure XP system), and library quality was assessed on the Agilent Bioanalyzer 2100 system. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia, San Diego, CA, USA) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Novaseq platform, and 150 bp paired-end reads were generated.

2.3. Bioinformatic Analyses

Raw data (raw reads) in the fastq format were firstly processed through fastp software. In this step, clean data (clean reads) were obtained by removing reads containing adapters, reads containing 1 ploy-N and low-quality reads from raw data. Moreover, the Q20, Q30 and GC content of the clean data were calculated. All the downstream analyses were based on clean data with high quality.
The reference genome and gene model annotation files were downloaded from the genome website directly. The index of the reference genome was built using Hisat2 v2.0.5, and paired-end clean reads were aligned to the reference genome using Hisat2 v2.0.5. We selected Hisat2 as the mapping tool because Hisat2 can generate a database of splice junctions based on the gene model annotation file and thus a better mapping result than those of other non-splice mapping tools.
The mapped reads of each sample were assembled via StringTie using a reference-based approach [34,35]. StringTie uses a novel network flow algorithm as well as an optional de novo assembly step to assemble and quantify full-length transcripts representing multiple splice variants for each gene locus.
Counts v1.5.0-p3 was used to count the read numbers mapped to each gene. Then, the FPKM of each gene was calculated based on the length of the gene and the read count mapped to each gene. The FPKM, the expected number of Fragments Per Kilobase of transcript sequence per million base pairs sequenced, considers the effect of the sequencing depth and gene length for the read count at the same time and is currently the most commonly used method for estimating gene expression levels.
The differential expression analysis of two conditions/groups (two biological replicates per condition) was performed using the DESeq2 R package. DESeq2 provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The resulting p-values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted p-value ≤ 0.05 found via DESeq2 were assigned as differentially expressed.
The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource used for understanding the high-level functions and utilities of the biological system, such as the cell, the organism and the ecosystem, from molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies (http://www.genome.jp/kegg/) (accessed on 25 December 2022) [36]. We used the cluster Profiler R package to test the statistical enrichment of differentially expressed genes in KEGG pathways.

2.4. Statistical Analysis

The summary statistics for data on the levels of chlorophyll, leaf water content and nitrogen in the leaves are shown as the mean and standard deviation for each day using the results from the four independent experiments. The Wilcoxon rank-sum test was used to determine the mean differences in levels across different groups because the dataset has a non-normal distribution based on the assessment of the Shapiro–Wilk test.We performed a post hoctest using the Bonferroni methodology to adjust p-values for pairwise comparisons. Significant differences were considered at an adjusted p-value of <0.05. All statistical analyses were carried out using SPSS V26.0.

3. Results and Discussion

3.1. Variations in the Physiological Parameters of Leaves

Figure 1 presents the variations in the chlorophyll value, water content and nitrogen in the leaves in three groups during the experiment. For the control and exposure groups, the mean levels of chlorophyll, water content and nitrogen showed a similar variation trend, which decreased in the first two days and gradually increased on the third and fourth days due to watering. Then, the mean levels of chlorophyll, water content and nitrogen lowered steadily to the normal levels in the leaves. This finding indicates that humidity had an essential role in promoting the levels of chlorophyll and nitrogen in the leaves [37]. Watering on the leaves had a one-day delay effect on the variations in the levels of chlorophyll and nitrogen in the leaves. There were no significant differences in the levels of chlorophyll, water content and nitrogen across the low-concentration exposure group, high-concentration exposure group and control group. It is worth noting that the mean levels of chlorophyll in the low-exposure group were observed to be slightly higher than those in the control group, which may be attributable to the hormesis effects from PAHs [38]. In comparison, the mean levels of chlorophyll were mildly lower than those of the control group, which may be ascribed to the inhibition effects of high levels of PAHs [38].

3.2. Transcriptome Statistics

After the original data filtering, sequencing error check and the GC content distribution check, the clean reads were obtained for a follow-up subsequent analysis, and the data are summarized in the table below. As shown in Table 1, the levels of clean reads varied from 38.5 to 46.8 Gb, whichwere obtained from the sequencing library, with Q30 values in the range of 91.0–93.2% [39].

3.3. Differentially Expressed Genes

We used the concentrations of individual PAHs at 5 μg L−1 to investigate the responses of leaves under exposure to PAHs at a high concentration. As shown in Figure 2A, significant differences in expressed genes were observed after 7 days of exposure between the high-concentration exposure group and the control group. The numbers of upregulated and downregulated genes of the high-concentration exposure group relative to the control group were 1150 and 5502, respectively. It is noteworthy that the concentrations of individual PAHs in the exposure group were much higher than the levels of individual PAHs at the environmentally relevant concentration in this study [23,24].
After exposure to individual PAHs at 0.1 μg L−1 for seven days, differences in expressed genes between the low-concentration exposure group and the control group were found to be insignificant. The numbers of upregulated and downregulated genes were 112 and 195, respectively. Although several studies have reported that exposure to PAHs could lead to changes in the physiological and biochemical processes of leaves, our results indicate that the changes in the expression of genes in leaves resulting from exposure to individual PAHs at environmentally relevant concentrations were insignificant.

3.4. KEGG Enrichment Analysis

Subsequently, we performed KEGG enrichment analysis to showwhich genes or pathways that the differentially expressed gene cluster was in. KEGGis a comprehensive database that integrates genome, chemistry and system function information [34]. A padj of less than 0.05 was used as the threshold for significant enrichment.For the comparison between the high-concentration exposure and the control group, significant differentially expressed genes were assigned into seven main pathways, which included flavone and flavonol biosynthesis, glyoxylate and dicarboxylate metabolism, RNA polymerase, ribosome biogenesis in eukaryotes, porphyrin metabolism, photosynthesis-antenna proteins and photosynthesis (Figure 3). The photosynthesis pathway includes 24 downregulated genes, and the photosynthesis-antenna proteins pathway has 15 downregulated genes. Moreover, the porphyrin metabolismpathway has 1 upregulated and 25 downregulated genes. Ribosome biogenesis in the eukaryotespathway has 2 upregulated and 43 downregulated genes. The RNA polymerasepathway has 3 upregulated and 26 downregulated genes. The glyoxylate and dicarboxylate metabolismpathway has 4 upregulated and 27 downregulated genes.
In contrast, there were no significant differentially expressed genes in the KEGG enrichment analysis between the low-concentration exposure group and the control group (Figure 4). In general, the gene transcription levels obtained from the KEGG enrichment analysis exhibited a similar pattern and degree of alterations in comparison to the physiological results, thereby demonstrating that the responses of leaves to exposure to PAHs at environmentally relevant concentrations were negligible.
In this study, the mixtures of 16 PAHs at environmental levels, not individual PAHs with high toxicity, were used to investigate the toxicological effects on the gene transcription levels of the leaves of Rosa chinensis Jacq. because these 16 types of PAHs have been frequently detected in the ambient environment [2,11,12,13,14,15,16]. Because individual PAHs with low toxicity could not influence the toxicity of other PAH species [1,7], the findings of this study could reflect the physiological conditions of leaves under exposure to PAHs at environmental levels in real situations. Future studies may focus on the physiological conditions of leaves under exposure to individual PAHs at environmental levels, which could aid in understanding the metabolism mechanism of leaves under exposure to individual PAHs with different toxicity.

4. Conclusions

In this study, we measured the physiological parameters and gene transcription levels of the leaves of Rosa chinensis Jacq. exposedtoPAHsin the environmental concentration exposure group (0.1 μg L−1) and high-concentration exposure group (5 μg L−1). After exposure toPAHs at 5 μg L−1, significant differences in the gene transcription levels of the leaves of Rosa chinensis Jacq. were observed relative to the control group. These significant differentially expressed genes pertain to seven main pathways, which are flavone and flavonol biosynthesis, glyoxylate and dicarboxylate metabolism, RNA polymerase, ribosome biogenesis in eukaryotes, porphyrin metabolism, photosynthesis-antenna proteins and photosynthesis. On the contrary, environmentally relevantconcentrations of PAHs had no significant effect on the physiological parameters and gene transcription levels of the leaves of Rosa chinensis Jacq. during the 7-day exposure period in this study. Our results illustrate that exposure to ambient PAHs may havea negligible impact on the leaves of terrestrial plants at gene transcription levels. The findings of our study highlight the important roles of terrestrial plants in regulating ambient PAHs.

Author Contributions

Q.L. and Y.L. conceived and designed the experiments. S.T., Q.M. and J.Q. performed the experiments. Q.L. and Y.L. analyzed the data. S.T. and Q.L. drafted the manuscript, J.L., M.Y. and P.S. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Beijing Natural Science Foundation (8232025 and 8222044) and the Beijing Academy of Science and Technology Innovation Project (23CA005-05).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beriro, D.J.; Cave, M.R.; Wragg, J.; Thomas, R.; Wills, G.; Evans, F. A review of the current state of the art of physiologically-based tests for measuring human dermal in vitro bioavailability of polycyclic aromatic hydrocarbons (PAH) in soil. J. Hazard. Mater. 2016, 305, 240–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Agus, B.A.P.; Rajentran, K.; Selamat, J.; Lestari, S.D.; Umar, N.B.; Hussain, N. Determination of 16 EPA PAHs in food using gas and liquid chromatography. J. Food Compos. Anal. 2023, 116, 105038. [Google Scholar] [CrossRef]
  3. Wang, X.; Meyer, C.P.; Reisen, F.; Keywood, M.; Thai, P.K.; Hawker, D.W.; Powell, J.; Mueller, J.F. Emission Factors for Selected Semivolatile Organic Chemicals from Burning of Tropical Biomass Fuels and Estimation of Annual Australian Emissions. Environ. Sci. Technol. 2017, 51, 9644–9652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Huang, W.; Huang, B.; Bi, X.; Lin, Q.; Liu, M.; Ren, Z.; Zhang, G.; Wang, X.; Sheng, G.; Fu, J. Emission of PAHs, NPAHs and OPAHs from residential honeycomb coal briquette combustion. Energy Fuels 2014, 28, 636–642. [Google Scholar] [CrossRef]
  5. Shen, G.; Tao, S.; Wei, S.; Chen, Y.; Zhang, Y.; Shen, H.; Huang, Y.; Zhu, D.; Yuan, C.; Wang, H.; et al. Field Measurement of Emission Factors of PM, EC, OC, Parent, Nitro-, and Oxy- Polycyclic Aromatic Hydrocarbons for Residential Briquette, Coal Cake, and Wood in Rural Shanxi, China. Environ. Sci. Technol. 2013, 47, 2998–3005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Liu, W.; Xu, Y.; Zhao, Y.; Liu, Q.; Yu, S.; Liu, Y.; Wang, X.; Liu, Y.; Tao, S.; Liu, W. Occurrence, source, and risk assessment of atmospheric parent polycyclic aromatic hydrocarbons in the coastal cities of the Bohai and Yellow Seas, China. Environ. Pollut. 2019, 254, 113046. [Google Scholar] [CrossRef]
  7. Liu, W.; Wang, J.; Li, W.; Lin, N.; Liu, Q.; Xu, M.; Du, W. Unexpected increase of PAH toxicity in ambient particulate matter under the implementation of clean air action: Evidence from two megacities in northern China. Air Qual. Atmos. Health 2022, 15, 749–760. [Google Scholar] [CrossRef]
  8. Shen, G.; Du, W.; Zhuo, S.; Yu, J.; Tao, S. Improving regulations on residential emissions and non-criteria hazardous contaminants—Insights from a field campaign on ambient PM and PAHs in North China Plain. Environ. Sci. Policy 2019, 92, 201–206. [Google Scholar] [CrossRef]
  9. Chen, J.; Song, Y.; Liu, Y.; Chen, W.; Cen, Y.; You, M.; Yang, G. DBP and BaP co-exposure induces kidney injury via promoting pyroptosis of renal tubular epithelial cells in rats. Chemosphere 2023, 314, 137714. [Google Scholar] [CrossRef]
  10. Song, S.; Chen, B.; Huang, T.; Ma, S.; Liu, L.; Luo, J.; Shen, H.; Wang, J.; Guo, L.; Wu, M.; et al. Assessing the contribution of global wildfire biomass burning to BaP contamination in the Arctic. Environ. Sci. Ecotechnol. 2023, 14, 100232. [Google Scholar] [CrossRef]
  11. Udom, G.J.; Frazzoli, C.; Ekhator, O.C.; Onyena, A.P.; Bocca, B.; Orisakwe, O.E. Pervasiveness, bioaccumulation and subduing environmental health challenges posed by polycyclic aromatic hydrocarbons (PAHs): A systematic review to settle a one health strategy in Niger Delta, Nigeria. Environ. Res. 2023, 226, 115620. [Google Scholar] [CrossRef] [PubMed]
  12. Klingberg, J.; Strandberg, B.; Sjöman, H.; Taube, M.; Wallin, G.; Pleijel, H. Polycyclic aromatic hydrocarbon (PAH) accumulation in Quercus palustris and Pinus nigra in the urban landscape of Gothenburg, Sweden. Sci. Total Environ. 2022, 805, 150163. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, M.; Tian, S.; Liu, Q.; Yang, Z.; Yang, Y.; Shao, P.; Liu, Y. Determination of 31 Polycyclic Aromatic Hydrocarbons in Plant Leaves Using Internal Standard Method with Ultrasonic Extraction-Gas Chromatography-Mass Spectrometry. Toxics 2022, 10, 634. [Google Scholar] [CrossRef] [PubMed]
  14. Sari, M.F.; Esen, F.; Tasdemir, Y. Characterization, source apportionment, air/plant partitioning and cancer risk assessment of atmospheric PAHs measured with tree components and passive air sampler. Environ. Res. 2021, 194, 110508. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Y.; Zhang, Z.; Xu, Y.; Rodgers, T.F.M.; Ablimit, M.; Li, J.; Tan, F. Identifying the contributions of root and foliage gaseous/particle uptakes to indoor plants for phthalates, OPFRs and PAHs. Sci. Total Environ. 2023, 883, 163644. [Google Scholar] [CrossRef] [PubMed]
  16. Giráldez, P.; Aboal, J.R.; Fernández, J.Á.; Di Guardo, A.; Terzaghi, E. Plant-air partition coefficients for thirteen urban conifer tree species: Estimating the best gas and particulate matter associated PAH removers. Environ. Pollut. 2022, 315, 120409. [Google Scholar] [CrossRef]
  17. Nowak, D.J.; Ellis, A.; Greenfield, E.J. The disparity in tree cover and ecosystem service values among redlining classes in the United States. Landsc. Urban Plan. 2022, 221, 104370. [Google Scholar] [CrossRef]
  18. Nowak, D.J.; Hirabayashi, S.; Bodine, A.; Greenfield, E. Tree and forest effects on air quality and human health in the United States. Environ. Pollut. 2014, 193, 119–129. [Google Scholar] [CrossRef] [Green Version]
  19. Oksanen, E.; Kontunen-Soppela, S. Plants have different strategies to defend against air pollutants. Curr. Opin. Environ. Sci. Health 2021, 19, 100222. [Google Scholar] [CrossRef]
  20. Han, D.; Shen, H.; Duan, W.; Chen, L. A review on particulate matter removal capacity by urban forests at different scales. Urban For. Urban Green. 2020, 48, 126565. [Google Scholar] [CrossRef]
  21. Corada, K.; Woodward, H.; Alaraj, H.; Collins, C.M.; de Nazelle, A. A systematic review of the leaf traits considered to contribute to removal of airborne particulate matter pollution in urban areas. Environ. Pollut. 2021, 269, 116104. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, B.; Liu, S.; Liu, Y.; Li, X.; Lin, X.; Liu, M.; Liu, X. PAHs uptake and translocation in Cinnamomum camphora leaves from Shanghai, China. Sci. Total Environ. 2017, 574, 358–368. [Google Scholar] [CrossRef] [PubMed]
  23. Fellet, G.; Pošćić, F.; Licen, S.; Marchiol, L.; Musetti, R.; Tolloi, A.; Barbieri, P.; Zerbi, G. PAHs accumulation on leaves of six evergreen urban shrubs: A field experiment. Atmos. Pollut. Res. 2016, 7, 915–924. [Google Scholar] [CrossRef]
  24. Tian, L.; Yin, S.; Ma, Y.; Kang, H.; Zhang, X.; Tan, H.; Meng, H.; Liu, C. Impact factor assessment of the uptake and accumulation of polycyclic aromatic hydrocarbons by plant leaves: Morphological characteristics have the greatest impact. Sci. Total Environ. 2019, 652, 1149–1155. [Google Scholar] [CrossRef]
  25. Prigioniero, A.; Zuzolo, D.; Niinemets, Ü.; Postiglione, A.; Mercurio, M.; Izzo, F.; Trifuoggi, M.; Toscanesi, M.; Scarano, P.; Tartaglia, M.; et al. Particulate matter and polycyclic aromatic hydrocarbon uptake in relation to leaf surface functional traits in Mediterranean evergreens: Potentials for air phytoremediation. J. Hazard. Mater. 2022, 435, 129029. [Google Scholar] [CrossRef]
  26. Yli-Pelkonen, V.; Viippola, V.; Rantalainen, A.-L.; Zheng, J.; Setälä, H. The impact of urban trees on concentrations of PAHs and other gaseous air pollutants in Yanji, northeast China. Atmos. Environ. 2018, 192, 151–159. [Google Scholar] [CrossRef]
  27. De Nicola, F.; Concha Graña, E.; López Mahía, P.; Muniategui Lorenzo, S.; Prada Rodríguez, D.; Retuerto, R.; Carballeira, A.; Aboal, J.R.; Fernández, J.Á. Evergreen or deciduous trees for capturing PAHs from ambient air? A case study. Environ. Pollut. 2017, 221, 276–284. [Google Scholar] [CrossRef]
  28. Anna, K.-I.; Emanuel, G.; Anna, S.-R.; Błońska, E.; Lasota, J.; Łagan, S. Linking the contents of hydrophobic PAHs with the canopy water storage capacity of coniferous trees. Environ. Pollut. 2018, 242, 1176–1184. [Google Scholar] [CrossRef]
  29. Nemecek-Marshall, M.; MacDonald, R.C.; Franzen, J.J.; Wojciechowski, C.L.; Fall, R. Methanol Emission from Leaves (Enzymatic Detection of Gas-Phase Methanol and Relation of Methanol Fluxes to Stomatal Conductance and Leaf Development). Plant Physiol. 1995, 108, 1359–1368. [Google Scholar] [CrossRef] [Green Version]
  30. Tao, M.; Xu, Y.; Liu, Q.; Liu, Y.; Tian, S.; Schauer, J.J. Penetration of submicron amino-functionalized graphene quantum dots in plant stomata, implication for the depollution of atmospheric soot particles. Environ. Chem. Lett. 2023, 21, 1281–1286. [Google Scholar] [CrossRef]
  31. Sabaratnam, S.; Gupta, G.; Mulchi, C. Effects of nitrogen dioxide on leaf chlorophyll and nitrogen content of soybean. Environ. Pollut. 1988, 51, 113–120. [Google Scholar] [CrossRef] [PubMed]
  32. Brown, L.A.; Williams, O.; Dash, J. Calibration and characterisation of four chlorophyll meters and transmittance spectroscopy for non-destructive estimation of forest leaf chlorophyll concentration. Agric. For. Meteorol. 2022, 323, 109059. [Google Scholar] [CrossRef]
  33. Terzaghi, E.; De Nicola, F.; Cerabolini, B.E.L.; Posada-Baquero, R.; Ortega-Calvo, J.-J.; Di Guardo, A. Role of photo- and biodegradation of two PAHs on leaves: Modelling the impact on air quality ecosystem services provided by urban trees. Sci. Total Environ. 2020, 739, 139893. [Google Scholar] [CrossRef]
  34. Sangwan, R.S.; Tripathi, S.; Singh, J.; Narnoliya, L.K.; Sangwan, N.S. De novo sequencing and assembly of Centella asiatica leaf transcriptome for mapping of structural, functional and regulatory genes with special reference to secondary metabolism. Gene 2013, 525, 58–76. [Google Scholar] [CrossRef] [PubMed]
  35. He, J.; Fu, T.; Zhang, L.; Wanrong Gao, L.; Rensel, M.; Remage-Healey, L.; White, S.A.; Gedman, G.; Whitelegge, J.; Xiao, X.; et al. Improved zebra finch brain transcriptome identifies novel proteins with sex differences. Gene 2022, 843, 146803. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, F.; Fu, M.; Li, J.; Chen, L.; Feng, K.; Huang, T.; Cai, Y.-D. Analysis and prediction of protein stability based on interaction network, gene ontology, and KEGG pathway enrichment scores. Biochim. Biophys. Acta (BBA)-Proteins Proteom. 2023, 1871, 140889. [Google Scholar] [CrossRef]
  37. Chen, Y.; Cai, X.; Tang, B.; Xie, Q.; Chen, G.; Chen, X.; Hu, Z. SlERF.J2 reduces chlorophyll accumulation and inhibits chloroplast biogenesis and development in tomato leaves. Plant Sci. 2023, 328, 111578. [Google Scholar] [CrossRef]
  38. Calabrese, E.J.; Blain, R.B. Hormesis and plant biology. Environ. Pollut. 2009, 157, 42–48. [Google Scholar] [CrossRef]
  39. Zhou, S.; Yang, Q.; Song, Y.; Cheng, B.; Ai, X. Effect of Copper Sulphate Exposure on the Oxidative Stress, Gill Transcriptome and External Microbiota of Yellow Catfish, Pelteobagrusfulvidraco. Antioxidants 2023, 12, 1288. [Google Scholar] [CrossRef]
Toxics 11 00610 g001 550
Figure 1. Effects of PAH stress on the mean levels of chlorophyll value (SPAD), water content (H2O) and nitrogen element (N) in the leaves.
Figure 1. Effects of PAH stress on the mean levels of chlorophyll value (SPAD), water content (H2O) and nitrogen element (N) in the leaves.
Toxics 11 00610 g001
Toxics 11 00610 g002a 550Toxics 11 00610 g002b 550
Figure 2. Number of upregulated and downregulated genes for high-concentration exposure group relative to the control (A) and low-concentration exposure group relative to the control (B).
Figure 2. Number of upregulated and downregulated genes for high-concentration exposure group relative to the control (A) and low-concentration exposure group relative to the control (B).
Toxics 11 00610 g002aToxics 11 00610 g002b
Toxics 11 00610 g003 550
Figure 3. KEGG enrichment analysis of high-concentration exposure group relative to the control group. Scatter plot of KEGG enrichment, showing all pathways significantly enriched (p < 0.05) in each module.
Figure 3. KEGG enrichment analysis of high-concentration exposure group relative to the control group. Scatter plot of KEGG enrichment, showing all pathways significantly enriched (p < 0.05) in each module.
Toxics 11 00610 g003
Toxics 11 00610 g004 550
Figure 4. KEGG enrichment analysis of low-concentration exposure group relative to the control group. Scatter plot of KEGG enrichment, showing all pathways significantly enriched (p < 0.05) in each module.
Figure 4. KEGG enrichment analysis of low-concentration exposure group relative to the control group. Scatter plot of KEGG enrichment, showing all pathways significantly enriched (p < 0.05) in each module.
Toxics 11 00610 g004
Table
Table 1. Summary of sample sequencing data quality.
Table 1. Summary of sample sequencing data quality.
SampleLibraryRaw
Reads		Raw
Bases		Clean
Reads		Clean
Bases		Error
Rate		Q20Q30GCpct
C1FRAS220304683-1r452492346.79 G444314766.66 G0.0397.4292.7645.9
C2FRAS220304686-1r478786407.18 G468456107.03 G0.0397.3692.748.22
C3FRAS220304684-1r441270626.62 G427985026.42 G0.0397.0391.9848.1
C4FRAS220304680-1r402731306.04 G394459305.92 G0.0397.4992.945.73
L1FRAS220304681-1r427915046.42 G420582366.31 G0.0397.5493.0245.14
L2FRAS220304682-1r463997206.96 G451278206.77 G0.0397.3392.6947.96
L3FRAS220304685-1r452095266.78 G432447886.49 G0.0397.0392.1248.26
L4FRAS220304687-1r453944866.81 G440534366.61 G0.0397.492.7647.11
H1FRAS220304677-1r461233006.92 G449092226.74 G0.0397.1892.448.88
H2FRAS220304679-1r391836285.88 G385606425.78 G0.0397.4492.7245.4
H3FRAS220304676-1r422994686.34 G413569726.2 G0.0397.2892.4545.95
H4FRAS220304678-1r394504725.92 G385463985.78 G0.0396.6591.0345.57
Raw reads: Number of reads in the original data. Raw bases: Number of bases of the raw data (raw base = raw reads × 150 bp). Clean reads: Number of reads filtered from raw data. Clean bases: Number of bases filtered from the raw data (clean base = clean reads × 150 bp). Error rate: Overall data sequencing error rate. Q20: Percentage of bases with a Phred value greater than 20 in the total. Q30: Percentage of bases with a Phred value greater than 30 in the total. GCpct: Percentages of G and C in the four bases in clean reads.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, S.; Liu, Q.; Qu, J.; Yang, M.; Ma, Q.; Liu, J.; Shao, P.; Liu, Y. Whole-Transcriptome Analysis on the Leaves of Rosa chinensis Jacq. under Exposure to Polycyclic Aromatic Hydrocarbons. Toxics 2023, 11, 610. https://doi.org/10.3390/toxics11070610

AMA Style

Tian S, Liu Q, Qu J, Yang M, Ma Q, Liu J, Shao P, Liu Y. Whole-Transcriptome Analysis on the Leaves of Rosa chinensis Jacq. under Exposure to Polycyclic Aromatic Hydrocarbons. Toxics. 2023; 11(7):610. https://doi.org/10.3390/toxics11070610

Chicago/Turabian Style

Tian, Shili, Qingyang Liu, Jingming Qu, Ming Yang, Qiaoyun Ma, Jia Liu, Peng Shao, and Yanju Liu. 2023. "Whole-Transcriptome Analysis on the Leaves of Rosa chinensis Jacq. under Exposure to Polycyclic Aromatic Hydrocarbons" Toxics 11, no. 7: 610. https://doi.org/10.3390/toxics11070610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop