Integrated Metabolomics Study of the Milk of Heat-stressed Lactating Dairy Cows

Heat stress (HS) damages the global dairy industry by reducing milk yields and quality, harming health, and damaging the reproduction of dairy cows, causing huge economic losses each year. However, an understanding of the physiological mechanism of HS lactating dairy cows remains elusive. Here, a metabolomics study using LC-MS and 1H NMR spectroscopy was performed to analyze the metabolomic differences in the milk between HS-free and HS dairy cows, and discover diagnostic biomarkers and changes in the metabolic pathway. A total of 53 discriminating metabolites were significantly up- or down-regulated in the HS group compared with the HS-free group (P < 0.05). These biomarkers were involved in pathways of carbohydrate, amino acid, lipid, and gut microbiome-derived metabolism. Comparing these potential biomarkers with previously identified HS candidate biomarkers in plasma, significant correlations between the levels of lactate, pyruvate, creatine, acetone, β-hydroxybutyrate, trimethylamine, oleic acid, linoleic acid, lysophosphatidylcholine 16:0, and phosphatidylcholine 42:2 in milk and plasma were found, indicating that the blood-milk barrier became leaky and the levels of these 10 biomarkers in milk can reflect HS-induced metabolomic alterations in blood. These novel findings can support more in-depth research to elucidate the milk-based changes in metabolic pathways in HS lactating dairy cows.

understanding of the pathophysiological alterations stimulated by environmental changes [12][13][14][15] . In dairy cattle, milk can be very conveniently collected, provides information on the changes in the lactation mechanisms of dairy cows, and can be directly analyzed to determine its nutritional quality. Therefore, milk is regarded as the ideal biological sample for monitoring physiological alterations. Compared to blood sampling, milk sampling is non-invasive and occurs on a daily basis, and thus, dairy cows' metabolic states can be observed in real time. As a result, cows can be monitored and their HS state determined, and prompt management strategies can be implemented to reduce the impact of HS. Hence, the present study used LC-MS and 1 H NMR spectroscopy to identify metabolic differences in the milk from mid-lactation cows with and without HS to identify HS biomarkers and explore the alterations in the metabolic pathways of lactating dairy cows in different HS states. LC-MS in multiple reaction monitoring (MRM) mode was also used to verify the reliabilities of the discriminating metabolites. Partial Pearson's correlation analyses of the candidate biomarkers found in milk with those in blood 16 were performed to track the causes of the disturbed metabolites in milk. The concentrations of cytokines, c-reactive protein, and cytochrome c were detected using enzyme-linked immunosorbent assays (ELISAs) to obtain insights into HS-induced inflammation and apoptosis. Overall, the results of this study improve our understanding of the metabolic alterations in dairy cows exposed to HS. An overview of the study design is shown in Fig. 1.

