977
Views
2
CrossRef citations to date
0
Altmetric
Original Article

Water holding capacity and microstructure of sturgeon (Acipenser gueldenstaedti) fillets as affected by low temperature vacuum heating

, ORCID Icon, ORCID Icon, , &
Pages 1061-1073 | Received 22 Mar 2021, Accepted 29 Jun 2021, Published online: 25 Jul 2021

ABSTRACT

This study was aimed to investigate the effect of low-temperature vacuum heating (LTVH; versus the traditional thermal method– TC - 100°C boiling water) on water content/distributions, microstructure, and tenderness of sturgeon fillets. Although all fillets were well done, those fillets heated at lower temperature/time showed higher moisture content (P< .05). Furthermore, microstructure by paraffin section and shear force results also indicated LTVH70-30 (LTVH method; 70°C, 30 min) shown the most similar to TC fillets, the less firm and less juicy. Oppositely, LTVH60-15/30 (LTVH method, 60°C, 15/30 min) showed the best tenderness/juiciness among all fillets. Therefore, the LTVH method is agood alternative to the TC for fish processing in the catering industry, with tenderness/juiciness gains.

Graphical abstract of sturgeon fillets by low temperature vacuum heating GA

Introduction

Fish are the main superb source of animal protein, with ahigh nutritive value owing to their favorable essential amino acid composition which is necessary in adiet .[Citation1,Citation2] Although fish are sometimes eaten in raw (e.g., sashimi), fish are usually consumed after cooking in boiling water, steaming, frying, or subjected to other heating treatments. The main advantages of heating treatment are sterilization, which ensures safe consumption, increased digestibility, [Citation3] enhanced flavor and taste, and the increase shelf life .[Citation4] As aculinary tradition, boiling is popular in China, the USA, Japan, and other countries and has accepted the most attention in the context of home cooking .[Citation5]

In addition to their valuable caviar, sturgeons are also popular for their boneless flesh with excellent taste; in fact the sturgeon flesh accounts for approximately 300 g/kg skinless fillets .[Citation6] However, in contrast to marine fish and animal meat, sturgeon, as afreshwater fish, has asmaller muscle fiber diameter, less connective tissue, and atender taste .[Citation7] This means that the sturgeon fillets have poor thermal tolerance and their fiber muscles are much more easily broken during thermal processing. The natural features of the protein molecules are denatured during thermal treatment, and the reactive groups of proteins are exposed, resulting in loss of water, soluble vitamins, and proteins, and ultimately leading to changes in the surface and taste of food .[Citation3] Anumber of studies exploring the effects of thermal treatments on the eating quality of cooked fish occurring by protein denaturation have been reported.Citation8, Citation9, Citation10 Of note, the changes during processing depended on the processing method used and the heating time and temperature. For instance,Citation11reported that the heat treatment at high temperatures causes protein oxidation overhand, which consequently leads to changes in the physicochemical and organoleptic characteristics of food, including texture deterioration induced by the loss of water-holding capacity, and alteration of flavor. However, thermal treatment of food at lower heating temperatures usually requires alonger heating time, which is generally not appropriate for industrial applications because of the lower productivity of this cooking method .[Citation12]

Low-temperature vacuum heating (LTVH) is similar to the traditional domestic boiling method. Samples are heated in boiling water, but at mild temperatures (under 100°C) vis the use of special self-made equipment that allows controllable vacuum conditions. Our previous study showed that the LTVH method can reduce the oxidation of protein and lipid in fillets, compared to those observed in the context of traditional domestic cooking (100°C boiling water) .[Citation13] Additionally, the LTVH fillets had agreat shelf life after 7 days storage, performing abetter storage quality than the control groups thermal treatment at 100°C boiling water .[Citation14] However, there is still limited data available on the water-holding capacity and tenderness of sturgeon fillets after processing via different heating time/temperature combinations.

Hence, the aim of this study was to evaluate the effect of heating time/temperature combinations by LTVH method (using once again the traditional heating method as the control condition), on water-holding capacity (WHC), water distributions, and the microstructure of sturgeon fillets. Abetter understanding of the water mobility and fiber structure changes induced in fillets during heat-mediated processing is essential to elucidate their role in the protein matrix structure and for the development of new healthy aquatic products.

Materials and methods

Sample preparation

