Introduction
Understanding the exposure dynamics of zoonotic vector-borne diseases is important to inform public health measures, such as the provision and planning of mosquito control programmes. However, such studies are difficult to conduct in human populations as seroprevalence data are often derived from clinically ill patients. Thus, subclinical infection rates in humans are not known and cannot be incorporated into infectious disease models. Sentinel animals, that is, naïve animals strategically placed for the monitoring of incursion and transmission of economically and medically important exotic and endemic diseases, are therefore often used in lieu [Reference Racloz, Griot and Stärk1]. Livestock, such as sheep and cattle, are routinely used as sentinel animals for the Australian National Arbovirus Monitoring Programme [2]. However, ruminants are not relevant to the transmission of Ross River virus (RRV), as they only express a variant of the alphavirus receptor, Mxra8, which cannot bind RRV [Reference Kim3], and hence they are not involved in the transmission dynamics of RRV [Reference Stephenson4].
Ross River virus is an alphavirus, belonging to the family Togaviridae, causing most arboviral infections in Australia, with an average of around 5,000 case notifications annually [Reference Yuen and Bielefeldt-Ohmann5]. RRV transmission is maintained between the reservoir host, such as marsupials (e.g., brushtail possums and wallabies), and mosquitoes [Reference Yuen and Bielefeldt-Ohmann5]. Impressively, more than 40 species of mosquitoes have been identified to be capable of transmitting RRV [Reference Russell6, Reference Harley, Sleigh and Ritchie7]. To date, only humans and horses have been confirmed as susceptible hosts of RRV disease. Infected humans are reported to experience fever, joint pain, and lethargy [Reference Condon and Rouse8–Reference Westley‐Wise11], with some patients experiencing chronic clinical signs of more than 12 months’ duration [Reference Condon and Rouse8, Reference Selden and Cameron10, Reference Westley‐Wise11]. Published case studies/series have consistently reported similar clinical signs in horses to those described in infected humans [Reference Azuolas12–Reference Kay15].
Previous research identified high seroprevalence of RRV in horses in a subtropical region (Lockyer Valley) of South-East Queensland (SE QLD), Australia [Reference Yuen16]. This region is characterized by intensive horticulture production and has extensive pastures that support a high density of horses. It is adjacent to large metropolitan areas, including Brisbane, the Gold Coast, and the Sunshine Coast, where millions of people reside, and large cohorts of tourists visit at all times of year. The year-round presence and interaction between marsupials (reservoir hosts), the mosquito population (vectors), and horses (susceptible hosts) in this region [Reference Yuen16] provide a unique opportunity for the study of exposure dynamics of RRV and for conducting an assessment of the potential usefulness of horses as sentinel animals for RRV transmission.
Vector-borne disease transmission cycles are depended on reservoir host and vector dynamics, which is influenced by meteorological conditions. These factors fluctuate in different parts of Australia depending on landscape features and local weather. Various types of predictive models have been developed to identify environmental risk factors associated with increased RRV transmission. It appears that rainfall, with a lag period of 0 to 3 months, is the most frequently identified meteorological risk factor for RRV outbreaks in humans across Australia [Reference Yuen and Bielefeldt-Ohmann5, Reference Kelly-Hope, Purdie and Kay17–Reference Hu24]. Temperatures (both minimum and maximum) and relative humidity, with 1–2 months lags, were also frequently identified as environmental drivers for increased RRV notification in humans [Reference Tong18–Reference Tong and Hu21, Reference Jacups23]. Significant RRV hotspots for humans have been identified in peri-urban suburbs in SE QLD, Australia, where residential areas, agricultural practice, and conserved natural landscapes intersect [Reference Murphy25].
The study presented here aimed to determine the exposure dynamics of RRV in horses over a 3.5-year period. Specifically, (i) estimation of the rate of waning of RRV-specific neutralizing maternal antibodies in foals, (ii) estimation of the rate of acquiring natural infection of RRV in naïve horses, and (iii) identification of meteorological variables (as a surrogate for vector activity) that are associated with RRV infection in naïve horses. This represents the first longitudinal study to identify environmental risk factors of RRV transmission in an agricultural region in Australia.
Material and methods
Study area
We conducted a prospective longitudinal study over a 3.5-year period, commencing 1 August 2020 and ending 31 June 2023. The study was conducted in a subtropical region of SE QLD, Australia, in the local administrative area of the Lockyer Valley (Figure 1). This region is characterized by intensive horticulture production, which include cucurbits, legumes, Brassica spp., corn, and sorghum. Fertile and well-irrigated soils in the region also provide excellent pastures for horses. The study area is adjacent to the administrative areas of Brisbane, the Gold Coast, and the Sunshine Cost. The rural University of Queensland (UQ) Gatton Campus is located in the heart of the study region. The campus covers an area of more than 1,000 ha, with approximately 100 horses at any one time within a year, and also incorporates crop and grazing fields, dairy and sheep production, a piggery, and a wildlife centre [26]. Marsupials such as red-necked wallabies (Notamacropus rufogriseus), Northern brown (Isoodon macrourus), and long-nosed (Perameles nasuta) bandicoots, common brushtail possums (Trichosurus vulpecula), rufous bettong (Aepyprymnus rufescens), and Eastern grey kangaroos (Macropus giganteus) are regularly reported in the study area [27].
Foals born and residing at the Equine Unit on the UQ Gatton Campus from 2020 to 2022 (n = 32) were enrolled in this study. Each cohort of horses born in the same year were agisted together in the same paddock, with regular paddocks rotations, throughout the study period. All paddocks were irrigated periodically as required depending on pasture availability and rainfall. Water troughs, with automatic filling system, were present in all paddocks and hay were supplemented as required depending on pasture quality and availability. Paddocks were free of trees with shelters provided.
Sample collection
Blood samples from foals were collected via jugular venipuncture using a 20G hypodermic needle pre-suckle, post-suckle (24 h post foaling), then approximately monthly until six months of age, after which foals were weaned and sampled opportunistically approximately every two to three months. Blood samples from mares were collected at the time of foaling. No clinical signs were recorded from any animals seroconverted to RRV. Animal ethics approval was obtained from The University of Queensland (UQ) Production Animal Ethics Committee (permit numbers SVS/344/18 and 2021/AE000763) prior to commencement of the study.
Cell culture and virus amplification
Cell culture and virus amplification were performed as described previously [Reference Yuen16]. Briefly, Vero cells (African green monkey (Chlorocebus sp.)-derived kidney epithelial cells) and C6/36 cells (Aedes albopictus larvae mid-gut cells) were cultured as previously described [Reference Prow28]. The virus seed stock used in this study was RRV isolate Cairns 2007 [Reference Jansen29]. C6/36 cells were used for virus stock amplification. Amplified virus stock was stored at −80 °C in 1-mL aliquots. Virus titre was determined by virus titration assay as previously described [Reference Prow28] and calculated as TCID50 infectious units/mL.
Virus neutralization test
Virus neutralization tests (VNT) were performed as previously described [Reference Yuen16]. Briefly, each heat-inactivated serum sample was tested for the presence of neutralizing antibodies to RRV in duplicates serially diluted from 1:20 to 1:160. Positive (no serum) and negative (no virus) controls were included in each microtitre plate.
Meteorological data
Half hourly weather data collected from the UQ Gatton Campus weather station (station no. 040082) between 1 August 2020, and 3 July 2023, were made available for analysis. The weather station is owned and maintained by the Bureau of Meteorology (Australian Government). Due to equipment failure, no data were available between 26 February 2022, and 27 May 2022. Time periods representing the highest mosquito activities were used, i.e. dawn and dusk [Reference Webb, Doggett and Russell30]. Air temperature, dew point temperature, relative humidity, wind speed, and wind gust speed between 5 am–8 am (dawn) and 4 pm–7 pm (dusk), and daily cumulative rainfall were extracted. The three-hourly data for dawn and dusk were averaged, respectively, to represent the morning and afternoon mosquito activity. The averaged dawn and dusk data were then further averaged to represent the daily average data. Missing data points were computed using a multiple imputation approach (see Section title “Statistical analysis”).
Statistical analysis
Seropositivity was defined as a RRV neutralizing antibody titre >1:20. Depending on the temporal stage of seropositivity, each sample was assigned to group ‘0’, ‘1’, or ‘2’. Foals with RRV maternal antibodies detected after colostrum intake were noted as ‘1’. The first sample that tested negative to RRV neutralizing antibody after colostrum intake was assigned as ‘0’. Once natural infection was acquired, foals were assigned as ‘2’. Therefore, for the purpose of maternal antibody waning analysis, the failure event was detected when sample group changed from ‘1’ to ‘0’; and for natural infection, when ‘0’ changed to ‘2’.
Days for the waning of protective RRV neutralizing maternal antibodies and days from seronegative till RRV seroconversion were visualized using Kaplan–Meier survival curve and median survival time with 95% confidence intervals were calculated for these time periods [Reference Cleves, Gould and Marchenko31]. Day at risks for RRV infections commenced from the date when the first RRV-seronegative sample was collected after colostrum intake. Log rank tests were performed to compare survival time in waning of maternal antibodies and time till RRV seroconversion across years [Reference Cleves, Gould and Marchenko31]. The correlation between days till RRV seroconversion and days until loss of maternal antibodies was displayed in a scatter plot and quantified by the calculation of the Pearson’s correlation coefficient.
Time series analysis was performed on the meteorological daily data. Missing data were computed using multiple imputation [32]. For this, a univariate predictive mean matching imputation method with nearest neighbours set at 21 observations and posterior estimates calculated from a bootstrap sample was applied. Time series decomposition [33] was performed for each meteorological variable by year using the unobserved-component smooth-trend model with a seasonal component. Post-estimation prediction of unobserved components of trend, seasonal, and residuals were calculated using all sample information.
Bivariate Pearson’s correlation coefficients were calculated to explore the correlation between any two meteorological variables. If high correlation was identified (r > 0.95), one of the two variables was omitted.
To investigate the impact of the meteorological variables on the time till seroconversion in foals, semi-parametric Cox proportion hazards regression was used [Reference Cleves, Gould and Marchenko31]. As seropositivity for RRV was determined approximately every month, daily averaged meteorological data were averaged per month and daily cumulative rainfall was summarized per month.
Three analytical approaches were used to investigate the impact of meteorological variables on time until seroconversion. In the first approach, individually monthly meteorological variables in their original units of measurements were considered as time-varying covariates with lag periods of up to two months. In the second approach, correlation between monthly meteorological variables was accounted for by developing principal components of meteorological variables and the estimated components scores were used as covariates in the analysis. In the third approach, meteorological variables were standardized by calculating the percent change of a monthly mean in comparison with monthly mean in the previous month, and the percent changes of individual covariates were used in the analysis.
The proportional hazard assumption was tested using Schoenfeld residuals and the link test [Reference Cleves, Gould and Marchenko31]. Covariates with hazard ratios at p < 0.2 in the univariate models were considered for inclusion in the multivariate Cox regression models using a stepwise forward model building strategy. Likelihood-ratio tests were performed to compare models. Confounding was assessed by individually adding each covariate back to the final model, retaining any variables that resulted in >20% change in hazard ratio of covariates already in the model. Goodness of fit of the final model was assessed using Cox–Snell residuals. Statistical analyses were performed in Stata BE 17.0 (StataCorp. 2021. Stata Statistical Software: Release 17. College Station, TX: StataCorp LLC).
Results
Time for waning of RRV maternal antibodies and time until RRV natural infection
All pre-suckle samples were negative for RRV-specific antibodies. A total of 27 foals had neutralizing antibodies against RRV after colostrum intake. Foals (n = 5) born to RRV-seronegative mares remained seronegative to RRV after colostrum intake. The overall median time for the waning of RRV-neutralizing maternal antibodies in foals (n = 27) was 137 days (95% CI: 113–154), which is equivalent to 4.6 months (Figure 2, Supplementary Table S1). The overall median time from seronegative till seroconversion was 429 days (95% CI: 294–582), with the median age at seroconversion being 69 weeks (95% CI: 53–75). The overall cumulative incidence risk of seroconversion was 56% (95% CI: 38–74) over the study period of 3.5 years, with a yearly incidence risk of 18% (95% CI: 5–40) and 36% (95% CI: 19–56) in 2021 and 2022, respectively, and 22% (95% CI: 6–48) for the 6-month observation period in 2023. All foals from the 2020 cohort (n = 10) had seroconverted to RRV by the conclusion of the study (31 June 2023); however, six and eight foals from the 2021 and 2022 cohorts, respectively, remained seronegative (right-censored) at the end of the study period (Figure 3, Supplementary Table S1). The days for maternal antibodies to wane and time until seroconversion were not correlated (r = −0.0062; p = 0.98; Supplementary Figure S1).
Meteorological data
The daily averaged values for each meteorological variable were averaged by month and were presented by years in Table 1, and results are visualized in the Supplementary Material (Supplementary Figures S2–S13). Overall, air temperatures follow a seasonal pattern: spring (September to November), summer (December to February), autumn (March to May), and winter (June to August). Increased daily cumulative rainfall was observed between October 2021 and July 2022. All seroconversion events occurred between December and March each year.
Note: Boxes coloured in green indicate events of seroconversion occurring in respective months; boxes in grey refers to no data available (i.e. months outside the study period or months with failure of data recording equipment).
As expected, all meteorological variables were significantly correlated. Notably, wind speed and wind gust were highly correlated (r = 0.987, p < 0.0001); and monthly cumulative rainfall and monthly average rainfall were highly correlated (r = 0.9707, p < 0.0001). As the maximum wind speed would have a greater effect on mosquito dispersion than averaged wind speed, wind gust was chosen as the variable to be retained. Cumulative monthly rainfall was considered more important as a proxy for mosquito-breeding grounds than average monthly rainfall and was retained. Therefore, five meteorological variables (monthly averaged) were considered as covariates for time until seroconversion, namely air temperature, dew point temperature, relative humidity, wind gust, and cumulative rainfall.
Association between change in meteorological variables and time until new RRV infection
Modelling approach 1: Impact of individual meteorological variables in original units of measurement on seroconversion
The univariate analysis revealed that monthly air temperature, dew point temperature, and relative humidity (time-lag 0) were significantly (p < 0.05) associated with time till RRV seroconversion (Table 2). Every unit increase in air temperature and dew point temperature slightly increased the hazard ratio, while relative humidity slightly decreased the hazard ratio. The meteorological variables in the previous 2 months (time-lag 1 and 2) were not significantly (p > 0.05) associated with RRV seroconversion (except time-lag 2 for relative humidity). Wind gust and cumulative rainfall were not significantly associated with RRV seroconversion (p > 0.05) under modelling approach 1. Using air temperature (without time lag) as the base model, the addition of each meteorological variable did not result in a significant (p < 0.05) improvement of the final model, nor was any confounder identified. Therefore, the univariate regression results represent the final models for this modelling approach (Table 2). Goodness-of-fit results are shown in Supplementary Figure S14.
Bold entries highlight p-values that are significant at p < 0.05.
Given the non-significant results in time-lag variables from the first model, time-lags were not considered in the remaining two modelling approaches.
Modelling approach 2: Impact of principal components of meteorological variables on seroconversion
Principal component analysis identified five components, with three components having Eigenvalues <1, therefore being excluded from further analysis (Supplementary Table S2). Temperature related variables and wind gust had high positive loadings in component one, while moisture related variables had high positive loadings in component two (Figure 4). In the univariate analysis, an overall increase in each unit of component one (temperature-driven) led to an increased hazard ratio of seroconversion by 2.3 times (p = 0.013; Table 3). Component two (moisture-driven) was not significantly associated with seroconversion (p = 0.086; Table 3). Results of the goodness-of-fit test are shown in Supplementary Figure S15.
Bold entries highlight p-values that are significant at p < 0.05.
Modelling approach 3: Impact of percent change of individual meteorological variables compared to previous month on seroconversion
The univariate analysis in the third modelling approach revealed that every percent change in air temperature, dew point temperature, and wind gust compared to the previous month resulted in significantly increased hazard ratios by 2.9, 1.4, and 2.0 times, respectively (Table 4). The final multivariate model only included air temperature and wind gust, with both increasing significantly (p < 0.05) the hazard ratio for seroconversion (Table 5). Results of the goodness-of-fit assessment are shown in Supplementary Figure S16.
Bold entries highlight p-values that are significant at p < 0.05.
Bold entries highlight p-values that are significant at p < 0.05.
Discussion
Horses constitute one of the very few domestic mammalian species that are naturally susceptible to many arbovirus infections of medical concern and they display clinical signs similar to those of infected humans. In Australia, RRV infection is a notifiable disease in humans. The last annual report (2014–2015) on human infection was published in 2019 indicating that RRV remains the most prevalent arbovirus infection [Reference Knope34]. Since RRV infection in animals is non-notifiable in Australia, the prevalence of arbovirus infection is only known from published serological surveys. Since an outbreak of neurological arboviral disease in horses in 2011, where RRV infection was also reported in South Australia [Reference Roche35], there have only been three serological surveys conducted in horses [Reference Yuen16, Reference Prow36, Reference Gummow37]. Currently an active Australian National Arbovirus Monitoring Programme is conducted in cattle, sheep, and goats to monitor the transmission of bluetongue virus, Akabane virus, and bovine ephemeral fever virus [38]. These viruses are primarily transmitted by Culicoides spp. and can have a significant impact on livestock production, but they do not infect humans and horses and hence are not considered of public health concern.
Maternal antibodies protect foals from infection for the first three to four months of life post-partum. After this time, maternal neutralizing antibody levels fall below the detectable limit but may still interfere with endogenous immunoglobulin production in the foal upon infection or vaccination until six months of age [Reference Wilson39, Reference Bielefeldt-Ohmann40]. In the present study, protective levels of maternal antibodies specific for RRV in foals born to RRV seropositive mares were detectable for 4.6 months, with a 95% confidence interval of 3.8 to 5.1 months. Once maternal antibodies declined to non-protective levels (< 1:20 in VNT), foals are susceptible to natural infection of RRV by the bite of an infected mosquito. Interestingly, the foals in the current study were at risk of infection for 14.3 months (with a 95% confidence interval of 9.8 to 19.4 months) and all seroconversion events occurred between December and April consistently across all three foal cohorts. This supports the suggestion that naïve horses may be suitable sentinel animals for RRV circulation considering that notification rates in humans in SE QLD consistently peak later (between March and May), as identified by several 10-year retrospective studies [Reference Yu41, Reference Gatton42].
Across all three modelling approaches, temperature-related variables, especially air temperature, were consistently associated with seroconversion. As seroconversion in horses occurs over summer to early autumn (December to March), it is not surprising that the increase in temperature increased the hazard ratio for seroconversion. It appears that seroconversion events occurred when the averaged air temperature was approximately between 18 °C and 25 °C. This is likely related to mosquito activities as larvae development and host-seeking activities are optimal between 15 °C and 32 °C [Reference Drakou43, Reference Bellone and Failloux44], and as mosquito activity ceases when temperature falls below the minimum metabolic requirements [Reference Baril45]. In addition, percent change in wind gust was significantly associated with seroconversion. While the hazard ratio for percent change in model approach 3 does not indicate whether a positive or negative percent change in wind gust increases the hazard ratio, it is likely that the decrease in wind gust increases the hazard ratio as strong wind would disperse mosquitoes [Reference Kay and Farrow46] and decrease air temperature, thereby reducing overall mosquito activities and thus limiting the chances for infection to occur.
Owing to the differences in landscape and climate features across Australia, the transmission of RRV varies across different regions even within a state. For example, transmission of RRV occurs year-round in northern QLD [Reference Tong and Hu21, Reference Yu41, Reference Mackenzie47], but it is only seasonal in SE QLD [Reference Murphy25]. Studies identified a seroprevalence of RRV in horses at 91% in northern QLD [Reference Gummow37], and approximately at 50% in SE QLD, with the location for the current study having the highest prevalence at 85% [Reference Yuen16]. Results from the present study confirmed that transmission of RRV in horses in SE QLD follows a seasonal pattern with seroconversion events only being detected between December and March consistently throughout a 3.5-year period and being associated with air temperature without time-lag. In contrast, epidemiological modelling using human notification data in QLD showed that human infection peaked between March and May with rainfall, relative humidity, and air temperature identified as risk factors with 0–3 months’ time-lags [Reference Yuen and Bielefeldt-Ohmann5, Reference Kelly-Hope, Purdie and Kay17–Reference Hu24]. Year-round irrigation usage and filled water troughs on farms and in paddocks in this subtropical study area would provide ideal breeding grounds for mosquitoes. Within urban areas, the presence of water puddles after rainfall, poorly maintained swimming pools, bird baths, and water accumulating in plant pots may provide more short-term breeding grounds for mosquitoes. Also, in regional and rural areas, reservoir hosts, such as marsupials, may be present in larger numbers than in build-up areas and would maintain the transmission cycle all year. With the close proximity of mosquito breeding grounds and the shared habitat between reservoir hosts and horses, it is not surprising that a time-lag effect was not observed in this study.
Overall, this further supports the conjecture that RRV infection in horses may provide a proxy for increased RRV transmission in the human population. Real time monitoring of seroconversion events in horses in peri-urban areas may be more effective and reliable than predictions from mathematical modelling which uses historical data. This is especially true in the era of climate change as meteorological factors change unpredictably, making disease transmission forecasting through mathematical models very challenging. While mosquito surveillance programmes provide information about arboviruses circulating in the mosquito population, they do not confirm active exposure dynamics between vector and susceptible hosts. A near real-time monitoring of active circulation of RRV in horses would support the implementation of time-sensitive public health measures. Moreover, the seroprevalence of 50% in horses residing in southern QLD means that the transmission of RRV in the horse population is high enough to provide confidence in detecting seroconversion events in any given year [Reference Yuen16].
There are a few considerations before conducting such real-time monitoring. Firstly, a young horse population is preferable, as horses older than 6 years of age are two times more likely to have already seroconverted than 2–6 years old [Reference Yuen16]. Secondly, strategic or risk-based sampling would lower the costs of these studies. Given transmission in SE QLD consistently starts in December, sampling could begin in mid-spring each year and end in autumn when human cases would have peaked. Thirdly, a simple ELISA test or point-of-care test (e.g., lateral flow assay) could be developed and used as a diagnostic tool. As the goal of such monitoring is to detect the presence of neutralizing antibodies, the labour-intensive and time-consuming gold standard virus neutralization assay is not required to determine antibody titre. This would increase the turn-around time and help lower the cost of laboratory testing. In addition, since Getah virus remains exotic to Australia [Reference Rawle48] and since no alphavirus vaccine is currently available in Australia, antibody cross-neutralization would be of minimal concern. A similar monitoring approach could be adapted across Australia, and in neighbouring countries where RRV is endemic, or spillover events have been sporadically detected.
The study presented here used the presence of neutralizing antibodies as an indicator of infection. However, detectable seroconversion only occurs approximately two weeks after infection. Therefore, it should be acknowledged that a two-week lag exists in the present data. However, since blood sample collection was conducted on a monthly basis and analyses of meteorological risk factors were summarized monthly, the two-week lag effect of the detection of seroconversion from initial infection would have little-to-no impact on the analysis results. While mosquito testing was not performed to confirm the presence of RRV in the vector, mosquito trappings revealed that the majority of the mosquito population on the UQ Gatton campus consists of Culex annulirostris (data not shown), a known competent, efficient, and important vector for RRV transmission [Reference Russell6].
Conclusion
This study represents the first prospective longitudinal study of RRV transmission in a region of Australia with high transmission rate. In summary, the increase in air temperature and probably decrease in wind gust are associated with seroconversion events in horses without time-lag effect. These results support the use of horses as sentinel animals for RRV transmission to inform public health measures and horses could potentially be integrated into an early warning detection system of RRV epidemics or outbreaks. Given the importance of RRV infection medically and economically in Australia, serosurveillance of horses in peri-urban areas as in parallel with real-time human notification data in surrounding regions (a One Health approach) is warranted to confirm the suitability and effectiveness of using horses as sentinels for RRV infection.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/S0950268824000554.
Data availability statement
All data have been presented in the manuscript and Supplementary Material.
Acknowledgements
The authors would like to thank all staff at the Equine Unit (UQ) for assisting with animal handling and sample collection; Dr. Natalie A. Prow (QIMR Berghofer Medical Research Institute, Brisbane) for providing the RRV seed stock; and Dr. Vincent Mellor (School of Agriculture and Food Sustainability, UQ) for providing the meteorological data.
Author contribution
Nicholas Yuen contributed to methodology, sample collection, investigation, analysis, writing original draft, review, and editing. Helle Bielefeldt-Ohmann participated in conceptualization, methodology, investigation, analysis, review and editing, and supervision. Mitchell Coyle coordinated and was involved in sample collection. Joerg Henning guided and contributed to analysis, and contributed to review and editing, and supervision. All authors have read and agreed to the published version of the manuscript.
Funding statement
The study was funded by a grant from the Australian Infectious Diseases Research Centre (to Helle Bielefeldt-Ohmann), and a donation was received from BIOHMPATHOLOGY. Nicholas Yuen is a recipient of a UQ Research Training Scholarship.
Competing interest
Helle Bielefeldt-Ohmann is the proprietor of the consultancy firm BIOHMPATHOLOGY. All other authors declare no conflict of interest.