Ecological and life history traits are associated with Ross River virus infection among sylvatic mammals in Australia

Background Ross River virus (RRV) is Australia’s most important arbovirus given its annual burden of disease and the relatively large number of Australians at risk for infection. This mosquito-borne arbovirus is also a zoonosis, making its epidemiology and infection ecology complex and cryptic. Our grasp of enzootic, epizootic, and zoonotic RRV transmission dynamics is imprecise largely due to a poor understanding of the role of wild mammalian hosts in the RRV system. Methods The current study applied a piecewise structural equation model (PSEM) toward an interspecific comparison of sylvatic Australian mammals to characterize the ecological and life history profile of species with a history of RRV infection relative to those species with no such history among all wild mammalian species surveyed for RRV infection. The effects of species traits were assessed through multiple causal pathways within the PSEM framework. Results Sylvatic mammalian species with a history of RRV infection tended to express dietary specialization and smaller population density. These species were also characterized by a longer gestation length. Conclusions This study provides the first interspecific comparison of wild mammals for RRV infection and identifies some potential targets for future wildlife surveys into the infection ecology of this important arbovirus. An applied RRV macroecology may prove invaluable to the epidemiological modeling of RRV epidemics across diverse sylvatic landscapes, as well as to the development of human and animal health surveillance systems. Electronic supplementary material The online version of this article (10.1186/s12898-019-0220-5) contains supplementary material, which is available to authorized users.


Background
Ross River virus (RRV) is simultaneously Australia's most significant vector-borne and zoonotic pathogen. There are approximately 5100 cases annually [1], which is tenfold higher than all other zoonoses [2]. Transmission of this alphavirus involves multiple mosquito vectors, which in turn exhibit heterogeneous host preferences in diverse dryland and wetland ecosystems [3]. Moreover, this heterogeneity modulates the dynamics of enzootic, epizootic and zoonotic RRV transmission, imparting a unique infection ecology to each. Furthermore, reservoir and amplification hosts among Australian wildlife have been shown to influence the landscape epidemiology of RRV, particularly with respect to human spillover [4,5]. Nevertheless, the contribution of specific mammalian hosts to RRV epidemiology remains ill-defined despite several decades of human surveillance and wildlife sampling across Australia [1].
Ross River virus infection has been identified definitively in 21 native mammalian hosts in several field surveys and experimental studies [6][7][8][9][10][11][12][13][14]. However, many of these surveys target only single species in highly localized areas. None consider species' ecology and how this may relate to (1) infection dynamics, such as the role of wildlife as maintenance or amplification hosts, or (2) the facilitation of human exposure in anthropogenically altered landscapes [15]. By contrast, identifying the biological and life-history characteristics of wildlife hosts has the potential to offer novel insight into the epidemiology and infection ecology of RRV and may suggest a specific wildlife-human interface as particularly vulnerable to spillover. This is a critical consideration as Australian sylvatic landscapes have undergone and continue to undergo rapid change due to habitat loss [16]. To date there has been no investigation of species-level biological and life history characteristics of RRV wildlife hosts and their comparison to non-hosts in endemic areas. Recent work in this domain has identified associations between specieslevel life history characteristics of wildlife reservoirs for other pathogens and human spillover [17]. This approach may prove useful in identifying traits associated with RRV infection status, thus enabling macroecology to inform sylvatic RRV epidemiology.
Recent studies have explored intriguing differences in life history characteristics related to the ways in which fast versus slow life history may mediate disease systems with high potential for spillover. Key delineators of fast living commonly include (1) high volume reproductive output with lower investment in early development as typified for example by short gestation, short inter-birth intervals, large litters, and early sexual maturity, and (2) a shorter lifespan. Conversely, slow living species would be characterized by lesser reproductive output with greater investment in early development as well as longer lifespans [18]. Some work has shown that fast life history is associated with greater reservoir competence [19,20], and may in turn drive zoonotic transmission [20,21], due to immunological mediation of transmission dynamics [22]. A comprehensive survey of all mammalian hostpathogen systems provided strong evidence confirming this relationship between reservoir hosts and living fast [17]. The current investigation explored whether the RRV system converged with, or diverged from, the typical fast life history of many mammalian host-pathogen systems that has been identified globally.
This study compared biological and life history traits of sylvatic mammalian species with demonstrated RRV infection to the traits of species without evidence of infection. Toward this end and as a causally sound approach to epidemiological inference, interspecific biological and life history characteristics were interrogated with respect to species' RRV infection status using a piecewise structural equation model (PSEM) with phylogenetic generalized linear models (PGLM) as its component structures to account for the taxonomic correlation.