Frozen Russian sturgeons, male, 3year old and 3 ± 1kg weight, provided by Quzhou Sturgeon Aquatic Food Technology Development Co. Ltd, Zhejiang, China, were transported to the laboratory within 2 h. The middle part of the sturgeon was cut into fillets (approximately 1.0cm thickness and 65 to 75 g weight) after thawing with running water for half anhour. The fillets were bathed without any package in the boiling water using the self-made equipment under controllable vacuum conditions and randomly divided into 5 groups: traditional thermal method (TC, fillets were treated at 100°C– boiling water– for 15 min under normal pressure), LTVH60 (including LTVH60-15 and LTVH60-30, fillets were treated at 60°C for 15 or 30 min, respectively, under a− 812 mbar pressure), and LTVH70 (including LTVH70-15 and LTVH70-30, fillets were treated at 70°C for 15 or 30 min, respectively, at a− 700 mbar pressure). After thermal treatment, the fillets were collected, placed into with plastic bags, and cooled to room temperature in an ice bath for subsequent quantitative and qualitative analyses.

Cooking loss rate

The cooking loss rate was calculated based on the weight loss (before versue) after heating treatment. Determinations were carried out in triplicate. The performed equation is shown to be the following equation:

Eq.1 Cookinglossrate(%)=[(W1W2)/W1]×100Eq.1

where W1 is the mass (g) of the sturgeon fillet before heating and W2 is the mass (g) of the sturgeon fillet after heating.

Moisture content

The moisture content was determined based on the weight loss after drying 5 g of heated fillet in an assisted air circulation oven (GZX-9140MBE, Boxun Industrial Co., Ltd., Shanghai, China) at 105°C to aconstant weight (differences of weight limited to 2mg). Determinations were duplicated in triplicate. The data were calculated with the following equation:

Eq.2 Moisturecontent(%)=[(m1m2)/m1]×100Eq.2

Where m1 is the mass (g) of the sturgeon fillet before drying and m2 is the mass (g) of the sturgeon fillet after drying to aconstant weight.

Water-holding capacity [WHC)

The WHC of fillets was determined according to the method described by,Citation15 using an high-speed refrigerated centrifuge (H2050R, Xiangyi Laboratory Instrument Development Co., Ltd, Hunan, China]. All determinations were carried out three times, independently. The following equation was used for the determination of the WHC:

Eq.3 WHC(%)=(C2/C1)×100Eq.3

Where C1 is the mass (g) of the sturgeon fillet before centrifugation (1000 × g for 10 min at 4°C) and C2 is the mass (g) of the sturgeon fillet after centrifugation (1000 × g for 10 min at 4°C), g.

Low-field nuclear magnetic resonance [LF-NMR)

The water distribution in sturgeon fillet subjected to thermal-treatment was determined via LF-NMR relaxation experiments according to the method of proposed by,Citation16Briefly, transverse relaxation time (T2] measurements were performed using aMeso MR23-060 V–I nuclear magnetic resonance analysis system (Niumag Electric Corporation, Shanghai, China). Astandard copper sulfate solution was placed in the RF coil (25mm diameter) center of the magnet case, under the following condition: magnet temperature of 32°C, the magnet strength of 0.5 T. Using the FID sequence, the center frequency SFO1 was 23.20 MHz, automatically find the 90° and 180° pulse widths, determined the 90° pulse time (P1) was 30.00 μs, and the 180° pulse time (P2) was 56.00 μs. T2 values were measured by collecting CPMG (Carr-Purcell-Meiboom-Gill) decay signals with the sampling parameters are set as follows: repeat sampling waiting time (Tw) of 2000 ms, analog gain (RG1) of 20, digital gain (DRG) of 2, pre-amplification gain (PRG) of 1, accumulation number (NS) of 8, echo number (NECH) of 1000, and delay (DL1) of 0.1 ms. The Multi-Exp Inv Analysis software (Niumag Electric Corporation, Shanghai, China) was used for the distributed multi-exponential fitting of the CPMG decay curves and visualization of the lateral relaxation spectrum of the fillets. Three independent scan repetitions were conducted for each group, and the average value was recorded.

Proton density imaging

