Development of Forensically Important Sarcophaga peregrina (Diptera: Sarcophagidae) and Intra-Puparial Age Estimation Utilizing Multiple Methods at Constant and Fluctuating Temperatures

Simple Summary Constant and fluctuating temperatures can affect the growth and development of insects, which will affect the accuracy of minimum postmortem interval estimation. Our study found that the development time of S. peregrina was prolonged under fluctuating temperatures. The estimation of the pupal stage has a significant effect on the estimation of PMI. We explored the feasibility of different methods for intra-puparial age estimation and found that the DEG expression profile, ATR-FTIR, and CHCs detection, combined with chemometrics, had the potential to estimate the intra-puparial age of S. peregrina under constant and fluctuating temperatures. Abstract Sarcophaga peregrina (Robineau-Desvoidy, 1830) has the potential to estimate the minimum postmortem interval (PMImin). Development data and intra-puparial age estimation are significant for PMImin estimation. Previous research has focused on constant temperatures, although fluctuating temperatures are a more real scenario at a crime scene. The current study examined the growth patterns of S. peregrina under constant (25.75 °C) and fluctuating temperatures (18–36 °C; 22–30 °C). Furthermore, differentially expressed genes, attenuated total reflectance Fourier-transform infrared spectroscopy, and cuticular hydrocarbons of S. peregrina during the intra-puparial period were used to estimate age. The results indicated that S. peregrina at fluctuating temperatures took longer to develop and had a lower pupariation rate, eclosion rate, and pupal weight than the group at constant temperatures did. Moreover, we found that six DEG expression profiles and ATR-FTIR technology, CHCs detection methods, and chemometrics can potentially estimate the intra-puparial age of S. peregrina at both constant and fluctuating temperatures. The findings of the study support the use of S. peregrina for PMImin estimation and encourage the use of entomological evidence in forensic practice.


Introduction
Entomology evidence could be the only evidence that can be used to estimate the minimum postmortem interval (PMI min ) in decaying corpses. After insects are collected from the corpse and identified, the PMI min is estimated using development data of necrophagous flies [1,2]. Therefore, development data collection and accurate age estimation of necrophagous flies are meaningful for PMI min estimation.
The intra-puparial period, which occupies half of the immature stage [19], is crucial for estimating PMI min . Many methods for estimating pupal age have been reported, including puparium color variation [20,21], intra-puparial morphological characteristics [22], x-irradiation [23], micro-computed tomography imaging [24], and optical coherence tomography [25]. However, morphological methods require morphological knowledge, and the application of development data in the intra-puparial period is limited. Micro-CT techniques, for example, are uncommon at forensic workstations and inconvenient to use at crime scenes. Therefore, it is crucial to find a fast, accurate, and convenient method for estimating the intra-puparial period age.
Studying differentially expressed genes (DEGs) as a method for age estimation has been used for blow flies and flesh flies [26,27]. Fourier-transform infrared (FTIR) spectroscopy can determine the metabolic properties of insects based on absorption characteristics, such as lipids, proteins, and cellular processes [28]. Attenuated total reflection Fourier-transform infrared (ATR-FTIR) has been used to evaluate the effects of the cause of death on PMI estimation [29], to discriminate between normal and pathological populations [30][31][32] and to detect SARS-CoV2 infections on the spot [33]. However, there are a limited number of reports of FTIR applications in forensic entomology. Cuticular hydrocarbons (CHCs), the major components of the waxy layer of the insect cuticle [34,35], which has been proved to be a reliable method of predicting the age of the fly [36][37][38][39]. As a result, further research into the value of DEGs, ATR-FTIR, and CHCs in intra-puparial age estimation is required.
In the current study, developmental data of S. peregrina were collected at constant and fluctuating temperatures for PMI min estimation. In addition, six DEG expression profiles, the ATR-FTIR spectroscopy characters, and the CHC profiles of S. peregrina were investigated and analyzed during the intra-puparial period. This study provides a basis for rapid and accurate PMI min estimation.

