{'fraction_non_alphanumeric': 0.037074042911232646, 'fraction_numerical': 0.022822885990744637, 'mean_word_length': 4.876699629171817, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Gender discrimination of Populus tomentosa barks by HPLC fingerprint combined with multivariate statistics

Abstract A high‐performance liquid chromatography (HPLC) fingerprint method with multivariate statistical analyses was applied to discriminate the male and female barks of Populus tomentosa for the first time. The samples of 11 male and 13 female barks of mature P. tomentosa were collected in Beijing. The chemical fingerprint of methanol extract was established by HPLC method with diode array detector (DAD). The principal component analysis (PCA), hierarchical clustering analysis (HCA), and supervised orthogonal partial least squares discriminant analysis (OPLS‐DA) were applied to discriminate male and female barks based on the area of common peaks identified in HPLC fingerprints. A clear grouping trend (R 2 X, 0.83; Q 2, 0.595) among the male and female samples was exhibited by PCA score plot. Two groups were clearly divided into male and female samples by HCA. Both male and female samples were well discriminated with OPLS‐DA (R 2 X, 0.775; Q 2, 0.795). Seven potential chemical markers were screened by variable importance in projection (VIP values >1.0) of OPLS‐DA model and four of them were identified as micranthoside, siebolside B, sakuranin, and isosakuranin. The HPLC fingerprint combined with multivariate statistical analyses could be used to discriminate the gender of barks of P. tomentosa and revealed the differences in chemical components, which enriched the basic studies on dioecious plant.

. It is a great value to discover the differences in physiology of dioecious plants, just as the differences of androgen and estrogen contents in human body (Cao et al., 2018).
The physicochemical properties of male and female individuals of dioecious plants may different. For example, the contents of polysaccharides and tichosanthin in male plant were higher than those in female plant, which provide reasonable explanation for that male plant was better . So, if there were differences between male and female barks, the pharmacological activities may different.
The plant of Populus tomentosa Carrière (Fam. Salicaceae) is an endemic afforestation tree in north China. Some researches in discrimination of male and female plant based on buds, leaves, and flowers were reported. The phenomena of higher content in volatile compounds of methyl benzoate, benzyl benzoate, and 2-cyclohexen-1-one in female flower buds and higher content of ethyl benzoate in male flower buds can discriminate the mature male and female buds of P. tomentosa (Xu, et al., 2019a). The polymorphism between male and female trees was tested by using the methods of RAPD and bulked segregate analysis and searched the primer S60 could be used for selecting the gender of P. tomentosa (Hou et al., 2009). The isoperoxidase in female plants could be separated into 10-12 bands but male plants only separated in to 4-5 bands by using isoelectric focusing in polyacrylamide gel electrophoresis. In addition, the protein content solvable, peroxidase activity, and the scanning curves of electrophoresis pattern of isoperoxidase between male and female plant of P. tomentosa were obviously different (Tian et al., 1993). Transcriptome analysis identified genes significantly differentially expressed between the male and female plant of P. tomentosa, including genes related to floral development, phytohormone synthesis and metabolism, and DNA methylation (Song et al., 2013). However, whether there are differences in chemical composition of the barks which are not the sexual organ of P. tomentosa, and whether the differences can be used to identify the male and female individual of dioecious plants have not been reported so far. The barks of P. tomentosa have the activities of clearing heat, removing dampness, relieving cough, and resolving phlegm (Editorial Board of Chinese Medica, 1999) for the treatment of hepatitis, dysentery, stranguria with turbid urine, cough, and phlegmatic asthma and prevention of infectious hepatitis (Jiangsu New Medical College, 1985;Lv & Tao, 1985). So, it would be valuable and instructive for the clinical use to differentiate male and female barks. In addition, the research on non-sexual organs would enrich the range of dioecious plant.
Previous researchers mostly studied the gender identification on the differences in the mature individuals in morphology , physiological and biochemical (Song et al., 2014), peroxidase isozyme (Tian et al., 1993), transcriptome, phytohormone, and DNA methylation analysis (Song et al., 2013). It is a complex system to discriminate the gender of dioecious plants. The gender identification based on gene expression products is inaccurate and unreliable, and cannot be primarily applied to the early gender identification of dioecious plants (Yin et al., 2003). The RAPD method has poor stability and reproducibility because the results are highly dependent on certain factors such as the concentrations of the template DNA and Mg ions (Du et al., 2019).
High-performance liquid chromatography (HPLC) is a widely used technique for species differentiation and quality control of herbal medicines. Fingerprinting strategies provide analytical signals related to the composition in a non-selective way such as by collecting a spectrum or a chromatogram (Cuadros-Rodríguez et al., 2016).
HPLC fingerprint is a type of quality assessment system with integrated, macroscopic, and "fuzzy" nonlinear characteristics that comprehensively analyzes the complex components of food, plant, and herbal medicines. It emphasizes systematic characterization of the analytical target that contains large amount of information and relates to the classification of samples based on integral chemical background (Wipawee et al., 2008). Identifying the remarkable differences in fingerprints only depend on content determination is hard through it could be an effective strategy to help assess the quality of herbs. It was real a challenge to find remarkable compound in the past, but the wide application of multivariate statistical analysis makes it easy nowadays (Zhang et al., 2015).
The multivariable statistical approach has become more important for analytical chemistry. Such approaches are used to predict a pure component feature, and its concentration profiles from a set of spectra, chromatographs, and any situation where multiple measurements are acquired (Beebe & Kowalski, 2012). Recently, fingerprint methods combined with multivariate analysis have been used to successfully classify and discriminate different complex herbal sources Qin et al., 2015;Tang et al., 2015).
In this study, a total of 24 of male and female barks of P. tomentosa was collected. A HPLC fingerprint method was established and multivariate statistical analyses including principal component analysis (PCA), hierarchical clustering analysis (HCA), and orthogonal partial least squares discriminant analysis (OPLS-DA) were applied to discriminate the male and female barks of P. tomentosa. Then, the potential chemical markers were screened for the discrimination of gender of P. tomentosa.