The proton density images of fillets subjected to thermal treatment were acquired by using aMINI MR-60 device [Niumag Electric Corporation, Shanghai, China) according to the method of,Citation17and with minor modifications. The device was operated at afrequency of 23 MHz and 32°C. The main parameters were set as follows: repetition time (TR] of 500 ms, viewing field of 50mm×50mm, and echo time (TE) of 19.0 ms. The final images were obtained by using the respective image processing software [Shanghai Niumag Corporation, Shanghai, China).

Extraction of myofibrillar proteins

The extraction and calculation of myofibrillar proteins was performed as described by,Citation13The determinations were repeated three times. The extraction rate of myofibrillar proteins was calculated as follows (Citation16):

Eq.4 Extractionrateofmyofibrillarprotein(%)=E1/E2×100Eq.4

where E1 is the quantity (mg) of myofibrillar proteins extracted from heated sturgeon fillet and E2 is the quantity [mg) of myofibrillar protein extracted from raw sturgeon fillet.

Light microscopy microstructural analysis

The microstructures of the heated fillet were investigated following the method of reported by.Citation18 Briefly, the fillet samples were fixed in fixative (including 4% paraformaldehyde], embedded in paraffin, and then cut (5–8 μm thick sections) using amicrotome (RM2016, Leica Instruments Co., Ltd., Shanghai, China). After de-paraffinization with xylene, the samples were stained with Ehrlich’s hematoxylin for 5 min and dehydrated in absolute ethanol for 5 min twice. The microstructures of the samples were observed and photographed using alight microscope (Nikon Eclipse E100, Nikon Co., Tokyo, Japan) coupled to an imaging system [Nikon DS-U3, Nikon Co., Tokyo, Japan).

Shear force analysis

According to the method reported by,Citation19but with minor modifications, the shear force analysis in the context of sturgeon fillets was performed using atexture analyzer (TA.XT Express, Stable Micro Systems, Surrey, UK] combined with aWarner Bratzler shearing device and aspeed of 5cm/min and adistance of 25cm. Six replicates were used per fillet.

Statistical analysis

All results are expressed as mean values ± standard deviation. Duncan’s multiple comparison test was used to determine statistical significance, set at P< .05. Pearson’s correlation analysis was performed to understand the relationships between the parameters.

Results and discussion

Extraction rate of myofibrillar proteins

Myofibrillar proteins are the most abundant proteins in the sturgeon (Acipenser gueldenstaedti) fillets, constituting approximately 50–55% of the total protein content .[Citation13] Therefore, changes in the myofibrillar proteins during thermal treatment play an important role in the alteration of the characteristics of sturgeon fillets. In daily life, we usually judge the readiness of cooked fillets based on their visual appearance and color. For instance, well-done fillets are pale, grayish-brown, and dry. Of note,,Citation16reported that the abalones were effectively cooked when the extraction rate of myofibrillar proteins was limited to less than 10%. As shown in , the extraction rates of myofibrillar proteins in all heat-treated fillets were below 10%. These results suggest that all fillets (LTVH and TC) were well done. Of note, the extraction rates of LTVH fillets decreased with the increase in the heating temperature (P< .05) and heating time (P> .05) (). Therefore, altogether, these results suggest that increased heating temperatures have agreater impact on decreased myofibrillar proteins than prolonged heating time.Citation20 reported that the myofibrillar proteins were sensitive to heating temperatures and get denatured during thermal treatment, resulting in protein aggregation and decreased extraction rate. Notably, as shown in , there were no significant differences between the extraction rates of the LTVH70-15, LTVH70-30, and TC groups. The denaturation of myofibrillar proteins may have been completed at the respective temperatures.,Citation21reported that proteolysis and aggregation occur together during thermal treatment. The insignificant difference in LTVH70 and TC was presumably due to the combination of protein aggregation and degradation at high temperatures.,Citation22,and Citation23also reported that protein digestibility decreased as the heating temperature increased, which was attributed to protein aggregation. Overall, in our study, LTVH60 fillets exhibited better protein retention than both LTVH70 fillets and TC fillets.

Table 1. Changes on the extraction rate of myosfibrillar protein, cooking loss, moisture content and shears force of sturgeon fillets by low temperature vacuum heating

Cooking loss rate, moisture content and WHC