Insect Rearing and Experimental Temperature Settings
Sarcophaga peregrina was captured using pig carcasses in Changsha (28 • 12 N, 112 • 58 E), Hunan Province, China, in July 2017 and reared in Guo's Lab (Changsha China). Traditional morphological characteristics [56] the long cytochrome oxidase subunit I (CO I) sequences were used to identify the species. Adults were fed in an artificial climate chamber (LRH-250-GSI, Taihong Co., Ltd., Shaoguan, China) with 25.0 ± 1.0 • C, 75% RH, a 12 h light and 12 h dark cycle, and the feeding conditions were fresh water and mixed milk powder and sugar.
The experimental constant and fluctuating temperatures ( Figure 1) were designed based on the temperature variations in Changsha, China, from June to October 2011 to 2021 ( Figure S1). The temperature groups were recorded as fluctuating model A (group A: 18-36 • C), fluctuating model B (group B: 22-30 • C), and constant model C (group C: 25.75 • C). The temperature tolerance of the artificial climate chamber was 1 • C with 75% RH and a 12 h: 12 h (L: D) cycle, as shown in Figure 1 and Table S1, respectively. the development of S. peregrina through commonly used literature search platforms (Pub Med and Web of Science).
In the current study, developmental data of S. peregrina were collected at constan and fluctuating temperatures for PMImin estimation. In addition, six DEG expression pro files, the ATR-FTIR spectroscopy characters, and the CHC profiles of S. peregrina were in vestigated and analyzed during the intra-puparial period. This study provides a basis for rapid and accurate PMImin estimation.

Insect Rearing and Experimental Temperature Settings
Sarcophaga peregrina was captured using pig carcasses in Changsha (28°12′ N, 112°58 E), Hunan Province, China, in July 2017 and reared in Guo's Lab (Changsha China). Tra ditional morphological characteristics [56] the long cytochrome oxidase subunit I (CO I sequences were used to identify the species. Adults were fed in an artificial climate cham ber (LRH-250-GSI, Taihong Co., Ltd., Shaoguan, China) with 25.0 ± 1.0 °C, 75% RH, a 12 h light and 12 h dark cycle, and the feeding conditions were fresh water and mixed milk powder and sugar.
The experimental constant and fluctuating temperatures ( Figure 1) were designed based on the temperature variations in Changsha, China, from June to October 2011 to 2021 ( Figure S1). The temperature groups were recorded as fluctuating model A (group A: 18-36 °C), fluctuating model B (group B: 22-30 °C), and constant model C (group C 25.75 °C). The temperature tolerance of the artificial climate chamber was 1 °C with 75% RH and a 12 h: 12 h (L: D) cycle, as shown in Figure 1 and Table S1, respectively.

Collection of Samples for Developmental Analysis and Multi-Method Analysis
A pig lung was used to induce the larvae to hatch. Larvae (ca. 1500) produced within 1 h were divided into three groups: A, B, and C. Ten larval samples were randomly taken every eight hours. Larval instar was identified using a Zeiss 2000-C stereomicroscope (Carl Zeiss, Jena, Germany). Larval length was measured using a digital caliper (Mineette, Shanghai, China). The developmental time, pupation rate, eclosion rate, and pupal weight were recorded during the experiment. Data were analyzed using Origin Pro 8.6 software (OriginLab, Northampton, MA, USA, SCR: 015636). Three replications were performed.
Larvae produced within 1 h were cultured at constant temperature (group C) and fluctuating temperature (group A). Thirty pupae were harvested at 24 h intervals until the emergence of adults. The samples collected each time were divided into three groups on average for DEGs, ATR-FTIR, and CHCs analysis, respectively. Three replications were performed. A total of 1620 pupae were collected: 810 in group A and 810 in group C.

DEGs Study
Total RNA from each intra-puparial tissue obtained by removing the puparium was extracted using the SteadyPure Quick RNA Extraction Kit (Code No. AG21017) (Accurate Biotechnology (Hunan) Co., Ltd. Changsha, China) and stored at −80 • C following the manufacturer's instructions. Total RNA quantification was performed using NanoDrop 2000 (Thermo Fisher Scientific, Wilmington, DE, USA). The ratios of optical density (OD) of 260/280 and 260/230 were selected as the evaluation indicators. cDNA was synthesized with the Evo M-MLV RT Mix kit (Code No. AG11728) (Accurate Biotechnology (Hunan) Co., Ltd. Changsha, China), and then qPCR was performed.
The candidate DEGs were identified according to the genome [51] (NCBI no. JABZEU000000000) and transcriptome data [57] (NCBI nos. PRJNA795032 and PR-JNA810308) of S. peregrina. Quantitative real-time PCR operation was performed according to the instructions of the SYBR ® Green Premix Pro Taq HS qPCR Kit (ROX Plus) (Code No. AG11718) (Accurate Biotechnology (Changsha, China) Co., Ltd.) on an ABI 7500 Real-Time PCR system (Applied Biosystem, Carlsbad, CA, USA). Primers for the DEGs and reference genes were designed using Primer Premier 5.0 (Biosoft Premier, Palo Alto, CA, USA) and are listed in Table S2. The 2 − Ct method was used to calculate the relative quantification (RQ) of the DEGs. All statistical analysis was conducted using the software Origin Pro 8.6.

