Abstract

In this study, a fast, simple, precise, and sensitive hydrophilic interaction liquid chromatography (HILIC) method was established for simultaneous determination of free amino acids in three different varieties of duckweed including Spirodela polyrhiza (L.) Schleid., Landoltia punctata (G. Mey.) Les & D. J. Crawford, and Lemna aequinoctialis Welwitsch by ultrahigh performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS/MS). Method validation was processed in terms of linearity, precision, stability, repeatability, and accuracy as well as limits of detection and quantification. The developed method was applied for quantification of 59 batches of samples. Then chemometric analysis was used to evaluate different duckweeds by principle component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA). The results demonstrated that there was no significant difference in FAAs’ profile among three varieties of duckweed.

1. Introduction

Amino acids are kinds of nitrogenous components that play vital and diverse roles in metabolism and have attracted significant attention in food, feedstuff, and alimentary supplements. Amino acids are basic units and important ingredients of proteins and involved in the progress of biosynthesis for glycoprotein, porphyrins, neurotransmitters, polyamines, and nitric oxide [15]. Nutritional studies show that amino acids can modulate gene expression and enhance the growth of skeletal muscle and small intestine [6].

The family of Lemnaceae colloquially known as duckweed has been consumed as human food since long [7]. Duckweeds have the striking capacity of explosive reproduction [8]. They can grow almost everywhere with appropriate temperature and nutrition in water. But their greatest potential is to produce large quantities of protein-rich biomass that is a promising food resource for humans [9] and is suitable for feeding a wide range of animals including fish, poultry, and cattle [1012]. The contents of free amino acids (FAAs) are very important for the evaluation of the protein-rich food. Fresh duckweed contains a large amount of water with moisture content of 86%–94%. But for dried samples, their protein concentration ranges from 25% to 40% with much lower fibrous material than that of land forage [13]. The habitat influences the content of proteins and fibers dramatically. Content of amino acids and protein always depend on their species pattern and growth environment, such as sunshine, nutrition media, and their growth pattern: gregarious or solitary. Under high nutritional conditions, more proteins and fewer fibers are accumulated in duckweed and vice versa. The amino acid type of duckweed is very close to that of animals and thus could be used by animals efficiently [7, 14]. The first limited amino acid (lysine) is similar to soybean and is more superior to that in other food such as sorghum and maize.

There are several methods reported in literature for analysis of free amino acids. In general, the FAAs are extracted with solvents water, formic acid, hydrochloric acid, or ethanol, followed by filtration or centrifugation before analysis [1517]. Several techniques have been described for detection of free amino acids including ion-exchange chromatography [18], precolumn derivatization followed by reversed high-performance liquid chromatography coupled with diode array or fluorescence detector [19, 20], gas chromatography, mass spectrometry, and capillary electrophoresis [21]. In recent years, determination of FAAs using ultrahigh performance liquid chromatography hyphenated to mass spectrometry has gained great notice because of the merit of selective separation for the polar underivatized analytes under HILIC condition with acetonitrile-dominated mobile phase [2225].

In present study, a feasible and reliable UPLC-QTRAP-MS/MS method was developed and validated for quantitative analysis of 24 underivatized FAAs simultaneously. Multivariate statistical analysis was employed to assess the differences in the profiles of FAAs among three species of duckweed. This simple and fast analytical method would be used for quality control of duckweeds.

2. Materials and Methods

2.1. Reagents and Materials

59 batches of duckweed were collected from different regions in China as shown in Table 1. The samples were identified and classified by PCR [2628], and the number of populations was as follows: 25 populations from species Spirodela polyrhiza, 16 populations from species Landoltia punctata (its another previously invalid name is Spirodela oligorrhiza (Kurz) Hegelm.), and 18 populations from Lemna aequinoctialis. Fresh materials were washed by tap water, dried in sunshine, and pulverized, and the powder was screened through 60-mesh sieve. Reference standards of tyrosine (Tyr), alanine (Ala), glutamic acid (Glu), phenylalanine (Phe), histidine (His), isoleucine (Ile), aspartic acid (Asp), lysine (Lys), glycine (Gly), valine (Val), citrulline (Cit), cystine (Cys2), proline (Pro), tryptophan (Trp), leucine (Leu), arginine (Arg), cysteine (Cys), hydroxylproline (Hpro), methionine (Met), threonine (Thr), asparagine (Asn), and serine (Ser) were purchased from National Institutes for Food and Drug Control (Beijing, China) with the purity more than 98%. Reference standards of γ-aminobutyric (GABA) and glutamine (Gln) with the purity more than 98% were obtained from Aladdin Industrial Corporation (Shanghai, China). Acetonitrile (HPLC grade) was purchased from Merk (Darmstadt, Germany). Formic acid and ammonium formate (HPLC grade) were provided by Mreda Technology Incorporation (USA). Deionized water used for analysis procedure was produced by a Milli-Q Academic ultrapure water system (Millipore, Bedford, MA, USA).