Results
Comparison of the metabolic profiles of HS-free and HS cows. The orthogonal partial least squares discriminate analysis (OPLS-DA) plots of the metabolomic data showed a clear separation between the HS-free and HS groups, without any overlap ( Fig. 2A,B), indicating that the metabolic profile was significantly altered in the HS milk samples. The 1 H NMR data of the milk serum and fat ( Fig. 2A) and LC-MS data of the milk metabolome (Fig. 2B) for the HS-free and HS groups identified one predictive component and two orthogonal components with satisfactory modelling and predictive abilities of 81.36% < R 2 (Y) < 88.2% and 63.9% < Q 2 (cum) < 86.1%. To avoid model overfitting, cross-validation across three components with 999 random permutation tests was performed and produced intercepts of 0.015 < R 2 < 0.239 and − 0.275 < Q 2 < − 0.190 for all of the data (Fig. 2C,D), indicating that the model was valid 16-18 . Identification of metabolic candidates and their structures. By analyzing the 1 H NMR data, we identified a total of 25 metabolic candidates, which are displayed in Table 1. These metabolites were carbohydrate, amino acid, lipid, and gut microbiome-derived metabolites, suggesting that these metabolic pathways were   altered in the HS group. Figure 3 shows an example of the 1 H NMR spectra of the metabolites identified from milk fat. Analyzing the LC-MS data revealed that 119 metabolites differed between the HS-free and HS groups, with 33 of these metabolites remaining after the elimination of redundant variables using a combination of extracted ion chromatograms and the R-package Collection of Algorithms for Metabolite Profile Annotation (CAMERA) 16,19 . To verify these 33 metabolites, a confirmatory detection was performed using LC-MS in MRM mode. Ultimately, 28 candidates were verified and identified ( Table 2). With the exception of urea, the other 27 compounds were lipid metabolites, indicating that lipid metabolism was disturbed in the HS group.
Metabolic alterations. The 1 H NMR-identified metabolites that differed between the HS and HS-free groups are shown in Table 1, and a heatmap was constructed based on the alterations in these potential candidates (Fig. 4A). HS-induced alterations in carbohydrate-related metabolites were observed for lactate, pyruvate, galactose-1-phosphate, and citrate, the concentrations of which were increased by 1.27-to 1.79-fold in the HS group compared with the HS-free group (P < 0.002), with the exception of fumarate decreased by 0.78 fold in the HS compared with the HS-free group (P < 0.001). The concentrations of isoleucine, proline, orotate, glycine, and phosphocreatine, which are metabolites related to protein or amino acids, were changed by 0.69-to 0.88-fold in the HS group compared with the HS-free group (P < 0.001). Creatine was increased by 1.59-fold in the HS group compared with the HS-free group (P < 0.001). The concentrations of butyrate, acetone, β -hydroxybutyrate (BHBA), FA, PUFA, and UFA, lipid-related metabolites, were increased by 1.08-to 1.38-fold in the HS group compared with the HS-free group (P < 0.02). The concentrations of N-acetylsugar A, N-acetylsugar B, N-acetylsugar C, N-acetylsugar D, and TC, were changed by 0.68-to 0.91-fold in the HS group compared with the HS-free group (P < 0.01). HS-induced perturbations of the gut microbiome-derived metabolism were also identified, with the concentrations of trimethylamine (TMA) changed by 0.52-fold in the HS group compared with the HS-free group (P < 0.001).
The LC-MS-identified metabolites changed by HS are shown in Table 2, and a heatmap was constructed based on the alterations in these potential candidates between the HS-free and HS groups (Fig. 4B). Most of these metabolites are lipids, with the concentrations of oleic acid, linoleic acid, and lysophosphatidylcholine (lysoPC) 16:0 being changed 1.28-to 1.41-fold in the HS group compared with the HS-free group (P < 0.005), and phosphatidylcholine (PC), sphingomyelin (SM), monoacylglycerol (MG), diacylglycerol (DG), and triradylglycerol (TG) changed 0.48-to 0.83-fold in the HS group compared with the HS-free group (P < 0.02). The concentration of urea, which is related to amino acids, was increased by 1.32-fold in the HS group compared with the HS-free group (P < 0.01). Partial Pearson's correlations between the candidate plasma and milk metabolites with corrections for the treatment groups (HS and HS-free). Table 3 lists the correlation coefficients of candidate metabolites in plasma and milk with corrections for the treatment groups (HS and HS-free). Significant correlations were observed between the lactate, pyruvate, creatine, acetone, BHBA, and TMA levels in plasma and in Scientific RepoRts | 6:24208 | DOI: 10.1038/srep24208 milk (P < 0.001). The significance of the correlations for oleic acid, linoleic acid, lysoPC 16:0, and PC 42:2 was less than 0.01. Moreover, significant correlations were also detected between the lactate and pyruvate levels, and acetone and BHBA levels in milk and plasma (P < 0.01). Table S6 lists changes in the concentrations of TNF-α , IL-2, IL-10, IL-12, IL-15, c-reactive protein, p53, Bax, Bcl-2, cytochrome c, caspase-3, caspase-8, and caspase-9 in plasma and milk from the HS and HS-free groups. The levels of all of these proteins were significantly higher in the HS group than in the HS-free group (P < 0.05), with the exception of Bcl-2 and IL-12 being significantly lower in the HS group than in the HS-free group (P < 0.05). Table S7-S8 list the correlation coefficients between the rectal temperature and candidate metabolite levels in plasma and milk, respectively. Significant correlations were observed between the lactate, pyruvate, creatine, proline, lysine, glycine, threonine, isoleucine, leucine, ornithine, citrulline, and arginine levels in plasma and the rectal temperature (P < 0.01), as well as between the lactate, pyruvate, creatine, and citrate levels in milk and the rectal temperature (P < 0.01). Supplementary Table S9 displays the correlation coefficients for the cytokine and candidate metabolite levels in plasma and milk. Significant correlations were noted between p53, Bax, Bcl-2, and cytochrome c and acetone and BHBA in both plasma and milk (P < 0.01). Supplementary Table S10 lists the significant correlations between the cytochrome c and c-reactive protein levels in plasma and milk (P < 0.01).

