Different Proteome Profiles between Male and Female Populus cathayana Exposed to UV-B Radiation

With increasing altitude, solar UV-B radiation is enhanced. Based on the phenomenon of male-biased sex ratio of Populus cathayana Rehder in high altitude alpine area, we hypothesized that males have a faster and more sophisticated responsive mechanism to high UV-B radiation than that of females. Our previous studies have shown sexually different responses to high UV-B radiation were existed in P. cathayana at the morphological, physiological, and transcriptomic levels. However, the responses at the proteomic level remain unclear. In this study, an isobaric tag for relative and absolute quantification (iTRAQ)-based quantitative proteome analysis was performed in P. cathayana females and males. A total of 2,405 proteins were identified, with 331 proteins defined as differentially expressed proteins (DEPs). Among of these, 79 and 138 DEPs were decreased and 47 and 107 DEPs were increased under high solar UV-B radiation in females and males, respectively. A bioinformatics analysis categorized the common responsive proteins in the sexes as related to carbohydrate and energy metabolism, translation/transcription/post-transcriptional modification, photosynthesis, and redox reactions. The responsive proteins that showed differences in sex were mainly those involved in amino acid metabolism, stress response, and translation/transcription/post-transcriptional modification. This study provides proteomic profiles that poplars responding to solar UV-B radiation, and it also provides new insights into differentially sex-related responses to UV-B radiation.


