- 1Department of Agricultural Biotechnology, Seoul National University, Seoul, Republic of Korea
- 2Department of Forest Environment Protection, Kangwon National University, Chuncheon, Republic of Korea
- 3National Institute of Biological Resources, Incheon, Republic of Korea
- 4Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea
Phasmatodea, commonly known as stick insects, are recognized as noteworthy pests globally, impacting agriculture and forest ecosystems. Among them, the outbreak of Ramulus mikado has emerged as a notable concern in East Asian forests. Recently, Metarhizium phasmatodeae has been identified as utilizing stick insects as hosts. We have observed evidence of this entomopathogenic fungus infecting stick insects. Given the increase in these occurrences during the rainy period, this study investigated the relationship between the survival of R. mikado and the M. phasmatodeae infection during the rainy seasons of 2022 and 2023. We collected stick insects in two representative forests of the Republic of Korea and examined insect survival, fungal infection, and various environmental factors. No infections were detected in specimens collected in June before the rainy season, but from July onwards, both the mortality of R. mikado and the fungal infection substantially increased. By the last sampling date of each year, 75% (2022), 71.4% (2023) of the specimens were infected, and over 90% of the total individuals succumbed as a result. Fungi isolated from deceased R. mikado were successfully identified as M. phasmatodeae using morphological and taxonomic approaches. Various statistical analyses, including principal component analysis and modeling, revealed a robust association between fungal infection and the survival of stick insects. The results highlight the correlation between mass deaths of stick insects and fungal infection, particularly during the summer rainy season. These findings offer valuable insights for forecasting R. mikado population in the upcoming year and developing effective pest control strategies.
1 Introduction
Stick insect, Ramulus mikado (Phasmatodea: Phasmatidae), is a wingless and parthenogenetic forest pest. R. mikado follows a univoltine reproductive cycle with obligatory embryonic diapause. Eggs are primarily laid in summer, and individuals overwinter as pharate first-instar nymphs (Yamaguchi and Nakamura, 2015). The stick insects pose a significant threat by causing serious defoliation in various deciduous trees. They exhibit polyphagous behavior by feeding on leaves. The long-term and repetitive defoliation caused by a high stick insect population can lead to premature seedling death and a reduction in the diameter growth of trees (Mazanec, 1967, 1968; Liu, 2021). In 2020, an outbreak of the stick insect in the urban forests near the metropolitan area in the Republic of Korea caused damage to 19 ha, escalating to 158 ha in 2021 and 981 ha in 2022 (Korean Forest Service, 2023). A comparable situation occurred in Japan, where R. mikado led an outbreak near residential areas over 2–3 years (Yano et al., 2021). Considering the economic and quality risks associated with the pest outbreak, various control methods, encompassing physical, chemical, and biological approaches, have been proposed (Ji et al., 2011; Klapwijk et al., 2016). Typically, synthetic pesticides such as fenitrothion and ethofenprox are applied to control R. mikado (Korean Forest Service, 2023). However, this method can affect non-target insects living in fields, leading to issues like bioaccumulation and the risk of side effects (Blacquière et al., 2012; Rahman et al., 2020; Yan et al., 2020).
Entomopathogenic fungi have emerged as highly promising alternative biocontrol agents to chemical insecticides (Erler and Ates, 2015; Perumal et al., 2023; Vivekanandhan et al., 2023; Perumal et al., 2024a,b). Metarhizium and Beauveria species are two main generalist entomopathogenic fungi widely employed in pest control across various insect species, including agricultural pests (Mascarin and Jaronski, 2016; Peng et al., 2021; Krutmuang et al., 2023; Vivekanandhan et al., 2024; Perumal et al., 2024a,b). Notably, the fungi of the genus Metarhizium have been recognized for their ability to infect and eliminate a broad range of arthropods (Zimmermann, 2007; Brunner-Mendoza et al., 2019; Vivekanandhan et al., 2022). This process involves the release of cuticle-degrading enzymes and the application of mechanical pressure exerted by the appressorium (Leger et al., 1986; Barelli et al., 2016). Upon successful penetration of the host insect, asexual spores (blastospores) are produced to facilitate dispersal in the insect hemocoel. Throughout the invasion of the insect body, depletion of nutrients in the hemocoel and the fat body occurs, eventually leading to the death of the host (Gao et al., 2011; Peng et al., 2022). In the case of the order Phasmatodea (Insecta), M. phasmatodeae has recently been reported as a novel species, utilizing stick insects as its hosts (Mongkolsamrit et al., 2020; Thanakitpipattana et al., 2020). In domestic cases, there have been reports of R. mikado being infected by M. anisopliae, and its biocontrol has been confirmed through artificial inoculation (Kim et al., 2010; Jung et al., 2023).
Before applying cultured entomopathogenic fungi to the field, it is necessary to conduct preliminary research on how these biocontrol agents react to environmental conditions during host infection. Their biocontrol activities may vary under actual field conditions compared to laboratory environments during inoculation (De La Rosa et al., 2000; Stafford and Allan, 2014). Particularly, the subtropical frontal zone of East Asia experiences heavy summer rainfall (Choi et al., 2020). The Korean rainy season, known as Changma, is characterized by a prominent increase in rainfall from late June to August, leading to elevated temperatures and relative humidity (Wang et al., 2007). Studies suggest that, entomopathogenic fungi, including the genus Metarhizium tend to enhance their virulence in response to increased temperature and humidity (Arthurs and Thomas, 2001; Goble et al., 2016). The genus Metarhizium is recognized for its presence in soil environments and its infection of insects inhabiting these areas (Mcguire and Northfield, 2020; Tamayo-Sánchez et al., 2022). However, research on Metarhizium species infections in stick insects, predominantly residing in the canopy layer throughout most of their life cycle, is limited. Taken together, we hypothesized that the survival of R. mikado could be affected by infection with entomopathogenic fungi, particularly with enhanced insecticidal performance during the rainy season.
The objectives of this study were 1) to assess the survival of R. mikado and the impact of M. phasmatodeae infection during the rainy season and 2) to explore the environmental factors affecting the biocontrol activity of M. phasmatodeae. Our investigation revealed that the survival of R. mikado is influenced by entomopathogenic fungi in the forest ecosystem, along with other environmental factors. We also observed a notable increase in the natural infection of M. phasmatodeae during the summer rainy season. These findings provide valuable insights for predicting the population dynamics of R. mikado in the upcoming years and for devising effective pest control strategies.
2 Materials and methods
2.1 Sample collection
Stick insects were identified morphologically by using taxonomic keys from National Institute of Biological Resources.1 Nymphs and adults of R. mikado were collected from Mt. Cheonggye (N 37° 24′ 21″ E 127° 00′ 02″) and Mt. Geumam (N 37° 30′ 42″ E 127° 10′ 52″) in Gyeonggi Province in 2022 and 2023, respectively. In each area, sampling was conducted by hand-picking and beating vegetation. The sampling was performed four times in 2022 and six times in 2023.
Approximately 40 to 50 individuals were collected per sampling, and 10 individuals were reared in cages (30 cm × 30 cm × 30 cm) at room temperature (RT) to minimize spatial stress. Specimens were chosen randomly and placed in separate cages for rearing. Cherry and oak branches with leaves (Prunus sargentii and Quercus mongolica) were supplied as food sources. Dead individuals were separated and individually transferred to petri dishes (90 × 15 mm) containing moistened filter paper to promote fungal sporulation. The dishes were incubated at 25 ± 1°C in darkness.
2.2 Data investigation
To evaluate survival and infection rates, dead individuals were counted every day and incubated to observe fungal conidia. Survival and infection rates were calculated as follows:
To obtain environmental information, data loggers, HOBO Pro v2 Temp / RH (Onset, Seoul, Republic of Korea) were mounted on trees for each sampling site and set to record temperature and relative humidity every hour. And then they were represented as average values using data from the 14 days preceding the sampling date. For precipitation, daily data obtained from the nearby Automatic Weather Station (AWS, https://data.kma.go.kr) at the sites were represented as cumulative precipitation data over the same period.
2.3 Fungal isolation
The strains used in this study were isolated from deceased R. mikado specimens were collected in the forest. We collected stick insects ranging from the 5th instar nymphs to adults. We randomly selected insects per sampling site and cultured fungi on potato dextrose agar (PDA: Difco, Detroit, United States) plates (90 × 15 mm) containing 100 μg/mL of penicillin and 100 μg/mL of streptomycin. The plates were sealed with parafilm and incubated at 25 ± 1°C in darkness. The emergence of mycelia was monitored daily, and the fungus was reisolated by single spore isolation to obtain a pure fungal culture. All pure isolates were stored in 20% glycerol at −80°C until use.
2.4 Morphological observation
To assess macro- and micro-morphological characters, each fungal isolate was grown on PDA at 25 ± 1°C for 14 days. The fungal spores were also obtained from PDA plates cultured for 14 days. Specimens were observed and photographed at 400 × magnification in a DM6 B microscope (Leica Microsystems, Wetzlar, Germany) equipped with a Leica DMC6200 camera.
2.5 Molecular identification and phylogenetic analyses
Fungal genomic DNA was extracted following the standard protocol (Specht et al., 1982) with minor modifications. Briefly, a 7-day-old fungal mass was placed in a 1.5 mL tube containing a mini plastic pestle after adding 450 μL of grinding buffer (20 mM Na4EDTA; 0.1 M Tris–HCl, pH 7.5; 1.4 M NaCl). Ground mycelia were mixed by inversion and centrifuged at 13,000 rpm for 10 min. The supernatant was transferred to a new tube, and the same volumes of isopropanol were added. Nucleic acids were precipitated by centrifugation at 13,000 rpm for 10 min at 4°C, dried for 15 min, and resuspended in 20 μL of distilled water. The extracted DNA was preserved at −20°C and used as a template for PCR amplifications. As previously described (White et al., 1990), PCR amplifications were carried out in 50 μL reactions using Taq DNA polymerase (Takara Shuzo, Japan), 20 pmol of each primer, and roughly 50 ng of template DNA for four loci of all strains; the internal transcribed spacer (ITS) region of 18 S-28 S nuclear ribosomal DNA, and 5′ intron-rich region translation elongation factor 1-α (5’tef). PCR primer pairs used to amplify the gene regions for this study were: ITS1 (GARTGYCCDGGDCAYTTYGG)/ITS4 (CCNGCDATNTCRTTRTCCATRTA) for ITS (White et al., 1990), EF1T (ATGGGTAAGGARGACAAGAC)/EF2T (GGAAGTACCAGTGATCATGTT) for 5’tef (Bischoff et al., 2006). Thermocycler conditions for amplification of the DNA regions followed previously described protocols (Rehner and Buckley, 2005). Amplified PCR products were purified using the MEGAquick-spinTM plus Fragment DNA Purification kit (Intron, Seongnam, Republic of Korea) and sequenced (Bioneer, Daejeon, Republic of Korea).
Analysis of DNA sequences was performed with Seqman Pro (DNAStar, Madison, WI, United States) to assemble and edit forward and reverse sequences. These sequences were aligned using Clustal W by MEGA-X (Kumar et al., 2018). Sequences of ITS and 5’tef from related species were retrieved from the NCBI Genebank to elucidate relationships in Metarhizium (Mongkolsamrit et al., 2020). Maximum likelihood (ML) analyses were performed with the general time reversible model with Invariant sites and Gamma distribution (GTR + I + G) with 1,000 bootstrap replicates using RAxML (Edler et al., 2021). This model was selected due to its comprehensive approach, accounting for base substitution rates, invariant sites, and rate variations among sites (Abadi et al., 2019). Bayesian posterior probabilities (BPP) were performed using Markov Chain Monte Carlo sampling (MCMC) in MrBayes v. 3.2.7 (Ronquist et al., 2012). Six simultaneous Markov Chains were run for 100,000 generations, and trees were sampled every 1,000th generation. The phylogenetic trees were visualized in Fig Tree v. 1.4.4 (Rambaut, 2009). The NCBI accession numbers of the three strains in this study are listed in Supplementary Table S2.
2.6 Insect bioassay
To verify the pathogenicity of M. phasmatodeae against R. mikado, adult R. mikado specimens were used in the bioassay as described previously (Kim et al., 2020). Three fungal isolates were randomly selected and cultured on PDA plates (90 × 15 mm) at 25 ± 1°C for 14 days. To induce exposure to conidia, individual adult R. mikado specimen was placed on the fungal culture. After 1 h, the infected adult was transferred to a plastic cup (20 mm × 90 mm × 81.6 mm) along with a fresh cherry branch (Prunus sargentii) as a food source. The stick insects were reared in the cups at 25 ± 1°C for 14 days, with one milliliter of sterile distilled water supplied every 2 days to maintain high humidity. The number of deceased adults was recorded daily. Once stick insects died, they were moved to petri dishes (90 × 15 mm) with moistened filter paper to assess infection by M. phasmatodeae. This bioassay was repeated three times.
2.7 Statistical analysis
The survival of R. mikado post-sampling was visualized using the Kaplan–Meier curve (Kaplan and Meier, 1958). Any insects surviving beyond this time were regarded as “censored.” Log-rank (Mantel-Cox) pair-wise comparisons were used to evaluate the difference in survivorship based on the sampling dates.
Additionally, principal component analysis (PCA) was used to determine the multiple correlations based on the sampling dates (Wold et al., 1987). In the PCA biplot, survival and infection rates, and environmental factors were represented by vectors, and sample IDs were represented by symbols. Angles between vectors reflect their correlation, which can be estimated as the cosine of the angle between any two vectors; vectors pointing in the same direction have a positive correlation, those pointing in opposite directions are negatively correlated, and vectors at right angles mean that vectors are not correlated.
The relationships between survival rate and each of the variables (infection rate, temperature, relative humidity, and precipitation) were investigated. In this study, nonlinear or linear models that accurately describe these relationships were fitted using the R packages nls (nonlinear model) or lm (linear model). The performance of each model was estimated using the coefficient of determination (R2) based on the observed value and the predicted value from the model. All data analyses were performed using the software R 4.3.2 and SPSS 25.
3 Results
3.1 Survey on the survival of R. mikado in the summers of 2022 and 2023
Diseased stick insects exhibiting paralysis-like symptoms and green spore formation on their body surfaces were observed in some mountains of the Republic of Korea (Figure 1A). Some of these infected stick insects were found deceased, with a notable prevalence during the summer season. Therefore, we aimed to monitor the survival of these insects over time and uncover potential causes of mortality during specific periods.
Figure 1. Disease symptoms of Ramulus mikado. (A) R. mikado infested on the body surface presenting a symptom of paralysis. Green spores (red arrows) were formed on the carcass of R. mikado. (B) Collected R. mikado samples. (C) Live (left) and dead (right) adults of R. mikado were observed. (D) Fungal spore masses covered the carcasses of R. mikado after incubation in humid conditions.
To assess the survival of R. mikado, we collected a total of 461 individuals through four samplings in Mt. Cheonggye (CM) and six samplings in Mt. Geumam (GM) in the summers of 2022 and 2023, respectively (Figures 1B,C; Table 1). The collected individuals underwent a two-week rearing process, revealing that dead insects were covered by fungal spore masses (Figure 1D). In the CM group, the survival rates for CM 3 and CM 4 were 2.7% and 7.5%, respectively, by day 14. In instances where spore masses were detected in dead insects, we classified them as “infected” and investigated infection rates among the deceased. Notably, infection rates exceeded 50% for both the CM and GM groups. In the GM group, survival rates for GM 5 and GM 6 were less than 50%, with fungal infection initiating at GM 3 and gradually escalated, reaching a peak infection rate of 70% by GM 6.
A comparison of survival rates based on the sampling dates (CM 1–4 and GM 1–6 groups) indicated a significant decrease in survival rates as the sampling date progressed for both CM and GM groups (Figure 2). Log-rank pair-wise comparisons of survival curves were conducted to assess the statistical significance, providing the chi-square values (X2) and significance levels (p) (Supplementary Table S3). A statistically significant difference in survival rates was observed depending on the sampling period. As an example, in both 2022 and 2023, the final sampling group (CM 4, GM 6) exhibited a rapid mortality rate compared to the first sampling group (CM 1, GM 1) (X2(CM) = 48.547, p < 0.001; X2(GM) = 81.694, p < 0.001). Additionally, the mortality due to fungal infection varied based on the sampling period, with over 70% of all specimens being infected in the final sampling group over the course of 2 years (Table 2). Collectively, our findings suggest a significant decline in the survival rates of R. mikado as the sampling period progresses. This prompts further exploration into whether the mortality of these insects is linked to fungal infection.
Figure 2. Survival curves according to Kaplan–Meier analysis of R. mikado. Survival rates were analyzed for R. mikado samples which were collected in (A) 2022 and (B) 2023. Mortality was examined daily by counting the number of dead insects for 14 days. Asterisk indicates a significant difference in survival among the groups (p < 0.001, log-rank test).
3.2 Fungal species identification
As described above, a significant number of fungus-infested insects were collected during the summer season. Notably, these insects exhibited green spore masses covering their bodies, a distinctive trait indicative of Metarhizium species infection (Thanakitpipattana et al., 2020). Thus, we isolated fungal strains from the carcasses of R. mikado to clearly identify the Metarhizium strains infecting R. mikado. Diseased stick insects were placed on PDA to isolate pure fungal cultures, and three putative Metarhizium strains were successfully obtained. Then we analyzed the phenotypic characteristics of vegetative growth and asexual reproduction. Colonies of all strains showed yellow with a white edge and robust sporulation when grown on PDA (Figure 3A). Conidia produced on PDA were single-celled and cylindrical shape (Figure 3B). Furthermore, a multigene alignment was conducted using DNA barcode genes, including the internal transcribed spacer (ITS) and 5′ intron-rich region of the translation elongation factor 1-alpha (5’tef). All three strains exhibited high similarity to the barcode gene sequence of the M. phasmatodeae strain BCC49272, with ITS showing 99.56–99.71% similarity and 5’tef showing 98.54–98.78%. In addition, phylogenetic analyses employing both maximum likelihood (ML) and Bayesian methods further confirmed a close relationship between the strains in this study and the M. phasmatodeae strain (BCC49272). As previously reported, M. phasmatodeae is a recently identified species known to infect stick insects (Phasmatodea) (Mongkolsamrit et al., 2020; Thanakitpipattana et al., 2020). In summary, based on comprehensive morphological and taxonomic characterization, the three strains were conclusively identified as belonging to M. phasmatodeae, suggesting their potential as biocontrol agents against Phasmatodea insects (Figure 3C).
Figure 3. Morphological characterization and phylogenetic reconstruction of isolates in this study. (A) Colonies on PDA after 14 days. (B) Conidia on PDA. Scale bar =10 μm. (C) Phylogenetic reconstruction of M. anisopliae species complex based on ITS and 5′tef sequences using ML and Bayesian analysis. Number on the nodes are ML bootstrap / Bayesian posterior probability values above 70% (MLBS) or 0.7 (BPP).
3.3 Insect bioassay
To assess the pathogenicity of isolated M. phasmatodeae against R. mikado, healthy R. mikado adults were directly exposed to the 14-day-old fungal culture for 1 h and subsequently transferred to a clean dish (Figure 4A). Stick insects exposed to a clean PDA culture served as the control group. Following the treatments, all individuals in the experimental group succumbed within 6 days, in contrast to the control group, which survived throughout the experimental period. The external emergence of mycelia from the carcasses and the formation of green conidia were observed from the third day onward (Figure 4B). Taken together, our findings affirm that the M. phasmatodeae strains examined in this study exhibit biocontrol activity against R. mikado.
Figure 4. Biocontrol effects of M. phasmatodeae on R. mikado. (A) Each insect sample was exposed to the culture of M. phasmatodeae. Pure PDA culture was used as a control. The treated insect was transferred to a plastic dish and checked every day until 100% mortality was achieved. (B) Representative images of R. mikado after treatment with M. phasmatodeae. Dead insects were removed from the dishes and incubated for an additional indicated period (1 or 3 days) under moist conditions until mycelial outgrowth was observed.
3.4 Association between R. mikado survival and environmental variables
Given the observed mortality of R. mikado during the summer season, our inquiry focused on understanding the factors contributing to the heightened virulence of M. phasmatodeae in forest ecosystems. As described earlier, an increase in rainfall was noted during the summer seasons of both 2022 and 2023, resulting in hot and humid conditions. We collected data on average temperature and relative humidity using data loggers at the sampling sites, and accumulated precipitation data were obtained from AWS (Table 2). Subsequently, each environmental factor was compared with survival and infection rates (Figure 5).
Figure 5. Investigation of the relationship between environmental factors and M. phasmatodeae infection against R. mikado. The values in the bar graph indicate various environmental factors in (A) 2022 and (B) 2023. Solid and dashed lines indicate survival and infection rates, respectively. Average temperature and Average relative humidity are represented as average values for 2 weeks derived from HOBO data loggers at each sampling site. Exceptionally, temperature and relative humidity values of GM1 were investigated using Automatic Weather Station (AWS) data. The values of precipitation represent the total accumulation for 2 weeks derived from the nearby AWS.
The analysis revealed a decline in survival rates with increasing temperature and humidity, accompanied by a rapid increase in infection rates. This implies an inverse relationship between the infection rate and survival rate, indicating a positive association between the infection rate and environmental factors. Notably, the high infection rate on the last sampling dates in both 2022 and 2023 was observed during periods of low precipitation. In this period, the warm and humid conditions following previous rainfall during the rainy season may have led to elevated infection rates, even as precipitation declined. In summary, our findings suggest fluctuations in the survival of R. mikado were influenced by the rainy season, which in turn affected temperature and humidity conditions. Our observation aligns with previous studies indicating that the insecticidal performance of entomopathogenic fungi could be influenced by environmental conditions such as humidity and temperature (De La Rosa et al., 2000; Arthurs and Thomas, 2001; Athanassiou et al., 2017). Therefore, we aimed to analyze their correlation through further statistical analysis.
3.5 Validation of the R. mikado survival model
To validate the R. mikado survival model, principal component analysis (PCA) was employed to understand relationships among various environmental factors (Table 2; Supplementary Table S4). This analysis aimed to identify correlated variables and ascertain which factors contributed most significantly to the variation in stick insect survival (Figure 6). The cumulative proportion of explanation for the first three dimensions (PC1, PC2, and PC3) was 75.0, 13.9, and 10.0%, respectively, with a total of 98.9% of the variances explaining the model. PC1 and PC2 with an eigen values of 3.75 and 0.69 revealed a notable negative correlation between R. mikado survival and M. phasmatodeae infection, while precipitation exhibited a somewhat weaker correlation with the survival rate and infection rate (Figure 6A). Conversely, PC1 and PC3 with an eigen values of 3.75 and 0.50 indicated a positive correlation between precipitation and the infection rate, with both factors negatively correlated with the survival rate (Figure 6B). To summarize, the survival rate displayed a significant negative correlation with infection rate, temperature, relative humidity, and precipitation. However, using the Spearman correlation coefficient, no statistically significant correlation was observed between precipitation and the infection rate, as well as precipitation and temperature (p > 0.05) (Supplementary Figure S1).
Figure 6. Principal component analysis (PCA) on the survival and infection rates of R. mikado, relative humidity, and temperature. Each biplot represents the combination of (A) PC1 and PC2, as well as (B) PC1 and PC3. The percentages of variance explained by axes or components (PC1-3) are shown in parentheses. Vectors indicate significant correlations between survival rate, infection rate, and environmental factors.
Further, considering the close relationship observed between insect survival and each variable (infection rate, temperature, relative humidity, and precipitation), we investigated whether these associations could be proposed as a host-pathogen interaction model. A significant curvilinear or linear relationship was identified between survival rate and infection rate (R2 = 0.95), temperature (R2 = 0.97), relative humidity (R2 = 0.68), and precipitation (R2 = 0.67) (Figure 7). This relationship exhibited a rapid decline in survival rate with low infection rates, followed by a gradual decrease as infection rates increased (Figure 7A). Temperatures ranging from 14°C to 23°C had minimal impact on the survival rate, while beyond 23°C, the survival rate decreased to less than 0.5 (Figure 7B). Additionally, as relative humidity increased, there was a progressively sharp decline in the survival rate, dropping below 50% at humidity levels exceeding approximately 85% (Figure 7C). The relationship between precipitation and survival rate showed a negative linear correlation (Figure 7D). Overall, these results suggest that the survival of R. mikado exhibits a nonlinear/linear relationship with various factors, including fungal infection, temperature, relative humidity, and precipitation.
Figure 7. Modeling between survival rate and other factors. (A) infection rate; (B) temperature; (C) relative humidity; (D) precipitation. Red line shows the relationship and the equation of each red line is described beside the red line. The performance of each model was examined using the coefficient of determination (R2). Gray dot (highlighted using red dash line) excluded when creating the model.
4 Discussion
The order Phasmatodea, encompassing stick insects, is known for harboring significant phytophagous pests affecting agriculture and forestry in various countries (Campbell, 1960; Prathapan et al., 2008). In the Republic of Korea, five species of insects, namely R. mikado, R. koreanus, Micadina phluctainoides, M. yasumatsui, and Phraortes illepidus, belong to Phasmatodea. Among these, R. mikado is recognized as a forest insect pest that causes severe defoliation by feeding on leaves from diverse deciduous trees (Jung et al., 2020).
To manage stick insects, a common approach involves the use of commercial sticky traps, which are encased around the tree trunks to capture hatched young nymphs climbing up from the ground (Korean Forest Service, 2023). Additionally, chemical control methods are employed, utilizing synthetic pesticides such as fenitrothion and ethofenprox. However, the application of sticky traps comes with inherent challenges, requiring individual installation in the forest and incurring substantial labor costs. Similar to chemical control methods, sticky traps also pose the disadvantage of impacting non-target organisms (Blacquière et al., 2012; Rahman et al., 2020). Moreover, the efficacy of chemical control may diminish over time due to the development of resistance in target insects (Hemingway and Ranson, 2000).
In this study, M. phasmatodeae, first reported in Korea and previously documented in Thailand and China, has demonstrated pathogenicity for stick insects as hosts (Thanakitpipattana et al., 2020; Zhao et al., 2023), suggesting its potential as a promising biocontrol agent. Our observation revealed that all insects inoculated with M. phasmatodeae succumbed within 5 days. Considering the rapid onset, M. phasmatodeae holds promise as a biocontrol agent against stick insect outbreaks in forest ecosystems (Sabbahi et al., 2022; Vivekanandhan et al., 2022; Perumal et al., 2023). However, for application in forest ecosystems, additional host specificity tests are necessary. To enhance the effectiveness of fungal infection, it is essential to develop spray techniques that consider the ideal timing in accordance with environmental conditions.
We established the correlations among environmental factors, insect survival, and fungal infection through various statistical analyses. PCA indicated a minimal correlation between precipitation and infection rate in PC1 and PC2, while a positive correlation was identified between precipitation and infection rates in PC1 and PC3. This disparity can be attributed to the characteristics of the rainy season, where precipitation undergoes rapid fluctuations, experiencing sudden increases followed by subsequent decreases (Wang et al., 2007). These changes, driven by the hot and humid conditions during the rainy season, persist even with reduced rainfall, enhancing the insecticidal performance of M. phasmatodeae. Consequently, the infection rate peaked on the last sampling date for each year, despite decreased precipitation, influencing the outcomes of PCA. Furthermore, modeling was conducted to examine the survival rates and various other factors. According to the infection rate-survival rate formula, a 16% or higher infection rate results in a 50% mortality rate among stick insects. Therefore, even a relatively low infection rate could induce a significant number of deaths in the stick insect.
The extensive mortality observed in R. mikado may be influenced by the population density of stick insects. A similar study conducted on locust populations reported a low infection rate (Shah et al., 1998), suggesting that population density could play a role in this disparity. The sampling sites in this study were forests experiencing significant stick insect outbreaks, resulting in much higher population densities than usual. The heightened density of hosts facilitates the acceleration of horizontal transmission among individuals (Roy and Pell, 2000). Additionally, the collective deaths of the stick insects could be related to the reproduction strategy. The formation of colonies through asexual parthenogenesis may indicate relatively equal susceptibility to a fungus (Steinkraus, 2006).
In conclusion, this study confirmed the natural infection of R. mikado by M. phasmatodeae and observed the mass mortality phenomenon in the stick insects due to the increased virulence of M. phasmatodeae in the hot and humid environment formed during the rainy season. Our findings will contribute to predicting population levels for the next year based on the rainy season. The adult female stick insects normally lay eggs per day in the summer (Kobayashi et al., 2014; Lee et al., 2018). This trait implies that the total number of eggs laid may vary depending on the onset of the rainy season, consequently impacting population levels in the subsequent year. The fecundity of insects is one of the crucial factors in predicting the size of the next population (Peters and Barbosa, 1977; Powell and Bentz, 2009). Forecasting the future population could potentially assist in formulating pest control strategies. To ascertain whether this phenomenon is exclusive to R. mikado or common among stick insects within the Phasmatodea order, further research will be necessary.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
DM: Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. SS: Formal analysis, Investigation, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. N-HL: Data curation, Methodology, Software, Writing – original draft. MB: Funding acquisition, Resources, Writing – original draft. S-JP: Funding acquisition, Resources, Writing – original draft. K-HK: Data curation, Methodology, Software, Writing – original draft. HS: Project administration, Supervision, Writing – original draft, Writing – review & editing. J-KJ: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by a grant from the National Institute of Biological Resources (NIBR), funded by the Ministry of Environment (MOE) of the Republic of Korea (NIBR202405104).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1383055/full#supplementary-material
Footnotes
References
Abadi, S., Azouri, D., Pupko, T., and Mayrose, I. (2019). Model selection may not be a mandatory step for phylogeny reconstruction. Nat. Commun. 10:934. doi: 10.1038/s41467-019-08822-w
Arthurs, S., and Thomas, M. B. (2001). Effects of temperature and relative humidity on sporulation of Metarhizium anisopliae var. acridum in mycosed cadavers of Schistocerca gregaria. J. Invertebr. Pathol. 78, 59–65. doi: 10.1006/jipa.2001.5050
Athanassiou, C. G., Kavallieratos, N. G., Rumbos, C. I., and Kontodimas, D. C. (2017). Influence of temperature and relative humidity on the insecticidal efficacy of Metarhizium anisopliae against larvae of Ephestia kuehniella (Lepidoptera: Pyralidae) on wheat. J. Insect Sci. 17:22. doi: 10.1093/jisesa/iew107
Barelli, L., Moonjely, S., Behie, S. W., and Bidochka, M. J. (2016). Fungi with multifunctional lifestyles: endophytic insect pathogenic fungi. Plant Mol. Biol. 90, 657–664. doi: 10.1007/s11103-015-0413-z
Bischoff, J. F., Rehner, S. A., and Humber, R. A. (2006). Metarhizium frigidum sp. nov.: a cryptic species of M. Anisopliae and a member of the M. Flavoviride complex. Mycologia 98, 737–745. doi: 10.1080/15572536.2006.11832645
Blacquière, T., Smagghe, G., Van Gestel, C. A. M., and Mommaerts, V. (2012). Neonicotinoids in bees: a review on concentrations, side-effects and risk assessment. Ecotoxicology 21, 973–992. doi: 10.1007/s10646-012-0863-x
Brunner-Mendoza, C., Reyes-Montes, M. D. R., Moonjely, S., Bidochka, M. J., and Toriello, C. (2019). A review on the genus Metarhizium as an entomopathogenic microbial biocontrol agent with emphasis on its use and utility in Mexico. Biocontrol Sci. Tech. 29, 83–102. doi: 10.1080/09583157.2018.1531111
Campbell, K. (1960). Preliminary studies in population estimation of two species of stick insects (Phasmatidae Phasmatodea) occurring in plague numbers in highland forest areas of South-Eastern Australia, Proceedings of the Linnean Society of new South Wales Volume. (Sydney proc. Linn. Soc. N. S. W.).
Choi, J. W., Kim, H. D., and Wang, B. (2020). Interdecadal variation of Changma (Korean summer monsoon rainy season) retreat date in Korea. Int. J. Climatol. 40, 1348–1360. doi: 10.1002/joc.6272
De La Rosa, W., Alatorre, R., Barrera, J. F., and Toriello, C. (2000). Effect of Beauveria bassiana and Metarhizium anisopliae (Deuteromycetes) upon the coffee berry borer (Coleoptera: Scolytidae) under field conditions. J. Econ. Entomol. 93, 1409–1414. doi: 10.1603/0022-0493-93.5.1409
Edler, D., Klein, J., Antonelli, A., and Silvestro, D. (2021). raxmlGUI 2.0: a graphical interface and toolkit for phylogenetic analyses using RAxML. Methods Ecol. Evol. 12, 373–377. doi: 10.1111/2041-210X.13512
Erler, F., and Ates, A. O. (2015). Potential of two entomopathogenic fungi, Beauveria bassiana and Metarhizium anisopliae (Coleoptera: Scarabaeidae), as biological control agents against the June beetle. J. Insect Sci. 15:44. doi: 10.1093/jisesa/iev029
Gao, Q., Jin, K., Ying, S.-H., Zhang, Y., Xiao, G., Shang, Y., et al. (2011). Genome sequencing and comparative transcriptomics of the model entomopathogenic fungi Metarhizium anisopliae and M. Acridum. PLoS Genet. 7:e1001264. doi: 10.1371/journal.pgen.1001264
Goble, T. A., Gardescu, S., Fisher, J. J., Jackson, M. A., and Hajek, A. E. (2016). Conidial production, persistence and pathogenicity of hydromulch formulations of Metarhizium brunneum F52 microsclerotia under forest conditions. Biol. Control 95, 83–93. doi: 10.1016/j.biocontrol.2016.01.003
Hemingway, J., and Ranson, H. (2000). Insecticide resistance in insect vectors of human disease. Annu. Rev. Entomol. 45, 371–391. doi: 10.1146/annurev.ento.45.1.371
Ji, L., Wang, Z., Wang, X., and An, L. (2011). Forest insect pest management and forest management in China: An overview. Environ. Manag. 48, 1107–1121. doi: 10.1007/s00267-011-9697-1
Jung, J.-K., Jung, B.-N., Lee, C. Y., Kim, K. E., Kim, J., Lee, Y. S., et al. (2023). Evaluation of insecticidal efficacy of six eco-friendly agricultural materials and Metarhizium anisopliae against Ramulus mikado. J. Korean Soc. For. Sci. 112, 117–125. doi: 10.14578/jkfs.2023.112.1.117
Jung, J.-K., Nam, Y., Kim, D., Lee, S.-H., Lim, J.-H., Choi, W. I., et al. (2020). Tree-crown defoliation caused by outbreak of forest insect pests in Korea during 2020. Korean J. Appl. Entomol. 59, 409–410. doi: 10.5656/KSAE.2020.10.0.054
Kaplan, E. L., and Meier, P. (1958). Nonparametric estimation from incomplete observations. J. Am. Stat. Assoc. 53, 457–481. doi: 10.1080/01621459.1958.10501452
Kim, J. C., Lee, S. J., Kim, S., Lee, M. R., Baek, S., Park, S. E., et al. (2020). Management of pine wilt disease vectoring Monochamus alternatus adults using spray and soil application of Metarhizium anisopliae JEF isolates. J. Asia Pac. Entomol. 23, 224–233. doi: 10.1016/j.aspen.2019.12.012
Kim, W.-G., Seok, S.-J., Weon, H.-Y., Lee, K.-H., Lee, C.-J., and Kim, Y.-S. (2010). Isolation and identification of entomopathogenic fungi collected from mountains and islands in Korea. Kor. J. Mycol. 38, 99–104. doi: 10.4489/KJM.2010.38.2.099
Klapwijk, M. J., Bylund, H., Schroeder, M., and Björkman, C. (2016). Forest management and natural biocontrol of insect pests. Forestry 89, 253–262. doi: 10.1093/forestry/cpw019
Kobayashi, S., Usui, R., Nomoto, K., Ushirokita, M., Denda, T., and Izawa, M. (2014). Does egg dispersal occur via the ocean in the stick insect Megacrania tsudai (Phasmida: Phasmatidae)? Ecol. Res. 29, 1025–1032. doi: 10.1007/s11284-014-1188-4
Korean Forest Service (2023). Forecasting and control plans for forest insect pests and diseases in 2023. Available at: https://www.forest.go.kr
Krutmuang, P., Rajula, J., Pittarate, S., Chanbang, Y., Perumal, V., Alford, L., et al. (2023). Biocontrol efficacy of Beauveria bassiana in combination with tobacco short stem and modified lure traps. Int. J. Trop. Insect Sci. 43, 1591–1600. doi: 10.1007/s42690-023-01063-x
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Lee, J., Baek, S., Kang, C., Lee, Y. S., Lee, Y., and Lee, J.-H. (2018). Temperature-dependent development and oviposition models of Ramulus irregulariterdentatus (Phasmida: Phasmatidae). J. Asia Pac. Entomol. 21, 903–913. doi: 10.1016/j.aspen.2018.07.003
Leger, R. S., Charnley, A., and Cooper, R. (1986). Cuticle-degrading enzymes of entomopathogenic fungi: synthesis in culture on cuticle. J. Invertebr. Pathol. 48, 85–95. doi: 10.1016/0022-2011(86)90146-1
Liu, H. (2021). Biology and ecology of the northern walkingstick, Diapheromera femorata (say)(Phasmatodea: Diapheromerinae): A review. J. Appl. Entomol. 145, 635–647. doi: 10.1111/jen.12902
Mascarin, G. M., and Jaronski, S. T. (2016). The production and uses of Beauveria bassiana as a microbial insecticide. World J. Microbiol. Biotechnol. 32:177. doi: 10.1007/s11274-016-2131-3
Mazanec, Z. (1967). Mortality and diameter growth in mountain ash defoliated by Phasmatids. Aust. For. 31, 221–223. doi: 10.1080/00049158.1967.10675444
Mazanec, Z. (1968). Influence of defoliation by the Phasmatid Didymuria violescens on seasonal diameter growth and the pattern of growth rings in alpine ash. Aust. For. 32, 3–14. doi: 10.1080/00049158.1968.10675456
Mcguire, A. V., and Northfield, T. D. (2020). Tropical occurrence and agricultural importance of Beauveria bassiana and Metarhizium anisopliae. Front. Sustain. Food Syst. 4:6. doi: 10.3389/fsufs.2020.00006
Mongkolsamrit, S., Khonsanit, A., Thanakitpipattana, D., Tasanathai, K., Noisripoom, W., Lamlertthon, S., et al. (2020). Revisiting Metarhizium and the description of new species from Thailand. Stud. Mycol. 95, 171–251. doi: 10.1016/j.simyco.2020.04.001
Peng, Z.-Y., Huang, S.-T., Chen, J.-T., Li, N., Wei, Y., Nawaz, A., et al. (2022). An update of a green pesticide: Metarhizium anisopliae. All Life 15, 1141–1159. doi: 10.1080/26895293.2022.2147224
Peng, G., Xie, J., Guo, R., Keyhani, N. O., Zeng, D., Yang, P., et al. (2021). Long-term field evaluation and large-scale application of a Metarhizium anisopliae strain for controlling major rice pests. J. Pest. Sci. 94, 969–980. doi: 10.1007/s10340-020-01313-8
Perumal, V., Kannan, S., Alford, L., Pittarate, S., and Krutmuang, P. (2024a). Study on the virulence of Metarhizium anisopliae against Spodoptera frugiperda (JE smith, 1797). J. Basic Microbiol. doi: 10.1002/jobm.202300599
Perumal, V., Kannan, S., Alford, L., Pittarate, S., Mekchay, S., Reddy, G. V., et al. (2023). Biocontrol effect of entomopathogenic fungi Metarhizium anisopliae ethyl acetate-derived chemical molecules: an eco-friendly anti-malarial drug and insecticide. Arch. Insect Biochem. Physiol. 114, 1–19. doi: 10.1002/arch.22037
Perumal, V., Kannan, S., Pittarate, S., and Krutmuang, P. (2024b). A review of entomopathogenic fungi as a potential tool for mosquito vector control: A cost-effective and environmentally friendly approach. Entomol. Res. 54:e12717. doi: 10.1111/1748-5967.12717
Peters, T. M., and Barbosa, P. (1977). Influence of population density on size, fecundity, and developmental rate of insects in culture. Annu. Rev. Entomol. 22, 431–450. doi: 10.1146/annurev.en.22.010177.002243
Powell, J. A., and Bentz, B. J. (2009). Connecting phenological predictions with population growth rates for mountain pine beetle, an outbreak insect. Landsc. Ecol. 24, 657–672. doi: 10.1007/s10980-009-9340-1
Prathapan, K., Anith, K., Faizal, M., Lekha, M., and Dhanya, M. (2008). A report on Sipyloidea stigmata Redtenbacher (Diapheromeridae: Necrosciinae) as the first phasmid crop pest in India and its redescription. Zootaxa 1959, 58–64. doi: 10.11646/ZOOTAXA.1959.1.3
Rahman, M. S., Sumon, K. A., Uddin, M. J., and Shahjahan, M. (2020). Toxic effects of fenitrothion on freshwater microcosms in Bangladesh. Toxicol. Rep. 7, 1622–1628. doi: 10.1016/j.toxrep.2020.12.002
Rambaut, A. (2009). FigTree. Tree figure drawing tool. Available at: http://tree.bio.ed.ac.uk/software/figtree
Rehner, S. A., and Buckley, E. (2005). A Beauveria phylogeny inferred from nuclear ITS and EF1-α sequences: evidence for cryptic diversification and links to Cordyceps teleomorphs. Mycologia 97, 84–98. doi: 10.1080/15572536.2006.11832842
Ronquist, F., Teslenko, M., Van Der Mark, P., Ayres, D. L., Darling, A., Höhna, S., et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542. doi: 10.1093/sysbio/sys029
Roy, H. E., and Pell, J. K. (2000). Interactions between entomopathogenic fungi and other natural enemies: implications for biological control. Biocontrol Sci. Tech. 10, 737–752. doi: 10.1080/09583150020011708
Sabbahi, R., Hock, V., Azzaoui, K., Saoiabi, S., and Hammouti, B. (2022). A global perspective of entomopathogens as microbial biocontrol agents of insect pests. J. Agric. Food Res. 10:100376. doi: 10.1016/j.jafr.2022.100376
Shah, P. A., Gbongboui, C., Godonou, I., Hossou, A., and Lomer, C. J. (1998). Natural incidence of Metarhizium flavoviride infection in two grasshopper communities in northern Benin. Biocontrol Sci. Tech. 8, 335–344. doi: 10.1080/09583159830144
Specht, C. A., Dirusso, C. C., Novotny, C. P., and Ullrich, R. C. (1982). A method for extracting high-molecular-weight deoxyribonucleic acid from fungi. Anal. Biochem. 119, 158–163. doi: 10.1016/0003-2697(82)90680-7
Stafford, K. C., and Allan, S. A. (2014). Field applications of entomopathogenic fungi Beauveria bassiana and Metarhizium anisopliae F52 (Hypocreales: Clavicipitaceae) for the control of Ixodes scapularis (Acari: Ixodidae). J. Med. Entomol. 47, 1107–1115. doi: 10.1603/me10019
Steinkraus, D. C. (2006). Factors affecting transmission of fungal pathogens of aphids. J. Invertebr. Pathol. 92, 125–131. doi: 10.1016/j.jip.2006.03.009
Tamayo-Sánchez, F., Tamayo-Mejía, F., Marín-Jarrillo, A., Santillán-Galicia, M. T., Pérez-Valdez, A., and Guzmán-Franco, A. W. (2022). Dynamics of natural infection of white grub larvae by Beauveria and Metarhizium in maize crops from Mexico. Biocontrol Sci. Tech. 32, 1177–1193. doi: 10.1080/09583157.2022.2098925
Thanakitpipattana, D., Tasanathai, K., Mongkolsamrit, S., Khonsanit, A., Lamlertthon, S., and Luangsa-Ard, J. J. (2020). Fungal pathogens occurring on Orthopterida in Thailand. Persoonia 44, 140–160. doi: 10.3767/persoonia.2020.44.06
Vivekanandhan, P., Alahmadi, T. A., and Ansari, M. J. (2024). Pathogenicity of Metarhizium rileyi (Hypocreales: Clavicipitaceae) against Tenebrio molitor (Coleoptera: Tenebrionidae). J. Basic Microbiol. 2300744:e2300744. doi: 10.1002/jobm.202300744
Vivekanandhan, P., Swathy, K., Alford, L., Pittarate, S., Subala, S. P. R. R., Mekchay, S., et al. (2022). Toxicity of Metarhizium flavoviride conidia virulence against Spodoptera litura (Lepidoptera: Noctuidae) and its impact on physiological and biochemical activities. Sci. Rep. 12:16775. doi: 10.1038/s41598-022-20426-x
Vivekanandhan, P., Swathy, K., Lucy, A., Sarayut, P., and Patcharin, K. (2023). Entomopathogenic fungi based microbial insecticides and their physiological and biochemical effects on Spodoptera frugiperda (JE smith). Front. Cell. Infect. Microbiol. 13:1254475. doi: 10.3389/fcimb.2023.1254475
Wang, B., Jhun, J.-G., and Moon, B.-K. (2007). Variability and singularity of Seoul, South Korea, rainy season (1778–2004). J. Clim. 20, 2572–2580. doi: 10.1175/JCLI4123.1
White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. Academic Press, Cambridge, MA.
Wold, S., Esbensen, K., and Geladi, P. (1987). Principal component analysis. Chemom. Intell. Lab 2, 37–52. doi: 10.1016/0169-7439(87)80084-9
Yamaguchi, H., and Nakamura, K. (2015). Effects of environmental factors on the regulation of egg diapause in the walking-stick insect, Ramulus irregulariterdentatus (Phasmatodea: Phasmatidae). Eur. J. Entomol. 112, 35–40. doi: 10.14411/EJE.2015.011
Yan, S., Meng, Z., Tian, S., Teng, M., Yan, J., Jia, M., et al. (2020). Neonicotinoid insecticides exposure cause amino acid metabolism disorders, lipid accumulation and oxidative stress in ICR mice. Chemosphere 246:125661. doi: 10.1016/j.chemosphere.2019.125661
Yano, K., Ozaki, T., Suzuki, T., Yamazaki, H., Nasuno, M., Degawa, Y., et al. (2021). Outbreak of the stick insect, Ramulus mikado (Phasmatodea, Phasmatidae), in the akashina area of Japan (Azumino City, Nagano prefecture). Entomol. Sci. 24, 196–200. doi: 10.1111/ens.12467
Zhao, J., Chen, Y., Keyhani, N. O., Wang, C., Li, Y., Pu, H., et al. (2023). Isolation of a highly virulent Metarhizium strain targeting the tea pest, Ectropis obliqua. Front. Microbiol. 14:1164511. doi: 10.3389/fmicb.2023.1164511
Keywords: Ramulus mikado, Metarhizium phasmatodeae, rainy season, biocontrol agent, entomopathogenic fungi
Citation: Min D, Shin S, Lee N-H, Baek MJ, Park S-J, Kim K-H, Son H and Jung J-K (2024) Survival dynamics of stick insect and the impact of environmental factors on natural fungal infection during the rainy season. Front. Microbiol. 15:1383055. doi: 10.3389/fmicb.2024.1383055
Edited by:
P. Vivekanandhan, Chiang Mai University, ThailandReviewed by:
Ma. Cristina Del Rincón-Castro, University of Guanajuato, MexicoYasir Iftikhar, University of Sargodha, Pakistan
K. Meenambigai, Kirupananada Variyar Engineering College, India
Copyright © 2024 Min, Shin, Lee, Baek, Park, Kim, Son and Jung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hokyoung Son, aG9nb25naTdAc251LmFjLmty; Jong-Kook Jung, amtqdW5nQGthbmd3b24uYWMua3I=
†These authors have contributed equally to this work and share first authorship