Results
The final PGLM of RRV infection history and species traits is presented in Table (Table 1), wherein increasing specialization and gestation length, and lower population density, were all associated with a positive infection history. The PSEM demonstrated that no additional associations with infection history were mediated through unanticipated indirect effects ( Fig. 1), suggesting that the single PGLM model adequately represented the relationships between species' traits and infection history. The model that emerged among these covariates showed that diet breadth, population density, and gestation length were influential to RRV infection history among sylvatic mammals that have been surveyed for this arbovirus. Additionally, reporting effort was not associated with infection history, thus suggesting that RRV-positive species did not appear to be over-represented in sampling effort among the sylvatic mammals surveyed.

Discussion
This investigation suggests that a narrower repertoire of food sources, indicating relative diet specialization, lower population density, and a tendency toward longer gestation characterize wild mammalian RRV infection history. The extent to which perturbations such has habitat loss or resource provisioning among competent hosts with this profile could enhance their potential to act as bridging species for zoonotic transmission is unknown from this, or previous, studies. Such phenomena were not testable in the current study because a clear delineation between competent and non-competent hosts was not possible. Nevertheless, this should be formally tested in future work through better virological surveillance of wildlife and by testing the epidemiological impact of habitat conservation on the prevention of zoonotic RRV transmission. Interestingly, slow-living, as indicated by longer gestation, emerged from the life history profile of RRV-positive species. However, this will also require validation with future work, particularly since only this single life history characteristic was relevant in the current study and the possibility that slow-living may simply reflect greater opportunity for exposure. The ecological context for the relationships between RRV infection history and diet breadth and population density is not clear from the current study, however it may be that animals with a narrower dietary breadth may expend greater time in foraging, and possibly over a greater area, thereby exposing them to mosquito vectors for longer periods of time [23]. Despite its broader ecological significance, the relationship with dietary specialization is interesting because specialists also tend to be more sensitive to habitat fragmentation [24,25]. While such changes can sometimes result in species' inability to adapt to the altered landscapes, species with low diet breadth have also been shown to exhibit the lowest levels of population decline following habitat fragmentation [26]. As such, habitat fragmentation or resource provisioning may induce population perturbations that lead to the persistence of important maintenance hosts, or increase contact with these hosts, both of which in turn could enhance the potential for enzootic, epizootic and zoonotic transmission if biodiversity is also lost [27]. The inverse association with population density is also interesting, but may seem counterintuitive given that vector-borne diseases are commonly assumed to demonstrate frequency-dependent transmission dynamics [28]. It is worth noting, however that this characterization is generally made without adequate parameterization of local spatial heterogeneity in rates of contact, and frequency-dependence can mask density-dependence unless contact networks are adequately defined [29]. The association between infection history and population density thus may reflect an underappreciation of the importance of local heterogeneity in contact rates for sylvatic RRV transmission, but it is also possible that the association with population density is an artifact of species with lower population density exhibiting greater dietary specialization [30], which exhibited a strong association with infection history in this study. As described above, the results from the current study are insufficient to make definitive claims given the inability to distinguish between competent and non-competent hosts, but rather suggest ecological contexts to explore in future RRV disease ecology fieldwork.
This study also contributes to a growing body of work endeavoring to define zoonotic infection ecology by classifying mammalian hosts according to species' life history. The extent to which species "live fast" or "live slow" may have important implications for immunogenicity and pathogenicity [22], thereby influencing pathogen persistence and shedding and subsequent spillover to novel hosts [19][20][21]. This study provided some evidence that a slower life history, as characterized by longer gestation, is characteristic of RRV infection history among sylvatic mammals. It would appear that the RRV system in the Australian mammalian landscape is more complex with respect to life history than other host-pathogen systems that demonstrate a higher proclivity to infection susceptibility among fast-living mammals [17,20,21]. This contrast notwithstanding, a similar pattern of slow-living life history was identified in another vector-borne arbovirus (Rift Valley fever virus) with significant zoonotic disease burden [31]. The importance of increased early life investment may reflect the significance of species with fewer but extended pulses in reproductive outputs to the seasonal cycling of RRV transmission. For example, increased time to reach maturity may correspond to a longer period of immunological naivety, which may extend opportunities for infectious exposure to vectors to within the critical temporal window that defines the transmission dynamics of the seasonal cycling of RRV. Moreover, it has been shown that life history modulates immune function with fast-living species investing more in the development of innate immunity, while slow-living species invest more in adaptive immunity [32,33], which may further affect the infection susceptibility of mammalian species [20][21][22]. Nevertheless, the association between infection history and slow life history may simply reflect a longer exposure opportunity for slow-life species. As the current study could not evaluate acute infection in these species, the findings require validation by studies that do.
The epidemiology of Australia's most important vectorborne and zoonotic pathogen has heretofore been poorly informed by the infection ecology of its wildlife hosts. While macropods are typically considered the primary drivers of RRV circulation and spillover to humans [3,7,8], the current interspecific interrogation defines an ecological profile of potential wildlife hosts according to their biological and ecological traits, and thus expands the scope for more targeted modeling of zoonotic transmission [15]. Moreover, the current study helps to synthesize an ecological model from previous wildlife surveys, which have frequently been conducted in isolation. Of the wildlife surveys conducted to date, only five were generalist surveys that sampled a broad selection of potential wildlife species for evidence of RRV infection [7,9,11,13,34]. The remaining studies were either part of larger outbreak investigations or targeted single or limited species for study. A recent systematic review of RRV virology parameters in vertebrate hosts identified the highest seroprevalence surveyed in macropods and high viremia in experimentally infected marsupials [5]. However, while high in marsupials, this review also showed no significant difference in viremia between marsupials, placentals, and birds during experimental infection. Moreover, the investigators make a sound argument that placental hosts may also contribute to RRV circulation [5]. The current study largely agrees with this systematic review, highlighting the potential additional contribution of non-native animals to RRV transmission dynamics and therefore cautions that a broader focus may yet be required to definitively articulate enzootic, epizootic, and zoonotic transmission. Nevertheless, the current study focused specifically on wildlife to identify how macroecology may inform the epidemiology of sylvatic RRV transmission. Wildlife species are an important focus in their own right given the rapid incursion of the human population into wildlife habitat and the largely unknown consequences of emerging wildlife-human interfaces on pathogen transmission dynamics. Wildlife species abundance is highly variable across diverse wetland and dryland ecosystems, is often sensitive to anthropogenic alteration of these landscapes, and can modify the presence, diversity, and frequency of vector mosquitoes, with further implications for zoonotic transmission risk [3,35,36].
Beyond what has already been described, this study has some additional limitations that are discussed in greater detail below. First, the number of wildlife species surveyed for RRV is relatively small (n = 46) and therefore some relationships may have been missed due to the diminished power of the sample. Second, as described in "Methods" section, there were considerable missing data for some species' biological and life history characteristics, which could not be included in these analyses. Data were imputed where appropriate (i.e. variables within a threshold of < 70% missingness) using methods previously employed [17]. As such, the imputed data represent a narrower spectrum of Australian mammal biology and life history traits, but it is expected that inferences based on these traits will be more reliable if also incomplete. Third, the classification of infectionpositive versus infection-negative species was based on serology in many instances, which only identifies these animals as susceptible to infection and does not indicate definitively whether these are competent reservoir hosts that may serve to maintain, or amplify, virus circulation. Fourth, while the inferences made here with respect to RRV host traits are ecologically relevant, if the ultimate goal is a more complete description of zoonotic RRV epidemiology then this will require the landscape analysis of these species, their mosquito vectors, and the coincident human cases in real time as a broad application of virology-based RRV surveillance in rural, urban, and peri-urban space. Valuable initial efforts have been made to define the landscape epidemiology of RRV, but only at coarse spatial [4,37] and temporal scale [38,39], and based on limited surveillance data.