Discussion
In the present study, 53 potential metabolite biomarkers that could be used to diagnose the HS states of lactating dairy cows were identified. These potential biomarkers were involved in pathways for carbohydrate, amino acid, lipid, and gut microbiome-derived metabolism, indicating that these metabolic pathways were influenced by HS. HS up-regulated the pyruvate and lactate levels in the milk from HS dairy cows in the present metabolomics study, as observed in our previous metabolomics study of the plasma from HS dairy cows 16 . Correlation analyses (Table 3) of the pyruvate or lactate levels in plasma and milk from HS and HS-free cows revealed that the plasma pyruvate and lactate concentrations were significantly correlated with their respective concentrations in milk, indicating that these two metabolites are directly secreted from the blood into the milk via the mammary gland 20 .
The milk citrate concentration was increased in the HS group compared with the HS-free group. This species has been used as a biomarker of energy balance in dairy cows and is correlated with the presence of ketone bodies in milk and de novo FA synthesis 21 . Because the mammary epithelium is impermeable to citrate in both directions, the milk citrate levels reflect an HS-induced disturbance in mammary function rather than a disturbance in general metabolism. The insignificant differences in the plasma citrate levels between the HS and HS-free groups   Table 3. Partial Pearson's correlations between the candidate metabolites in plasma and milk with corrections for the treatment groups (HS and HS-free). *P < 0.01. **P < 0.001. NS, no significant correlations. noted in our previous study 16 further confirm the hypothesis that the milk citrate levels reflect an HS-induced disturbance in mammary function rather than alterations in the blood citrate metabolism. The reduced fumarate concentration in the HS milk samples may have been caused by a disturbed energy metabolism and impaired function of the tricarboxylic acid cycle 22 .
Lactose is the predominant sugar in most species' milk. The role of glucose metabolism in milk production is to transport glucose from the blood into the mammary gland for lactose synthesis, and has a vital function in milk yield relating to the osmoregulation of milk 23,24 . Although the present study did not identify significant differences in the milk lactose concentrations between the HS and HS-free groups, the milk yields of the HS group were greatly reduced (Supplementary Table S4) relative to the HS-free group (P < 0.01). As a result, lactose yields were reduced, suggesting that less blood glucose was available to the mammary glands or lactose synthesis was disrupted.
The increased galactose-1-phosphate concentration in milk was likely attributable to leakage of this component from the mammary epithelial cells into milk, possibly through cell apoptosis, and seems to be strongly dependent on the energy balance states 25 . A negative energy balance (NEB) can greatly increase the apoptotic index in the mammary gland 25 . HS dairy cows are in a NEB state, which may have caused the increased galactose-1-phosphate levels in milk in the present study. The decreased N-acetylsugar A, N-acetylsugar B, N-acetylsugar C, and N-acetyl sugar D levels in the HS group indicated that HS altered their synthesis by reducing the levels of their precursors or influencing the related enzymes.
Decreased concentrations of isoleucine, proline, orotate, glycine, and phosphocreatine were detected in the milk samples from the HS group compared with those from the HS-free group. However, the concentrations of these metabolites were up-regulated in the plasma of the HS dairy cows in our previous metabolomics study 16 . These amino acids may be the main precursors for glucose production via gluconeogenesis or energy production through deamination and oxidation 26 , possibly decreasing their distribution into milk via the mammary gland. Inside the cell, amino acids are involved in metabolic reactions yielding, inter alia, CO 2 , urea, polyamines, and nonessential amino acids (NEAAs). The increased urea concentration in the plasma samples from HS dairy cows indicates that more amino acids were metabolized into urea than into milk proteins inside the mammary gland cells. The rectal temperature was significantly correlated with the proline, lysine, creatinine, threonine, isoleucine, leucine, ornithine, citrulline, arginine, and creatine levels in plasma (Supplementary Table S7), suggesting that HS increased the catabolism of the skeletal muscle protein into amino acids in plasma 27 . Significant correlations were observed between the rectal temperature and lactate and pyruvate levels in both plasma and milk (Supplementary  Table S8), indicating enhanced anaerobic glycolysis.
In lactating cows, NEB is associated with the mobilization of the body's energy reserves, leading to extended lipolysis of adipose tissue and partitioning of non-esterified FA into the bloodstream 25 . This partitioning will influence the FA profile in milk serum, such as by increasing the unsaturated FA (UFA) concentrations, mainly C18:1 (oleic acid) and C18:2 (linoleic acid). The same trend was also observed in our results for the milk fat composition, which showed a higher concentration of UFA and polyunsaturated FA (PUFA) in the milk of the cows in the HS group. In cows, the presence of UFA has been related to inflammatory diseases, such as mastitis and metritis. Cholesterol is transported throughout the body via lipoprotein particles. A previous study on the plasma of lactating dairy cows under HS showed that HS reduced the concentrations of high-density lipoprotein (HDL) and very-low-density lipoprotein/low-density lipoprotein (VLDL/LDL), decreasing their supply to the mammary cells 16 . The concentration of cholesterol in milk fat was lower in the HS group than in the HS-free group, indicating that the synthesis and transport of cholesterol into milk is limited under HS.
The BHBA and acetone concentrations were significantly increased in the HS group compared with the HS-free group. These changes agreed with our previous plasma metabolomics study of HS dairy cows 16 . The partial correlation analysis revealed that the BHBA and acetone concentrations in the milk of lactating dairy cows in the HS and HS-free states were significantly correlated with their plasma concentrations reported in a previous study (P < 0.01) 16 , indicating that the BHBA and acetone in the blood passed unchanged into the milk via the mammary cells. Thus, elevated concentrations of BHBA and acetone in milk can indicate increased blood BHBA and acetone levels, which are well-known markers of the energy status in dairy cows and reflect excessive protein mobilization and insufficient glucose supply 21 .
Increased P53 protein, Bax, cytochrome c, caspase-3, caspase-8, caspase-9, and decreased Bcl-2, all of which are apoptosis-related molecules 28 , were detected in both plasma and milk of HS group. Increased P53 protein can facilitate the biosynthesis of routing-proteins for apoptosis. The balance of Bcl-2 and Bax ratio acts as an important regulator of the mitochondrial response to apoptotic signals, and its imbalance can lead to cytochrome c release, predicting cell apoptosis. The shear activation of caspase-8 can initiate caspase cascade reaction. Caspase-3 and caspase-9 are regarded as apoptosis executioner and initiator, respectively. The BHBA and acetone concentrations were significantly correlated with P53 protein, Bax, Bcl-2 and cytochrome c in both plasma and milk (Supplementary Table S9), suggesting that the NEB increased apoptosis 29 .
The concentrations of PC, including PC 40:1, PC 40:0, PC 42:2, PC 44:3, and PC 44:1, were lower in the HS group than in the HS-free group, whereas the concentrations of lysoPC 16:0 were higher in the HS group than in the HS-free group. LysoPCs are metabolites of PC catabolism, which is regulated by phospholipases A 1 , A 2 , and D 30,31 . The observed alterations in the PC and lysoPC concentrations were consistent with our previously plasma metabolomics study of HS dairy cows 16 . Other lipid-related metabolites in milk, such as SM, MG, DG, and TG were also down-regulated in the HS group. Taken together, these findings clearly show that HS induces overall alterations in the lipid metabolism in both the blood and milk of lactating dairy cows.
There was a strong correlation between the TMA levels in milk and plasma, indicating that alterations in the plasma TMA levels 16 influenced its concentration in milk (Table 3). We found that HS caused a marked reduction in the milk TMA levels in HS cows compared with the HS-free cows (Table 1). This leads us to believe that HS induced variations in the numbers and/or activities of intestinal microbes.
Scientific RepoRts | 6:24208 | DOI: 10.1038/srep24208 The immune function of cows can be affected by HS 32 . Increased c-reactive protein, TNF-α , IL-2, IL-10, IL-15, and decreased IL-12 in both the plasma and milk of HS group are regarded as a response to HS-induced inflammation 29,33,34 , confirming the impairment of the immune system in the HS group. The plasma levels of cytochrome c and c-reactive protein were significantly correlated with their respective concentrations in milk (Supplementary Table S10), suggesting that HS induced an inflammatory response of the whole body of dairy cows 33 . During a mammary immune response, the blood-milk barrier is compromised 35 . Therefore, our previously reported up-or down-regulation of the lactate, pyruvate, creatine, acetone, BHBA, TMA, oleic acid, linoleic acid, lysoPC 16:0, and PC 42:2 levels in the plasma 16 were also observed in milk in the present study (Table 3). Although the levels of these 10 candidate biomarkers in milk were significantly correlated with their respective concentrations in blood, the remaining 43 potential biomarkers (Tables 1 and 2) in milk did not show similar correlations, indicating that the former 10 biomarkers can reflect HS-induced perturbations in blood metabolites, whereas the remaining 43 reflect a disturbed mammary function rather than general metabolic disturbances.

