3.1. Fungal bioaerosol concentrations
The mean concentration and variation of fungal bioaerosol at various sampling sites in the Mansurpur sugar mill is provided in Figure 2. The highest mean fungal concentration was observed at the cutter site (4022 ± 321 cfu/m3), followed by the bagasse (2026 ± 157 cfu/m3) and mill sites (1897± 152 cfu/m3). At packing, market, and residential sites, obtained concentrations were 1116 ± 117, 1411 ± 125, and 1313 ± 115 cfu/m3 respectively. Whereas the lowest mean fungal concentration was observed at the storage site, i.e., 832 ± 85 cfu/m3, during the sampling period. Additionally, the cutter site contributed 1.12, 0.99, 2.60, 3.83, 1.85, and 2.06 times more load to the fungal bioaerosol as compared with that of the mill, bagasse, packing, storage, market, and residential sites, respectively.
One of the reasons for maximum concentration at cutter site could be due to fact that it is crowded with humans, tractors, bullock carts, etc. They usually keep standing in queue to unload their sugarcane. Pushing tons of sugarcane could increase the emission of bioaerosol at this site. Furthermore, the persistent human activity that releases the enormous number of microbes along with organic droplets while talking, wheezing, and coughing, possess protection to the microbes present in the air and enhance their survival at that particular site for a longer duration (Balyan et al., 2020). Additionally, at the mill site, the high concentration could be due to the presence of huge installed machines, conveyor belts, rake carrier (carrying fiberized sugarcane) and tunnels, which might have influenced/enhanced the fungal concentration. Whereas, the high bioaerosol concentration at the bagasse site is likely due to the generation of a huge amount of bagasse, i.e., 2320 tons per day or 97 tons per hour. This site consists of dry fibrous organic waste material from the sugarcane and decomposed organic litter that may act as a host for various pathogenic fungal bioaerosol. The previous studies conducted at different sites have reported the availability of nutrients and other environmental factors that affect the growth of fungi (Balyan et al., 2019). The low fungal bioaerosol loading at the storage site (where 10,000 Tons of sugar bags are stored) could be on account of a clean, open, and wide area. This site includes an active ventilation system. The past study reveals that the availability of a proper ventilation system provides an adequate amount of air exchange and is reason for the low bioaerosol concentration (Faridi et al., 2015).
3.2. Variation of fungal bioaerosol during the study period
Spatial variation of fungal bioaerosol concentration over the different months at various sampling sites is given in Figure 3.The obtained concentrations were high in January and low in March month at all the sampling sites. The maximum fungal bioaerosol concentration was observed at the cutter site (6834 ± 241 cfu/m3), followed by the mill (3918 ± 189 cfu/m3), and bagasse sites (3365 ± 185 cfu/m3). Whereas the minimum concentration was obtained during March month at the residential site (333 ± 39 cfu/m3). A low concentration was observed at the storage site during all the months. The spatial variation of fungal bioaerosol concentration levels were in following order: Cutter Site > Mill Site > Bagasse Site > Residential Site > Market Site > Packing Site > Storage Site. Further, the cutter site, which is highly polluted, contributes 2.19 times more than the less polluted site, such as the storage site. The detail of the monthly mean concentration of fungal bioaerosol at all the seven sites is listed in Table 1.
As, winter season is categorized by low temperatures, moderate relative humidity, less rainfall (Nayak et al., 1998), and usually lowest mixing height which favors the agglomeration of suspended particles present in the ambient air and thus the higher fungal load (Lal et al., 2017; Jothish & Nayar, 2004). Further, the absence of rain washout due to less rainfall throughout the winter months intensifies the growth of fungi (Nayak et al., 1998). The above-mentioned reasons might be the cause of the higher number of fungal bioaerosol concentrations at all the sampling sites during January month. The present study’s results are similar to the work performed by Agarwal et al. 2016, Huang et al. 2002 and Wu et al. 2000,according to the reported data, the highest fungal bioaerosol concentrations were observed during the winter season, while a study carried out at the historical museum, Egyptby Awad et al. 2020reported a high concentration of fungal bioaerosol during summer and lowest in fall. Another study was done in a largely industrialized region of Iran reported that maximum fungal microorganisms occurred in June and lowest during March month (Hosseini et al., 2021). Hence, geographical, meteorological, and other favourable conditions of a specific location promote the survival and growth of microbial organisms suspended in the air (Srivastava et al., 2021; Awad et al., 2020).
In addition, the concentration observed during January month is 5 times more than that of March month. The present study revealed that the highest fungal concentrations at all the sites were observed in January and February month which is in close proximity with the results of the previous study conducted in a coir, sawmill, chlor-alkali and flour mill, that had obtained the maximum spore count in winter months due to low rainfall, moderate temperature, and relatively high humidity, favouring fungal sporulation (Nayar et al., 2007; Jothish & Nayar, 2004; Nayak et al., 1998; Misra and Jamil, 1991). In this study, the meteorological parameters, suitable growth substrate, and crushed sugarcane, especially at the cutter site might have provided the growth of fungi.
3.3. Distribution of size-segregated fungal bioaerosol concentration
The bioaerosol stages are divided into two broad categories, fine (3.3 to < 0.65 microns (stage 4-6)) and coarse (3.3 to >7 micron (stage 1-3)) bioaerosols. The idea was basically to understand the variation (monthly and spatial) of bioaerosol that may perhaps differ on changing the size range. Figure 4 and Figure 5reveal that the maximum concentration of fungal bioaerosol was still high in January month at site 1 (Cutter), followed by site 2 (Mill) and site 3 (Bagasse). Further, minimum concentrations in all months were observed at site 5 (storage) and site 7 (residential) during March for both fine and coarse size fractions. Hence, the site didn't play any major role in the distribution of size segregated fungal bioaerosol.
To understand the dependence of fine and coarse size fungal bioaerosol with each other, regression analyses were carried out for each month. Considering fine as an independent and coarse as a dependent variable. Regression analyses results illustrate that fine bioaerosol controls the variation of coarse bioaerosol as shown in Supplementary Figure S2. It may be observed from the plots that there is a good regression between fine and coarse size, in all four months at a level of significance p ≤ 0.05. The R2 and p values are 0.8029 and 0.0063, 0.9587 and 0.0001, 0.845 and 0.0034 and, 0.6779 and 0.0228 for December, January, February, and March months respectively.
3.4. Statistical analysis of fungal bioaerosol
In the present work, the box plot has also been plotted as box plot in (Figure 6)represents the fungal bioaerosol concentration (cfu/m3) at various sampling sites. Further, the Fligner- Killeen test was performed to test the homogeneity of variances that showed a chi-square value of 17.739, and p-value of 0.006918, a significant difference was observed in the variance of all sampling sites (p < 0.05). Kruskal Wallis test was also performed for fungal bioaerosol concentrations at various sampling sites, which showed the statistic value of the chi-square test as 8.0172. However, the obtained p values (p = 0.2368) were not significant. Hence, no multivariate analysis was performed. Thus, Parametric One-Way ANOVA (Analysis of Variance) tests were applied for further analysis, followed by Tukey's multiple comparison test (Tukey's HSD) in order to acquire the comparative analysis of the fungal mean concentration in various sampling sites. One-Way ANOVA was carried out for fungal bioaerosol concentration, which showed a p-value of 0.0098 as provided in Supplementary Table S1. A significant difference was seen in the mean fungal bioaerosol concentration between other sampling sites at p < 0.05.
Thebox plot (Figure 7)represents the fungal concentration (cfu/m3) in different months. The Fligner- Killeen test was done to test the homogeneity of variances which represents a chi-square test statistic as 6.4232, and a p-value as p = 0.09274. Additionally, the Kruskal Wallis test was performed on fungal bioaerosol concentration in different sampling months, and the results of the test showed a chi-squared test statistic as 15.473, which was found to be significant at p = 0.001454 (p < 0.05). Therefore, Parametric One-Way ANOVA tests were applied for further analysis, followed by Tukey's multiple comparison test (Tukey's HSD) to compare the mean concentration of fungal bioaerosol in various sampling months. The result of the One-Way ANOVA test, showed a p-value of 0.006665 for fungal bioaerosol concentration has been shown in Supplementary Table S2. A significant difference was observed between the mean fungal bioaerosol concentration and different sampling months at p < 0.05. The outliers can be explained owing to the presence of relatively high relative humidity and low temperature throughout the sampling period at the cutter site which favours the growth of fungi with respect to other sites. As far as concentration relation with meteorological parameters is concerned, it may be noted that sites (cutter, mill, and bagasse) with higher relative humidity (as shown in Table 2) have more fungal bioaerosol concentration. Conversely, storage sites that have the minimum concentration with lowest level relative humidity amongst all. In an earlier study carried out by Baghani et al. 2020 similar results were found.
3.5. Inter-relationship between fungal bioaerosol concentration and meteorological parameters
Microclimatic conditions widely vary at different sites. The presence or absence of a ventilation system, i.e., the number of windows and doors, may affect bioaerosol load at a particular site. In the present study, the lowest temperature and high relative humidity were both associated with the cutter site, i.e., 19 ± 5 °C and 73 ± 21%, respectively, as compared to other sites. To validate whether there is any significant relationship among the meteorological parameters with the mean fungal bioaerosol concentration (cfu/m3), the Pearson correlation coefficient (r) test was performed. The p-value was calculated, and the correlation between the mean fungal bioaerosol concentration (cfu/m3), and wind speed (m/s), relative humidity (%), and temperature (°C) were estimated at the significance level of p < 0.05 for the sampling period. The results analyzed indicated that wind speed (m/s) and the concentration of fungal bioaerosol are positively correlated at p < 0.05 (p = 0.003 and r = 0.921). Additionally, the results of correlation analysis calculated for relative humidity (%) and concentration of fungal bioaerosol showed a positive correlation at p < 0.05 (p = 0.004 and r = 0.912) whereas a negative correlation for temperature was observed at p < 0.05 (p = 0.002 and r = - 0.927). This suggests that all the meteorological factors are significantly related to airborne fungi and play a major role in their occurrence and sustenance.
According to Awad et al., 2019, the interaction between daily variation in temperature, relative humidity, sunlight, wind speed, air stream passage, geographical factors and anthropogenic activities may synergistically affect fungal bioaerosol concentration over a particular location. Low temperature and relatively high relative humidity appear to favour the growth of fungal bioaerosol. Our result conforms with the reports of similar strong correlations between various meteorological parameters and concentrations of fungal bioaerosol, which includes relative humidity and temperature in a cardboard and waste paper recycling factory in Iran (Baghani et al., 2020). It is also in close similarity with some other works, done by Awad et al. 2019 on indoor fungal pollution in a historical museum in Egypt, Patil and Kakde in 2017 on a landfill site in Mumbai as well as on landfill site in Barranquilla by Gamero et al. in 2018. Whereas, a study conducted at various sites in Delhi (Lal et al., 2017) reported no consequent association between relative humidity and temperature with fungal bioaerosol.
3.6. Correlation matrix between different stages
Correlation matrix of fungal bioaerosol concentration was developed to test the degree of association between different stages during the sampling period. The p-value was calculated, and the significance level was set at p < 0.01 and p < 0.05. The results were analyzed using a 2-tailed test. Table 2 shows the strong correlation between different stages of fungal bioaerosol with each other, where color codes varying from yellow to dark green represent a strong correlation. The strong relationship between different stages means that change in fungal bioaerosol concentration in one stage is strongly correlated with the changes in the fungal bioaerosol concentration in other stages.
Therefore, in accordance with statistics, all the stages are significantly correlated with each other, except stage 6, which is not significantly correlated with other stages. Table 2 showed thata strong correlation exists between different stages during the sampling period. The correlation and r values between the stages are as follows: stage 1 and 2 (r = 0.953); stage 1 and 3 (r = 0.912); stage 1 and 5 (r = 0.937); stage 4 and 5 (r = 0.900), whereas good correlation exists among stage 1 and 4 (r = 0.811); stage 2 and 5 (r = 0.818); stage 2 and 4 (r = 0.875); stage 2 and 3 (r = 0.863); stage 3 and 4 (r = 0.894); and lastly, stage 5 and 6 (r = 0.849). Hence, a conclusion can be drawn from the correlations as mentioned above that stage 1 has the maximum number of correlations with other stages followed by stages 2, 3, and 4. However, in contrast to this, stages 5 and 6 do not have a significant correlation with other stages. The possible reason for such a correlation pattern might be due to the fact that the spores of fungi are mostly kept attached with various size fractions of particulate matter. Thereby, they stay separated as per their size fraction in various stages of the air sampler. The present study confers that the most dominant spores of fungi exist in size range varying from stages 1 to 5, and a low concentration of spores was detected in stage 6.
This goes well with the fungal bioaerosol study conducted by Maharia and Srivastava, 2015 and Kumar et al. 2021 that reported a similar correlation and distribution of fungal bioaerosol in the range of stages 1 to 5.
3.7. Health risk assessment
The risk to human health owing to the exposure of fungal bioaerosol generated in and around the sugar mill was evaluated as per the USEPA guidelines. The calculations were performed using the mean fungal bioaerosol concentration during the sampling period. The health index (HI) for HIinhalation and HIdermal values were calculated as 4.83 x 10-12 and 3.82 x 10-7, respectively. Further, the observed HIinhalation value was 105 higher than that of the HIdermal value. The higher HIinhalation indicated that the primary route of bioaerosols intake by workers was through the inhalation route. The findings of some other works are similar to the present work. For instance, the study conducted by Li et al. 2013 and Wang et al. 2018 at the wastewater treatment plant in Xi'an and Tianjin, China, respectively, observed that on-site workers are mainly exposed through inhalation. However, the present study results are inconsistent with the study conducted at the irrigation of municipal wastewater site (Carlander et al., 2009). According to their results, ingestion plays a crucial and predominant role when compared to another exposure pathways. The comparison of findings of this study with the other carried out studies representing Mean Health Quotient (HQ) for dermal and inhalation are shown in Table 3.
It was also found that the lifetime average daily dose (LADDdermal) was highest at the cutter site (1.52 x 10-4 cfu (kg d)-1), whereas lowest at the storage site (3.15 x 10-5 cfu (kg d)-1). Furthermore, the maximum exposure due to lifetime average daily dose (LADDinhalation) was observed at the cutter site (192.482 cfu (kg d)-1) and minimum at the storage site (39.827 cfu (kg d)-1). Thus, it can be the inference that LADDinhalation is ~5 times over LADDdermal. As shown in Figure S3 and S4,the trend of both LADDdermal and LADDinhalation values for the sampling sites are as follows: Cutter site > Bagasse site > Mill site > Market site > Residential site > Packing site > Storage site.
According to the previous studies, HI >1 or HQ>1, indicate the potential risk or harmful health effects such as carcinogenic effects. While, HI <1 or HQ <1, suggests acceptable hazard risk levels such as non-carcinogenic effects, which are not of great concern since the dose level is lesser than the reference dose (RfD) (Dehghani et al., 2018; Baghani et al., 2020; USEPA, 2011b, Deng et al., 2014).In the present study, both HQ and HI calculated were observed to be less than unity (HQ and H >1), which suggests an "acceptable hazard risk level".
The risk assessment carried out at the sugar mill provides a valuable perspective in order to understand the potential risks of bioaerosol emitted from the sugar mill. Although the calculation of risk indices is based on the absolute data, there could be a possibility that present results might underestimate the risk of the workers. However, the significant exposure time for a longer duration, less awareness about safety measures and improper food nutrition may develop various visible health symptoms among workers.
Lastly, the "Heat Map" was also generated to represent the variation of extent of health risk quotient factor related to fungal bioaerosol at different sampling sites during the whole sampling period. It is illustrated in Figure 8. The color code expressed the extent of risk exposure as the red color shows a prominent and maximum index, whereas the green color shows fewer prevailing effects. This includes lifetime exposure of fungal bioaerosol through LADDinhalation + LADDdermal pathways. The different colors were generated on the basis of the observed Z-score of the fungal bioaerosols at different sampling sites in the mill among all the sampling periods.
3.8. Fungal identification
Total 8 airborne fungal species from sampled bioaerosol were identified from all the sites of sugar mill during the sampling period. They are given in supplementary Table S3. Cladosporium one of the genera found is commonly observed in food products and dead plants. It may easily grow in cold weather conditions and damped areas. It is mostly known to cause asthma and allergic reactions in human beings. Cladosporium genera also secrete secondary metabolites, i.e., mycotoxins. Emodin and Cladosporin are the most common among all the types of secreted mycotoxin by Cladosporium species. The most common infections caused by Cladosporium species in humans are chromoblastomycosis which is a chronic infection of the skin and subcutaneous tissue i.e., the deepest layer of the skin (Ogórek et al., 2012).
Further, Aspergillus niger is a common fungus known to produce aflatoxins which is a classic mycotoxin example. It is known as highly toxic to many animals and humans. In humans, ingestion of aflatoxin has been associated with liver cancer along with Hepatitis B as a risk factor, and inhalation causes lung cancer and other lung-related diseases (Diba et al., 2007). Penicillium is predominantly found in food, grains, soil, indoor house dust, and damped buildings.Previous studies reported that numerous hypersensitivity pneumonitis epidemics are caused by Penicillium species (Larone et al., 1995). Aspergillus oryzae does not produce aflatoxins or any other toxic metabolites (Gomi, 2014). Fusarium causes keratomycosis (infection of the cornea) (Balajee et al., 2009). Rhizopus is an agent of mucormycosis that includes rhinocerebral mucormycosis, gastrointestinal, mucocutaneous, pulmonary & disseminated infections in immunocompetent patients (Hoog et al., 2000). Curvularia has been widely known to cause allergic fungal sinusitis (Hoog et al., 2000). Trichoderma inhalation causes lung infection, allergic reactions and asthma (Misra and Jamil, 1991). Alternaria infects crops worldwide and have been known to cause allergic rhinitis when inhaled (Thomma et al., 2003).
Penicillium, Alternaria and Aspergillus were the most abundant inside the sugar mill, whereas the Cladosporium was found outdoor. The presence of Penicillium, Aspergillus, and Cladosporium, in abundance, indicates an allergic environment for humans.