| MATERIAL S AND ME THODS
The samples of 11 male barks and 13 female barks of mature P. tomentosa were collected on February 14, 2020, near the Sanyuan Bridge which was located at 39°57'37.31" north latitude, 116°27'28.53" east Chromatography grade methanol was purchased from Fisher Scientific (Fair Lawn). Ultrapure water was purchased from Hangzhou Wahaha Group Co., Ltd. Reference standards of micranthoside, siebolside B, sakuranin, and isosakuranin were isolated and purified from male barks of P. tomentosa by conventional column chromatography, their structures were identified by ESI-MS, 1 H-NMR, and 13 C-NMR in comparison with the data from literatures, and related data have been published (Liu et al., 2012b;Xu, et al., 2019b). Their purities were determined to be higher than 95% based on HPLC analysis with a method of peak area normalization.
All of barks were dried by airing and crushed into a powder with a grinder for 1 min and passed through a 24-mesh sieve. The extraction method is modified from a previously reported method (Wang et al., 2005). Powder samples of 3.0 g were weighed accurately into 150-mL flat-bottomed flasks and 50-mL aliquots of methanol were added. Samples were then weighed and refluxed for 120 min, and were weighed again after cooling. Weight losses from reflux of extracts were replaced with methanol. The filtrates were passed through 0.45μm membrane filter as sample solutions.
All sample solutions were prepared as described in preparation of sample solution. The precision of the chromatographic method was established by analyzing the same sample solution six times within 1 day. The stability was determined by analyzing a single-sample solution that was stored at room temperature at 25℃ for 0, 3, 6, 9, 12, 24, and 48 hr. The repeatability was determined by analyzing six-sample solutions prepared from the same sample. The similarity between each of the sample fingerprints and a reference fingerprint were calculated to estimate precision, stability, and repeatability.
The raw HPLC chromatographic data of the 24 samples were converted into AIA format files for further processing. All of these files were imported into the Similarity Evaluation System for Traditional Chinese Medicine Chromatographic Fingerprinting (Version 2004 A; Committee for the Pharmacopoeia of PR China). One sample was randomly selected to generate the template. Other samples were performed multipoint correction and automatic match to generate a fingerprint that represented the characteristic mode. The similarities of every HPLC fingerprint were calculated based on the reference fingerprint. For further analysis, the peak areas of all common peaks were calculated simultaneously and exported as an excel file for further statistical analysis including PCA, HCA, and OPLS-DA.
Multivariate analyses of PCA, HCA, and OPLS-DA were performed using the SIMCA-P version 14.1 software (Umetrics) to evaluate relationships on the basis of similarities or differences between groups of fingerprint data. Unit-variance scaling was applied in data preprocessing for both PCA and OPLS-DA.
PCA, an unsupervised pattern recognition pattern, reduced the dimension of the data matrix and converted the original variables into new independent variable called principal components (PCs) (Albu et al., 2015). It revealed the interrelationships between different variables and interpret sample patterns, groupings, similarities, and differences (Xu et al., 2019). R 2 X and Q 2 values were the important parameters of the modeling in PCA. R 2 was a measure of how well the model fitted the data. R 2 X showed how much of the variation in the X variable could be explained by the selected components. Q 2 indicated how well the model predicts new data.
PCA result was displayed in the form of score plots. The HCA was performed with SIMCA-P version 14.1 software (Umetrics, Malmö, Sweden) based on the PCs scores of samples, and distances between samples were calculated using Ward's minimum-variance method, with results presented as a dendrogram.
The OPLS-DA extended a regression of the PCA, used the class membership to maximize the variation, introduced an orthogonal signal correction (OSC) filter to separately handle the systematic variation correlated with or uncorrelated with the Y variable, and therefore had better discriminant ability for the samples with larger within-class divergence than PCA (Bylesjoe et al., 2008). The quality of the model was estimated by parameters of Q 2 , R 2 X, and R 2 Y which described how well the model was fitted. The values of R 2 and Q 2 ranged from 0 to 1, where 1 indicated perfect fitness and predictivity (Huang et al., 2018). In order to identify the components that contribute greatly to the sample classification, the variable importance in projection (VIP) scores were used to screen differential markers. VIP scores were a weighted sum of squares of PLS weights, with scores larger than 1 indicating variables were important in the sample classification (Zheng et al., 2018).
The student's t-test was performed on the SPSS 20.0 software for windows next to confirm the significance of differences in differential markers screened from OPLS-DA between male and female barks. Sample M11 was chosen to identify the chromatographic peaks and analyzed using an Acquity UPLC H-Class from Waters (Milford, MA, America), equipped with a binary solvent delivery pump, an autosampler, a column compartment, and a PDA detector. A reverse-phase Acquity UPLC BEH C 18 column (2.1 × 50 mm, 1.7 μm) was used. Flow rate was 0.35 ml/min; injection volume, 1 μL; and column temperature, 35°C. Mobile phases were methanol (A) and water (B) and the elution conditions applied were as follows: 0-3 min, linear gradient 10%-30% A; 3-8 min, linear gradient 30%-35% A; 8-10 min, linear gradient 35%-55% A; and 10-12 min, linear gradient 55%-100% A. UV-visible spectra were recorded at 230 nm.
All MS data acquisitions were performed on a Vion IMS quadrupole time-of-flight (QToF) configuration (Waters) equipped with an electrospray ionization (ESI) source. The mass spectrometer detector was operated in positive mode, with the capillary voltage at 3.0 kV, 30 V cone voltage, 450 ℃ desolvation temperature, and 120 ℃ source temperature. Nitrogen was used as the desolvation and cone gas at flow rates of 1,000 L/hr and 50 L/hr, respectively. Data acquisition was recorded in the mass range m/z 50-1000 with a scan time of 0.2 s. Waters MassLynx V4.1 software was used for data acquisition and processing.