Conclusion
Because the levels of the candidate biomarkers lactate, pyruvate, creatine, acetone, BHBA, TMA, oleic acid, linoleic acid, lysoPC 16:0, and PC 42:2 in milk showed strong correlations with their respective concentrations in the blood of HS-free and HS dairy cows, the detection of these milk metabolites can diagnose HS states and HS-induced perturbations in general metabolism. Other potential biomarkers can reflect HS-induced alterations in the mammary function of HS dairy cows. The results identify potential biomarkers that could be used to accurately monitor the HS states of lactating dairy cows and merit further investigation to elucidate the physiological mechanisms underlying the HS-induced changes in metabolic pathways. This could lead to better management practices for dairy cows exposed to HS.

Materials and Methods
All experiments involving animals were conducted according to the principles of the Chinese Academy of Agricultural Sciences Animal Care and Use Committee (Beijing, China), which approved the study protocols. , where T db is the dry-bulb temperature (°C) and RH is the relative humidity (%). Forty-four Holstein cows in second parity at mid-lactation that were fed the same diet were used in this trial. The HS-free group consisted of 22 cows, with milk samples obtained in the spring season under a THI of 50-55 for one month. The HS group consisted of 22 cows, with samples obtained in the summer season, after THI gradually increased from 68 to 80 over 1 month and remained stable at 80 for 1 week. Morning milk samples were collected before feeding, and stored at − 80 °C until further use. The detailed characteristics of the selected dairy cows, feed composition, temperature, humidity, rectal temperature, respiration rate, and production characteristics of the two groups of dairy cows are shown in Supplementary Table S1-S4 referenced to the reported paper 16 . NMR spectroscopic analysis of milk serum and fat. Milk samples were thawed at room temperature, homogenized, and centrifuged at 20000 g for 30 min at 4 °C. All the samples were analyzed at 298 K using a VARIAN VNMRS 600 MHz NMR SPECTROMETER (Varian Inc., Palo Alto, CA) operating at 599.871 MHz using a 5-mm inverse-proton (HX) triple resonance probe with z-axis gradient coil. The 1 H NMR spectra of the milk serum were recorded using the water-suppressed standard 1D CPMG pulse sequence (RD-90°-(τ -180°-τ ) n-ACQ), where a fixed total spin-spin relaxation delay 2nτ of 320 ms was applied to attenuate the broad NMR signals from slowly tumbling molecules (such as proteins) and retain those from low-molecular weight compounds and some lipid components 36 . The free induction decays (FIDs) were collected via 64 K data points with a spectral width of 12000 Hz and 128 scans. The FIDs were zero-filled to double size and multiplied by an exponential line-broadening factor of 0.5 Hz before Fourier transformation (FT). The 1 H NMR spectra of milk fat do not require presaturation of the water resonance. Here, we used a simple 901 pulse-acquire sequence because we were measuring fully relaxed spectra. We collected 64 transients into 32,768 data points; the other parameters were as follows: relaxation delay = 2 s, spectral width = 20 ppm, and acquisition time = 1.36 s. The metabolites were identified by inserting the experimental spectra into the Chenomx spectral database (Edmonton, AB, Canada), and comparing them with the spectra of standard compounds. were as follows: ESI source voltage 5.5 kV; vaporizer temperature, 550 °C; drying gas (N 2 ) pressure, 60 psi; nebulizer gas (N 2 ) pressure, 60 psi; curtain gas (N 2 ) pressure, 30 psi; and declustering potential, 50 V. The scan range was m/z 60-1,000. Data acquisition and processing were performed using Analyst ® TF 1.6 Software (Applied Biosystems/MDS Sciex). Auto-calibration was performed using the Calibrant Delivery System in both positive and negative ion modes. The metabolites were identified using the information-dependent acquisition mode in the MS/MS analyses. The collision energy was 35 eV. High-resolution MS, isotope abundance ratios, MS/MS, the Human Metabolome database, the METLIN database 37 , a literature search, and comparisons with standards were employed to identify the ion structures. Quality control (QC) samples were detected once every 10 LC-MS runs to monitor the reproducibility of the instrument. The plasma samples from the different groups were randomly alternated during the analysis. Auto-calibration was performed using the Calibrant Delivery System in both positive and negative ion modes.