ATR-FTIR Study
The intra-puparial tissue was ground using liquid nitrogen, and later, transferred into a 1.5 mL EP tube and shaken using a vortex mixer (Guangzhou GenXion Biotechnology Co. LTD, Guangzhou, China). An 850 Fourier-transform infrared spectrometer equipped with an ATR attachment was used (Tianjin Gangdong Technology Development Co. LTD, Tianjin, China) to collect spectra. The ATR attachment was cleaned with absolute ethanol before each measurement, and background spectra were acquired before each operation. Spectra were recorded at a resolution of 4 cm −1 in the frequency range from 4000 to 900 cm −1 , with 32 scans per sample. OMNIC, version 9.2 (Thermo Fisher Scientific, Wilmington, MA, USA), was used to process the spectral data. Unnecessary background information was reduced by spectral preprocessing. The standard normal variate (SNV) was used to decrease the impact of light scattering and spectral size due to sample thickness [58]. Furthermore, all spectral data were smoothed and denoised (Savitzky-Golay convolution algorithm; number of smoothing points: 15). Finally, the spectral range of 1800-900 cm −1 was used for further analysis.
In this study, Principal Component Analysis (PCA), Orthogonal Partial Least Squares Discriminant Analysis (OPLS-DA), and Partial Least Squares Discriminant Analysis (PLS) were performed using SIMCA 14.1 software to observe the spectral characters of the intrapuparial tissue at different times of sampling. PCA can reduce the dimensions of the data and retain original data features. OPLS-DA is a supervised statistical method of discriminant analysis [59,60]. For OPLS-DA, values of R 2 (cum) and Q 2 (cum) close to 1.0 indicate an excellent model in theory, >0.5 indicates a better model, and >0.4 indicates an acceptable model. The permutation test was used to examine if the model was overfitted. PLS regression analysis can decompose X variable and extract latent variable under the guidance of Y variable [61,62]. R 2 and the root mean square error (RMSE) were used to evaluate the PLS model; values close to 1 for R2 and 0 for RMSE indicate a better model. The variable importance in projection (VIP) was used to assess the contribution of each variable to the OPLA-DA and PLS models, and a VIP value of >1.0 was considered to be significant.