Conclusions
The findings described here highlight dietary specialization, population density, and gestation length as characteristic of RRV infection in wildlife, which provides the first such interspecific comparison and offers a richer understanding of RRV infection ecology. This applied RRV macroecology may add substantive value to the epidemiological modeling of RRV epidemics across varied wetland and dryland habitat. However, the value added will depend entirely on the commitment to active monitoring of wildlife host movement and their biotic and abiotic interactions coupled with ongoing vector, animal, and human RRV surveillance. Nevertheless, the contextualization of RRV epidemiology using macroecology may suggest potential One Health initiatives that integrate wildlife conservation with human public health for the maximal benefit of both.

Methods
Forty-six sylvatic mammalian species surveyed for RRV infection were identified in the literature. Species reported as positive for infection by serology, polymerase chain reaction amplification, or isolation of virus were classified as RRV positive species [6-14, 34, 40-43]. The remaining sylvatic mammal species sampled for RRV screening, but which tested negative, were classified as RRV negative species. It should be noted that the broad classification scheme based on the three modalities described above is not able to delineate RRV positive species as viral reservoirs or competent hosts. Rather, it simply marks these species as susceptible to infection. Twenty-one native mammalian species, comprising diprotodonts (11 macropods and 1 possum species), peramelemorphs (3 species), megabats (2 species), and rodents (4 species), were classified as RRV-positive, while 5 introduced mammalian species (brown rats, rabbits, pigs, dogs, and cats) were classified as positive. The full list of species classified as RRV hosts and non-hosts based on their documented infection status in the literature is presented in Additional file 2. Only wild/feral species were considered in this analysis because this study sought to describe the macroecology of sylvatic RRV infection history, which may offer unique insight into the epidemiological context of wilderness landscapes. Moreover, contextualizing zoonotic risk with macroecology may also present opportunities to integrate conservation and public health initiatives.
Species-level biological and life history characteristics were obtained from the PanTHERIA mammalian dataset [44]. The following life history traits were selected as representative metrics of the fast-slow continuum concept of life history [45], and which exhibited low levels of missing data (< 70%) in the PanTHERIA database for the host taxonomies under study: gestation length (months), sexual maturity age (days), weaning age (months), litter size (number of animals per litter), inter-birth interval (months) and maximum longevity (months), i.e. lifespan; whereas adult body mass (kg), population density (animals/km 2 ), home range (km 2 ), and diet breadth were selected as additional ecologically relevant species traits that also exhibited low levels of missing data (Additional file 3). Species-level diet breadth describes the mean number of dietary categories consumed throughout the year, wherein food categories were classified as vertebrate, invertebrate, fruit, flowers/nectar/pollen, leaves/ branches/bark, seeds, grass, and roots/tubers. Any remaining species-level missing data for these characteristics were imputed using a random forest machine learning algorithm, which has been shown to be a robust approach to this application [17,46]. The algorithm was implemented using the rfImpute function in the random-Forest package [47].

