Introduction
In recent decades, the global incidence of dengue fever has surged significantly, largely attributed to inadequate surveillance of Aedes mosquito breeding sites (Mehlhorn et al., Reference Mehlhorn, Schmahl and Schmidt2005; Buliva et al., Reference Buliva, Chasela, Elnossery, Tayyab, Mahboob, Marinda, Kinge, Okwarah, Brennan and Abubakar2023). The World Health Organization (WHO) has alarmingly estimated that more than 2.5 billion individuals are at risk of contracting dengue viruses worldwide (Piot et al., Reference Piot, Larson, O'Brien, N'kengasong, Ng, Sow and Kampmann2019). This upsurge in dengue cases has been accompanied by an extensive spread of the virus and its primary mosquito vector, Aedes species, which are now prevalent in regions between latitudes 35°N and 35°S across the world (Wichmann et al., Reference Wichmann, Gascon, Schunk, Puente, Siikamaki, Gjø, Lopez-Velez, Clerinx, Peyerl-Hoffmann and Sundøy2007; Qualls et al., Reference Qualls, Moser, Pereira and Koehler2022). Currently, dengue fever is endemic in over 100 countries, spanning tropical and non-tropical regions. Additionally, instances of imported cases have been reported in non-endemic nations (Domingo et al., Reference Domingo, Alves, de Ory, Teichmann, Schmitz, Müller and Niedrig2015). The most substantial burden of this disease is borne by regions such as Southeast Asia, the Americas, and the Western Pacific. The first dengue cases and outbreaks in Pakistan were reported in 1960 and 1994, respectively (Rasheed et al., Reference Rasheed, Butlin and Boots2013; McKenzie et al., Reference McKenzie, Stevens, McKenzie, Bozic, Mathias and Zohdy2019). In 2011 and 2013, the largest (Khan et al., Reference Khan, Khan and Amin2016) and second-largest dengue epidemics worldwide occurred in Karachi, Lahore, and Islamabad Regions (Wesolowski et al., Reference Wesolowski, Qureshi, Boni, Sundsøy, Johansson, Rasheed, Engø-Monsen and Buckee2015; Saeed et al., Reference Saeed, Ali, Khan, Muhammad, Reboita, Khan, Goheer, Khan, Kumar and Ikram2023). The country's climate, urbanisation, and poor sanitation infrastructure made it an ideal site for Aedes mosquitoes, leading to a high burden of dengue fever (Muhammad et al., Reference Muhammad, Ahmad, Abbas, Arif, Jan, Sarwar, Aziz, Alam, Sumbal and Hameed2023; Saeed et al., Reference Saeed, Ali, Khan, Muhammad, Reboita, Khan, Goheer, Khan, Kumar and Ikram2023). Analysis of these outbreaks showed an increasing trend of dengue occurrence in other zones (Bhakkar, Kallurkot, Darya Khan, and Mankera) with varying degrees of intensity (Khan et al., Reference Khan, Khan, Ghaffar and Khalid2018). In response to these outbreaks, vector control measures, such as larval surveillance and control, have been implemented (Bakhsh et al., Reference Bakhsh, Sana and Ahmad2018). The effectiveness of these measures in reducing dengue transmission has yet to be assessed in some areas (Kazi et al., Reference Kazi, Qazi, Ahsan, Khawaja, Sameen, Saqib, Khan Mughal, Wajidali, Ali and Ahmed2020). The escalating incidence of dengue across different regions of Pakistan, combined with the conducive environmental conditions for Aedes mosquito proliferation, highlights the urgent need for effective vector control measures.
Weather conditions, such as temperature, relative humidity, and rainfall, significantly influence local dengue transmission (Custódio et al., Reference Custódio, Nogueira, Souza, Fernandes, Oshiro, Oliveira, Piranda and Oliveira2019). Warm, wet climates provide ideal habitats for Aedes mosquitoes (Brady et al., Reference Brady, Johansson, Guerra, Bhatt, Golding, Pigott, Delatte, Grech, Leisnham and Maciel-de-Freitas2013). Defining the relationship between surveillance and weather conditions provides an effective model for devising control strategies (Colón-González et al., Reference Colón-González, Lake and Bentham2011). Implementing these strategies, we can effectively reduce the larval population and generate a climatogram that serves as an early warning message. This integrated approach plays a pivotal role in effectively managing dengue transmission and mitigating its impact on public health.
Dengue virus (DENV) is transmitted in Pakistan and worldwide primarily by Aedes albopictus (Skuse) and Aedes aegypti (L.) (Khan et al., Reference Khan, Khan, Ghaffar and Khalid2018; Brady and Hay, Reference Brady and Hay2020). Ae. albopictus also serves as a secondary vector for viruses (Gardner et al., Reference Gardner, Chen and Sarkar2017), such as West Nile virus, yellow fever virus, and Zika virus (Benelli and Romano Reference Benelli and Romano2017; Leta et al., Reference Leta, Beyene, De Clercq, Amenu, Kraemer and Revie2018). Aedes larval surveillance and control, which involves enhanced monitoring of mosquito development sites and the application of larvicides to control mosquito populations, have reduced dengue fever incidence in other countries (Nguyen et al., Reference Nguyen, Tuyet-Hanh, Mulhall, Minh, Duong, Chien, Nhung, Lan, Minh and Cuong2022). However, the implementation and effectiveness of these measures in Pakistan have not been evaluated (Khan et al., Reference Khan, Khan, Ghaffar and Khalid2018). In 2019–2020, the Punjab government developed an online electronic-based Dengue Tracking System (DTS) (Mahroof and Haider, Reference Mahroof and Haider2022). Using this data management system, we collected data for analysing the relationship between environmental factors and dengue fever occurrence. The aim of this study was to implement an integrated approach to prevent the onset of a dengue epidemic by effectively managing the distribution and abundance of vectors. The specific objectives were as follows: (1) Reduction of vectors population across the entire area by identifying and analysing the abundance of Aedes species and eliminating their immature developmental sites (both indoor and outdoor). (2) Cross-validate and eliminate immature developmental sites through Third-Party Validation (TPV), employing both mechanical and chemical larviciding methods. (3) Investigate the association between Aedes population incidence and key weather factors. (4) Assess the effectiveness of an Anti-dengue mobile application in enhancing Aedes surveillance and monitoring dengue occurrence trends.
Materials and methods
Aedes surveillance
The entomological survey was conducted in four towns (locally called tehsil) in Bhakkar District (Bhakkar, Kallurkot, Darya Khan, and Mankera Tehsil) during 2019 and 2020. Every town was surveyed weekly, as it took a week to cover each district. A total 48 surveys were conducted each year. The district's health department assigned two indoor teams and one outdoor team to each of the district's 42 Union Councils (UC) to conduct door-to-door surveillance for the presence of Aedes larvae within natural and artificial containers (Table 1). Every surveillance team had an electronic device containing the Anti-dengue Android mobile application to report containers with confirmed Aedes larvae. Surveillance was monitored by the online Dengue Tracking System (DTS) developed by PITB (Punjab Information Technology Board). The teams uploaded data to the online dengue dashboard that tracked all the reported Aedes positive containers, enabling the implementation of control measures and the monitoring of positive containers.
Most venerable containers for mosquito immature development (considered containers with larvae only).
Responsibilities of third-party validation and his role in dengue vector management
Larval Aedes and DENV case surveillance were conducted by district health department teams. TPV team's role was the identification and elimination of larval species through field inspection to ensure the precise allocation of resources and to avoid wasting time and effort. TPV teams conducted field inspections to validate data collection, assure accuracy, and evaluate the authenticity of the mosquito control activities and dengue cases after laboratory confirmation until sites were completely clear from the vectors. The third party here consisted of a paid independent analyst, which governmental organisations rely on to confirm performance (Abramson et al., Reference Abramson, Aldrich, Paolino and Rohde1995). TPV confirmed the location, larval presence, and species identification and suggested control strategy according to the situation (Mintz, Reference Mintz2014; Rentz et al., Reference Rentz, Ulman and Galster2021). Subsequent visits were continued until the area was entirely free from Aedes larvae.
The Punjab government selected a scientific officer from the Arid Zone Research Institute, Bhakkar (ARZI), as the TPV agency to determine the accuracy of anti-dengue activity data within district Bhakkar. TPV of dengue vector surveillance (larval surveys) was performed through an Android mobile phone Anti-dengue TPV application, which was registered by PITB to use for the data management system DTS (fig. 1). Once a potential positive container was identified by Health Department surveillance teams, three types of surveillance were done on the android-based application; namely, audit larvae, audit patients, and audit activities (fig. 2) were uploaded via the Anti-dengue TPV application by the ARZI teams. At the same time, a photograph was taken of each larval habitat and uploaded. ARZI teams then performed the following tasks: (1) Checked the containers based on audit larvae and the photograph. (2) Validated previous and current dengue fever cases based on blood test reports. (3) Suggest chemical control measures (larvicides: Temephos 50EC @ 1 ml/liter water; adulticides: Pirimiphos-methyl CS applied IRS @ 1000 mg/m2 (Actellic® 300CS, Syngenta) with hand compression sprayer and Deltamethrin and Cypermethrin (2 ml/Liter water) spray applied at 7-day intervals with ULV fogger machine sprayer (ULV application rate: 0.5-1Liter/Hectare; Wall spray: 75 ml/10-Liter water). (4) Aedes larvae were collected from each sampling site and identified for analysing Aedes species abundance – larvae were morphologically identified using taxonomic keys (Rueda, Reference Rueda2004). (5) Finally, the teams validated the door-to-door visits based on audit activities and the GPS location of past activities.
Association of environmental data with vector population
Environmental data on temperature, rainfall, and relative humidity were obtained from regional weather stations (Arid Zone Research Institute, Bhakkar) and evaluated for their role in mosquito population dynamics. The average distance between the surveillance sites and weather stations was approximately 10–20 kilometres. Data were analysed using yearly data combined from all containers.
Statistical analysis
Data on total visited spots (sites visited for surveillance), number of positive containers (confirmed Aedes-larvae number ≥1 for 30 days), active patients (number of DEN positive fever cases at each location), incidence per village (new DEN cases per unit of population per time interval), and weather were derived through the DTS system and processed in two ways. (1) Correlation analysis: The incidence of positive containers per month was correlated with temperature (°C), rainfall (mm), and relative humidity (%) on monthly basis (all the variables were assumed to be random and distributed normally). (2) Multiple linear regression analysis: This was used to calculate the relationship in the form of an equation (the incidence of positive containers was considered random and distributed normally). Descriptive statistics were used to analyse Aedes species abundance and the indoor and outdoor container contribution to Aedes immature development, i.e. the number of positive containers (n) and percentage (%). The coefficient of determination (R2) for the regression models and impact (%) were calculated using statistical software Origin 2023b. All graphical representations were made in Origin 2023b.
The prevalence of vectors in different towns was compared using the following indices.
Results
Indoor surveillance indicated that evaporative air coolers (air water coolers), discarded tyres, tyre shops, discarded house freezer trays, bird baths (pans containing drinking water for flying birds), and water storage drums were the containers found most frequently positive for larvae in 2019 and 2020. Containers most frequently associated with evaporative air conditioners contributed 30% and 33% of positive larvae in 2019 and 2020, respectively. Outdoor surveillance found containers positive for larvae at tyre shops, nurseries, and cemeteries, as well as standing water in tube wells in 2019 and 2020. Tyre shops were the most frequently positive containers, contributing 30% and 32% in 2019 and 2020, respectively (Table 2).
N, number of positive Aedes containers; % = Container contributed to Aedes immature development.
The Aedes species abundance from the mean number of immatures collected from indoor and outdoor positive containers in each town of Bhakkar district was determined by identifying larval populations in a laboratory (Table 3). The average number of Aedes larvae that were collected during surveillance revealed that Ae. albopictus was more abundant species as compared to Ae. aegypti from 2019 to 2020. Specifically, more Ae. albopictus larvae were collected from the field than Ae. aegypti. Among the sampling locations, Bhakkar and Darya Khan had the highest abundance of Ae. albopictus larvae confirmed after identification.
Mean values with standard error of immature collected from indoor and outdoor positive Aedes containers.
In areas where surveillance was increased to control the development of Aedes by eliminating stagnant water and other potential developmental sites, there was a 50.20% reduction in the number of positive mosquito larval sites following an increased TPV (the follow-up surveillance of the visited spots by TPV) of 22.64% (Table 4). Additionally, there was a 70.47% decrease in the number of active febrile patients, with the rate of TPV satisfactorily (successfully treated spots that had been identified to have Aedes mosquito) conducted at 9.09% higher in 2020 compared to 2019. These mosquito control activities were carried out throughout the province, which has 36 districts. This 70.47% reduction in patients reflects the overall impact of all surveillance and control activities in district Bhakkar only (Table 4). All cases were individuals with a confirmed history of travel outside of the province.
Visited containers: Sites visited to identify Aedes larvae.
Positive spot: Visited containers where Aedes larvae were confirmed.
TPV conducted: The follow-up surveillance of the visited containers by TPV.
Active febrile patients: Confirmed patients who were carrying the dengue virus.
TPV satisfactory (%): Successfully treated containers that had been identified to have Aedes mosquito.
The number of containers positive was highest in September 2019, and August 2020 was significantly associated with temperature, whereas there was no significant correlation with relative humidity and rainfall. During the peak months of April, May, August, and September in 2019 and 2020, positive containers included newly positive and repeated positives in surveys to eradicate the Aedes immatures populations until these containers had 0% infestation after repeated surveys (Table 5). A regression equation showed the impact of these factors on per-unit change in container frequency. Temperature was accounted for 43.3 and 48.9% of the variance in the number of Aedes positive developmental sites recorded during the years 2019–20, respectively, whereas relative humidity was responsible for 0.1% and 0.8% during 2019 and 2020, respectively (Table 6).
Note: Containers = numbers of positive containers × numbers of valid investigation; valid investigation means investigation of at least one dengue larva being identified.
X, Containers; Y, Average temperature; Z, Relative humidity; R, Coefficient of determination.
* = Significant, ns = non-significant.
The prevalence of Aedes presence in towns was compared by calculating the House (HI), Container (CI), and Breteau Indices (BI) (fig. 3) from indoor and outdoor positive containers. The Bhakkar township had the highest HI and BI indexes, followed by Kallurkot, Darya Khan, and Mankera in both 2019 and 2020. Different types of standing water had varied impacts on container prevalence and index calculations. The prevalence index refers to a measure of the abundance of Ae. albopictus and Ae. aegypti mosquitos in a given area. Evaporative air coolers, discarded house refrigerator trays, bird baths, water drums, tyre shops, nurseries, standing water in tube wells, and cemeteries all contributed to Aedes developmental sites. Maximum HI (%) index of 15.2 and 14.9%, CI (%) index of 25.3 and 21.1%, and BI (%) index of 10.8 and 7.5% during 2019–20 were recorded at tehsil Bhakkar while the lowest HI (%) index, CI (%) index, and BI (%) index were recording at Mankera, respectively.
Discussion and conclusion
Dengue fever is a serious problem in Pakistan but can be controlled through the validation of potential developmental sites and well-managed surveillance of Aedes-positive containers, which can result in a decrease in the number of mosquito developmental sites (Mohamud et al., Reference Mohamud, Qazi, Latif, Khan and Anwar2020). Surveillance on the cross-validation and management of developmental sites of Ae. albopictus and Ae. aegypti mosquitoes using Anti-dengue mobile application technology was implemented, coupled with TPV control operations (natural, mechanical, and chemical control), to stop the local transmission of the dengue virus from imported cases in the Bhakkar district (Table 2; Fig. 2). People traveling to other areas, including those within Pakistan, were the main reason for the dengue cases detected in the province. If local Aedes mosquitoes bite these individuals, they may become virus carriers, leading to locally-acquired dengue cases (Rimal et al., Reference Rimal, Shrestha, Pandey, Nguyen, Bhandari, Shah, Acharya, Adhikari, Rijal and Ghimire2023). Studies have shown that most of Pakistan's affected areas include cases with a travel history, which aligns with our findings (Ahmad et al., Reference Ahmad, Khan and Jamal2020; Rehman et al., Reference Rehman, Haq, Asghar, Afridi and Faisal2020). The reduction of a significant number of positive containers of Aedes mosquito control resulted in no local transmission from imported cases from 2019 to 2020.
Ae. albopictus predominated over Ae. aegypti in all the zones of District Bhakkar, although declining gradually in 2020. In agreement, immatures of the were more abundant in Punjab than those of the Ae. aegypti (Saleem et al., Reference Saleem, Ghouse, Hussain, Saleem and Abbas2014; Akhtar et al., Reference Akhtar, Gupta and Singh2022). These findings were in accordance with Ae. albopictus immatures found more abundant species in Punjab province as compared to Ae. aegypti (Wint et al., Reference Wint, Jones, Kraemer, Alexander and Schaffner2022; Ullah et al., Reference Ullah, Hafeez, Ali, Arshad, Akram, Ali, Saba, Abbas, Iftikhar and Jafir2023).
Our enhanced surveillance for Aedes developmental sites identified peak periods in April, May, August, and September for positive indoor and outdoor containers, indicating increased immature Aedes incidence due to favourable environmental conditions. The months of August, September, and October in Pakistan have the most significant number of Aedes larval developmental sites (Table 5). The implications of identifying peak periods for Aedes larval activity were profound, as they hold crucial significance for public health and vector control efforts (Patil et al., Reference Patil, Yadav, Dasgupta, Zehr, Barwale and Char2021). Many studies are completely aligned with our findings that continuous anti-dengue activities (surveillance and control measures) posed direct reflection on the reduction of vector population and induced firm grip on elimination of local transmission (Sharmin et al., Reference Sharmin, Glass, Viennet and Harley2018; Xia et al., Reference Xia, Dweck, Lutomiah, Sang, McBride, Rose, Ayala and Powell2021). Our study supports continuous vectors management strategies using mobile technology (reduction in larval developmental sites) during the identified peak months, as this strategy disrupts the mosquito developmental cycle (egg, larvae, and pupae) and inhibits the Aedes adult emergence, thereby reduces local transmission (Hafeez et al., Reference Hafeez, Akram and Shaalan2011; Faraji and Unlu, Reference Faraji and Unlu2016). Integrating our findings into existing public health frameworks enhances preparedness and responsiveness to potential Aedes developmental sites. As we move forward, continued research and surveillance efforts using mobile-based technology will be instrumental in refining vector control strategies and safeguarding communities against the threats posed by Aedes mosquitoes.
The temperature was responsible for the variation in the number of positive Aedes developmental sites recorded during the years 2019 and 2020, respectively (Table 6). Aedes mosquitoes thrive in warm and humid environments, with temperatures ranging from 25–28 °C being optimal for their growth and reproduction (Aryaprema et al., Reference Aryaprema, Steck, Peper, Xue and Qualls2023). When temperatures rise above this range, mosquito populations can decrease rapidly because of the reduction of Aedes developmental sites. Temperature is a fundamental driver for the Aedes physiology and developmental periods (Ware-Gilmore et al., Reference Ware-Gilmore, Sgrò, Xi, Dutra, Jones, Shea, Hall, Thomas and McGraw2021). Our results agree with the study that temperature exerted a significant impact, while relative humidity and rainfall had limited contributions to the transmission of dengue fever (Ware-Gilmore et al., Reference Ware-Gilmore, Sgrò, Xi, Dutra, Jones, Shea, Hall, Thomas and McGraw2021). Our results on rainfall agree with the study that rainfall had a positive but non-significant correlation with the positive containers, whereas average temperature exhibited considerable influence on containers (Naqvi et al., Reference Naqvi, Sajjad, Waseem, Khalid, Shaikh and Kazmi2021; de Souza et al., Reference de Souza, de Camargo Guaraldo, Honório, Câmara, Sukow, Machado, Duarte dos Santos and da Costa-Ribeiro2022). The incidence and spatial distribution of Aedes larvae were also found to be influenced by temperature (Couret et al., Reference Couret, Dotson and Benedict2014; de Oliveira Lemos et al., Reference de Oliveira Lemos, Júnior and de Assis Mendonça2021). Surveillance of Aedes developmental sites during these particular environmental conditions will be effective.
Different sources of standing water had varying impacts on container and frequency index calculations here and in previous studies. Standing water sources such as evaporative air coolers, discarded house refrigerator trays, bird pots, water drums, tyre shops, nurseries, tube wells, and cemeteries had maximum containers and contribution for the Aedes developmental sites (Table 2) with the average house index (HI), container index (CI), and breteau index (BI) of 8.75, 14.0 and 5.6% in different towns of Bhakkar (fig. 3). This agreed with previous findings that water tanks, evaporative air coolers, discarded tyres, drums, tree holes, pots, flower vessels, and discarded cups were the most common Aedes developmental sites (Kesetyaningsih et al., Reference Kesetyaningsih, Andarini, Sudarto and Pramoedyo2018; Fatima, Reference Fatima2019; Khan Reference Khan2021; Hafeez et al., Reference Hafeez, Naeem-Ullah, Akram, Arshad, Iftikhar, Naeem and Saleem2022). A study found the highest number of dengue mosquito larvae and pupa found in plastic boxes, grinding stones, discarded cups, earthen pots, and discarded tyres, with nearly identical trends of HI, CI, and BI, as reported in our study (Saleem et al., Reference Saleem, Ghouse, Hussain, Saleem and Abbas2014). Interestingly, few researchers reported high values of HI, CI, and BI because the Bhakkar district largely consists of urban areas with relatively low intensity of dengue virus as compared to other larger, more-populated cities. In our study, dengue fever cases were successfully prevented from local transmission from imported cases as a result of vector prevention efforts, and there were no cases of locally transmitted dengue fever reported because of the elimination of developmental sites, insecticide spraying, and public education campaigns. The management of any disease and its vector, particularly dengue and Aedes mosquitoes, depends critically on public awareness and education of Aedes developmental sites.
The study on enhanced larval surveillance of Aedes mosquitoes conducted in Bhakkar provided positive results which demonstrated the potential of technology in improving dengue control efforts. Anti-dengue mobile application and online integrated PITB-DTS have been successfully implemented in district Bhakkar, demonstrated the significance of technology in advancing public health (fig. 2). The dynamic and complex nature of dengue transmission and control was enhanced through the use of smartphones and mobile apps (Regis et al., Reference Regis, Monteiro, Melo-Santos, Silveira, Furtado, Acioli, Santos, Nakazawa, Carvalho and Ribeiro2008). These innovative ways provided real-time information on dengue cases, vector presence, and control activities, enabling users to take preventive measures and avoid potential exposure to the disease (Orong et al., Reference Orong, Durias, Paredes and Pasay2015; Razaq et al., Reference Razaq, Batool, Ali, Khalid, Saif and Naseem2016). It facilitated coordination and collaboration between various stakeholders involved in dengue control, including public health authorities, healthcare providers, community leaders, and the general public (Pley et al., Reference Pley, Evans, Lowe, Montgomery and Yacoub2021). Further research and implementation of these technologies to combat other vector-borne diseases and public health issues are crucial.
In conclusion, this study provided valuable insights into the factors contributing to the proliferation of Aedes developmental sites and highlighted the potential of technology-based solutions for disease surveillance and control. This study has identified optimal temperature conditions (23–28℃) and key months (April, May, August, and September) conducive to the development of the dengue vector. By emphasising the critical need to eliminate standing water sources such as evaporative air coolers and tyre shops, our findings offer enhanced strategies for more effective dengue prevention and control. The effectiveness demonstrated by the online PITB-DTS configured system and Anti-Dengue mobile application in District Bhakkar highlights the significance of harnessing technology to enhance public health outcomes. This is particularly crucial in resource-limited settings, where traditional disease surveillance and control methods may prove insufficient. The findings of this study can serve as a model for other regions with similar characteristics, providing valuable insights into the development of effective and sustainable dengue control programs. Moreover, this study contributes to the broader field of vector-borne disease control by confirming the importance of taking a comprehensive approach to disease prevention and control. Overall, this study contributes significantly to our comprehension of dengue transmission and control, providing valuable insights that can guide the formulation of more efficacious prevention and control strategies in regions affected by the disease.
Acknowledgements
We thanks for the support of Government of Punjab, Pakistan, Ministry of Science and Technology of China (2022YFD1500701), and the Italian National Biodiversity Future Center [CN00000033].
Availability of data and materials
Not applicable
Author contributions
Sohail Abbas: Conceptualisation, Methodology, Writing-Original Draft, Writing-Review and Editing, Data Curation, Introduction, Results and References; Muneer Abbas: Methodology, Data collection, Data Analysis, Discussions, Writing-Original Draft, Writing-Review and Editing, Supervision; Aleena Alam: Schematic Diagram and Data Curation, Writing-Review and Editing, Niaz Hussain: Writing-Review and Editing, Muhammad Irshad: Writing-Review and Editing; Mudassar Khaliq: Resources; Xiao Han: Software; Faisal Hafeez: Writing-Review and Editing; Donato Romano: Methodology and Writing-Review and Editing, Supervision; Chen Ri Zhao: Supervision, Methodology, Funding, Resources, Writing-Review and Editing.
Financial support
We thank the Ministry of Science and Technology of China (2022YFD1500701) for the funding
Competing interests
None.
Ethical approval
Not applicable