INTRODUCTION
Natural levels of ultraviolet-B (UV-B: 280-320 nm) radiation act as an environmental regulatory factor in plants and control gene expression, growth and development (Jenkins, 2009;Hideg et al., 2013;Czégény et al., 2016). In alpine areas, with increasing altitude, UV-B radiation is enhanced. High doses of UV-B radiation negatively affect cellular/subcellular and macromolecular structures in vegetative plant tissues, including DNA, RNA, and proteins, modulate metabolites and reduce photosynthesis, biomass, and seed production (Jansen et al., 1998;Ries et al., 2000). Thus, plants must sense natural UV-B intensity and quickly respond to high levels of UV-B radiation as an acclimation process (Jenkins, 2009). A number of studies have reported that the responsive and adaptive mechanisms of plants to high altitude with high UV-B radiation differ among species, populations and genotypes (Ren et al., 2006(Ren et al., , 2007Wang et al., 2016). However, few studies have focused on the sexual differences in the responses.
In nature, the sex ratio of Populus is consistent with a 1:1 equilibrium in low-altitude habitats (1,600-2,000 m). However, at high altitude (2,000-2,600 m), there is a significant male-biased sex ratio (Wang et al., 2011;Lei et al., 2017). Therefore, we hypothesized that males are more resistant and have a more sophisticated response to high UV-B radiation than that of females, thus leading to a male-biased sex ratio at high altitude. To test this hypothesis, based on previous work, we exposed Populus cathayana Rehder to different UV-B radiation intensities and showed differences in the sexual responses to high UV-B radiation (Xu et al., 2010;Feng et al., 2014;Jiang et al., 2015). At the physiological level, males appear to be more resistant to high UV-B radiation. Observations of the organelle ultrastructure, photosynthetic rate and biomass accumulation indicated that these parameters were less affected by increased UV-B intensity in males than in females (Xu et al., 2010;Feng et al., 2014), and UV-B exposure enhanced bud break in male clones (Stromme et al., 2015). A comparative transcriptome analysis indicated that in P. cathayana males, sex-related transcriptional reprogramming occurred in certain important metabolic processes, and showed that sex-biased gene regulation under solar UV-B radiation was more responsive in males than in females (Jiang et al., 2015). At the metabolic level, UV-B radiation has been reported to mediate metabolic rearrangement in gray poplars (P. canescens syn.; Kaling et al., 2015), and increased UV radiation has been shown to promote the emission of more volatile organic compounds in females of European aspen (P. tremula L.) than in male plants (Randriamanana et al., 2015;Maja et al., 2016). However, sexual differences in the response to UV-B radiation at the proteomic level have not been reported. Because of the complicated and sophisticated modification and regulation of gene expression processes in plant cells, UV-B radiation may promote different patterns of change in the pathways at the mRNA, protein and metabolite levels. Therefore, determining how plants respond to high UV-B radiation at the leaf proteomic level will provide new insights into UV-B responsive processes and ultimately increase our understanding of the sex-ratio variation in high-altitude alpine areas. In this study, an iTRAQ-based quantitative proteome analysis was performed. The objectives were to (1) understand how UV-B radiation triggers overall proteomic changes and (2) assess whether males are more responsive than females to high UV-B radiation.

Plant Materials and Experimental Design
The material treatment and experimental design were as described by Jiang et al. (2015). Briefly, male and female cuttings were produced from F 1 individuals derived from a controlled cross between two P. cathayana genotypes with divergent phenotypes. Sixty male and female cuttings were planted in 10 L plastic pots (one plant per pot) filled with 8 kg of homogenized soil and 8 g of slow-release fertilizer (13% N, 10% P, and 14% K) in a greenhouse located in the Wanglang National Nature Reserve (32 • 98 ′ N and 104 • 08 ′ E, Supplementary Figure 1). The altitude, mean annual rainfall, annual temperature, and annual summer temperature in the area are 2,600 m, 801 mm, 2.9 and 12.7 • C, respectively.
Two UV-B radiation intensities were used: 2.1 and 8.5 kJ m −2 day −1 of average daily biologically effective UV-B (UV-BBE) radiation. In the low solar UV-B radiation treatment (2.1 kJ m −2 day −1 ), a polyester film (0.13 mm, Shanghai HiTeC Plastics Co. Ltd., Shanghai, China) was used to selectively exclude solar transmission, including UV-B and UV-C radiation. In the high solar UV-B radiation treatment (8.5 kJ m −2 day −1 ), 0.13 mm cellulose diacetate film (Qingzhou Yi-Run Agricultural Film Factory, Qingzhou, Shandong, China) was employed, which allowed the transmission of both UV-A and UV-B radiation (wavelength ≥290 nm). The spectral irradiance was weighted using the generalized plant response function normalized at 300 nm to obtain the UV-BBE radiation (Caldwell, 1971). The polyester and cellulose diacetate films were both replaced weekly. The spectral irradiance of the solar radiation at the plant level was determined using a USB2000 Fiber Optic Spectrometer (Ocean Optics, Inc., Dunedin, New Zealand) with a CC-3-UV Cosine Corrector. Prior to the measurements, the spectrometer was calibrated using a DH2000-CAL Radiometric Calibrated Deuterium Tungsten Source (210-1,050 nm, OceanOptics, Inc., Minneola, Florida, US). Figure 1 shows the high and low solar UV-B radiation intensities from June to October 2012.
The experimental layout was completely randomized with two factors: sex and UV-B radiation. Therefore, there were four treatments: (1) males exposed to low solar UV-B radiation (MC), (2) females exposed to low solar UV-B radiation (FC), (3) males exposed to high solar UV-B radiation (MU), and (4) females exposed to high solar UV-B radiation (FU). The duration of the experiments was 100 days (from 1 June to 10 August).

Gas Exchange and Chlorophyll Fluorescence Measurements
Five cuttings of the 4th fully expanded leaves from each sex per treatment were randomly selected for the gas exchange and chlorophyll fluorescence measurements. The net photosynthesis rate (P n ), stomatal conductance (g s ), intercellular CO 2 concentration (C i ), and transpiration rate (E) were measured using a LI-COR 6400 portable photosynthesis measuring system (LI-COR, Lincoln, NE, USA) between the hours of 08:30 and 11:30. Prior to conducting the measurements, the samples were illuminated with saturated PPFD provided by the LED light source of the equipment for 10 min to achieve full photosynthetic induction. The saturated photosynthetic photon flux density (PPFD) was determined by preliminary experiments, and a carbon dioxide gas cylinder (LI-COR, Lincoln, NE, USA) was used to provide constant and stable CO 2 during the measurements. A standard LI-COR leaf chamber (2 × 3 cm) was used, and the optimal parameters were as follows: leaf temperature of 25 • C, leaf-air vapor pressure deficit of 1.5 ± 0.5 kPa, CO 2 concentration of 400 ± 5 µmol mol −1 , and photosynthetic photon flux density of 1,400 mmol m −2 s −1 .
Chlorophyll fluorescence was measured using a PAM chlorophyll fluorometer (PAM 2100, Walz, Effeltrich, Germany). The leaf samples were placed in the dark for 30 min using an aluminum foil cover, followed by measurement of the minimal fluorescence (F o ) and the maximal fluorescence (F m ). Then, the leaves were illuminated with actinic light at an intensity of 250 µmol m −2 s −1 . The actinic light was removed and the minimal fluorescence (F ′ o ) and the maximal fluorescence (F ′ m ) were measured by illuminating the leaves with far-red light for 3 s. Finally, a saturating white light with a pulse of 8,000 µmol m −2 s −1 was applied for 0.8 s to measure the F m and F ′ m values. The measurements were conducted from 07:30 to 09:00. The chlorophyll fluorescence kinetics parameters (Fv/Fm, maximum efficiency of PSII; Yield, the effective quantum yield of PSII; qP, photochemical quenching coefficient; qN, non-photochemical quenching coefficient) were measured and calculated as described by Vankooten and Snel (1990).

Leaf Transmission Electron Microscopy
For transmission electron microscopy analysis, small leaf sections (2 mm in length; avoiding the midrib) were selected from the middle position of the leaves. Five biological replicates were used. Leaf sections were fixed in 2.5% (v/v) glutaralpentanedial in 0.2 M phosphate buffered saline (pH 7.0) for 3 h at 22 • C and postfixed in 2% osmium tetroxide (OsO 4 ) for 2 h. Then, the leaves were sequentially dehydrated in 30, 50, 70, and 90% acetone, respectively, and embedded in Epon 812 for ∼2 h. Ultra-thin sections (80 nm) were sliced, stained with uranyl acetate and lead citrate, and mounted on copper grids for viewing in a H-600IV TEM (Hitachi, Tokyo, Japan) at an accelerating voltage of 60.0 kV.

Hydrolyzed Amino Acid Compound Measurements
Powdered dried leaves (0.1 g) were transferred into a 20-ml hydrolysis tube with 10 ml of 6 mol l −1 hydrochloric acid. The hydrolysis tube was sealed under vacuum and transferred to a constant-temperature drier at 110 ± 2 • C for 22 h. After the tubes were removed from the dryer and cooled, the hydrolysis liquid was filtered, transferred into a 50-ml volumetric flask, and diluted with deionized water to scale. Then, 1 ml of the diluted hydrolysis liquid was withdrawn and dried in a vacuum drier at 40-50 • C. The residue was dissolved in 1 ml of sodium citrate-hydrochloric acid buffer solution (pH 2.2) for analysis. The concentration and composition of amino acids were determined by an L-8800 automatic amino acid analyzer (Hitachi, Tokyo, Japan).

Protein Extraction, Digestion, and iTRAQ Labeling
Total proteins were extracted from three biological replicates (fresh leaves from tree cuttings) in each treatment using acetone methods as previously described . The leaves were ground to a fine powder and suspended in 0.5 M triethylammonium bicarbonate (TEAB) buffer with 1 mM phenylmethyl sulfonyl fluoride and 0.1% SDS (w/v), and the samples were then sonicated for 5 min and centrifuged at 25,000 g for 20 min. The supernatant was transferred to another tube, 0.5 M TEAB buffer was added to the pellet to repeat the protein extraction, and the sample was centrifuged at 25,000 g for 20 min. The proteins in the combined supernatant were reduced (10 mM DTT, 56 • C for 60 min), alkylated (55 mM iodoacetamide, dark room temperature for 45 min), precipitated by pre-cooled acetone at −20 • C for 2 h, and then centrifuged at 25,000 g for 20 min. The pellet was washed twice with acetone, and the final pellet was dissolved in 0.5 M TEAB buffer with 0.1% SDS, sonicated for 15 min, and centrifuged at 25,000 g for 20 min. The supernatant was used for liquid digestion, and the protein concentration was determined using the Bradford assay.
The processed protein (100 µg) was removed from each sample solution and digested with Trypsin Gold (Promega, Madison, WI, USA) at a protein: trypsin ratio of 20:1 at 37 • C for 12 h. After trypsin digestion, the peptides were dried by vacuum centrifugation, followed by reconstitution in 0.5 M TEAB and processing according to the manufacturer's protocol for 8-plex iTRAQ (Applied Biosystems). In this work, three biological replicates from each treatment and sex were analyzed. Thus, 12 samples were divided into two sets for iTRAQ. Males exposed to high (MU) and low (MC) solar UV-B radiation were in set 1, and females exposed to high (FU) and low (FC) solar UV-B radiation were in set 2. The MC1 samples were labeled in both sets as the control group. The detailed labeling conditions are listed in Table 1.

Separation of Peptides by Strong Cation Exchange (SCX) and ESI Mass Spectrometric Analysis
The labeled samples were fractionated using a LC-20AB highperformance liquid chromatography (HPLC) system (Shimadzu, Kyoto, Japan) and a 4.6 × 250 mm Ultremex strong cation exchange (SCX) column (Phenomenex, Torrance, CA, USA). After reconstitution of the labeled peptide mixtures with 4 ml of buffer A [10 mM NaH 2 PO 4 in 25% acetonitrile (CAN, v/v), pH 2.6], SCX separation was performed at a flow rate of 1 ml min −1 using elution buffer A for 10 min, followed by a linear gradient of 5-60% buffer B (25 mM NaH 2 PO 4 , 1 M KCl in 25% ACN, pH 2.7) for 20 min and 100% buffer B for 2 min. The column was equilibrated with buffer A for 10 min prior to the next injection. The eluted fractions were monitored by measuring the absorbance at 214 nm, desalted with a Strata X C18 column (Phenomenex, Torrance, CA, USA), and vacuum dried.
MS analysis was performed using a Triple TOF 5600 System (AB SCIEX, Concord, ON) fitted with a Nanospray III source (AB SCIEX, Concord, ON) and a pulled quartz tip as the emitter (New Objectives, Woburn, MA). The data were acquired using an ion spray voltage of 2.5 kV, curtain gas of 30 psi, nebulizer gas of 15 psi, and an interface heater temperature of 150 • C. The MS was operated with an RP ≥30,000 FWHM for the TOF MS scans. For the IDA, survey scans were acquired in 250 ms, and as many as 30 product ion scans were collected if a threshold of 120 counts per second (counts s −1 ) was exceeded at a 2+ to 5+ charge state. The total cycle time was fixed to 3.3 s. The Q2 transmission window was 100 Da for 100%. Four time bins were summed for each scan at a pulse frequency value of 11 kHz via monitoring of the 40 GHz multichannel TDC detector with a four-anode channel detection ion. A sweeping collision energy setting of 35 ± 5 eV coupled with iTRAQ-adjusted rolling collision energy was applied to all pre-cursor ions for collisioninduced dissociation. In this experiment, full MS scans were acquired in a mass range of m/z 350 to 1,500 in a scan time of 250 ms. Fragment ion spectra were acquired in the mass range of m/z 100-2,000 and excluded for further fragmentation over 15 s.

Protein Identification and Quantification
The raw data files were converted into MGF files and then searched against a local poplar database (ftp://ftp.jgi-psf.org/ pub/compgen/phytozome/v9.0/Ptrichocarpa/annotation/) using the Mascot server (version 2.3.02, Matrix Science, Boston, MA). The data were downloaded from Populus trichocarpa V 3.0 (73013 protein-coding transcripts; https://phytozome.jgi.doe. gov/pz/portal.html#!info?alias=Org_Ptrichocarpa). The Mascot search settings were as follows: one missed cleavage site by trypsin was allowed with a fixed modification of carbamidomethyl (C), iTRAQ8plex (N-term) and iTRAQ8plex (K) and variable modifications of Gln-> pyro-Glu (N-term Q), oxidation (M), and iTRAQ8plex (Y). The fragment mass tolerance was ±0.1 Da, and the peptide mass tolerance was ±0.05 Da. The max missed cleavages value was 1. The specified false discovery rate (FDR) was 1% when automatic searching using software of pFind 2.0 (Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China). The detailed information of pFind can be referenced to Wang et al. (2007) and Li et al. (2005). For protein identification, at least 2 unique peptides identified were considered. For proteins that were not identified in the local poplar database, SwissProt database (August, 2016, including 553,474 sequences) was used for re-searching using the same parameter settings. The ratio of treatments to controls over ±1.5 and P ≤ 0.05 (t-test) were considered to indicate the differentially expressed proteins (DEPs). The proteins identified within a family were grouped in the Mascot protein family summary. The sample labeled 113 (MC1) was used as the reference (the control group) based on the weighted average of the intensity of report ions in each identified peptide. The final ratios of protein were then normalized by the median average protein ratio for the mixes of different labeled samples. This normalization corrects the systematic error. For each sex, only DEPs identified in all biological replicates (including the three controls and three UV-B treatments) were further used in the analysis of sexually differential expression. DEPs that were detected in both sexes were defined as commonly changed DEPs, and those detected in only one gender were defined as sexual DEPs. All of the DEPs were functionally categorized using Blast2go, a web-based bioinformatics tool that groups proteins based on their GO annotations. According to the molecular functions listed on the UniProt and Gene Ontology websites, the DEPs were classified into different functional categories in both sexes.

Statistical Analysis
For the physiological parameters, the effects of UV-B radiation, sex and their interaction were analyzed by an analysis of variance (ANOVA) using a randomized complete-block design in SPSS 16.0 (SPSS, Chicago, IL, USA). Prior to the analysis, the data were checked for normality and homogeneity of variance. Post-hoc comparisons were tested using Tukey's test at a significance level of P < 0.05. The mean values and standard errors were determined for each variable. For protein identification, at least 2 unique peptides identified were considered. The ratio of treatments to controls over ±1.5 and P ≤ 0.05 were considered to indicate the DEPs.

Gas Exchange and Chlorophyll Fluorescence Parameter Changes in Both Sexes
Overall, compared with the low solar UV-B radiation treatment, the high solar UV-B radiation significantly reduced the P n , g s , C i , E, qP, and chlorophyll fluorescence Yield values in females but not in males. For all parameters, the interactive effects of sex and UV-B were significant (Table 2). Additionally, the g s , C i , E, qP, and Yield values showed significant sexual differences under the high UV-B radiation treatment but not under the low UV-B radiation treatment.

Ultrastructural Morphological Changes in Both Sexes
As shown in Figure 2, under low UV-B radiation, the chloroplasts of both sexes had a typical ellipsoidal shape and typical structure, and 15-30 well-arranged thylakoids (on average 20) were included in each granum (Figures 2A,C). Under high UV-B radiation, the number of starch grains was lower (Figures 2B,D) than that in the individuals under the low UV-B conditions, and additional plastoglobuli were present in the chloroplasts. The accumulation of plastoglobuli is observed as lipid droplets derived from thylakoid degradation, and it was greater in females ( Figure 2D) than in males ( Figure 2B). However, the high intensity of UV-B radiation did not cause significant changes in the chloroplast granum number, mitochondria, or cellular membrane in the ultra-structures of the mesophyll cells.

Hydrolyzed Amino Acid Compound Changes in Both Sexes
In this study, a total of 17 types of amino acids were detected in both males and females (Table 3). High UV-B radiation significantly decreased the aspartic acid, leucine, arginine, and proline concentrations, whereas cystine exhibited smaller changes in both sexes. Interestingly, the other amino acids were significantly decreased in males but not in females under high UV-B radiation.

Protein Expression Profiles in Response to Solar UV-B Radiation in Both Sexes
To investigate the changes in the protein profiles in response to solar UV-B radiation, proteins were extracted from the poplar leaves and analyzed using an iTRAQ-based shotgun proteomics strategy. Using the Mascot search engine, a total of 6,091 peptides from the trypsin-digested proteins were identified; of these, 5,210 peptides were unique. A total of 2,405 proteins were identified. The protein information is listed in Supplementary Table 1. According to the selection criteria of expression abundance ≥1.5 (or ≤ 0.67) and P ≤ 0.05, 331 proteins were defined as DEPs. In females, high UV-B radiation was associated with the decreased abundance of 79 proteins and the increased abundance of 47 proteins (Figure 3 and Supplementary Table 2). In males, high UV-B radiation was associated with the decreased abundance of 138 proteins and the increased abundance of 107 proteins (Figure 3 and Supplementary Table 3). According to the identified protein IDs, decreases and increases were observed in 22 and 15 proteins common in both sexes, respectively (Figure 3 and Table 4).

Functional Classification of DEPs
All of the DEPs were functionally categorized using Blast2go, a web-based bioinformatics tool that groups proteins based on their GO annotations. According to the molecular functions listed on the UniProt and Gene Ontology websites, the DEPs were classified into 11 functional categories in both sexes (Figure 4). Among these, the common DEPs in both sexes were involved in carbohydrate and energy metabolism, translation/transcription/post-transcriptional modification, photosynthesis, and redox reactions ( Table 4,  Supplementary Tables 2, 3 and Figure 5). The sex-specific DEPs were mainly involved in translation/transcription/posttranscriptional modification, amino acid metabolism, and stress responses (Supplementary Tables 2, 3). Additionally, a greater number of DEPs was observed in males than in females in all the categories, particularly in the translation/transcription/post-transcriptional modification proteins.

Common Changes between the Sexes in the Proteome in Response to UV-B Radiation
The commonly changed proteins in both sexes under high UV-B radiation were mainly involved in translation/transcription/posttranscriptional modification, carbohydrate and energy metabolism, photosynthesis, and redox reactions (Supplementary Tables 2, 3). Of these, RNA-binding proteins play roles in the transport, stability, and translation of mRNA (Sahi et al., 2007). Several RNA-binding proteins (e.g., glycinerich RNA-binding protein and glycine-rich RNA-binding protein 2) were up-regulated by UV-B radiation. Glycine-rich RNA-binding protein (GRP) is involved in pre-mRNA splicing, nucleocytoplasmic mRNA transport, mRNA stability and decay, and translation (Dreyfuss et al., 1993;Simpson and Filipowicz, 1996). It plays an important role in the process of regulating gene expression, principally at the post-transcriptional level. Glycine-rich RNA-binding protein 2 may play a role in RNA transcription, which is likely because of its RNA chaperone activity during the stress adaptation process (Kim et al., 2007;Zhang et al., 2012). Numerous peptidyl-prolyl cis-trans isomerases were up-regulated in both sexes, particularly in males under UV-B radiation. The stromal 70 kDa heat shock-related proteins were also increased in abundance in both sexes. Heat shock proteins perform a chaperone function by stabilizing new proteins to ensure correct folding or by aiding the refolding of damaged proteins. The up-regulation of peptidyl-prolyl cis-trans isomerases and stromal 70 kDa heat shock-related proteins may be an adaption to avoiding protein misfolding and promoting the formation of disulfide bonds in nascent proteins under stressed conditions (Marshall and Keegstra, 1992;Schönbrunner and Schmid, 1992). Additionally, in conjunction with rRNA, ribosomal proteins compose the ribosomal subunits involved in the cellular process of translation. The evidence showed that UV-B radiation damaged ribosomes by crosslinking cytosolic and chloroplast ribosomal proteins to RNA. Ribosomal damage that accumulated during the initial stage of UV-B radiation correlated with a progressive decrease in new protein production; however, de novo synthesis of some ribosome proteins increases in the later stage of UV-B radiation (Casati and Walbot, 2004). In our study, a set of cytosolic and chloroplast ribosomal proteins (40S and 60S) was down-regulated in both sexes under UV-B radiation. This decrease may be attributable to ribosome rearrangement. Thus, these changes may enhance the efficiency of protein translation under UV-B stress. Ribosome rearrangement has been reported in many stress studies and is thought to be essential in maintaining efficient translation under stress. In addition, the down-regulation of elongation factors (e.g., elongation factor 1-alpha, beta and Ts, and elongation factor 2) that participate in protein synthesis in the cell cycle and facilitate translational elongation in the ribosome suggests that protein synthesis was hindered under high UV-B radiation. Furthermore, UV-B radiation might change protein phosphorylation levels and accelerate the degradation of misfolded proteins as suggested by the observation that LRR receptor-like serine/threonine-protein kinase At5g45780 was down-regulated and several proteases (e.g., thiol protease aleurain and protease Do-like 8) were up-regulated. The former catalyzes protein phosphorylation at serine/threonine (Gao et al., 2010), and the latter plays a role in proteolysis and leads to the mobilization of nitrogen under stressed conditions. Therefore, the results indicated that UV-B radiation might have a complex effect on gene expression regulation in P. cathayana, and a large variety of changes in gene expression processes, particularly posttranscriptional regulation and post-translational modifications, occurs in poplars. However, further study is needed to test this hypothesis.
Under high UV-B radiation conditions, many biological processes, including ion transport, ROS scavenging, and UVabsorbing compound synthesis, require extra energy. ATP is produced via light-driven photosynthetic reactions in the thylakoid membrane, and it is primarily utilized in carbon fixation reactions (Yin et al., 2010). In this study, ATP synthase subunits were down-regulated in females but up-regulated in males, indicating that more ATP might be produced in males than in females under UV-B radiation. However, several ATP carrier-and exchange-related proteins were decreased in abundance in both sexes, e.g., 3-isopropylmalate dehydrogenase, ATPase and ADP/ATP carrier proteins, suggesting that energy metabolism was hindered by high UV-B radiation in both sexes. However, this result was contrary to other reports that these proteins were up-regulated in plants suffering from short-term UV-B radiation, such as Lonicera japonica and Arabidopsis (Yin et al., 2010;Zhang et al., 2013). One reason might be that these changes in the proteins involved in energy metabolism varied with plant species, stress duration and development. In another study, ATP was produced by light-driven photosynthetic reactions and carbohydrate metabolism. However, in this study, several key enzymes involved in these pathways were down-regulated by UV-B radiation, e.g., PSI, PSII, cytochrome b6, glyceraldehyde-3phosphate dehydrogenase and glucan endo-1,3-beta-glucosidase. Additionally, the high UV-B treatment lasted for a long time (one growth season), and excessive energy production will lead to accumulation of ROS and plant cell damage. Therefore, the decrease in ATP production may be an adaption to high UV-B conditions in nature, especially under higher intensities.
The photosynthetic system and subcellular organelles are vulnerable to UV-B damage. Previous studies have shown that high UV-B radiation lowers chlorophyll pigment content, damages the photosystem, and distorts chloroplast structure in P. cathayana, thereby leading to a decline in P n (Ren et al., 2006(Ren et al., , 2010Xu et al., 2010;Feng et al., 2014). The decrease in chlorophyll content also has a very distinct effect on chlorophyll-binding proteins and protein complexes, thus further hindering light capture and delivery (Varsano et al., 2006;. In this study, the abundance of many photosynthetic proteins, especially light-dependent reaction proteins, was changed in both sexes (Figure 5). For example, many chlorophyll a/b binding proteins (CABs) and PSI and PSII reaction center proteins were lower in abundance under solar UV-B radiation, whereas some thylakoid lumen proteins and PsbP domain-containing proteins were increased. CABs are photosynthetic proteins that are typically downregulated in plants in response to enhanced UV irradiation (Taylor et al., 1996;Zhang et al., 2013). These proteins are responsible for binding to chlorophyll pigment molecules in light-harvesting complexes and transferring excitation energy during photosynthesis (Bassi et al., 1990;Jespersen and Huang, 2015). In plants, the degradation of the photosystem is linked to the remodeling of light-harvesting antenna, e.g., chlorophyll a/b binding proteins (Moseley et al., 2002;. The lower expression of photosystem reaction center proteins implies that UV-B radiation causes electron transport discordance and redox homeostasis imbalance between PSI and PSII (Joshi et al., 2011). The cytochrome b6f complex, which mediates electron transport between PSI and PSII, is also a thylakoid-bound protein that is an essential component of the electron transport chain and serves as the final electron acceptor in the cyclical electron flow pathway of photosynthetic light reactions (Munekaga et al., 2004;Cramer et al., 2006;. The PsbP and thylakoid lumen proteins were up-regulated in both sexes, indicating that solar UV-B radiation may result in instability of the photosystem (Jiang et al., 2015). PsbP proteins are essential for the regulation and stabilization of PSII (Ifuku et al., 2005). Recent reports have indicated that thylakoid lumen proteins play roles in regulating thylakoid biogenesis and the activity and turnover of photosynthetic protein complexes, especially the PSII and NAD(P)H dehydrogenase-like complexes. However, the function of the majority of luminal proteins in Populus remains unknown (Järvi et al., 2013;. Additionally, the abundance of OEE and Calvin cycle proteins was increased in males only, indicating sexual differences in the response of the photosynthetic system (Figure 5). OEE and Calvin cycle proteins participate in the light-dependent and light-independent reactions of photosynthesis. The upregulation of these proteins might increase the photosynthetic capacity of males, which could partially explain the higher photosynthetic rate in males than in females under high solar UV-B radiation.

Sex-Specific Changes in the Proteome under UV-B Radiation
Although a common set of proteins was changed in both sexes, there were many of sex-specific DEPs (Figure 5). For example, the 26S proteasome non-ATPase regulatory subunit 1 was up-regulated in females, whereas other proteasomes and their subunits were up-regulated in males (e.g., cysteine proteinase 2 and 3, proteasome subunit beta type-2-A and -6, and subtilisin-like proteases). Calnexin homolog 1 and luminalbinding protein 5, which control the apparatus of incorrectly folded proteins in the endoplasmic reticulum, were downregulated in males only. The sex-specific changes in these proteins may be attributed to two reasons: (1) sex-specific gene expression and regulation may lead to sex-specific protein abundance and (2) limited proteomic detection technology (e.g., the number of identified peptides per proteins is still low) increases the difficulty of identifying all of the protein profiles in the two sexes, especially for the alkaline proteins and membrane proteins at low abundance.
A considerable number of proteomic profiles indicated that more proteins involved in amino acid metabolism were upregulated to a greater degree in males than in females under high UV-B radiation (Figure 5). This result is consistent with that of the transcript profiles (Jiang et al., 2015). Alanine aminotransferase 2 and glycine dehydrogenase [decarboxylating] participate in the degradation of L-alanine and glycine, and isoaspartyl peptidase/L-asparaginase 1 are involved in aspartate production. The up-regulation of these proteins indicates that alanine, glycine, and aspartate catabolic processes changed more in males upon UV-B radiation. The changes in hydrolyzed amino acid compound observed in this study also confirmed this phenomenon. Additionally, several lyases, e.g., glycerate dehydrogenase and LL-diaminopimelate aminotransferase, were down-regulated in females, and 3-dehydroquinate synthase, diaminopimelate decarboxylase 2 and N-carbamoyl-L-amino acid hydrolase were down-regulated in males. These lyases catalyze the different chemical reactions that utilize and produce special amino acids. Glycerate dehydrogenase participates in glycine, serine, and threonine metabolism and glyoxylate and dicarboxylate metabolism, and their down-regulation only in females indicates that these metabolic processes are regulated more stringently in females than in males. Our results provide further evidence that the female strategy of gene regulation differs from that of males at the proteomic level.
For the females, except for 2 proteins that changed in common in both sexes, most of the stress responsive proteins were down-regulated. For example, 12-oxophytodienoate reductase 3 catalyzes the biosynthesis of jasmonic acid (JA), which is important as a gene regulator for development and defense (Schaller et al., 2000;Tani et al., 2008). 1-aminocyclopropane-1-carboxylate oxidase 3 (ACO3) is an essential gene during plant senescence and development and may play a role in ethylene synthesis (Hunter et al., 1999). The auto-inhibited plasma membrane P-type H+ ATPases were activated by 14-3-3 proteins (Jahn et al., 2002). Additionally, actin-7, phosphate carrier protein and protein canopy-1 were down-regulated in females, and several proteins were also down-regulated in males including allene oxide cyclase 3 (AOC3), probable linoleate 9Slipoxygenase 5 (LOX1.5) and two protein aspartic proteases in guard cell 1 (ASPG1). AOC3 is involved in the production of 12-oxo-phytodienoic acid (OPDA), a pre-cursor of JA. ASPG1 play essential roles in restricting bacterial growth, plant defense and drought avoidance (Willmann et al., 2011;Yao et al., 2012). Plant lipoxygenases may be involved in a number of diverse aspects of plant physiology, including growth and development, pest resistance, senescence, and the wound response (Kolomiets et al., 2000). The down-regulation of these proteins indicates that higher intensity UV-B radiation may induce a sex-specific  decrease in the defensive ability of poplars. Interestingly, several stress-related proteins were up-regulated only in males, such as heavy metal-associated isoprenylated plant protein 26 (HIPP26), heme-binding-like protein At3g10130, pathogenesisrelated protein 1A (PRP1A), and stable protein 1 (SP1). PRP1A are plant disease-resistance proteins and play roles in plant defense. HIPP26 and At3g10130 are metal-binding proteins and play roles in metal detoxification (Barth et al., 2009;Vanhee et al., 2011). High UV-B radiation resulted in an increase in these stress-response proteins in P. cathayana males but not in females. Thus, these proteins might play roles in enhancing UV-B resistance to a greater extent in P. cathayana males than females. Our previous parallel studies have shown that at the physiological level, solar UV-B radiation significantly decreases the photosynthetic capacity of P. cathayana, with a greater decrease in females than males (Xu et al., 2010;Feng et al., 2014). At the transcriptome level, solar UV-B radiation induced changes in the expression of a set of genes, and sexual differences were observed. For example, many differentially expressed genes involved in amino acid metabolism are up-regulated in P. cathayana males but down-regulated in females (Jiang et al., 2015). However, a full understanding of poplar regulatory mechanisms under UV-B solar radiation remains to be attained. Sex-related adaption to solar UV-B radiation in P. cathayana is such a sophisticated process that we cannot identify all of the regulating pathways that play dominant roles. For instance, flavonoids were the most frequently reported secondary metabolites responding to UV-B radiation in plants (Bassman, 2004;Schreiner et al., 2012). Our previous transcript study detected two flavonoid genes in P. cathayana males (Jiang et al., 2015). However, we did not identify any protein related to flavonoid metabolism at the proteomic level in this study (Figure 5). Numerous studies have indicated that the correlation between protein abundance and mRNA transcript levels is limited in woody plants (Lippert et al., 2009;Zhang et al., 2010;Dong et al., 2016). It is widely assumed that post-transcriptional regulation and post-translational modification play important roles in translational efficiency, which leads to discordance between mRNA expression and protein abundance. Therefore, to further understand the evolutionary and adaption mechanisms of poplars to UV-B radiation, differences in the translation, posttranscriptional modification, and metabolism responses should be investigated.

CONCLUSIONS
Our results showed that P. cathayana females and males under high UV-B radiation exhibited common and sexspecific responses in terms of physiology and proteome dynamics. These common responsive proteins in both sexes were mainly categorized into carbohydrate and energy metabolism, translation/transcription/post-transcriptional modification, photosynthesis and redox reactions. The sexually different responsive proteins were involved in translation/transcription/post-transcriptional modification, amino acid metabolism, and stress responses. P. cathayana males showed a greater scope of protein changes than females under high UV-B radiation, suggesting that males are more sophisticated in their response to UV-B radiation compared with females. Although our study is limited by the low number of replicates (three replicates each treatment), and the results may vary depending on plant growth conditions (e.g., temperature, condensation and available light spectrum), it reveals some dynamic and sex-specific changes in poplars at the proteomic level in response to high UV-B radiation and provides new insights into the mechanism of the malebiased sex ratio at high altitude, thus complementing existing knowledge.

AUTHOR CONTRIBUTIONS
YunxiangZ was responsible for the analysis of proteomic data and parts of manuscript writing. LF did the field work, and she was also responsible for obtaining proteomic data and analysis of physiological parameters. HJ and YuanbinZ did much of the field work and the measurements. SZ had the initial research idea, was responsible for parts of manuscript writing and acquired the funding for the project which was done in his laboratory.