The cooking loss rate, playing an important role in food processing, is responsible for the juiciness and yield of final products .[Citation24] Changes in the cooking loss rates and moisture contents of the different groups of fillets in this study are presented in . The heating time played acritical role in cooking loss rates. When the fillets were heated at 60°C for 15 min (LTVH60-15), the cooking loss rate was 27.56% ± 0.53%; while significantly increased to 31.26%±1.01% at LTVH60-30 (P< .05). However, there is no significant difference between the cooking loss rate of LTVH70-15 and LTVH70-30 (P> .05). Notably, the highest value of cooking loss rate was observed at TC fillets, and significantly higher than those heated by LTVH method (P< .05, ). Some water-soluble proteins were lost with the water during cooking process, [Citation25] which enhanced the cooking loss in the TC fillets. The results illustrate the advantages of the LTVH method compared with TC in serving better juiciness and yield of final products.

As presented in , the moisture content of the heated samples tended to decline with the heating time (P< .05) and heating temperature (P> .05) in fillets processed via the LTVH method. However, there was no significant difference between the WHC of LTVH and TC fillets (P> .05). There are three main states of water in muscles: bound water, immobilized water, and free water [introduced in moisture distribution).,Citation26reported that immobilized water and free water, entrapped in the myofibrillar network or within the space between the thick and thin filaments, are sensitive to thermal treatment. It is, therefore, not surprising that increased the heating temperature or prolonging the heating time leads to water evaporation, and consequently, shrinking of myofibrils proteins (due to the weakened-binding forces between proteins and water] .[Citation27] In line with this notion, here, the conducting the boiling process at low temperature or for less time persevered more water in the fillet structures; the LTVH60-15 samples showed the highest moisture content and WHC (68. ± 0.19%, 75. ± 1.66%, respectively; ). We hypothesize that the capillary water loss during thermal treatment was one of the reasons for the reduction of free water in fillets, as reported by others .[Citation28] Overall, our results suggest that LTVH fillets treated at lower temperature had batter moisture content than TC fillets, probably due to the preservation of protein solubility.

Moisture distribution

Moisture distribution and migration are influenced by the heating temperature and heating time during thermal treatment. The different of T2 relaxation spectra obtained in this study are shown in ; of note, these data exhibit astrong correlation with water-holding capacity .[Citation29] As shown in , three peaks, T21, T22, and T23, are observed in the sturgeon fillets processed using low-temperature vacuum heating, indicating three distinct water states; T21 (ranging from 0 to 10 ms), T22 (ranging from 25 to 35 ms) and T23 (ranging from 140 to 1000 ms) were assigned to bound water, immobile water, and free water, respectively (). These results are in line with previous studies reporting that the existence of three water states in heated muscle tissues, [Citation30,Citation31] with apredominance of immobile water.

Figure 1. Changes on NMR T2 relaxation time of sturgeon fillets by low temperature vacuum heating

LTVH60-15/LTVH60-30, fillets were treated at 60°C for 15/30 min under −812 mbar vacuum pressure; LTVH70-15/LTVH70-30, fillets were treated at 70°C for 15/30 min under −700 mbar vacuum pressure; TC, fillets were treated at 100°C boiling water for 15 min with normal pressure. Aand Cwere enlarged plot of T21 peak (green dotted rectangle) and T23 peak (black dotted rectangle) in B, respectively;
Figure 1. Changes on NMR T2 relaxation time of sturgeon fillets by low temperature vacuum heating

As shown in , the relaxation time T21 shifted downward with the heating temperature and heating time, and its amplitude of T21 declined. However, there was no significant difference in T21, except for LTVH60-30 (, P> .05). As mentioned above, the relation time T21 represents bound water components, which bind to protein molecules via hydrogen bonds, having the lowest mobility in fillet muscles, [Citation32] which is responsible for its resistance to thermal treatment .[Citation26] Interestingly, Pearson’s correlation analysis showed that there was avery significant positive correlation (R = 0.910, P< .01, ) between the relaxation time T21 and the extraction rate of myofibrillar proteins in heated fillets. As shown in , A21, corresponding to the peak area of T21, increased with the heating temperature and heating time [P> .05), suggesting that thermal treatment indeed influences protein–water interactions.,Citation33also reported that the reduction in protein solubility was responsible for the decreased water-holding capacity due to weakening of the protein–water interactions and strengthening of the protein–protein interactions. In line with this nation, T21 of the LTVH70 fillets and TC fillets changed very little while that of LTVH60 fillets changed much more as myofibrillar proteins varied.

