The COVID-19 pandemic has impacted the health and well-being of the global population and poses a significant threat to the food security and nutrition of millions of people across the world(1). Different preventive measures taken by governments to contain COVID-19 (e.g. staying at home orders, closure of industries, border restrictions, etc.) may pose additional barriers to food systems, including marketing, logistics and trading, which could potentially impact food security(1). Moreover, the economic turmoil as a result of the pandemic has jeopardised physical and economic access to sufficient, and sufficiently nutritious, food(2). Recent evidence shows that prior to the end of the year 2020, food insecurity could nearly double to 265 million people worldwide due to the COVID-19 pandemic(3). Consequently, malnutrition, hunger, low dietary diversity and other nutrition-related problems are likely to increase due to COVID-19(2).
Household food security (HFS) is defined as when all members of a given household have physical and economic access to sufficient, safe and nutritious food at all times that can meet their dietary needs and food preferences for an active and healthy life(3). HFS is underpinned by food availability, accessibility, utilisation and stability(4). The novel COVID-19 pandemic has the ability to impact HFS by not only disrupting food systems but also jeopardising household incomes due to unemployment and physical access to food due to home quarantine orders(Reference Devereux, Béné and Hoddinott5). Several additional factors have been associated with HFS such as socio-economic status, family income, education level and family size(Reference Ahmed, Mozahid and Dhar6,Reference Saha, Frongillo and Alam7) . To tackle the potential impact of COVID-19, information is needed on how factors that influence HFS are currently being affected during the COVID-19 pandemic.
Another key indicator of a household’s ability to not only provide enough food, but a diverse range of foods for good nutrition, is dietary diversity. HDD is a proxy indicator measuring the economic ability of a household to access a variety of foods during a determined period and is widely used as an indicator of food security as a whole(Reference Huluka and Wondimagegnhu8). Previous studies revealed that higher HDD is associated with socio-economic status and HFS(Reference Huluka and Wondimagegnhu8). Further determinants, such as type of residency (urban v. rural area), type of family (joint family v. nuclear family) and occupation have a significant association with adequate HDD(Reference Mukherjee, Paul and Saha9).
In light of the COVID-19 pandemic, a new set of challenges for households to maintain a healthy and varied diet have arisen(Reference Naja and Hamadeh10). Only a handful of studies have explored food security, food system disruption and dietary diversity in the USA (11), China (12), India (29) and Bangladesh (13), during the COVID-19 pandemic. A recent study by Adams et al. noted that the number of families from the USA reporting ‘very low’ scores for food security increased by 20 % from data prior to COVID-19(Reference Adams, Caccavale and Smith11). Another recent study on dietary diversity among Chinese residents concluded that people living in areas with a high number of confirmed COVID-19 cases had a lower score of HDD(Reference Zhao, Li and Ke12). Lastly, one Bangladeshi study reported that families’ food insecurity increased by 51·7 % during the COVID-19 home quarantine (lockdown) period(Reference Hamadani, Hasan and Baldi13).
Although Bangladesh has achieved substantial macro-economic growth in the past decades, the COVID-19 pandemic has made a significant impact on households nationwide. A recent nationwide survey by BRAC, a development organisation in Bangladesh, showed that 93 % of participants have suffered a loss of earnings, with 54 % reporting no income in March 2020(14). Moreover, disruption in transportation systems has caused the dumping of perishable food products and dramatic price reductions, affecting food security for rural producers(15). Thus, despite several steps taken by the Bangladeshi government to address food insecurity during the COVID-19, such as relief, distributing cash aid, etc(Reference Chakrobarty, Rasheduzzaman and Basunia16), increased poverty and economic turmoil of this magnitude will likely drive many people to short-term food insecurity and poor dietary diversity.
Collectively, there is a lack of information about HFS and HDD amid the COVID-19 pandemic, and an exploration of these issues is urgently needed in Bangladesh. A better understanding of key determinants during the COVID-19 pandemic will help government agencies in the development and implementation of policies and practices that can positively impact food security and improve dietary diversity, in this and future shocks. Thus, this study aimed to explore HFS and household dietary diversity (HDD) status, and its associated factors among Bangladeshi households during the COVID-19 pandemic.
Methods
Study setting, design and participants
A community-based cross-sectional survey was conducted to assess the HFS and HDD of the Bangladeshi population from 7 September to 15 September 2020, approximately 5 months after the start of lockdown (e.g. staying at home, closure of industries, border restriction, etc.)(Reference Hamadani, Hasan and Baldi13). The Government of Bangladesh decided in early June to enforce zone-wise lockdowns across the region, dividing areas into three zones (red, yellow and green) based on the number of COVID-19 positive patients. Red, yellow and green were marked as high-, moderate- and low-risk zones, respectively(17). Considering these restrictions, data were collected using online platforms and face-to-face interviews. Online data were collected through social media (e.g. Facebook and WhatsApp) across the country using the snowball sampling method. Face-to-face data were collected from green and yellow zones only. At least one rural and one urban area (green or yellow zones only) from all nine divisions/states were chosen by the authors for face-to-face data collection, based on convenience of access. In total, 18 research staff were involved in face-to-face interviews. Prior to data collection, research staff was trained by the field research supervisor. Considering the limited number of studies on the topic, a 50 % response rate, 95 % confidence level and 2·5 % margin of error were used to calculate the needed sample size of 1537 to achieve 80 % power. In total, 1876 participants replied during the survey duration, and all participants were included in the final analyses. Data from 1164 (62 %) participants were collected randomly via face-to-face interviews, and data from 712 (38 %) participants were collected using online platforms. The eligibility criteria for the study was being a Bangladeshi resident, being an adult (≥ 18 years) and residing in Bangladesh during COVID-19.
Procedures
The pre-study questionnaire was written in English, and two bilingual researchers translated it into the local language (Bengali). Another independent bilingual expert back-translated the questionnaire to check for consistencies and to prevent any bias. The electronic forms were shared in the local language (Bengali). The survey questionnaire was pilot-tested with a small group (100 households) in both online and face-to-face surveys to ensure its clarity and to avoid any unnecessary/repeated questions. Face-to-face data were collected from the head of the household, whereas online data were collected from any adult member of the household who used internet. The authors distributed the survey link in all nine divisions of Bangladesh. By clicking on the link, the participant was automatically directed to study information and an informed consent page. If they agreed to participate in the survey, the participant had to fill in demographic information, after which, a set of questions followed. To capture households in divisions/states designated as red zone areas, social media avenues were targeted. Face-to-face data collection from door-to-door using appropriate personal protective equipment and maintaining social distance practices was employed in only yellow and green zones. A total of nine rural and nine urban areas were selected based on transport facilities and distances from respective divisional headquarters. All households were selected randomly from the selected areas. The study questionnaire included a short overview of the study context, purpose, procedures, confidentiality agreement and informed consent. The privacy of the electronic data was ensured using an anonymous form and storing data in a password-protected folder.
Survey contents
The survey consisted of 38 close-ended items, which took respondents approximately 20–25 min to complete. The questionnaire was split into four sections: participant characteristics (10 items), effects of COVID-19 on participants’ day-to-day life (5 items), HFS (11 items) and HDD (12 items). Socio-demographic data were collected on age, gender, educational status, occupation of household head, monthly family income, location of residence, and family size and type. Participants were asked about the effects of COVID-19 on income, occupation, food quality and food quantity. The questions to assess the COVID-19 effects on these variables were closed-ended and coded with numerical scores.
Measures
Household Food Security Scale
We evaluated HFS using the Household Food Security Scale (HHFS), a measure that reflects a household’s food security for the previous month(Reference Frongillo, Chowdhury and Ekström18). The HHFS consists of 11 items covering the topics of purchased rice and perishable food, cooking frequency, consumption of snacks and management strategies, with a score assigned for each item response(Reference Frongillo, Chowdhury and Ekström18). Higher and lower scores were assigned for responses which indicated a higher (i.e. more favourable) v. lower (i.e. less favourable) food security status, respectively. For example, for the frequency of daily cooking, a score of 1 was assigned to those who responded they typically never cook on a daily basis, which indicated that the household had little-to-no store of food to cook with. For those respondents who reported cooking four or more times per day, a score of 5 was assigned, which indicated that the household had an adequate supply of food to use for cooking(Reference Frongillo, Chowdhury and Ekström18). A higher score indicated a more favourable HHFS. The reliability of the HHFS was high (Cronbach’s α = 0·81).
Household Dietary Diversity Score
Household Dietary Diversity Score (HDDS) was calculated by summing up the number of food or food groups eaten over the past 24 h by any member of the household(Reference Swindale and Bilinsky19). In total, the 12 food groups (FG) were as follows: (FG1) cereals; (FG2) tubers and roots; (FG3) vegetables; (FG4) fruits; (FG5) meat, poultry, organ, etc.; (FG6) eggs; (FG7) fish and others seafood; (FG8) pulses, legumes and nuts; (FG9) milk and other dairy products; (FG10) oils and fats and butter; (FG11) sugar and honey; and (FG12) miscellaneous foods such as condiments and processed foods like snacks, and beverages. The authors’ assigned values for each group as ‘0’ for the negative answer (not consumed), or ‘1’ for the positive answer (consumed). Higher dietary diversity was indicated by a higher score, ranging from 0 to 12. We used the HDDS scale both as a continuous and categorical variable for analyses. For categorical HDD, scores were divided into three categories: high dietary diversity (7–12), medium dietary diversity (4–6) and low dietary diversity (0–3)(Reference Ahmed, Mozahid and Dhar6). The internal consistency of the HDDS scales was adequate (Cronbach’s α = 0·77).
Statistical analyses
Descriptive statistics were calculated for socio-demographic characteristics of the participants and variables assessing the impact of COVID-19 on respondents’ day-to-day life. Bivariate linear regression was employed to assess the associations between socio-demographic characteristics, the variables related to the impact of COVID-19 effect on daily life, HDD and HFS. Multicollinearity was checked via the variance inflation factor. The final model selection was computed using Hosmer and Lemeshow goodness-of-fit test, and the significance of variables was assessed using the Wald test. The estimates of the strengths of associations were demonstrated by the β (co-efficient) with a 95 % CI. All tests were two-tailed, and a P-value of <0·05 was set to determine statistical significance. Statistical analyses were performed using SPSS version 23.0 (IBM SPSS Statistics).
Results
The majority of the family/household heads were male (93·6 %), and above 40 years old (64·7 %). Nearly, half (45·9 %) of household heads were educated below the secondary level and 20·3 % were self-employed. The majority of the households were from rural areas (59·1 %), nuclear family (81·1 %), with a family size of ≤ 5 (65·5 %) and a monthly income of ≤ 20 000 BDT (238$) (65·4 %). Overall, 41·2 % of households had no refrigerator and 36·1 % of household heads reported getting their dietary/nutrition information from traditional media sources, such as television and radio.
Household dietary diversity score and household food security scores of participants
In the present study, the mean HDD score during the COVID-19 pandemic was 6·22 (sd 2·52). When HDD scores were split into the three categories (high/medium/low dietary diversity), 779 (41·5 %) households reported high dietary diversity, 842 (44·9 %) households had medium dietary diversity and 255 (13·6 %) households had low dietary diversity. The mean HFS score during the COVID-19 pandemic was 31·86 (sd 5·49). Table 1 represents the proportions of participants who consumed different food groups in the previous 24 h. The cereals group (100 %) had the highest consumption rate, followed by the vegetable group (79·9 %), with pulses, legumes and nuts (25·4 %), and oils and fats (26·9 %), the least reported consumption groups (Table 1).
Food groups consumption by Bangladeshi households during COVID-19 pandemic in the last 24 h (n 1876)
* Miscellaneous foods include condiments and most processed foods like snacks, beverages, tea, coffee, etc.
Association of socio-demographic variables with household dietary diversity and household food security score
Households from rural areas had lower HHD scores (β = −2·08, 95 % CI: −2·29, −1·87) and HFS scores (β = −4·68, 95 % CI: −5·14, −4·22) compared to households from urban areas. Household heads’ higher educational levels (below secondary, above secondary and a university graduate or above) were associated with increased scores on both the HDD and HFS, compared to the head of a household having no formal schooling. Households with a family head or main household earner who was involved in a private job, business, day labouring, farming or other occupation were significantly associated with lower HDD scores and HFS scores, compared to the families having a family head or main household earner who was involved in a government job. Households with a monthly income of ≤ 20 000 BDT had lower HDD scores (β = −3·04, 95 % CI: −3·24, −2·85) and HFS scores (β = −5·61, 95 % CI: −6·07, −5·16) compared to households with a monthly income of above 20 000 BDT. Family members getting dietary/nutrition information from online resources had significantly higher HDD scores (β = 1·57, 95 % CI: 1·07, 2·07), and families who did not know any dietary/nutrition information had lower HDD scores (β = −1·63, 95 % CI: −2·11, −1·16), compared to families who got dietary/nutrition information from health professionals (Table 2).
Associations between socio-demographic variables, household dietary diversity scores and household food security scores of Bangladesh households (n 1876)
R 2 = R-squared, AR 2 = adjusted R-squared.
* Others included outsourcing, rickshaw pulling, etc.
† Health professional included physician, dietitian, nutritionist, etc.
‡ Traditional media included TV, radio, newspaper, etc.
§ Others included friends, family members, books, etc.
** P < 0·01.
*** P < 0·001.
Variables related to COVID-19 impact on daily life, and household dietary diversity and household food security
A household head or main earner for the household who lost their main occupation due to COVID-19 had a significant negative association with HDD scores (β = −1·78, 95 % CI: −2·16, −1·39) and HFS scores (β = −4·61, 95 % CI: −5·43, −3·79). Similar results were seen for those individuals who had to change their occupation; lower HDD (β = −0·87, 95 % CI: −1·24, −0·50) and HFS (β = −2·83, 95 % CI: −3·62, −2·04) scores were observed compared to those who had retained their previous occupation during the COVID-19 pandemic. The households whose monthly income decreased as a result of the COVID-19 pandemic had lower HDD scores (β = −2·50, 95 % CI: −2·74, −2·27) and HFS scores (β = −4·54, 95 % CI: −5·07, −4·01) compared to those with same income before the pandemic. Households who reported that they faced increased food prices due to COVID-19 had significant negative associations with HDD (β = −0·95, 95 % CI: −1·43, −0·48) and HFS (β = −1·77, 95 % CI: −2·82, −0·73) scores. Families that did not get the same amount of food as before COVID-19 had a significant negative association with HDD score (β = −2·25, 95 % CI: −2·45, −2·04) and HFS score (β = −5·41, 95 % CI: −5·85, −4·98) compared to those that got the same amount of food as before COVID-19. Similarly, families that did not get the same type of food as they did prior to COVID-19 had a significant negative association with HDD (β = −2·74, 95 % CI: −2·94, −2·54) and HFS (β = −5·52, 95 % CI: −5·97, −5·08) scores compared to those households/families that were able to purchase the same types of food as they did prior to COVID-19 (Table 3).
Association between the impact of COVID-19 pandemic, the household dietary diversity scores and household food security scores of Bangladesh households (n 1876)
R 2 = R-squared; AR 2 = adjusted R-squared; HH = household head.
** P < 0·01.
*** P < 0·001.
Discussion
This study explored the determinants of HFS and HDD in Bangladesh during the COVID-19 pandemic. This is one of the first studies that has provided the data on HFS and HDD status during the COVID-19 pandemic. We found that both household socio-demographic variables and COVID-19-mediated factors were associated with HFS and HDD. The current COVID-19 pandemic has negatively impacted the food security and diet quality of millions of people around the globe(Reference Devereux, Béné and Hoddinott5,20) . The present study provides evidence of COVID-19’s impact on HDD and food security in a low-resource setting such as Bangladesh. Specifically, results from the current study show that the COVID-19 has impacted the economic status of many families, which in turn negatively affects their food security and dietary diversity. A recent Bangladeshi study provided evidence on the reduction of food insecurity during the COVID-19 lockdown compared to the pre-lockdown period(Reference Hamadani, Hasan and Baldi13).
Demographic determinants of household food security
One of the main findings from our study was the association between rural residence and lower food security scores. A recent report revealed a 15 % decrease in rural populations consuming three meals per day during COVID-19 compared to the pre-COVID-19 period(Reference Rahman and Matin21). Comparatively, respondents from lower socio-economic status groups (e.g. day labourers, farmers) and households that do not have a refrigerator had lower food security scores compared to their counterparts in the present study. Occupations other than government jobs do not have or have fewer unemployment benefit or wage compensation services and are unable to do their job from home, so their occupation can be affected by the lockdown, that may result in food insecurity(Reference Devereux, Béné and Hoddinott5). This finding was similar to prior studies of Bangladeshi residents before the pandemic that has shown households with lower socio-economic status have fewer capabilities and resources to produce their food, which can lead to an increase in food insecurity(Reference Saha, Frongillo and Alam7). A report during COVID-19 showed a 27 % decrease in per capita food expenditure in lower poverty regions(Reference Rahman and Matin21).
Higher educational status was positively associated with food security, which might be attributable to the fact that better-educated people can diversify their work environments and, through higher earning potential, have the ability to maintain their families’ food availability(Reference Ahmed, Mozahid and Dhar6). However, families with older family heads had more financial security, which might be explained by the fact that older people are more experienced in terms of working, and thus, have the means to ensure their families having stable access to a variety of food. Moreover, nuclear families had significantly higher food security. This finding suggests that higher household size contributes to lower HFS, which is similar to a previous Bangladeshi study(Reference Ahmed, Mozahid and Dhar6).
Demographic determinants of household dietary diversity
The mean HDD score of participants was 6·22 (sd 2·52), but was lower for households in rural areas and with less secure incomes. which is noticeably lower than a recent study in China during COVID-19 (9·7 (sd 2·1))(Reference Zhao, Li and Ke12), but somewhat similar to a study during the pre-COVID-19 period in India (6·28 (sd 1·3))(Reference Mukherjee, Paul and Saha9). Our data show that economic stresses imposed by COVID-19 play a significant role in reducing dietary diversity for the households reporting that the income of their household head was affected by COVID-19. This persistent low dietary diversity of many Bangladeshi households during COVID-19 implies the importance of tracking diet-related policies and practices (e.g. access to a variety of foods) during times where households are under additional strain, such as the COVID-19 pandemic. Our study revealed lower dietary diversity among households whose family heads were day labourers, private jobholders, businessman, farmers and others indicating a deficiency in dietary diversity among precariously employed households(Reference Rabbani22). Individuals who are self-employed, such as day labourers, do not have access to unemployment benefit or wage compensation systems, and their work is not conducted from home; therefore, due to restrictions on residents travel during COVID-19, their employment situation can be impacted resulting in less income to purchase foods(Reference Devereux, Béné and Hoddinott5). Respondents having less family income and whose family income was affected by COVID-19 experienced lower dietary diversity. This finding is supported by previous evidence that spending on food consumption is positively linked to greater diversity(Reference Rabbani22).
Furthermore, rural residence was associated with lower HDD scores in our study population, corresponding to other previous studies(Reference Zhao, Li and Ke12,Reference Andrissi, Mottini and Sebastiani23) . It is plausible that lower nutritional knowledge and awareness might lead rural people to less dietary diversity than urban people(Reference Saeidlou, Babaei and Ayremlou24). Furthermore, online ordering and delivery services, as well as consumption of junk foods, are common practices in urban areas which can increase dietary diversity albeit with nutritionally negative foods, among urban people: It has been shown before that dietary diversity increases with increased intake of junk food(Reference Bezerra and Sichieri25).
Higher educational status of the household head was significantly associated with higher dietary diversity scores, which may be due to improved nutritional awareness in the household(Reference Samuel, Darebo and Desta26). Farmer as family heads’ occupation had a significant association with low HDD scores. It is worth to mention that over 80 % of Bangladesh’s total population lives in rural areas and relies directly or indirectly on agriculture(27). In rural areas, there is potential for farming households to lose access to food from their production if they fall sick from COVID-19, are unable to access manpower and materials (e.g. fertilisers and seeds), or lose access to markets to sell their produce due to restrictions on trade and mobility(Reference Devereux, Béné and Hoddinott5). In our study, we found that the nuclear family had more dietary diversity compared to the joint family which is opposite to a previous finding(Reference Mukherjee, Paul and Saha9) and so likely depends on context.
Impacts of COVID-19 on household food security
In the present study, 71·8 % of respondents claimed that their income decreased due to COVID-19, and this was positively associated with household food insecurity. This finding aligns with another study indicating that approximately half of a families’ income fell below $1·90 (an indicator of extreme poverty) per day during COVID-19, which was associated with food insecurity(Reference Hamadani, Hasan and Baldi13). Another recent study revealed that 65·35 % and 74·68 % of wage earners had no earning source and did not receive a salary during COVID-19(Reference Sultana, Khan and Hossain28). A study in India during the COVID-19 reported that 60 % of surveyed farms’ income dropped by half(Reference Harris, Depenbusch and Pal29). This precarious situation might play a significant role behind food insecurity in a low-resource, high-agriculture country like Bangladesh.
Moreover, unemployment due to COVID-19 was significantly associated with lower food security scores, again suggesting that reduced income and poverty increased food insecurity during COVID-19(Reference Hamadani, Hasan and Baldi13). Collectively, these findings underscore the impact of economic stress on food access during shocks such as the pandemic(30). The World Food Program predicts that the worldwide pandemic has the capability of increasing food insecurity to 265 million people(3), thus, careful consideration of this situation is required as household food insecurity can also lead to negative health outcomes, such as maternal anxiety and depression(Reference Hadley and Patil31,Reference Taaffe McLearn, Minkovitz and Strobino32) .
Impacts of COVID-19 on household dietary diversity
Notably, 45·3 % and 61 % of the respondents reported that they did not get the same quantity or type of food as they did prior to COVID-19, respectively. These findings are similar to a study in India where 62 % of respondents reported that their diet was disrupted due to COVID-19(Reference Harris, Depenbusch and Pal29). Disrupted production and distribution of food due to COVID-19 have clear consequences such as decreased food diversity(Reference Harris, Depenbusch and Pal29,33) . In addition, increased food prices (94 %) due to COVID-19 decrease dietary diversity to a large extent, which has been found in previous studies pre-COVID-19(Reference Devereux, Béné and Hoddinott5).
Strengths and limitations
This study had several strengths. To the best of our knowledge, it was the first study among Bangladeshi residents that measured HFS and HDD during the COVID-19 pandemic. In addition, this study employed both face-to-face and online surveys for recruitment and sampling to ensure a broad reach of survey responses given the restrictions placed on the country due to the pandemic. Lastly, piloting the surveys prior to data collection added appropriateness and clarity to the survey tools for the context of the study and setting.
However, this study has some limitations. This study mainly includes the online users which may suggest sampling bias by unintentionally excluding those who do not have internet access. This may reduce the generalisability of the findings. In addition, due to the cross-sectional nature of the study, we were unable to evaluate any seasonal variation of HFS and HDD during the COVID-19 pandemic. Although the HDD scores can help determine food accessibility, they do not capture the amount of actual food consumption by households, and nor do they capture changes or reductions in diversity within food groups, which is important for nutrition. Lastly, all information was self-reported by the participants and may be open to reporting bias.
Conclusion
This study explored the socio-demographic and pandemic-mediated determinants of HFS and dietary diversity during the COVID-19 pandemic in Bangladesh. Rural residence, no formal education, occupation other than government job and low monthly income, loss of occupation and income due to the COVID-19 pandemic, and increased food prices are the major risk factors for low HFS and HDD scores. There is a need for more attention towards HFS and HDD during the COVID-19 pandemic, to maintain food security and health. Authorities should work to protect supply chains of food items in all parts of the community. Food and nutrition assistance needs to be a fundamental part of social protection programmes to ensure diverse food access for the most vulnerable populations by protecting their purchasing power and by directly providing food where necessary. COVID-19 awareness programmes should include HFS and dietary diversity during the pandemic.
Acknowledgements
Acknowledgements: The authors are grateful to all of the respondents who participated in this study voluntarily. Furthermore, the authors highly appreciate the contribution of Nomita Das (Research Associate, plant breeding division, Bangladesh Rice Research Institute, Gazipur, Bangladesh) and Abu Saim Mohammad Saikat during data collection period. Financial support: This research did not receive any specific grant from funding agencies in the public, commercial or not-for-profit sectors. Conflict of interest: The authors declare no conflicts of interest. Authorship: S.K.: conceptualisation, formal analysis and writing – original draft; M.H.A.B.: conceptualisation, data curation and writing – original draft; A.S.: conceptualisation, methodology, writing – original draft, and final approval; M.S.S.: writing – original draft and visualisation; K.B.: writing – review and editing; J.H.: writing – review and editing; M.M., I.J., M.T. A.: data curation, writing – original draft; M.S.I.K.: visualisation, review and editing, and validation. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects/patients were approved by the Research Ethical Committee (REC) of Department of Food Microbiology, Patuakhali Science and Technology University, Bangladesh (Approval No: FMB: 29/05/2020:07). Written informed consent or Electronic consent (e-consent) was obtained from all subjects/patients.
Open access