| RE SULTS
In order to evaluate the feasibility of the HPLC method, the similarities between each of the samples fingerprint and reference fingerprint were calculated to estimate precision, stability, and repeatability. The results of the similarities of precision, stability, and repeatability were found to be higher than 0.998, higher than 0.975, and above 0.998, respectively, indicating that the HPLC measurements were stable and under control.  (Tables S1 and S2), which demonstrated that the male and female samples were highly similar.
To evaluate the differences between male and female barks of P. tomentosa, PCA was carried out with the areas of 13 common peaks (Tables S3 and S4). The scores plot obtained by PCA was shown in Figure 3. The first contribution value of two principal components was accounted for 78.3% (R 2 X) of the total variance (PC1 described 62.4% and PC2 described 15.9% of the sample variability). The predictive ability of the model (Q 2 ) was 0.595, which demonstrated that it was a good model (Q 2 ≥ 0.50) (Triba et al., 2015;Zhang et al., 2018;Zheng et al., 2018). There was a clear separation between male and female barks of P. tomentosa. These results indicated that there were differences in some chemical compounds between male and female barks of P. tomentosa. HCA was performed based on two PCs from the above PCA model, and displayed relationships between male and female samples in the form of a dendrogram. As shown in Figure   HPLC fingerprint were higher than 1.0. It was indicated that these seven peaks possessed the most influence on the discrimination between male and female samples. The student's t-test was performed based on the peak areas of seven compounds screened in OPLS-DA model and the result showed that the differences of these seven peaks between male and female samples were significant (p <.05).