Table 2. Changes on NMR parameters obtained from sturgeon fillets by low temperature vacuum heating

Table 3. Pearson’s correlation measurement of indicators from sturgeon fillets by low temperature vacuum heating

As shown by the NMR date, immobilized water, accounting for more than 90% of moisture, was the highest in fillets after thermal treatment (]. Of note, the relaxation time T22 and T23 decreased, as the relaxation components shifted to the left in , indicating that water mobility weakened as the heating temperature/time increased. Notably, the peak area proportion of immobilized water (A22) was inclined, while that of free water [A23) increased with both the heating temperature and time.,Citation34reported that changes in water may be correlated with the shrinkage and toughening of the fillet structure and cooking loss due to heating. In our study, the correlation analysis revealed avery significant positive correlation (R = 0.779, P< .01] between the peak area proportion A23 and the cooking loss rate in heated fillets. These results suggest that immobilized water was transformed into free water, and then released into the broil during thermal treatment, causing an increase in the cooking loss rate.

Magnetic resonance imaging (MRI) is anoninvasive technique that has been used to analyze the water distribution in food .[Citation35] In this study, the distribution and migration patterns of water in fillets are also presented in the form of MRI-derived proton density images (; the high proton densities are highlighted in red, while low proton densities are highlighted in blue). From the red zone ratios of these images, it is clear that the thermal treatment had asignificant influence on the content and distribution of water in the fillets. Notably, there were some red lines obviously present in fillets (), demarking the connective tissue wrapping the myofibrils in fillet structures.,Citation36speculated that the proteins in the connective tissues can restore main water loss during thermal treatment. Our data suggests that the water, existing between muscle filaments, myofibrils, and membranes, did not easily flow because of the strong interactions with proteins. As shown in , the moisture content decreased when the heating temperatures increased. However, no significant difference was observed in fillets heated via the LTVH method at the same temperature but during different time, which is in line with the results of T2 relaxation. Of note, TC fillets showed fewer red zones than LTVH fillets, which may be due to the changes in the muscle tissues in the context of high-temperature boiling water.

Figure 2. Changes on proton density images of sturgeon fillets by low temperature vacuum heating

LTVH60-15/LTVH60-30, fillets were treated at 60°C for 15/30 min under −812 mbar vacuum pressure; LTVH70-15/LTVH70-30, fillets were treated at 70°C for 15/30 min under −700 mbar vacuum pressure; TC, fillets were treated at 100°C boiling water for 15 min with normal pressure.
Figure 2. Changes on proton density images of sturgeon fillets by low temperature vacuum heating

Microstructure

The microstructural changes of fillets subjected to different thermal treatments were observed on transverse and longitudinal sections of paraffin-embedded tissues via light microscopy are shown in . The tissue structure of the fillets was significantly affected by thermal treatment. Denaturation and shrinking of myofibrillar proteins and connective tissues were occurred during the heating process, causing achange in the tissue microstructure; gaps were clearly observed in fillets, as reported elsewhere .[Citation37] However, fillets showed better tissue integrity after low-temperature heating processing (LTVH method). As shown in , the LTVH60-15 fillets showed the most intact tissue architecture; in contrast, wider intermyofibrillar spaces and expansion of muscle fibers, were clearly observed in the LTVH60-30 fillets. The same trend was observed in the LTVH70-15 fillets versus LTVH60-30 fillets. However, when the heating time was increased to 30 minin the context of the 70°C LTVH processing, some breaks were observed in the muscle fibers [LTVH70-30). Notably, the myofibril bundles of the TC fillets were severely broken with adissolved myofibrillar package and larger gaps. These results are in line with the results of the study by,Citation38which reported that treatment at high temperatures caused shrinking and dissolution of the myofibrillar fibers, resulting in the re-arrangement of broken myofibril proteins. Of note, large intercellular spaces were observed in LTVH70-30 fillets and TC fillets and to agreater extent in the latter, in where some solutes produced after the dissolution/degradation of cellular networks were fulfilled. Altogether, these results indicate that fillets treated using the LTVH method have better tenderness and WHC ( and4].

Figure 3. Changes on microstructure graphs of transverse (A) and longitudinal (a) sections in sturgeon fillets by low temperature vacuum heating (400ⅹ)

