Corrigendum: Structural and functional changes of bioactive proteins in donor human milk treated by vat-pasteurization, retort sterilization, ultra-high-temperature sterilization, freeze-thawing and homogenization
- 1Emma Children's Hospital, Amsterdam UMC, Amsterdam, Netherlands
- 2Food Quality and Design Group, Wageningen University and Research, Wageningen, Netherlands
- 3Laboratory of Biochemistry, Wageningen University, Wageningen, Netherlands
- 4Department of Chemical Biology and Drug Discovery, Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Utrecht, Netherlands
- 5Danone Nutricia Research, Utrecht, Netherlands
Donor human milk is generally processed by holder pasteurization (HoP) at 62. 5◦C for 30 min. This temperature-time combination is sufficient for eliminating pathogens in donor milk, but also negatively affects several bioactive milk components. Long heating up times may further affect the bioactive properties of pasteurized milk. High-Temperature-Short-Time (HTST), a treatment with shorter processing times (72◦C for 15 sec), was investigated as a suitable alternative to HoP. In addition, pasteurization methods that follow the same temperature regime but with varying heating up times were compared. Human milk samples from four different donors were combined into one pool, which was then used to perform all analyses. The effects of these methods on the levels and functionality of immunoglobulin A, lactoferrin, lysozyme and bile salt-stimulated lipase, were evaluated with LC-MS/MS-based proteomics and activity assays, while the pasteurization efficacy was evaluated with an alkaline phosphatase test. HoP, a treatment with long processing times, times, caused the highest reduction in all proteins studied (reduced by 50–98%). Compounds such as lactoferrin and bile salt-stimulated lipase that are more sensitive to heat treatments were better retained with HTST, but their levels and functionality were still significantly lower than those of untreated donor milk (52 and 81% reduction of lactoferrin and bile salt-stimulated lipase activity, respectively). Our findings showed that a treatment with considerably shorter processing times, such as HTST, may reduce the thermal damage caused to the bioactive proteins compared to HoP, without affecting pasteurization efficacy. Since the vast majority of the donor human milk banks that are currently operating on a global level apply HoP to donor milk, our findings may provide relevant information for the optimization of donor milk processing.
Introduction
Tailored to each infant's need, mother's own milk represents the optimal source of neonatal nutrition. The unique composition of human milk (HM) promotes healthy infant development and growth, while the numerous bioactive factors it contains, protect infants against infections and various diseases (1, 2). HM has been shown for example to be bacteriostatic against a number or bacteria, which is mainly attributed to the functionality of immunoglobulin A (IgA), lactoferrin (LTF) and lysozyme (LYZ). LTF, amongst others, limits the availability of the iron required for the growth of iron-dependent pathogens, LYZ disrupts cell walls in gram-positive bacteria, and secretory IgA is generally considered as the main antibody system in HM. Synergistic effects of these proteins have also been reported (3–5). Especially for premature infants, when mother's own milk is not available, donor human milk (DHM) represents the best alternative form of nutrition. Current evidence suggests that DHM protects against necrotizing enterocolitis when compared to infant formula, while its provision is correlated with improved long-term outcomes in preterm infants (2, 6). DHM should be provided though established human milk banks (HMBs) that can ensure its safety (7).
Before provision, HMBs subject DHM to holder pasteurization (HoP), a low-temperature long-time heat treatment at 62.5◦C for 30 min (8). HoP effectively inactivates potential viral and bacterial agents, and it is known as the method recommended for DHM treatment in all international human milk banking guidelines (8, 9). However, this method has been shown to cause substantial losses in various bioactive milk components, due to the thermally induced denaturation occurring in these components (8, 10). In fact, as previously reviewed, losses from 20 to 90% in the levels and functionality of immunoglobulin A (IgA), lactoferrin (LTF), and lysozyme (LYZ) have been reported after this treatment (10). In addition, HoP results in the total inactivation of bile salt stimulated lipase (BSSL), a heat-labile enzyme involved in fat absorption and infant metabolism (10). The reported lower growth rates in preterm infants fed with DHM that underwent HoP, may be the result of this detrimental effect (11).
Different devices are currently being employed by HMBs to perform HoP. Standard pasteurizers are most commonly used, with water as the heating medium (8, 12). The ideal HoP process should be comprised of a rapid heating up phase to 62.5◦C, a phase of constant temperature (30 min) and finally a rapid cooling down phase to <10◦C (8, 9). A recent study investigating human milk banking practices in Europe, revealed that in contrast to the recommendations on pasteurization performance, DHM is currently exposed to long processing times (12). This could result in higher losses of DHM components. More specifically, losses of 1.6, 1.7, and 2.4% were documented for IgA, LYZ and LTF, respectively, for every minute spent at 62.5◦C (13). Therefore, since the duration of the heat treatment and the temperature at which DHM is exposed highly affects the preservation of its bioactive components, ensuring short heating up times seems essential when using thermal techniques.
Various studies report that high-temperature-short-time (HTST) treatment may be a suitable alternative to HoP, as it was found to provide similar microbial reduction (e.g., in E.coli, S.aureus, S. epidermidis, E. faecalis, P. aeruginosa, L. monocytogenes, S. agalactiae and C. sakazakii counts) and to better retain the DHM bioactive proteins (14–19). HTST is usually performed by heating the milk at 72◦C for a duration of 15 sec. This method is well-established in the dairy industry, and it usually involves the rapid heating of a thin layer of milk in a continuous flow system (20). The shorter treatment time as well as the shorter exposure time at the processing temperature may be the reason of the promising results reported after this treatment (8). However, substantial information on whether shorter heating up times during pasteurization can positively influence the retention of the DHM bioactive components, is at present lacking.
Our aim was to determine 1. the effect of shorter heating up times and 2. the effect of shorter pasteurization duration at a higher temperature on the preservation of DHM quality, with emphasis on the IgA, LTF, LYZ and BSSL levels and bioactivity by comparing pasteurization methods that follow the same temperature regime.
Materials and Methods
Milk Samples
DHM samples were collected from the Dutch Human Milk Bank (located at Amsterdam UMC, Amsterdam, The Netherlands). Written informed consent was received from all donors before recruitment. Donor screening, milk expression and collection was conducted as per standardized protocols that adhere to internationally published guidelines (9). The samples were expressed with a breast pump and were subsequently collected in disposable bisphenol A-free bottles (Sterifeed, Medicare Colgate Ltd, Devon, England). The samples were immediately placed in a freezer at −20◦C and were transported to the HMB, at a temperature of −20◦C (21). All samples were placed in a refrigerator (4◦C) overnight, before analysis. Next, milk samples from four different donors were combined into a pool, to ensure sufficient DHM amounts for all analysis. This pool was then divided into two aliquots (600 mL); one that was centrifuged at 6,500 × g for 30 min at 4◦C (with rotor 16.250, Avanti Centrifuge J-26 XP, Beckman Coulter, USA) to remove the fat, and one that remained unskimmed. Then, both the unskimmed and the skimmed samples underwent heat treatment (with the methods mentioned in section 2.2), except for the untreated control samples. Skimming before heat treatment was done for the assays focusing on milk serum proteins, to avoid interference from the milk fat globule membrane. After treatment, all samples were cooled in an stationary ice-water bath and stored at −20◦C. Figure 1 illustrates the experimental approach used in this study. The heating procedures were performed as two independent experiments (biological replicates).
Figure 1. Schematic workflow indicating the experimental approach used. HTST, FH, T99, T74, HoP and HoP_r stand for high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively. For each heating method, DHM treatment was done in biological duplicate. Each sample was then aliquoted for the subsequent analyses, which were performed in technical duplicate.
Heat Treatments
To compare the effects of processing times and temperatures on the DHM bioactive proteins, six pasteurization methods with different heating profiles were conducted; (1) HoP, (2) HoP with rapid heating up-times (HoP_r), (3) HTST, (4) flash heating (FH), (5) HTST with a thermomixer preheated at 99◦C (T99), and (6) HTST with a thermomixer preheated at 74◦C (T74). An overview of the different time-temperature profiles is shown as Supplementary Material (Supplementary Figures 1, 2).
Holder Pasteurization
A 130 mL single-use polypropylene bottle (Beldico SA, Marche-en-Famenne, BE) was filled with 100 mL of DHM and was pasteurized at 62.5◦C for 30 min, in a shaking water bath (SW22, Julabo GmbH, DE) at 150 rpm. To achieve a rapid cooling phase, the sample was placed in an ice-water bath after treatment, until reaching a core temperature of 4◦C. The temperature was recorded by a temperature data logger RS PRO 1,384 (RS Components B.V., The Netherlands). Three thermocouples were placed to monitor the temperature of the milk and of the two baths, during the whole process.
Holder Pasteurization With Rapid Heating Up-Times
The HTST system described in section High-Temperature Short-Time was used to achieve rapid heating up to the pasteurization temperature (62.5◦C). The milk (50 mL) was pumped to the heating section and once its temperature reached 62.5◦C, it was transferred in a shaking water bath for 30 min, as described in Holder Pasteurization.
High-Temperature Short-Time
A laboratory scale pasteurizer was built to simulate continuous HTST pasteurization. The system included a peristaltic pump (Watson Marlow 505S, Hudson, MA, USA), as well as a heating, a holding, and a cooling section (Figure 2). All sections were connected to an RS PRO 1,384 temperature data logger. The pump, which ran at a speed of 35 rpm, was connected with a plastic tube (ø 4 mm) to a copper heating coil (810 mm, ø 4 mm), that was fully submerged into a water bath (heating section). A thermocouple was placed at the end of the coil to determine whether the milk (50 mL) leaving the heating section indeed reached a temperature of 72◦C. The coil was then connected to a plastic tube (ø 4 mm) that remained submerged into a second water bath (holding section, 15 sec). A second thermocouple was placed at the end of the holding section. The milk passed then through a copper coil (1,395 mm, ø 4 mm) which was submerged in an ice-water bath (2◦C). Finally, a third thermocouple was linked with the coil to monitor the temperature at the end of the cooling phase. The cooled milk was then dispensed into a sterile bottle (Beldico SA, Marche-en-Famenne, BE).
Flash Heating
The set up used for this treatment was based on previous reports (22, 23). A 250 mL Duran bottle (GL 45, DWK Life Sciences GmbH, DE) was filled with 100 mL of DHM and was placed into an aluminum pan (2 L) containing 450 mL of water. The water and the submerged DHM bottle were heated simultaneously using an electric hot plate, until the water reached 100◦C. The bottle was removed once the milk's temperature reached 71.5◦C and was transferred into an ice-water bath (2◦C). The temperatures of both milk and water were recorded by a temperature data logger (RS PRO 1,384).
Thermomixer Preheated at 99◦C and 74◦C
A shaking thermomixer (12 ml block, Eppendorf Thermomixer R Mixer, Hamburg, DE) was preheated to either 99◦C or 74◦C. Eight 12 mL Greiner tubes (Greiner Bio-One International GmbH, DE), each filled with 10 mL of DHM, were placed in the thermomixer once the desired temperature was reached (99◦C or 74◦C). The samples were shaken during the whole treatment (800 rpm), while the temperature was monitored with a temperature data logger. When the core temperature of the milk reached 72◦C, the samples were removed from the thermomixer and after 15 s, they were then immediately transferred in an ice-water bath (2◦C) for a rapid cooling down.
Native Milk Serum and Total Protein Content
To obtain the native serum protein fraction from the skimmed DHM samples, caseins and denatured proteins were precipitated by acidifying the samples with HCl (1 mol/L) under stirring, until a pH of 4.6 was reached (24). After leaving it for 30 min at 4◦C to equilibrate, the samples were ultracentrifuged for 90 min at 100,000 × g and at 30◦C (Optima L-80, Beckman Coulter, USA). Next, the supernatant with the native serum protein fraction was collected and the pellet containing micellar casein and denatured proteins was discarded. Finally, total protein content was determined with the bicinchoninic acid (BCA) assay kit (Thermo Fisher Scientific, USA), as specified in the manufacturer's instructions.
Proteomics by Liquid Chromatography With Tandem Mass Spectrometry (LC-MS/MS)
Filter Aided Sample Preparation for Proteomics
Before LC-MS/MS measurements, the samples were prepared with the FASP method, as previously described (25). Briefly, after the samples (1.0 μg/μL) were reduced with 15 mM dithiothreitol, they were first diluted by urea (8 M) in 100 mM Tris/HCl (0.1 M, pH 8.0) and then alkylated with the addition of 20 mM of acrylamide (0.2 M). Next, 100 μL of the alkylated samples were transferred to a Pall 3 K omega filter (10–20 kDa cut off, OD003C34, Pall corporation, Port Washington NY, USA) and were centrifuged for 30 min at 16,900 × g. The samples were then washed with 50 mmol/L NH4HCO3 (ABC) and were centrifuged again at the same conditions. The filter units were placed into new low-binding tubes (2 mL) and the samples were subjected to overnight digestion by the addition of 100 μL trypsin in NH4HCO3 solution (5 ng/μL). Next, a centrifugion step at 16,900 × g for 30 min followed, and another one at the same conditions after the addition of 100 μL 1 mL/L HCOOH in water on top of the filter unit. Finally, 3 μL of 10% trifluoroacetic acid was added to the filtrate to adjust sample pH to around 3. Before injection into the LC-MS/MS system, all samples were stored at −20◦C.
LC-MS/MS Proteomics
All analyses were carried out by the department of Biochemistry at Wageningen University and Research. The parameters used were the same as previously reported (25). The samples were directly injected on a 0.10*250 mm ReproSil-Pur 120 C18-AQ 1.9 μm beads analytical column prepared in house, at 800 bar. Elution of the peptides was done at a flow of 0.5 μL/min, using an acetonitrile gradient (9–34% acetonitrile in water with 1 ml/L formic acid in 50 min). The eluent was then ejected trough the tip of a needle with an electrospray potential of 3.5 kV. Full scan FTMS in positive mode between m/z 380 and 1,400 were measured using Q Exactive HF-X mass spectrometer (Thermo Electron, San Jose, CA, USA). MS/MS scans of the twenty most abundant multiply charged peaks, were measured in data-dependent mode. MaxQuant software (v1.6.3.4) was used to analyze the obtained MS data, against the Uniprot human protein database and a database containing the sequences of common contaminants (26). Protein modifications were set as propionamide (C) (fixed) and oxidation (M) (variable), while enzyme specificity was set for trypsin and a maximum of two missed cleavages, 20 ppm peptide tolerance first search, 4.5 ppm main search and 20 ppm MS/MS fragment match tolerance. Requirement for further analysis was the protein identification by a minimum of two peptides that had at least one unique and one unmodified peptide. Proteins detected in less than half of the samples as well as keratins and trypsin were removed from the final list of identified proteins.
BSSL Activity
The assay used to determine BSSL activity was based on a previously published method, with minor modification (27). DHM lipase activity is determined fluorometrically through the utilization of two synthetic substrates; 4-methylumbelliferyl butyrate (4-MUB) and 4-methylumbelliferyl laurate (4-MUL). Defatted DHM samples were preincubated at 40◦C for 3 min, under 800 rpm, in a ThermoMixer (SmartBlock 1.5 ml, Eppendorf, Hamburg, DE). A stop solution of GuHCl (8 M) and HCl (1 M) in water was then used to stop the conversion of the added substrate and a neutralizing solution with Bis-tris (1M), NaOH (0.85 M) and EDTA (0.25 M) in water was added next to clarify the samples. The fluorescence released was measured by using a fluorimeter (SpectraMax ID3, Molecular Devices, San Jose, CA, USA) at an excitation of 355 nm and an emission of 460 nm.
Alkaline Phosphatase
The method used for the detection of ALP was according to an international standard protocol (ISO/TS 6090|IDF/RM 82A:2004). Finally, ALP activity was measured in a p-nitrophenol calorimeter (Lovibond APTW/7, Tintometer GmbH, Dortmund, DE).
Bacteriostatic Properties
In order to assess the effect of the different methods on the functionality of the three major HM antimicrobial proteins (IgA, LTF, LYZ), we evaluated the growth rate of two bacterial strains known to be inhibited by these proteins. Fresh cultures of Escherichia coli K12 (DSM 498, DSMZ, Braunschweig, Germany) and Staphylococcus aureus (ATCC 6,538, American Type Culture Collection, Manassas, USA) were prepared from frozen stocks in nutrient broth overnight at 37◦C (CM0001, Thermo Fisher Scientific, Massachusetts, USA). Bacterial pellets were obtained after a centrifugation step of 10 min at 4,000 × g (Microcentrifuge 5890R, Eppendorf, Hamburg, Germany) and were subsequently dissolved in peptone physiological salt solutions (PFZ; Tritium Microbiology, The Netherlands). Optical density was determined by using a spectrophotometer (Cary 50 UV–Visible Spectrophotometer, Agilent Technologies, USA). Next, E. coli and S. aureus cultures were inoculated into untreated samples and samples that were first subjected to heat treatment with the different methods, to a concentration of ~103 colony forming units (CFU)/mL. DHM samples inoculated with E. coli were incubated at 37◦C for 2 h and the samples inoculated with S. aureus for 4 h, at the same temperature. All samples were plated in duplicate onto selective media; violet red bile glucose agar for E. coli (CM0107B, Thermo Fisher Scientific, Massachusetts, USA) and mannitol salt agar for S. aureus (CM0085B, Thermo Fisher Scientific, Massachusetts, USA) and were then incubated overnight at 37◦C. Bacterial counts were determined by colony counting (CFU/mL) while the growth rare per hour was measured as , were Nt = bacterial counts after either 2 h or 4 h of incubation, N0 = bacterial counts immediately after incubation and t = incubation time.
Data Analysis
GraphPad Prism software 8.0 (GraphPad Inc., La Jolla, CA) was used for data analysis and visualization. The effects of the different treatments were compared by ANOVA and Tukey's HSD for post-hoc tests. Protein retentions (%) were determined after dividing the concentrations after treatment by the concentration of untreated samples, multiplied by 100. The intensity based absolute quantification (iBAQ) values obtained with MaxQuant, were analyzed in Perseus software (v.1.6.2.1, Martinsreid, Germany). The iBAQ values are considered as suitable indicators for absolute protein concentrations, as the values refer to the sum of all peptide intensities divided by the number of theoretically generated tryptic peptides (28). Perseus was used to estimate significant differences in the protein pattern after treatment, by Student's t-tests with permutation-based false discovery rate (FDR) correction. The correlation between the levels of IgA, LTF and LYZ retained and the bacterial growth rate was also determined, by creating a correlation matrix with R version 3.4.0 (29). A p-value < 0.05 was used to indicate significant differences among the compared groups. The analyses were performed in duplicate for each sample and all data are shown as mean ± standard deviation of two independent experiments.
Results
Temperature Profiles
The temperature profiles of the different treatments were broken down into three sections: the heating up time to the pasteurization temperature (referred to as heating up time), the time that DHM was held at this temperature, and the time required for DHM to cool down to 4◦C (Table 1). The time DHM spent above 55◦C was used as an indicator of the thermally induced protein denaturation that usually occurs above this temperature (30). During HoP, the samples were exposed to temperatures above 55◦C for about 44 min, which was the longest exposure observed among the different treatments. In contrast, HTST-treated DHM was exposed above this temperature for only 33 sec.
Native Milk Serum Protein Concentration and Quantitative Analysis of the Milk Serum Proteome
A combination of acidification and ultracentrifugation was applied to specifically isolate the native milk serum proteins. The total native milk serum protein concentration of the untreated and the treated samples is shown in Figure 3. Of the treatments tested in this study, only HoP caused a significant decrease in native protein concentration (p < 0.05), when compared to the untreated samples. In addition to the total native protein content, the native protein profile was assessed as well by LC-MS/MS. The impact of the different heat treatments on the DHM native protein profile was then visualized by a clustered heat map of the obtained iBAQ values (Figure 4). According to the clustering pattern (Figure 4), the native protein profile of the samples that were the longest exposed to temperatures >55◦C (samples treated with T74, HoP or HoP_r) differed the most from the protein profile of the untreated samples.
Figure 3. Total native milk serum protein concentration (mg/ml) of untreated and the differently heat-treated DHM samples. The analyses were performed in duplicate for each sample and all values are presented as mean ± standard deviation of two independent experiments. *Expresses statistically significant difference to untreated samples (p < 0.05). UN, HTST, FH, T99, T74, HoP and HoP_r stand for untreated DHM, and DHM treated with high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively.
Figure 4. Heat map indicating differences in protein profile of the differentially heated DHM samples. The color scale is based on z-score normalized iBAQ values and each row represents individual proteins. Hierarchical clustering was performed using a Euclidean distance metric. The analyses were performed in duplicate for each sample and all values are presented as mean ± standard deviation of two independent experiments UN, HTST, FH, T99, T74, HoP and HoP_r stand for untreated DHM, and DHM treated with high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively.
IgA, LTF and LYZ Levels After Processing
Overall, the retentions of the three proteins showed a decreasing tendency with increasing exposure time above 55◦C. Compared to untreated DHM, FH, T99, HoP_r, T74 and HoP significantly reduced the IgA, LTF and LYZ levels (p < 0.05), with average retention rates between 19 and 64% (Figure 5). HoP preserved IgA, LTF and LYZ levels the least, but HoP with rapid heating up times was shown to improve their retention, although the differences observed between the two methods were non-significant (mean ± SD retention rates of IgA, LTF and LYZ after HoP_r and HoP; 50 ± 5% vs. 44 ± 4%, 26 ± 18% vs. 19 ± 4%, 60 ± 18% vs. 50 ± 6%, respectively, p > 0.05). HTST, the treatment with the shortest processing times, was found to better retain the levels of the three proteins; the concentrations of IgA and LYZ were not significantly different from those of untreated samples (74 ± 9 and 82 ± 19%, respectively) but the concentration of LTF was significantly reduced (48 ± 13% of LTF was retained, p < 0.05). However, the concentration of LTF after HTST was significantly higher than after HoP and T74 (p < 0.05). Our data also showed that the IgA levels were significantly higher after HTST than after HoP (p < 0.05), but for LYZ, no significant differences were observed between the two treatments.
Figure 5. Effect of HTST, FH, T99, HoP_r, T74 and HoP on the IgA, LTF and LYZ content. Protein retention (%) was calculated based on the iBAQ intensities obtained by LC-MS/MS. The dotted line represents the untreated values (100%). The analyses were performed in duplicate for each sample and all values are presented as mean ± standard deviation of two independent experiments. *Expresses statistically significant differences to untreated samples (p < 0.05). UN, HTST, FH, T99, T74, HoP and HoP_r stand for untreated DHM, and DHM treated with high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively.
BSSL Level and Activity After Processing
The effects of the different heat treatments on BSSL level and activity were determined by means of LC-MS/MS and an activity assay. All treatments caused a major reduction on the enzyme's level and activity, with respect to untreated DHM (p < 0.05, Figure 6). BSSL was affected the most by HoP (LC-MS/MS, 2% and activity assay, 4%), but the values obtained after FH, T99, HoP_r and T74 were comparable to those of HoP (p > 0.05). HTST retained significantly higher BSSL level and activity than the other treatments (LC-MS/MS, 9% and activity assay, 19%, p < 0.05).
Figure 6. Effect of HTST, FH, T99, HoP_r, T74 and HoP on BSSL levels. BSSL retention (%) was calculated based on the iBAQ intensities and a lipase activity assay. Untreated values were set at 100% (dotted line). The analyses were performed in duplicate for each sample and all values are presented as mean ± standard deviation of two independent experiments. *Expresses statistically significant differences to untreated samples (p < 0.05). UN, HTST, FH, T99, T74, HoP and HoP_r stand for untreated DHM, and DHM treated with high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively.
ALP Activity After Processing
ALP is a very heat-sensitive enzyme and it is expected to be completely inactivated when the pasteurization is adequate (19). Untreated DHM samples exhibited a mean ALP activity of 0.257 ± 0.049 U/ml, whereas all heat-treated samples were below the detection limit.
Retention of Bacteriostatic Properties After Processing
To assess the impact of the different heat treatments on the DHM bacteriostatic capacity, the growth rates of S. aureus and E. coli were evaluated, in untreated and in heat-treated samples (Figure 7). The lowest growth rate for both strains was documented in untreated DHM (1.7 ± 0.36 and 3.6 ± 0.04-fold per hour, for S. aureus and E. coli, respectively), which indicates that untreated DHM samples exhibited the highest bacteriostatic capacity among all samples (p < 0.05). In contrast, HoP caused the highest reduction in bacteriostatic capacity (S. aureus and E. coli growth rate, 2.9 ± 0.01 and 5.53 ± 0.09-fold per hour, respectively, p < 0.05). Compared to the untreated samples, HTST resulted in a comparable S. aureus growth rate (2.0 ± 0.17-fold per hour, p > 0.05) but the E. coli growth rate increased significantly after this treatment (4.2 ± 0.07-fold per hour, p < 0.05). Similarly, the bacterial growth rate was significantly increased after FH, T99, HoP_r and T74 (p < 0.05). When compared to HoP, S. aureus growth was significantly lower after HTST, FH, T99, HoP_r and T74, but when E. coli growth was assessed, that was the case only for the samples after HTST, FH and T99 (p < 0.05).
Figure 7. Growth rate per hour of S. aureus and E. coli in untreated DHM samples and after HTST, FH, T99, HoP_r, T74 and HoP. The analyses were performed in duplicate for each sample and all values are presented as mean ± standard deviation of two independent experiments. *Expresses statistically significant differences to untreated samples (p < 0.05). UN, HTST, FH, T99, T74, HoP and HoP_r stand for untreated DHM, and DHM treated with high-temperature short-time, flash heating, high-temperature short-time with a thermomixer preheated at 99◦C, high-temperature short-time with a thermomixer preheated at 74◦C, holder pasteurization and holder pasteurization with rapid heating up-times, respectively.
Considering that the IgA, LTF and LYZ levels decreased while the S. aureus and E. coli growth rate increased, a negative correlation between bacterial growth rate and the retention of the three major antimicrobial proteins is expected, as was indeed found (Figure 8). The strongest negative correlation was observed between bacterial growth rate and the levels of LTF and LYZ (LTF; r = −0.91 and r = −0.96, LYZ; r = −0.81 and r = −0.80, for S.aureus and E.coli respectively, p < 0.05). The correlation between IgA levels and bacterial growth rate was weaker but still significant (S.aureus, r = −0.62 and E.coli, r = −0.51, p < 0.05).
Figure 8. Correlation matrix of bacterial growth rates and IgA, LTF and LYZ iBAQ intensities. Each box includes a Pearson correlation coefficient (r value).
Discussion
The current study illustrates that both shorter heating up times and shorter duration of pasteurization at a higher temperature preserved the levels and bioactivity of key DHM bioactive components better than HoP. In addition, the ALP assay showed that all tested time-temperature combinations resulted in the complete inactivation of the enzyme, indicating that sufficient heat load was applied.
Effects of Processing Times and Temperatures on IgA, LTF and LYZ Levels and Activity
Our results confirmed the major impact of HoP on the DHM protein profile (16). HTST, FH and T99, the treatments with the shortest processing times, seemed to better preserve the DHM protein profile when compared to HoP_r, T74 and HoP. Among the methods tested, HTST showed the least reduction in the levels of the studied bioactive proteins, while HoP showed the highest reduction in their levels. More specifically, IgA and LYZ levels after HTST were not significantly different to those of untreated samples, while LTF was significantly reduced after all thermal treatments. These differences were to be expected as DHM was exposed above 55◦C the shortest after HTST (0.33 min) and the longest after HoP (44 min).
Our findings on the advantages of HTST over HoP in the retention of key DHM bioactive components are in line with the literature (14, 16, 19, 31–33). However, considerable variations in the losses of IgA, LTF and LYZ after HTST are reported; 0–60% of IgA and 0–85% of LTF, while for LYZ, losses up to 40% and increases up to 28% were reported (32, 34–37). Several reasons may explain these variations, including the fact that the extent of protein denaturation depends both on their physicochemical characteristics and the nature of the thermal treatment (38). The different HTST devices (e.g., laboratory or industrial heat exchangers, other benchtop devices or by immersion in thermostatically controlled water baths in bulk processes), the differences in holding times (5–25 s) and temperatures (62–87◦C) as well as differences in the methods of analysis in these studies (e.g., ELISAs, radial immunodiffusion assays, enzymatic activity assays, mass spectrometric methods) may have also contributed to different protein retentions documented (15–19, 32, 33, 35, 36, 39).
In respect to the impact of the different heat treatment parameters on the studied DHM bioactive components, among the treatments following the same temperature regime as HTST (15 s at 72◦C), T74 caused the highest protein loss. This could be attributed to the longer heating up time documented during T74 (40 min), which was the result of the small heat exchanging surface area and the small temperature differences between the heating medium and the desired pasteurization temperature. Similarly, HoP_r performed slightly better, although not significantly, than HoP. The possible explanation could be that HoP and HoP_r both follow the same holding regime (30 min at 62.5◦C), but the heating up time is much shorter for HoP_r (0.10 min) than during HoP (27 min). When evaluating the performance of the two treatments with the shortest heating up times, HTST and HoP_r, it was clear that the considerably longer holding time during HoP_r was the reason of the higher protein damage caused. These observations suggest that the combination of processing parameters that leads to a prolonged exposure above temperatures of 55◦C is of crucial importance for the retention of the three studied proteins.
In accordance with our results, Buffin et al. also showed that an optimized HoP, with a mean plateau temperature of 1.5◦C lower and duration of 11 min shorter than HoP, preserved higher amounts of IgA, LTF and LYZ (40). In addition, Escuder-Vieco et al. found a 30% decrease in IgA concentration, regardless of the temperature-time combination used for HTST (5–25 s at 70–75◦C), while Mayayo et al. showed that the largest reductions in the IgA and LTF levels during HoP were documented during the first 5 min of treatment (45 and 70% for IgA and LTF, respectively), with the remaining 25 min of treatment causing <10% reductions (35, 41, 42). For LYZ, studies showed contradictory results in the effect of thermal treatments due to its stable structure, which may be explained by the different analytical approaches used to measure its activity (17, 35, 36, 43, 44).
When the effect of the different heat treatment parameters on the DHM bacteriostatic capacity was assessed, a decrease of the bacteriostatic capacity with increasing exposure times above 55◦C was observed. Furthermore, the correlations observed between IgA, LTF and LYZ levels and the growth rate of S. aureus and E. coli, which are sensitive to these proteins, indicate that these proteins may have a significant role in retarding their growth. Of the treatments performed in the current study, HoP-treated DHM was found to exert the lowest bacteriostatic capacity, while the bacteriostatic capacity of HTST-treated DHM was significantly decreased against E. coli but unaffected against S. aureus. These results are in good agreement with the IgA, LTF and LYZ losses documented after these treatments. Especially for LTF, heat treatments have been shown to reduce its iron-binding capacity (45), which may have contributed to the significant increase in the E. coli growth rate in all heat-treated DHM samples. Other studies investigating the effect of heat treatment on the DHM bacteriostatic capacity, found a similar decrease after HoP (5, 46, 47). The heat-induced denaturation and aggregation during HoP could further explain the loss of protein functionality (13, 48). In contrast to our findings, Silvestre et al. reported a higher decrease in the DHM bacteriostatic capacity after pasteurization at 75◦C for 15 s than after 63◦C for 30 min (49). These differences may be attributed to the different pasteurization designs and the different bacterial strains used. Taken together, treatments with longer processing times, such as HoP, have a significantly larger impact on the DHM bacteriostatic capacity.
Effects of Processing Times and Temperatures on BSSL Level and Activity
Since BSSL is a heat labile enzyme that starts inactivating at temperatures of 45◦C (30), the great loss documented after all thermal treatments was to be expected. Wardell et al. showed that even a short exposure at 55◦C can inactivate the enzyme, which explains the <20% retention that was documented after HTST. BSSL level and activity were almost completely abolished after HoP, as previously reported (10, 19, 50, 51). FH, T99, HoP_r and T74 affected BSSL in a similar manner, independently of the different heating up and holding times applied. When comparing HTST to HoP, the significantly higher BSSL activity detected after HTST, may be attributed to the considerably shorter DHM exposure time over 55◦C (33 s vs. 44 min). Similar observations have been previously documented (19, 51). These findings suggest that non-thermal processing methods, such as ultraviolet-C irradiation or high-pressure processing, may offer substantially better results (50, 52).
Conclusion
Heat treatments, such as HTST, with considerably shorter processing times than the currently recommended HoP, were found to improve the retention of key DHM bioactive components. Our findings suggest that both reduced heating up and holding times are an essential factor for pasteurization optimization, and lead to improved DHM quality. Since the recipients of pasteurized DHM are high-risk infants, these outcomes are of crucial importance. The pasteurization treatments used in this study are all assumed to result in a safe product based on the inactivation of alkaline phosphatase, but whether inactivation of spore-forming pathogens such as B. cereus is achieved with such treatments is yet unknown. Moreover, considering that all thermal treatments caused a major reduction in BSSL levels and activity, in LTF and in bacteriostatic capacity, future studies should additionally investigate the effects of non-thermal methods on these components.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
EK performed the experimental work, analyzed the obtained data, and wrote the manuscript. BS, JvG, RvE, and KH contributed to the conception, the supervision of the study, and critically reviewed the manuscript. SB performed the LC-MS/MS measurements and critically reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of Interest
JvG is the founder and director of the Dutch National Human Milk Bank and member of the Dutch National Health Council. BS is as Science Director of Human Milk Research and Analytical Sciences an employee of Danone Nutricia Research, Utrecht, Netherlands.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to express their gratitude to the Dutch Human Milk Bank for providing DHM samples. We also gratefully acknowledge the support of Dewien Kalpoe in the experimental execution and equipment designing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2021.798609/full#supplementary-material
References
1. Ballard AM. Human milk composition: nutrients and bioactive factors. Pediatr Clin North Am. (2013) 60:49–74. doi: 10.1016/j.pcl.2012.10.002
2. Arslanoglu S, Corpeleijn W, Moro G, Braegger C, Campoy C, Colomb V, et al. Donor human milk for preterm infants : current evidence and research directions. J Pediatr Gastroenterol Nutr. (2013) 57:535–42. doi: 10.1097/MPG.0b013e3182a3af0a
3. Monaco MH, Kim J, Donovan SM. Human Milk: Composition and Nutritional Value. 1st ed. Encyclopedia of Food and Health. Elsevier Ltd (2015). p. 357–62.
4. Lönnerdal B. Nutritional and physiologic significance of human milk proteins. Am J Clin Nutr. (2003) 77:1537S–43S. doi: 10.1093/ajcn/77.6.1537S
5. Christen L, Lai CT, Hartmann B, Hartmann PE, Geddes DT. The effect of UV-C pasteurization on bacteriostatic properties and immunological proteins of donor human milk. PLoS ONE. (2013) 8:1–9. doi: 10.1371/journal.pone.0085867
6. Moro GE, Arslanoglu S. Heat treatment of human milk. J Pediatr Gastroenterol Nutr. (2012) 54:165–6. doi: 10.1097/MPG.0b013e318235d629
7. World Health Organization. Guidelines on Optimal Feeding of Low Birth-Weight Infants In Low-And Middle-Income Countries. Geneva: WHO. (2011). p. 16–45.
8. Moro GE, Billeaud C, Rachel B, Calvo J, Cavallarin L, Christen L, et al. Processing of donor human milk: update and recommendations from the European milk bank association (EMBA). Front Pediatr. (2019) 7:1–10. doi: 10.3389/fped.2019.00049
9. Arslanoglu S, Bertino E, Tonetto P, De Nisi G, Ambruzzi AM, Biasini A, et al. Guidelines for the establishment and operation of a donor human milk bank. J Matern Neonatal Med. (2010) 23(Suppl. 2):1–20. doi: 10.3109/14767058.2010.512414
10. Peila C, Moro GE, Bertino E, Cavallarin L, Giribaldi M, Giuliani F, et al. The effect of holder pasteurization on nutrients and biologically-active components in donor human milk: a review. Nutrients. (2016) 8:1–19. doi: 10.3390/nu8080477
11. Andersson Y, Sävman K, Bläckberg L, Hernell O. Pasteurization of mother's own milk reduces fat absorption and growth in preterm infants. Acta Paediatr Int J Paediatr. (2007) 96:1445–9. doi: 10.1111/j.1651-2227.2007.00450.x
12. Kontopodi E, Arslanoglu S, Bernatowicz-Lojko U, Bertino E, Bettinelli ME, Buffin R, et al. Donor milk banking: Improving the future. A survey on the operation of the European donor human milk banks. PLoS ONE. (2021) 16:e0256435. doi: 10.1371/journal.pone.0256435
13. Czank C, Prime DK, Hartmann BEN, Simmer K, Hartmann PE, Biochemistry D, et al. Retention of the immunological proteins of pasteurized human milk in relation to pasteurizer design and practice. Pediatr Res. (2009) 66:374–9. doi: 10.1203/PDR.0b013e3181b4554a
14. Escuder-Vieco D, Espinosa-Martos I, Rodríguez JM, Fernández L, Pallás-Alonso CR. Effect of HTST and holder pasteurization on the concentration of immunoglobulins, growth factors, and hormones in donor human milk. Front Immunol. (2018) 9:2222. doi: 10.3389/fimmu.2018.02222
15. Dhar J, Fichtali J, Skura BJ, Nakai S, Davidson AGF. Pasteurization efficiency of a HTST system for human milk. J Food Sci. (1996) 61:569–73. doi: 10.1111/j.1365-2621.1996.tb13160.x
16. Baro C, Giribaldi M, Arslanoglu S, Giuffrida MG, Dellavalle G, Conti A, et al. Effect of two pasteurization methods on the protein content of human milk. Front Biosci Elit. (2011) 3:818–29. doi: 10.2741/289
17. Klotz D, Schreiner M, Falcone V, Jonas D, Kunze M, Weber A, et al. High-temperature short-time treatment of human milk for bacterial count reduction. Front Pediatr. (2018) 6:1–8. doi: 10.3389/fped.2018.00359
18. Terpstra FG, Rechtman DJ, Lee ML, Van Hoeij K, Berg H, Van Engelenberg FAC, et al. Antimicrobial and antiviral effect of high-temperature short-time (HTST) pasteurization apllied to human milk. Breastfeed Med. (2007) 2:27–33. doi: 10.1089/bfm.2006.0015
19. Giribaldi M, Coscia A, Peila C, Antoniazzi S, Lamberti C, Ortoffi M, et al. Pasteurization of human milk by a benchtop high-temperature short-time device. Innov Food Sci Emerg Technol. (2016) 36:228–33. doi: 10.1016/j.ifset.2016.07.004
20. Peila C, Emmerik NE, Giribaldi M, Stahl B, Ruitenberg JE, Van Elburg RM, et al. Human milk processing: a systematic review of innovative techniques to ensure the safety and quality of donor milk. J Pediatr Gastroenterol Nutr. (2017) 64:353–61. doi: 10.1097/MPG.0000000000001435
21. de Waard M, Mank E, van Dijk K, Schoonderwoerd A, van Goudoever JB. Holder-pasteurized human donor milk : how long can it be preserved? J Pediatr Gastroenterol Nutr. (2018) 66:10–2. doi: 10.1097/MPG.0000000000001782
22. Daniels B, Reimers P, King T, Schmidt S, Coutsoudis A. Effect of the PiAstra benchtop flash-heating pasteurizer on immune factors of donor human milk. Breastfeed Med. (2018) 13:281–5. doi: 10.1089/bfm.2018.0017
23. Chantry CJ, Israel-Ballard K, Moldoveanu Z, Peerson J, Coutsoudis A, Sibeko L, et al. Effect of flash-heat treatment on immunoglobulins in breast milk. J Acquir Immune Defic Syndr. (2009) 51:264–7. doi: 10.1097/QAI.0b013e3181aa12f2
24. Law AJR, Leaver J. Effect of pH on the thermal denaturation of whey proteins in milk. J Agric Food Chem. (2000) 48:672–9. doi: 10.1021/jf981302b
25. Liu Y, Xiong L, Kontopodi E, Boeren S, Zhang L, Zhou P, et al. Changes in the milk serum proteome after thermal and non-thermal treatment. Innov Food Sci Emerg Technol. (2020) 66:102544. doi: 10.1016/j.ifset.2020.102544
26. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. (2008) 26:1367–72. doi: 10.1038/nbt.1511
27. Krewinkel M, Baur C, Kranz B, von Neubeck M, Wenning M, Scherer S, et al. A sensitive and robust method for direct determination of lipolytic activity in natural milk environment. Food Anal Methods. (2016) 9:646–55. doi: 10.1007/s12161-015-0233-4
28. Hettinga K, Reina FM, Boeren S, Zhang L, Koppelman GH, Postma DS, et al. Difference in the breast milk proteome between allergic and non-allergic mothers. PLoS ONE. (2015) 10:1–11. doi: 10.1371/journal.pone.0122234
29. R Core Team,. A Language Environment for Statistical Computing; R Foundation for Statistical Computing. Vienna (2018). Available online at: https://www.R-project.org/
30. Wardell JM, Wright AJ, Bardsley WG, D'Souza SW. Bile salt–stimulated lipase and esterase activity in human milk after collection, storage, and heating: nutritional implications. Pediatr Res. (1984) 18:382–6. doi: 10.1203/00006450-198404000-00017
31. Moro GE, Stellwagen LM, Thome UH, Lembo D, Donalisio M, Rittà M, et al. High temperature-short time pasteurization has a lower impact on the antiviral properties of human milk than holder pasteurization. Front Pediatr. (2018) 6:304. doi: 10.3389/fped.2018.00304
32. Goldblum RM, Dill CW, Albrecht TB, Alford ES, Garza C, Goldman AS. Rapid high-temperature treatment of human milk. J Pediatr. (1984) 104:380–5. doi: 10.1016/S0022-3476(84)81099-9
33. Hamprecht K, Maschmann J, Müller D, Dietz K, Besenthal I, Goelz R, et al. Cytomegalovirus (CMV) inactivation in breast milk: reassessment of pasteurization and freeze-thawing. Pediatr Res. (2004) 56:529–35. doi: 10.1203/01.PDR.0000139483.35087.BE
34. Billeaud C. High hydrostatic pressure treatment ensures the microbiological safety of human milk including bacillus cereus and preservation of bioactive proteins including lipase and immuno-proteins: a narrative review. Foods. (2021) 10:1327. doi: 10.3390/foods10061327
35. Mayayo C, Ramos SJ, Martínez MJ, Miguel L, Lourdes C, María S. Effect of high pressure and heat treatments on IgA immunoreactivity and lysozyme activity in human milk. Eur Food Res Technol. (2016) 242:891–8. doi: 10.1007/s00217-015-2595-7
36. Picaud J-C, Buffin R. Human milk-treatment and quality of banked human milk. Clin Perinatol. (2017) 44:95–119. doi: 10.1016/j.clp.2016.11.003
37. Wesolowska A, Sinkiewicz-Darol E, Barbarska O, Bernatowicz-Lojko U, Borszewska-Kornacka MK, van Goudoever JB. Innovative techniques of processing human milk to preserve key components. Nutrients. (2019) 11:1–17. doi: 10.3390/nu11051169
38. Bogahawaththa D, Chandrapala J, Vasiljevic T. Thermal denaturation of bovine immunoglobulin G and its association with other whey proteins. Food Hydrocolloids. (2017) 72:350–7. doi: 10.1016/j.foodhyd.2017.06.017
39. Goelz R, Hihn E, Hamprecht K, Dietz K, Jahn G, Poets C, et al. Effects of different CMV-heat-inactivation-methods on growth factors in human breast milk. Pediatr Res. (2009) 65:458–61. doi: 10.1203/PDR.0b013e3181991f18
40. Buffin R, Hays S, Drai J, Sarda MN, Picaud JC. Better control of holder pasteurization results in higher retention of human milk lactoferrin, IgA, and lysozyme. Front Pediatr. (2018) 6:381. doi: 10.3389/fped.2018.00381
41. Mayayo C, Montserrat M, Ramos SJ, Martínez-Lorenzo MJ, Calvo M, Sánchez L, et al. Kinetic parameters for high-pressure-induced denaturation of lactoferrin in human milk. Int Dairy J. (2014) 39:246–52. doi: 10.1016/j.idairyj.2014.07.001
42. Escuder-Vieco D, Espinosa-Martos I, Rodríguez JM, Corzo N, Montilla A, Siegfried P, et al. High-temperature short-time pasteurization system for donor milk in a human milk bank setting. Front Microbiol. (2018) 9:926. doi: 10.3389/fmicb.2018.00926
43. Guerra AF, Mellinger-Silva C, Rosenthal A, Luchese RH. Hot topic: holder pasteurization of human milk affects some bioactive proteins. J Dairy Sci. (2018) 101:2814–8. doi: 10.3168/jds.2017-13789
44. Sousa SG, Santos MD, Fidalgo LG, Delgadillo I, Saraiva JA. Effect of thermal pasteurisation and high-pressure processing on immunoglobulin content and lysozyme and lactoperoxidase activity in human colostrum. Food Chem. (2014) 151:79–85. doi: 10.1016/j.foodchem.2013.11.024
45. Mata L, Sánchez L, Headon DR, Calvo M. Thermal denaturation of human lactoferrin and its effect on the ability to bind iron. J Agric Food Chem. (1998) 46:3964–70. doi: 10.1021/jf980266d
46. Van Gysel M, Cossey V, Fieuws S, Schuermans A. Impact of pasteurization on the antibacterial properties of human milk. Eur J Pediatr. (2012) 171:1231–7. doi: 10.1007/s00431-012-1750-4
47. Akinbi H, Meinzen-Derr J, Auer C, Ma Y, Pullum D, Kusano R, et al. Alterations in the host defense properties of human milk following prolonged storage or pasteurization. J Pediatr Gastroenterol Nutr. (2010) 51:347–52. doi: 10.1097/MPG.0b013e3181e07f0a
48. Xiong L, Li C, Boeren S, Vervoort J, Hettinga K. Effect of heat treatment on bacteriostatic activity and protein profile of bovine whey proteins. Food Res Int J. (2020) 127:108688. doi: 10.1016/j.foodres.2019.108688
49. Silvestre D, Ruiz P, Martínez-Costa C, Plaza A, López MC, Silvestre D, et al. Effect of pasteurization on the bactericidal capacity of human milk. J Hum Lact. (2008) 24:371–6. doi: 10.1177/0890334408319158
50. Pitino MA, Unger S, Doyen A, Pouliot Y, Aufreiter S, Stone D, et al. High hydrostatic pressure processing better preserves the nutrient and bioactive compound composition of human donor milk. J Nutr. (2019) 149:497–504. doi: 10.1093/jn/nxy302
51. Escuder-Vieco D, Rodríguez JM, Espinosa-Martos I, Corzo N, Montilla A, García-Serrano A, et al. High-temperature short-time and holder pasteurization of donor milk: impact on milk composition. Life. (2021) 11:1–16. doi: 10.3390/life11020114
Keywords: protein functionality, donor human milk, milk processing, holder pasteurization, HTST
Citation: Kontopodi E, Boeren S, Stahl B, van Goudoever JB, van Elburg RM and Hettinga K (2022) High-Temperature Short-Time Preserves Human Milk's Bioactive Proteins and Their Function Better Than Pasteurization Techniques With Long Processing Times. Front. Pediatr. 9:798609. doi: 10.3389/fped.2021.798609
Received: 20 October 2021; Accepted: 29 December 2021;
Published: 20 January 2022.
Edited by:
Ardythe Morrow, University of Cincinnati, United StatesReviewed by:
Javier Calvo, Fundació Banc Sang i Teixits de les Illes Balears, SpainMaryanne Perrin, University of North Carolina at Greensboro, United States
Copyright © 2022 Kontopodi, Boeren, Stahl, van Goudoever, van Elburg and Hettinga. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kasper Hettinga, a2FzcGVyLmhldHRpbmdhJiN4MDAwNDA7d3VyLm5s