2.2. Instrumentation and Chromatographic and Mass Conditions

Chromatographic analysis was performed on a Shimadzu LC-30 ultrahigh performance liquid chromatography system (Shimadzu, Japan), which consisted of a communication bus module (CMB-20A), a vacuum degasser, binary gradient pumps (LC-30AD), an autosampler (SIL-30A), and a column oven (CTO-30A) coupled with an AB Sciex QTRAP 5500 (AB SCIEX, USA). A Waters XBridge Amide column (2.1 mm × 150 mm, 2.5 μm, Waters, USA) was used for chromatographic separation. The binary mobile phase was composed of acetonitrile (A) and 10 mM ammonium formate with 0.2% (v/v) formic acid (B) at a flow rate of 0.4 mL/min. The linear gradient elution was carried out as follows: 0–5 min, 90%–85.5% A; 5–11 min, 85.5%–54% A; 11-12 min, 54%–10% A; 12–14 min, 10% A; 14-15 min, 10%–90% A; 15–20 min, and 90% A. The reequilibration time was 4 min with a total running time of 20 min. The column compartment was kept at 20°C, while the autosampler trial was maintained at 15°C, and injection volume was 2 μl. The needle was washed with mixtures of acetonitrile and water.

The mass spectrometry assay was performed on a triple quadruple mass spectrometer equipped with an electrospray ionization source. The spectra were recorded under positive-ion type with the multiple reaction monitoring (MRM) mode. The capillary voltage was 5500 V, and desolvation gas temperature was 550°C. The curtain gas was 35 psi with both nebulizer and drying gas of 55 psi. Nitrogen was used as source gas with purity more than 95% and collision gas with purity over 99.999%. The 24 components of amino acids were optimized in the tuning mode with mass only to obtain the mass of the precursor and product ion, the best declustering potential (DP), and collision energy (CE), respectively. The details are shown in Table 2. The entrance potential (EP) and collision cell exit potential (CXP) were set as default value. The typical LC-QTRAP-MS/MS chromatograms are shown in Figure 1.

2.3. Preparation of Reference Compound Solution

Stock solutions of individual standards were made by dissolving accurately weighed components in 0.5% formic acid at a concentration of 0.2 mg/mL approximately and stored at 4°C before use. The mixed working solutions of all the standards were diluted with 0.5% formic acid solution to a series of appropriate concentrations immediately before analysis and used to attain the calibration curves.

2.4. Preparation of Sample Solution

Each sample powder was weighed 0.2 g accurately and transferred into a 25 mL stainless tube with 6 stainless beads with a diameter of 0.6 mm. Then, the sample was extracted in a high-throughput tissue-grinding apparatus with vibrating frequency of 70 kHz for 120 s with 15 mL 0.5% formic acid. The sample solution was centrifuged at 15000 rpm for 10 min, and then the supernatant was screened through a 0.22 μm polytetrafluoroethylene membrane filter and stored in a glass bottle at 4°C for later LC-MS/MS analysis.

2.5. Method Validation

For each reference substance, the calibration curve was confirmed by linear regression of the peak area versus concentration. The limits of detection (LOD) and limits of quantification (LOQ) were determined by diluting the standard solution till the signal-to-noise ratios were about 3 and 10, respectively. Precision was determined by intraday and interday variability. The intraday variability was conducted by determining the same standard solution in six replicates on the same day. The interday variability was conducted by determining the same solution for three consecutive days. The relative standard deviation (RSD) values were calculated to denote the precision. Stability of the sample solution was analysed by peak areas of analytes at 0, 2, 4, 8, 12, and 24 h at room temperature. To evaluate the repeatability of the developed method, the same sample was analysed in six replicates and variations were expressed as RSD. The recovery test was performed to assess the accuracy of the method. For the test, a known amount of the 24 standard components was spiked into a certain amount (0.1 g) of sample. Then, the spiked sample was extracted and analysed as described. The recovery rate was calculated by using the following formula: recovery (%) = (measured amount − original amount)/spiked amount  × 100%. Six replicates were performed for the recovery test.

2.6. Data Analysis