LTVH60-15/LTVH60-30, fillets were treated at 60°C for 15/30 min under −812 mbar vacuum pressure; LTVH70-15/LTVH70-30, fillets were treated at 70°C for 15/30 min under −700 mbar vacuum pressure; TC, fillets were treated at 100°C boiling water for 15 min with normal pressure.
Figure 3. Changes on microstructure graphs of transverse (A) and longitudinal (a) sections in sturgeon fillets by low temperature vacuum heating (400ⅹ)

Shear force

Tenderness is one of the most critical characteristics of fillets after thermal treatment. As shown in , thermal treatment, especially the heating temperature had asignificant influence on the shear force value of the fillets. With respect to the LTVH60-15 and LTVH60-30 fillets (treated at 60°C with −812 mbar), the low heating temperature led to low shear force values, which consequently leading to abetter tenderness of the fillets. Structural alterations of proteins are associated with simultaneous changes in the physicochemical properties of the water within the fillet structure .[Citation39] An earlier study had reported that changes in shear force values were mainly related to the changes in the tissue strength and water-holding capacity, which was responsible by the changes in myofibrillar .[Citation40] As mentioned previously, fillets were treated at low heating temperatures, which led to the better preservation of water and the tissue structures, resulting in the tenderization of LTVH60-15 and LTVH60-30 fillets.,Citation41investigated the tenderness of meat after thermal processing (also using the Warner-Bratzler shearing device) and their results suggested that the tenderness of meat heated at 60°C was better than that of meat heated at 50°C. However, fillets treated at 70°C under −700 mbar or at 100°C under normal pressure, showed higher shear force values (non-statistically difference among the LTVH70-15, LTVH70-30, and TC groups). These results are in line with aforementioned results that the muscle tissues were damaged severely occurring to the decline in water-holding capacity during high processing temperature.

Figure 4. Changes on shear force of sturgeon fillets by low temperature vacuum heating

LTVH60-15/LTVH60-30, fillets were treated at 60°C for 15/30 min under −812 mbar vacuum pressure; LTVH70-15/LTVH70-30, fillets were treated at 70°C for 15/30 min under −700 mbar vacuum pressure; TC, fillets were treated at 100°C boiling water for 15 min with normal pressure. Different small letters indicate significant difference at P< .05.
Figure 4. Changes on shear force of sturgeon fillets by low temperature vacuum heating

Conclusion

Although all fillets investigated in this study were well-done fillets (extraction rate of myofibrillar proteins was below 10%), different heating time/temperature combinations, resulted in different characteristics in sturgeon fillets. Due to treated at lower temperature and vacuum condition with boiling water, LTVH fillets showed better results in protecting the thermal loss in myofibrillar proteins and moisture. As the results of integrity of muscle tissue and moisture retention capacity, LTVH60 fillet performed best tenderness values out of all of the fillets via thermal processing. TC fillets, heated at 100°C boiling water under normal pressure, showed apredominant loss of structure at the microstructural level. This was closely linked to lower firmness and juiciness, which consumers perceived as undesirable. Therefore, the application of the LTVH method is agood alternative to the conventional processing fish (below 100°C boiling water) in the catering industry; notably, the final product is expected to be higher quality, with better tenderness and lower cooking loss. This would benefit both customers and the food processing industry.

Highlights

  1. The comparison between LTVH method and domestic heating at 100°C boiling water.

  2. The effects on water-holding capacity and microstructure of fillets investigated.

  3. The water distribution in fillet was observed by LF-NMR relaxation experiment.

  4. All fillets were well done and LTVH60 samples preserved more moisture.

  5. LTVH method can be served to achieve better tenderness and juicy fillets.

Acknowledgments

This work was supported by National Key R&D Program of China (Grant No.2018YFD0400600), Department of Education of Zhejiang Province (Grant No. Y201924135), and Zhejiang Provincial Natural Science Foundation (Grant No. LQ21C200004). And thanks for Quzhou Sturgeon Aquatic Food Technology Development Co. Ltd, Zhejiang China supported sturgeons for this study. The authors declare no conflict of interest.

Additional information

Funding

This work was supported by National Key R&D Program of China (Grant No. 2018YFD0400600), Department of Education of Zhejiang Province (Grant No. Y201924135), and Zhejiang Provincial Natural Science Foundation (Grant No.  LQ21C200004). 

References