CHCs Study
Referring to Zhang et al.'s [49,57] operation, the intra-puparial tissue was immersed in a 2 mL glass vial containing 1 mL of redistilled hexane for 30 min at room temperature. An alkane mixture of heptane and tetracontane (C7-C40, 1 µg/mL, O2SI) dissolved in 1 mL of hexane was selected as an external standard. Using a syringe with a nylon membrane and filter, the soaking solution was transferred into a new 2 mL glass vial. Following that, the samples were vacuum concentrated and dried before being dissolved in 200 µL redistilled hexane for GC-MS analysis. GC-MS (Agilent Technologies, 7890B-5977A GC/MSD) fitted with a DB-5MS capillary column was used to analyze the CHC. Operations of sample loading were performed and instrument condition settings were used according to references [49] that have already been reported. Due to the length limitation, specific setting contents are not described here.
MSD ChemStation Data Analysis F.01.03 was used to integrated the peak areas, and the qualitative identification of CHC was carried out by comparing the chromatographic retention time, characteristic fragment ions, and abundance ratio with those of standard values. The Kováts retention index, NIST 11 library, and literature [49,[63][64][65][66][67][68][69][70] were selected to assist identification, which we called the CHCs database. PCA, OPLS-DA, and PLS models were used to visualize the filtered data.

Development Analysis
The correlation between larval body length and development time under different temperatures is shown in Figure 2. The model of fourth-order polynomials was used to establish nonlinear regression. Table 1 shows the equations and R 2 values of the model. The results show that the larval body length of group C was longer than those of both group A and group B before 80-90 h, and the longest larvae can be observed in group C.
S. peregrina could complete their life history under constant and fluctuating temperatures, and the development duration is listed in Table 2. Group C spent the minimum amount of time reaching full development. The development rate of group A was similar to that of group B and slower than that of group C. As shown in Figure 3, the pupariation rate, eclosion rate, and pupal weight were lower and statistically significant in groups A and B than they were in group C. These results indicate that fluctuating temperatures can affect the development of S. peregrina.

DEGs Analysis
The expression profiles of six DEGs (circRNA_2143, circRNA_3489, circRNA_2847, fln, UQCRFS1, and COX5A) changed with the increase in intra-puparial age ( Figure 4). Polynomial regression analysis described the relationship between the level of DEG expression of S. peregrina and developmental age (days) ( Table 3). The increasing expression tendencies of the six DEGs, which are similar under constant and fluctuating temperatures, indicate that these six DEGs can be gene marker candidates to estimate the intra-puparial age of S. peregrina at both constant and fluctuating temperatures.   Table 2. Group C spent the minimum amount of time reaching full development. The development rate of group A was similar to that of group B and slower than that of group C. As shown in Figure 3, the pupariation rate, eclosion rate, and pupal weight were lower and statistically significant in groups A and B than they were in group C. These results indicate that fluctuating temperatures can affect the development of S. peregrina.

DEGs Analysis
The expression profiles of six DEGs (circRNA_2143, circRNA_3489, circRNA_2 fln, UQCRFS1, and COX5A) changed with the increase in intra-puparial age ( Figure   Figure 3. The pupariation rate, eclosion rate, and pupal weight were lower and statistically significant in groups A and B than they were in group C. (a) Pupariation rate. (b) Eclosion rate. (c) Pupal weight (* represents p < 0.05, ** represents p < 0.01).
Polynomial regression analysis described the relationship between the level of DE pression of S. peregrina and developmental age (days) ( Table 3). The increasing expr tendencies of the six DEGs, which are similar under constant and fluctuating tem tures, indicate that these six DEGs can be gene marker candidates to estimate the puparial age of S. peregrina at both constant and fluctuating temperatures.     Figure 5 shows the average FTIR spectra of different intra-puparial development ages (days) in the range of 1800-900 cm −1 . Based on previous studies [71][72][73], the main infrared absorption band identification is shown in Table 4. The OPLS-DA classification model ( Figure 6) showed that the same day's spectra were well clustered and distinguished from those of the other days for both groups, except for a few data overlaps between day 5 and day 9 for group C. Common wavenumber bands with VIP > 1 (Table S3) Table 5: R 2 > 0.9 and RMSE of 0.4271 day and 0.7683 day for groups A and C, respectively. Figure 6 shows the regression model to predict the intra-puparial age. Wavenumbers with VIP > 1 are shown in Table S3. The common wavenumber bands were These results suggest that the spectral differences observed in this study were caused by changes in proteins, lipids, and other intra-puparial tissue components. The error of intra-puparial age inference is smaller in a fluctuating temperature environment.

CHC Analysis
A total of 41 hydrocarbons were identified, including 21 n-alkanes, 17 branched alkanes, and 3 alkenes with a carbon chain length between C9 and C32 (Table 6). According to the mini bar chart, branched alkanes were the most abundant compound from day 1 to day 8 at constant or fluctuating temperatures; on day 9, it was n-alkanes. The CHCs profiles (Table S4) show that the most abundant compounds in group A were 4-Methyl C17 Figure 6. Chemometrics was used to analyze the data that the ATR−FTIR study generated. (a,b) OPLS-DA analysis of group A and group C, respectively. (c,d) PLS analysis of group A and group C, respectively. Table 5. Model evaluation parameters and regression equations.

CHC Analysis
A total of 41 hydrocarbons were identified, including 21 n-alkanes, 17 branched alkanes, and 3 alkenes with a carbon chain length between C9 and C32 (Table 6). According to the mini bar chart, branched alkanes were the most abundant compound from day 1 to day 8 at constant or fluctuating temperatures; on day 9, it was n-alkanes. The CHCs profiles (Table S4) show that the most abundant compounds in group A were 4-Methyl C17 on days 1-3 and DiMethyl C29 on days 6-9. Additionally, a transitional stage could be observed on days 4-5. Except for a slightly different transitional stage on days 4-5, group C showed nearly the same change as group A did. In addition, dominant CHCs in the intra-puparial tissues showed a shift from a low to a high molecular weight with the age of the pupae. The OPLS-DA plot (Figure 7) shows that the clustering effect of group A was better than that of group C. Common variables (Table S6) of groups A and C with VIP > 1 (Table S5) were selected. 4-Methyl C10, Methyl C29, C19, C21, and 7-Methyl C27 were identified. The parameters of PLS analysis are listed in Table 5: R 2 > 0.9, and RMSEs of 0.4837 day and 0.7653 day for groups A and C, respectively. C27, C25, C23, and DiMethyl C29 were identified. The CHCs identified above suggest a possible link between changes in the most abundant CHCs and increasing intra-puparial age, particularly DiMethyl C29, which could be useful in estimating S. peregrina intra-puparial age at constant and fluctuating temperatures. Additionally, the error of intra-puparial age inference is smaller in a fluctuating temperature environment.

Discussion
Fluctuating temperatures can influence the development of flies. In studies, temperature fluctuations have been shown to either accelerate or slow the development of insects [13,14] or have no significant difference in development [74]. This study found that fluctuating temperatures influenced the feeding stage of S. peregrina. It led to smaller lengths and weights, longer developmental times, and a lower survival rate for larvae and pupae. The thermal mean and its proximity to the developmental threshold [17,18] could be the

Discussion
Fluctuating temperatures can influence the development of flies. In studies, temperature fluctuations have been shown to either accelerate or slow the development of insects [13,14] or have no significant difference in development [74]. This study found that fluctuating temperatures influenced the feeding stage of S. peregrina. It led to smaller lengths and weights, longer developmental times, and a lower survival rate for larvae and pupae. The thermal mean and its proximity to the developmental threshold [17,18] could be the main reasons for this phenomenon. The developmental rate difference between fluctuating and constant temperatures grows as the temperature fluctuation range increases [4]. When the range is within the suitable range for insect growth, the effect of constant and fluctuating temperatures will not be significant. However, when the range exceeds the suitable range, such as a lower or higher temperature, the effects of fluctuating temperature are relatively large [75]. Changes in the development time at fluctuating temperatures can lead to incorrect age estimations, which in turn leads to incorrect PMI min estimations.
More development data under fluctuating temperatures are required to accurately determine the intra-puparial age of necrophagous flies. However, this can greatly increase forensic entomologists' workload. Therefore, it is crucial to develop an approach that is faster, simpler, more accurate, and universally applicable.
The value of DEGs for inferring the age of flies has been confirmed [27,57,76]. Studies on Sarcophaga dux [6], Aldrichina graham [77], Lucilia illustris [78], Glossina morsitans [79], and Calliphora vicina [1,80,81] showed that DEGs can be used to determine the age of necrophagous flies. However, the expression trends of the DEGs in the above research were not singular, and most of them were performed under constant temperatures. DEGs that can stably express at constant and fluctuating temperatures and show a single temporal expression trend are expected to be genetic markers for age estimation. For example, a study on the expression of bcd, sll, and cs shows a significant linear trend during the development of blowfly eggs [27]. Candidate DEGs for S. peregrina were chosen from previously published results in our previous study [26]. Their expression trends were not singular, and they were demonstrated under constant temperatures. The candidate DEGs in this study were chosen from S. peregrina's reference genome and transcriptome data [49,51]. We found that the six DEGs showed a single up-regulated temporal expression, with the intra-puparial period increasing under constant and fluctuating temperatures. This result suggested that the six DEGs could become insect gene markers for estimating the intra-puparial age of S. peregrina. Notably, three circular RNAs (circRNAs), which are more stable than linear RNAs are [82,83], were introduced as candidate DEGs for the first time in this paper. Although DEGs can be used to estimate the PMI min , the number of studies on the profiles of necrophagous flies is still limited. More research should be conducted.
The amount of research on the application of vibrational spectroscopy in biochemistry and forensic science is increasing [32,84,85]. Based on vibrational bands, the biochemical information associated with the fingerprint region (1800-900 cm −1 ) depicts the structure and function of cell samples [86][87][88][89]. With the enormous quantity of data produced by IR spectroscopy, chemometrics can help with data pre-processing, sample classification, and reducing the dimensionality of chemical information [84]. Perez-Mendoza [90] and Pickering et al. [91] successfully identified the age of different species by using infrared spectroscopy. Our study found that the proteins, lipids, and other components of the intrapuparial tissues, represented by the common bands, were affected by the development of flies. These substances are important for the model's establishment and age estimation. The model parameters suggested that better results were obtained at fluctuating temperatures, demonstrating that temperature settings that match the real environment are crucial for accurate age estimation. Moreover, the automation and convenience of the test instruments make immediate analysis at crime scenes possible [86]. The use of spectroscopic methods in forensic science could improve the use of entomological evidence. Cuticular hydrocarbons are the main component of the waxy layer of insect cuticles [35] and play a key role in insect survival. They can act as a barrier for insects to prevent water loss, avoid mechanical damage, isolate microorganisms and chemicals, and they also act as pheromones for insects [35,[92][93][94]. Sarcophaga bullata CHC profiles during development have been published [95]. Studies on Chrysomya albiceps [38], Chrysomya rufifacies [37], and Sarcophaga peregrina [49] also confirmed that the dynamics of CHC profiles during the development of flies can be used to determine the age of necrophagous flies [96]. The carbon chain length of the CHCs identified in this study was between C9 and C32, which was within the range of C8-C36 discovered by Zhang et al. [49]. However, the most abundant CHC differed from that in Zhang et al.'s study [49]. This is because the tissue used in the study, intra-puparial tissue, differs from the pupal shell used in their study. Meanwhile, as the intra-puparial age increased, the dominant CHCs shifted from low-to high-molecular-weight hydrocarbons. This is similar to the study of Chrysomya rufifacies [37]. The next stage of S. peregrina after completing metamorphic development is the adult stage. The heavier, long-chain hydrocarbons, which are the most important compounds in the insect cuticle [48,97], are essential for waterproofing [98][99][100] and help adults to survive in complex external environments. Additionally, the model parameters showed that better results could be obtained from group A, which was placed under fluctuating temperatures. These results indicate that more accurate data can be obtained for precise age estimation in a more realistic simulated situation. Additionally, CHC's profiles can be used to determine the intra-puparial age of S. peregrina, as well as the PMI min .

Conclusions
This study investigated methods for estimating the intra-puparial age of Sarcophaga peregrina under constant and fluctuating temperatures, such as DEGs, ATR-FTIR, and CHC profiles. We found that the development of S. peregrina was different under constant and fluctuating temperatures and demonstrated that the three methods are meaningful for pupal age estimation.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ani13101607/s1, Figure S1: The temperature variation from June to October 2011 to 2021 in Changsha city, China; Figure S2: The visualization results of PCA analysis and permutation test of the data that the ATR-FTIR study generated. a,b: PCA analysis of group A and group C, respectively. c, d: Permutation test result of group A and group C, respectively; Figure  S3: The visualization results of PCA analysis and permutation test of the data that the CHCs study generated. a,b: PCA analysis of group A and group C, respectively. c,d: Permutation test result of group A and group C, respectively; Table S1: The temperature setting of development of Sarcophaga peregrina at constant and fluctuating temperatures; Table S2: Primers used in RT-qPCR analysis; Table  S3: The ATR-FTIR spectral with VIP>1 of PLS and OPLS-DA model; Table S4: CHCs profile of S. peregrina intra-puparial tissue at constant and fluctuating temperatures; Table S5: The CHCs data with VIP > 1 of PLS and OPLS-DA model; Table S6: Common VIP recognition.
Author Contributions: All authors contributed to the study's conception and design. Material preparation and data analysis were performed by Y.S.; The first draft of the manuscript was written by F.Y.; F.J.N., X.Z. and Y.F. commented on previous versions of the manuscript; Supervision: L.R. and Y.G. All authors have read and agreed to the published version of the manuscript. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: All raw data involved in this study are kept in Guo's laboratory. To request the data, please contact the corresponding author of the article.