The quantitative analyses of acquired data were processed by MultiQuant 3.0.2 (AB SCIEX, USA). Multivariate statistical analysis including PCA and OPLS-DA was used to classify the sample. SIMCA-P 14.1 (Umetrics AB, Sweden) was used to perform statistical analysis for the data of 59 batches of duckweed.

2.7. Analysis of Minerals and Total Proteins

Twenty-four minerals of each environment water where duckweed lived were quantified by ICP-MS. Total nitrogen and phosphorus of the water body were measured according to the alkaline potassium persulfate digestion method GB11894-89 and the ammonium molybdate spectrophotometric method GB11893-89 published by the State Bureau of Technological Supervision of China. Total proteins of each duckweed were determined by the Kjeldahl method. Each sample of duckweeds was hydrolyzed in acidic conditions, and then total amino acids were quantified by the HPLC-MSMS method mentioned above.

3. Results and Discussion

3.1. Optimization of Extraction Procedure

Free amino acids can dissolve in methanol, ethanol, water, diluted acid, and so on. In order to acquire an efficient extraction strategy of amino acids from duckweed, variables of the extraction process were investigated. The solvents (methanol, ethanol, water, 0.1% formic acid, and 0.5% formic acid), methods (heating reflux, ultrasonication, and tissue-grinding), and time (20, 30, and 40 min for ultrasonication; 60, 90, and 120 s for tissue-grinding) were studied. The results showed that the amount of amino acids had no significant difference between the ultrasonication time of 30 min in 0.5% formic acid and tissue-grinding time of 120 s in 0.5% formic acid. The amounts of extraction by other methods were not more than that by these two modalities. From the viewpoint of working efficiency, the tissue-grinding method was much more time-saving than the ultrasonication method by the extracting time of 120 s versus 30 min.

3.2. Optimization of UPLC Condition

Chromatographic conditions were optimized to attain a satisfactory separation for amino acid especially for isomeric molecules. Columns (Waters XBridge Amide, 2.5 μm, 2.1 × 150 mm and Agilent Poroshell 120 HILIC, 2.7 μm, 2.1 × 150 mm), column temperature (20, 25, 30, 35, and 40°C), mobile phase (acetonitrile-0.1% formic acid, acetonitrile-0.2% formic acid, acetonitrile with 0.2% formic acid-0.2% formic acid, and acetonitrile-10 mM ammonium formate with 0.2% formic acid), and different kinds of gradient elution were investigated.

The results demonstrated that the peak shape of amino acids became good when the concentration of formic acid rose, but retention time would be shortened significantly especially if some formic acid was complemented into the acetonitrile. The buffers such as ammonium formate increased the resolution degrees of amino acids. Better retention times and separations were observed by conducting analysis on Waters Amide column than on Agilent Poroshell 120 HILIC column, and the best chromatograms were gained with temperature controlled at 20°C. Especially for the isomers of leucine and isoleucine, baseline separation was successfully achieved when using the mobile phase consisting of acetonitrile-10 mM ammonium formate with 0.2% formic acid. Under other abovementioned conditions, these two isomeric amino acids were not well separated reciprocally.

3.3. Optimization of MS/MS Conditions

Primary experiments were conducted to acquire the best mass spectrometry conditions of amino acids individually. For this purpose, each reference component of amino acids with the concentration of approximate 100 ng/mL was delivered into the mass spectrometer directly via a syringe pump with the velocity of 7 μL/min. While tuning the standard dissolved in water, most response of amino acids was very low when the sodium-adducted ion, [M+Na]+, was relatively high. On the contrary, when dissolving the standard in 0.1% or higher concentration of formic acid solution, the response of protonated molecule, [M+H]+, was dramatically higher and the metal-adducted ion was barely observed. Consequently, precursor ion (Q1), product ion (Q3), declustering potential (DP), and collision energy (CE) of each amino acid were acquired and optimized.

3.4. Method Validation

As shown in Table 3, the calibration curves of twenty-four amino acids showed good linear regression (all correlation coefficients > 0.995) within a wide range of concentration. The LOD and LOQ of each analyte with the signal-to-noise ratio (S/N) of 3 and 10 were in the range of 0.37–1.94 ng/mL and 0.74–4.49 ng/mL, respectively. The RSD values of intraday and interday precision are summarized in Table 4. From the results, the present method was found to be precise with intraday variability less than 3.14% and interday variability less than 4.27%. The sample solutions were stable within 24 h with stability RSD < 3.82%. The RSD values of repeatability were 0.52%–3.86%. The recoveries of 24 amino acids ranged from 95.27%–107.51% with RSD < 4.03%.