Statistical analysis
A piecewise structural equation model (PSEM) was used to identify associations between species' infection history and species' ecology and life history [48]. This approach is considered an optimal application to causal inference since it evaluates the covariance among all independent variables as part of the modeling process. Within the PSEM framework, relationships between covariates and between covariates and the outcome (RRV infection history) are estimated as individual structured equations to quantify and map direct and indirect associations. The implementation of local estimation rather than global estimation for these models also makes this approach more suited to small sample sizes. The PSEM thus comprises a set of equations quantifying causal paths in the system, and has become a popular approach to estimating causality in epidemiological relationships. Phylogenetic generalized linear models (PGLMs) were used as the set of equations for these causal paths in this analysis to account for the phylogenetic correlation [48]. To correct for potential bias introduced by differences in reporting effort across species, the number of individuals sampled per species surveyed (across all surveys) was used to quantify reporting effort and this was included as an additional covariate in these PGLMs [49,50]. Furthermore, because of the relatively low number of species available with RRV infection history (n = 46), the maximum number of variables included in the final PGLM was constrained to ensure efficient parameter estimation and prevent overfitting. As such, 10 bivariate models were fitted wherein each model comprised a simple PGLM assessing the crude association between infection history and each species trait (although these crude models did adjust for reporting effort as described above): gestation length, sexual maturity age, weaning age, litter size, inter-birth interval, maximum longevity, population density, home range, diet breadth, and body mass. Using the Akaike information criterion (AIC), the four best-fitting trait variables were then selected and fitted together to assess the best fitting multivariable PGLM. Body mass was also included by default as this trait has the potential to modify many other ecological traits and is generally considered essential to account for when comparing species traits across orders [45]. Initially variables were log transformed for normal approximation, however this did not change the model associations and so the nontransformed data were retained and presented here for easier interpretation of the results.
The PSEM was then fitted to RRV infection status based on the best-fitting PGLM described above. Fisher's C was used to test for conditional independence (p-value > 0.05 indicates there are no missing paths among variables) and the AIC was used to evaluate goodness-of-fit [51,52]. The piecewiseSEM package [52] in the R statistical environment, v. 3.1.3 [53] was used to implement the PSEM. The ape package was used to quantify the phylogenetic correlation [54,55], and the phylolm package was used to fit the PGLM used in the PSEM [56].