F I G U R E 5
Orthogonal partial least-squares discrimination analysis (OPLS-DA) models for the classification of male and female barks of P. tomentosa. (a) Score plot of the OPLS-DA results. Triangles represent female samples, and diamonds represent male samples. X axis represents the scores of PCs during OSC; Y represents the scores of orthogonal components during OSC; the differences between groups and within groups can be seen from X and Y axes, respectively. (b) Plot of the coefficients for the OPLS-DA. The X axis represents 1-13 common peaks; the Y axis represents VIP value The analysis of MS data was showed in Table 1. The MS and MS/ MS spectrograms contained the quasi-molecular ion peak and major fragments, respectively, were displayed in Figures S2-S8. The peaks of No.6,No.7,No.8,and No.9 in HPLC fingerprint were identified by comparing the retention time of chromatographic peaks ( Figure S1) and the corresponding MS and MS/MS spectrograms of each chromatographic peak with reference standards (Figures S2-S5) obtained from male barks of P. tomentosa. The results showed that the peaks of No.6,No.7,No.8,and No.9 were micranthoside, siebolside B, sakuranin, and isosakuranin ( Figures S2-S5). The peak No.10 was identified as isograndidentatin A by comparing the MS and MS/MS data ( Figure S6) with literatures (Xu, et al., 2019b;Yang et al., 2013).
The peak No.11 was the isomer of sakuranin (No.8) and isosakuranin (No.9) for their same molecular weight and major MS fragments ( Figure S7). The molecular formula of peak No.3 was calculated as C 21 H 28 O 10 ( Figure S8).

| D ISCUSS I ON
In the presented study, an accurate and feasible analytical method Among the seven compounds screened based on HPLC fingerprint combined with multivariate statistics, four main compounds were identified. The pharmacological activities of sakuranin are allergypreventive (Yuko et al., 2006), antihyperlipidemic (Weber et al., 2008), anti-inflammatory activities (Zhang et al., 2006), and it can be used as an acetylcholinesterase inhibitor (Remya et al., 2014). Siebolside B and isograndidentatin A have allergy-preventive (Yuko et al., 2006) and antioxidant activity (Si et al., 2009). It is demonstrated that these compounds may be the bioactive components in barks of P. tomentosa.
The peak areas of these four compounds in male barks were higher than those in female barks, which might affect the effect of barks of P.
tomentosa of different sex. This had important guiding significance for the use and further study of barks of P. tomentosa.
The flowers as the sexual organ for their different chemical components could be used to discriminate the gender. Similarly, the barks, the asexual organ, also can be used to discriminate the gender based on our research, which provide a new idea for researching the dioecious plants.
Some dioecious plants are medicinal in traditional Chinese medicine, such as the barks of mature Eucommia ulmoides Oliv.
(Eucommiaceae) and Ilex rotunda Thunb. (Aquifoliaceae). However, there are no clear regulation about using male barks or female barks.
The present study demonstrates the differences between male and female barks, which implies the differences in clinical efficacy. The study of the gender discrimination of dioecious plants has referential significance for their medicinal value.

ACK N OWLED G M ENT
We thank the support by China Agriculture Research System (No: CARS-21).