3.5. Quantitative Analysis for Samples

The developed UPLC-MS/MS method was applied for simultaneous determination of 24 free amino acids in 59 batches of three species of duckweed. All the contents are displayed in Table 5. The results showed that the total amount of amino acids varied dramatically among each sample from different regions and species. Total content of AAs in SP25 from Baoding City of Hebei Province reached as high as 1899.50 μg/g, which was much more than that in SP17 of 166.05 μg/g from Lianyungang City of Jiangsu Province. For each amino acid, the content also varied significantly even in the same sample. Coincidentally, the contents of free amino acids that include sulfur in their molecules were much less than others, including cysteine, cystine, and methionine with the concentration of 0.24, 4.25, and 4.20 μg/g in maximum, respectively. Especially for cysteine, it was not found in all of Spirodela polyrhiza and only detected in several samples of other two species. The putative reason was that the thiol in cysteine was not stable and easily oxidized. But in general, there was no significant difference observed in the total content of all FAAs among different species.

Twenty-six minerals, total proteins, and the content of amino acids after hydrolysis are shown in Tables 5 and 6. From the tables, it could be found that the total protein contents in 59 samples ranged from 13.61% to 21.46%. The relative contents of amino acids such as Lys, Ala, Ser, Thr, GABA, Gln, and Asn were all a little higher than other components based on the average value of each species among three different varieties of duckweed. The contents of Cys, Cys2, and Met were higher after hydrolysis. Obviously, it could be seen that the contents of free amino acids in duckweed were much lower than that of protein-bound amino acids. Moreover, the data showed that the contents of free amino acids in duckweed were negatively correlated with binding amino acids.

The elements with higher content in the samples were sodium, calcium, potassium, iron, magnesium, and aluminum. Their concentration ranges were 366808.463–897950.747, 373586.252–799614.153, 169080.807–737471.818, 35506.241–439955.497, 145491.528–308322.609, and 8383.317–34467.794 μg/L, respectively. The concentration ranges of strontium, barium, zinc, manganese and titanium were 154.053–279.990, 32.520–142.280, 22.616–95.305, 2.110–64.537, 16.833–45.839 μg/L. The concentrations of vanadium, chromium, cobalt, nickel, copper, arsenic, selenium, molybdenum, cadmium, tin, antimony, thallium, and lead are very low. It ranged from 0.0 to 3.492 μg/L. Evidently, the contents of macroelements in the samples were higher while they contained much lower trace elements.

3.6. Statistical Analysis

As one of the most important multivariate analysis techniques, PCA was employed to evaluate the variation among three species of duckweed. PCA is an unsupervised pattern recognition method without prior information of the data set and retained maximum variance of multidimensional data while reducing its dimensionality into two- or three-dimension [2931]. Furthermore, all of the data information about amino acids was subjected to a supervised method, OPLS-DA [32], to get more information about different duckweeds. All the raw data were processed by normalization and Pareto scaling (Par) before modeling. Three dimensions were established in both the PCA and OPLS-DA models based on the content of free amino acids and the ratio of free to total amino acids with R2X (cum) value of 0.628, 0.596 and 0.758, 0.676, respectively. As shown in Figures 25, both PCA and OPLS-DA methods demonstrated that Spirodela polyrhiza, Landoltia punctata, and Lemna aequinoctialis were clustered together with each other in the score plot. There was no remarkable difference in different duckweeds.

4. Conclusions

In this study, a simple, fast, and convenient method with UPLC-QTRAP-MS/MS was established for the identification and quantification control of duckweeds by HILIC separation without derivatization. The results of method validation suggested that the developed method was sensitive, accurate, and precise for determination of 24 free amino acids. 59 batches of sample in three species were quantified, and then the data were statistically analysed by using PCA and OPLS-DA. The multivariate analysis illustrated that the content of free amino acids was not significantly different among Spirodela polyrhiza, Landoltia punctata, and Lemna aequinoctialis in pattern recognition. But as a food, more amino acids mean more nutrition it has and more valuable it is. In other words, all of the three species of duckweed could be used as food equally if the absolute amount of AAs was more or less the same. In conclusion, a reliable and feasible method was developed for the nutritional evaluation of duckweeds.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the Chinese National Natural Science Foundation (81073002), the Program of Collaborative Innovation Center of Chinese Medicinal Material Resources Industrialization of Jiangsu Province (2016), and the National “Twentieth Five-Year” Technology Support Program (2011BAI04B06).