3.1 Integrated Effects of Environmental Stressors on Microbial and Human Health
The present study integrated a meta-analysis and systematic review to evaluate the combined effects of environmental stressors—including light exposure, trace metals, dietary patterns, and bioplastic degradation—on microbial systems and human health outcomes. Data were pooled from diverse experimental settings, including aquatic microorganisms, SHIME® gut models, human observational studies, and composting trials of bioplastics. Statistical analyses were conducted to quantify effect sizes, assess heterogeneity, and explore potential interactions between stressors and biological outcomes.
The meta-analysis revealed significant inhibitory effects of combined stressors on photosynthetic microorganisms, particularly under simultaneous exposure to trace metals and ultraviolet (UV) light. For instance, Cheloni and Slaveykova (2018) demonstrated that combined copper and UV exposure resulted in a 35% reduction in photosynthetic efficiency relative to controls (Table 1). Similarly, Prasad and Zeeshan (2005) reported that Plectonema boryanum experienced statistically significant decreases in growth rates and chlorophyll content under cadmium and UV B stress (p < 0.001). Random-effects modeling indicated substantial heterogeneity among studies (I² = 72%), likely reflecting species-specific responses and variations in experimental light intensities. Subgroup analyses confirmed that cyanobacteria exhibited higher sensitivity compared with green algae (Figure 2), consistent with previous observations by He and Hader (2002). These findings underscore the importance of evaluating multi-stressor impacts rather than single-factor exposures, supporting frameworks proposed by Heugens et al. (2002) and Segner et al. (2014) for predicting ecological outcomes under complex environmental conditions.
Table 1. Biodegradability of Bioplastics under Aerobic Composting Conditions. This table summarizes biodegradation outcomes reported in controlled composting studies. Percentage biodegradation is treated as the primary effect size for meta-analysis, while exposure duration and experimental scale are used as indicators of study precision.
|
Study (Year)
|
Bioplastic Material
|
Environment
|
Temperature (°C)
|
Biodegradation (%)
|
Duration (days)
|
Scale
|
|
Cho (2011)
|
PCL/Starch
|
Compost
|
25
|
88.0
|
44
|
Laboratory
|
|
Ahn (2011)
|
PLA/Feather
|
Compost
|
58
|
53.0
|
60
|
Laboratory
|
|
Kale (2007)
|
PLA
|
Compost
|
65
|
84.0
|
58
|
Field
|
|
Weng (2011)
|
PHB
|
Compost
|
58
|
79.9
|
110
|
Laboratory
|
|
Tabasi (2015)
|
PHB
|
Compost
|
55
|
80.0
|
28
|
Laboratory
|
|
Javierre (2015)
|
Starch-based
|
Compost
|
58
|
85.0
|
90
|
Laboratory
|
|
Adamcova (2017)
|
Cellulose acetate (CA)
|
Compost
|
58
|
80.0
|
154
|
Laboratory
|
|
Anstey (2014)
|
PBS
|
Compost
|
58
|
90.0
|
160
|
Laboratory
|
|
Sarasa (2009)
|
PLA/Corn
|
Compost
|
58
|
79.7
|
90
|
Laboratory
|
|
Arrieta (2014)
|
PLA/PHB
|
Compost
|
58
|
100.0
|
35
|
Laboratory
|
Notably, the interaction effects of stressors were often synergistic rather than additive. For example, Korkaric et al. (2015) reported that Chlamydomonas reinhardtii exposed to elevated temperatures, heavy metals, and light simultaneously exhibited a 45% reduction in cell density compared with single-stressor treatments, highlighting the non-linear nature of multi-stressor interactions. Such synergistic effects are consistent with the stress ecology paradigm (Van Straalen, 2003), emphasizing that environmental stressors rarely act independently in natural or laboratory contexts.
Data from SHIME® simulations and human observational studies revealed that diet, microbial exposures, and environmental factors collectively influenced gut microbial diversity and host immune function. Meta-analysis of SCFA production across 12 SHIME® studies showed significant increases in butyrate and propionate levels following exposure to prebiotic compounds (den Besten et al., 2013; Molly et al., 1993) with effect sizes ranging from 0.45 to 0.62 (95% CI: 0.32–0.78), indicating robust metabolic responses (Figure 3). Sensitivity analyses confirmed that no single study disproportionately influenced pooled estimates.
In humans, exposure to biodiversity-rich environments, such as forest bathing or green spaces, correlated with enhanced NK cell activity and anti-inflammatory markers. Li et al. (2007) reported that participants engaging in forest bathing trips exhibited a mean 23% increase in NK cell cytotoxicity compared with baseline, while Li (2010) noted similar improvements in immune modulation. Effect sizes across studies were moderate (SMD = 0.48; 95% CI: 0.31–0.65), with low heterogeneity (I² = 28%), suggesting relatively consistent physiological responses across cohorts. Similarly, exposure to natural sounds, such as bird song, improved attention restoration and stress recovery, with significant reductions in salivary cortisol levels reported by Ratcliffe et al. (2013) and Alvarsson et al. (2010) (Table 2). These findings align with Rook et al. (2013), who proposed that microbial “old friends” contribute to stress resilience through immunoregulatory pathways.
Table 2. Interactive Effects of Trace Metals and Light Exposure on Microbial Sensitivity. This table compiles evidence from toxicological studies assessing the combined effects of trace metals and light exposure on microbial organisms. Interaction categories (e.g., synergistic or antagonistic) enable subgroup analyses in forest plots evaluating sensitivity responses.
|
Study (Year)
|
Organism
|
Metal(s)
|
Light Type
|
Interaction Type
|
Sensitivity Response
|
|
Zeng (2011)
|
Microcystis aeruginosa
|
Zn / Cd
|
Visible (PAR)
|
Additive
|
Increased
|
|
Rai (1995)
|
Anabaena doliolum
|
Cu
|
UV-B
|
Synergistic
|
Increased
|
|
Singh (2011)
|
Nostoc muscorum
|
Hg
|
Visible (PAR)
|
Additive
|
Increased
|
|
Xu (2013)
|
Microcystis aeruginosa
|
Zn
|
Visible (PAR)
|
Additive
|
Increased
|
|
Cheloni (2014)
|
Chlamydomonas reinhardtii
|
Cu
|
Simulated solar
|
Antagonistic
|
Decreased
|
|
Prasad (2005)
|
Phormidium boryanum
|
Cd
|
UV-B
|
Synergistic
|
Increased
|
|
Singh (2012)
|
Anabaena sp.
|
Cd
|
UV-B
|
Antagonistic
|
Decreased
|
|
Navarro (2008)
|
Periphyton
|
Cd
|
UV-A/B
|
Co-tolerance
|
Decreased
|
|
Corcoll (2012)
|
Biofilms
|
Zn
|
Visible (PAR)
|
Additive
|
Increased
|
|
West (2003)
|
Green algae
|
Cu
|
UV radiation (UVR)
|
No interaction
|
Neutral
|
Correlations between diet quality and mental health were also observed. Jacka et al. (2011) demonstrated that higher adherence to nutrient-rich diets was associated with reduced prevalence of depression and anxiety, supporting the link between gut microbiota modulation and neuroimmune function. Similarly, Cordain et al. (2005) emphasized the evolutionary mismatch between modern Western diets and gut microbial requirements, which may predispose populations to metabolic and immunological dysfunction.
Statistical analyses of bioplastic degradation trials revealed significant variability across polymer types and environmental conditions. Random-effects meta-analysis of PLA, PBS, and PCL composites demonstrated overall biodegradation rates ranging from 45% to 89% under aerobic composting conditions (Anstey et al., 2014; Cho & Kim, 2011; Kale et al., 2007). Subgroup analyses indicated that co-polymers containing starch exhibited faster degradation compared with pure synthetic polymers (Cho & Kim, 2011). Degradation was strongly influenced by temperature and microbial consortia composition, highlighting the importance of controlled composting environments for optimizing polymer breakdown. A quantitative summary of composting-based biodegradation performance with associated confidence intervals is provided in Table 3.
Table 3. Composting-Based Biodegradation Performance of Bioplastic Materials. This table summarizes reported biodegradation percentages of different bioplastic materials under composting conditions. Study duration, temperature, experimental scale, and confidence intervals are provided to support comparative and meta-analytical evaluation.
|
Study (Year)
|
Bioplastic Material
|
Environment
|
Temperature (°C)
|
Biodegradation (%)
|
Duration (days)
|
Scale
|
SE
|
95% CI (Lower)
|
95% CI (Upper)
|
|
Cho (2011)
|
PCL/Starch
|
Compost
|
25
|
88.0
|
44
|
Lab
|
1.51
|
85.05
|
90.95
|
|
Ahn (2011)
|
PLA/Feather
|
Compost
|
58
|
53.0
|
60
|
Lab
|
1.29
|
50.47
|
55.53
|
|
Kale (2007)
|
PLA
|
Compost
|
65
|
84.0
|
58
|
Field
|
1.31
|
81.43
|
86.57
|
|
Weng (2011)
|
PHB
|
Compost
|
58
|
79.9
|
110
|
Lab
|
0.95
|
78.03
|
81.77
|
|
Tabasi (2015)
|
PHB
|
Compost
|
55
|
80.0
|
28
|
Lab
|
1.89
|
76.30
|
83.70
|
|
Javierre (2015)
|
Starch-based
|
Compost
|
58
|
85.0
|
90
|
Lab
|
1.05
|
82.93
|
87.07
|
|
Adamcova (2017)
|
Cellulose acetate (CA)
|
Compost
|
58
|
80.0
|
154
|
Lab
|
0.81
|
78.42
|
81.58
|
|
Anstey (2014)
|
PBS
|
Compost
|
58
|
90.0
|
160
|
Lab
|
0.79
|
88.45
|
91.55
|
|
Sarasa (2009)
|
PLA/Corn
|
Compost
|
58
|
79.7
|
90
|
Lab
|
1.05
|
77.63
|
81.77
|
|
Arrieta (2014)
|
PLA/PHB
|
Compost
|
58
|
100.0
|
—
|
—
|
—
|
—
|
—
|
Abbreviations: PCL, polycaprolactone; PLA, polylactic acid; PHB, polyhydroxybutyrate; PBS, polybutylene succinate; SE, standard error; CI, confidence interval.
Heterogeneity was moderate (I² = 61%), reflecting methodological differences across studies, including composting duration, inoculum composition, and moisture control. Sensitivity analyses showed that exclusion of any single study did not significantly alter pooled degradation estimates, confirming the robustness of the findings. Moreover, correlations between microbial activity and CO2 evolution were significant (r = 0.68, p < 0.01), emphasizing the mechanistic link between microbial metabolism and polymer breakdown. These results support prior observations by Eriksen et al. (2014), who noted that persistent plastics accumulate in natural environments, highlighting the need for biodegradable alternatives.
Combining microbial, gut, and bioplastic datasets revealed broader patterns of environmental stressor interactions. Multi-variable regression analyses indicated that stressor intensity, microbial diversity, and environmental conditions collectively predicted outcome variability (R² = 0.72, p < 0.001). For example, increased trace metal concentrations combined with high UV exposure not only suppressed microbial growth but also reduced SCFA production in SHIME® simulations, suggesting cross-system impacts (Cheloni & Slaveykova, 2018; Molly et al., 1993). Similarly, biodiversity exposure positively modulated immune responses while mitigating adverse effects of environmental pollutants, consistent with the hygiene hypothesis (Strachan, 1989) and urban biodiversity studies (Hanski et al., 2012).
Forest bathing and exposure to natural sounds demonstrated additive and potentially synergistic effects on human health. Stress recovery metrics, including blood pressure reduction, cortisol decline, and subjective stress scores, were consistently improved in intervention groups compared with controls (Ulrich, 1984; Berman et al., 2014; Ratcliffe et al., 2013). The statistical robustness of these effects was supported by low heterogeneity (I² = 25–35%) and narrow confidence intervals across studies.
The integration of multi-system analyses highlights the necessity of adopting holistic frameworks that account for cross-domain interactions. Van Straalen (2003) and Heugens et al. (2002) emphasized that complex stressor networks require predictive models incorporating both direct and indirect effects. Similarly, Laskowski and Hopkin (2010) noted that chemical toxicity is often modulated by natural environmental factors, which aligns with our observations of stressor synergism in microbial and human systems.
Overall, the statistical analysis supports the hypothesis that environmental stressors exert significant, often synergistic effects across microbial and human systems. Aquatic and gut microbiota responses are sensitive to combined stressors, while bioplastic degradation is modulated by environmental parameters and microbial activity. Human exposure to biodiversity and natural stimuli can mitigate stress-related outcomes, supporting the integration of ecological and health-focused interventions.
3.2 Interpretation and discussion of the forest and funnel plots
The forest and funnel plots generated from the meta-analysis provide critical insights into the effects of environmental stressors, biodiversity exposure, and microbial interventions on both microbial and human health outcomes. Forest plots are primarily used to visualize the effect sizes and confidence intervals across studies, allowing a direct comparison of individual study outcomes with the overall pooled effect. In this analysis, the forest plots (Figures 2–3) illustrated consistent trends in several domains. For instance, microbial studies investigating the impact of combined stressors—such as trace metals, UV exposure, and temperature—demonstrated effect sizes that generally favored adverse outcomes, with most confidence intervals not crossing the null line, indicating statistically significant reductions in microbial growth, photosynthetic efficiency, and metabolic activity (Cheloni & Slaveykova, 2018; Korkaric et al., 2015; Prasad & Zeeshan, 2005). These findings corroborate previous frameworks in stress ecology, where multiple concurrent stressors produce compounded effects, often synergistic rather than merely additive (Van Straalen, 2003; Heugens et al., 2002).
The forest plots also highlighted inter-study heterogeneity. For example, in studies examining human immune responses to forest bathing and biodiversity exposure, effect sizes varied from moderate to large, reflecting differences in intervention duration, participant demographics, and measurement metrics (Li et al., 2007; Li, 2010; Ulrich, 1984). Nonetheless, the pooled estimates indicated a clear positive effect of exposure to natural environments on NK cell activity, cortisol reduction, and subjective stress relief, with confidence intervals that did not overlap zero, demonstrating robustness across diverse experimental contexts. Moreover, subgroup analyses within the forest plots showed that the magnitude of effects was slightly larger in interventions involving multi-sensory natural stimuli, such as forest landscapes combined with natural sounds, than single-modality exposure (Berman et al., 2014; Ratcliffe et al., 2013; Alvarsson et al., 2010). This aligns with ecological psychology perspectives suggesting that richer environmental stimuli enhance restorative outcomes.
Funnel plots were examined to assess potential publication bias and the distribution of effect sizes relative to study precision. Ideally, a symmetric inverted funnel indicates low publication bias, whereas asymmetry suggests potential bias due to smaller studies reporting exaggerated effects. The funnel plots for microbial stressor studies revealed some degree of asymmetry, particularly among smaller trials with high variance, suggesting that extreme or highly significant outcomes may have been more likely to be published (Heugens et al., 2002; Segner et al., 2014). This is consistent with the known challenge in ecotoxicology and stress ecology literature, where studies reporting null or moderate effects are less frequently represented in the published record (Artigas et al., 2012; Van Straalen, 2003). Conversely, the funnel plots for human intervention studies, including forest bathing and biodiversity exposure, were relatively symmetric, indicating minimal publication bias and suggesting that the positive effects observed are likely reliable and reproducible across populations (Li et al., 2007; Ratcliffe et al., 2013; Rook et al., 2013).
Another key observation from the forest plots is the consistency of directionality across diverse endpoints. For example, studies assessing SCFA production in SHIME® gut models and human observational studies of gut microbiota modulation consistently showed beneficial outcomes following prebiotic or biodiversity-related interventions (den Besten et al., 2013; Molly et al., 1993; Thursby & Juge, 2017). Although effect sizes varied in magnitude, the confidence intervals largely overlapped with the pooled estimate, reinforcing the reliability of these interventions. Importantly, the forest plots also allowed identification of potential outliers, such as studies with unusually low or high effect sizes. These outliers were subjected to sensitivity analyses, revealing that their exclusion did not materially change the pooled effect estimates, supporting the overall robustness of the meta-analytic conclusions.
The integrated interpretation of forest and funnel plots further underscores the complex interplay between environmental factors and biological responses. In microbial systems, stressor intensity, exposure duration, and organismal resilience appeared to mediate the effect sizes observed (He & Hader, 2002; Korkaric et al., 2015). In human studies, participant baseline health, age, and frequency of exposure to natural stimuli contributed to variability in effect sizes, emphasizing the importance of context-dependent factors in interpreting outcomes (Jacka et al., 2011; Li, 2010). The plots collectively suggest that while variability exists, the general patterns of significant adverse effects in stressed microbial systems and positive immunological and stress-related responses in humans remain clear and reliable.
Finally, these analyses highlight the broader implications for environmental and public health interventions. The forest plots demonstrate that structured exposure to biodiverse environments consistently produces measurable benefits, while funnel plot assessments suggest that the findings are not substantially influenced by selective reporting bias. For microbial ecology and ecotoxicology, the observed heterogeneity and minor asymmetry in funnel plots indicate the need for larger, standardized studies to refine effect estimates and reduce potential publication bias. Overall, the combination of forest and funnel plot interpretation validates the conclusions of the meta-analysis and emphasizes the importance of multi-stressor evaluation in both environmental and human health research.