LC-MS analysis of milk serum.
Multivariate data analysis. The raw 1 H NMR spectra were manually corrected for phase and baseline distortions using TopSpin software (version 3.0; BrukerBiospin) and were referenced to the anomeric doublet resonance of lactose (δ 5.233 ppm) for the milk serum spectra and the 3-(trimethylsilyl) propionic-2,2,3,3,d4 acid sodium salt (δ 0.0 ppm) signal for the milk fat spectra. The 1 H NMR spectra of all specimens were binned into 0.001 ppm integral regions and integrated over the 0.5-7.0 ppm region with the exclusion of water (δ 4.7-5.12) using the AMIX software package (version 3.8.3, BrukerBiospin). The spectra were normalized to the total sum of the spectral integrals to compensate for differences in sample concentrations. The two-dimensional data matrices from the milk serum and fat data were integrated into one Excel file for the subsequent SIMCA-P analysis.
The raw LC-MS data files (.wiff) were converted into mzXML format using ProteoWizard (http://metlin. scripps.edu/xcms/download/pwiz/pwiz.zip). The files were processed using an open-source XCMS package (version 1.20.1) in R statistical software (version 2.10.0) for peak discrimination, filtering and alignment 38 . During XCMS implementation, the R-package CAMERA was used to annotate the isotope, adduct, and product ion peaks. The resulting two-dimensional matrices, including the observations (sample names) in columns, the variables (m/z-retention time pairs) in rows, and the peak areas, were saved as CSV files. The IS was used for data quality control (reproducibility) and normalization.
The resulting data sets were then imported into the SIMCA-P 13.0 software package (Umetrics AB, Umeå, Sweden); pareto and centre scaling were applied to the LC-MS and 1 H NMR data, respectively, to reduce the noise and artifacts in the models. Principal component analysis (PCA) was performed on the scaled data to visualize the global clustering and separation trends or outliers. PLS-DA models were applied to validate the model against overfitting through 999 random permutation tests. Discriminating variables were selected according to the S-plots, variable importance in projection values (VIP > 1), and raw data plots in the OPLS-DA models [16][17][18] . Furthermore, independent t-tests (P < 0.05) (SPSS version 13.0) were used to determine the significance of each metabolite in discriminating the HS group from the HS-free groups.
To analyze the correlations between the levels of all the potential biomarkers discovered in milk with their respective concentrations in plasma reported previously from specimens collected at the same time, partial Pearson's correlation coefficients were calculated at two confidence levels (0.01 and 0.001) by correcting for the HS group in the first order using the ParCorA software program, which can be downloaded at http://mendes.vbi. vt.edu/tiki-index.php?page. Cluster 3.0 and Treeview were separately used for cluster analysis and visualization, respectively.

UFLC-MS/MS-based verification test.
The chromatographic gradient elution conditions were unchanged. The eluent was injected into a triple quadrupole-trap mass spectrometer equipped with an ESI source (QTRAP 5500, Applied Biosystems/MDS Sciex). All of the potential candidates discovered by LC-MS were detected by UFLC-MS/MS in MRM mode. For each analyte of interest, the collision energies and precursor/ fragment ion pairs were pre-optimized to generate an optimal signal-to-noise ratio. The MRM transitions are listed in Supplementary Table S1.
ELISAs for cytokines, c-reactive protein, and cytochrome c. The cytokines, c-reactive protein, and cytochrome c concentrations in both plasma and milk were analyzed using commercially available bovine-specific ELISA kits (Abcam ® , Cambridge, MA, USA), according to the manufacturers' instructions.