Current and future burden of Ross River virus infection attributable to increasing temperature in Australia: a population-based study.

Background: Ross River virus (RRV), Australia's most notifiable vector-borne disease transmitted through mosquito bites, has seen increased transmission due to rising temperatures. Quantifying the burden of RRV infection attributable to increasing temperatures (both current and future) is pivotal to inform prevention strategies in the context of climate change. Methods: As RRV-related deaths are rare in Australia, we utilised years lived with disability (YLDs) associated with RRV infection data from the Australian Institute of Health and Welfare (AIHW) Burden of Disease database between 2003 and 2018. We obtained relative risks per 1 °C temperature increase in RRV infection from a previous meta-analysis. Exposure distributions for each Köppen-Geiger climate zone were calculated separately and compared with the theoretical-minimum-risk exposure distribution to calculate RRV burden attributable to increasing temperatures during the baseline period (2003-2018), and projected future burdens for the 2030s and 2050s under two greenhouse gas emission scenarios (Representative Concentration Pathways, RCP 4.5 and RCP 8.5), two adaptation scenarios, and different population growth series. Findings: During the baseline period (2003-2018), increasing mean temperatures contributed to 35.8 (±0.5) YLDs (19.1%) of the observed RRV burden in Australia. The mean temperature attributable RRV burden varied across climate zones and jurisdictions. Under both RCP scenarios, the projected RRV burden is estimated to increase in the future despite adaptation scenarios. By the 2050s, without adaptation, the RRV burden could reach 45.8 YLDs under RCP4.5 and 51.1 YLDs under RCP8.5. Implementing a 10% adaptation strategy could reduce RRV burden to 41.8 and 46.4 YLDs, respectively. Interpretation: These findings provide scientific evidence for informing policy decisions and guiding resource allocation for mitigating the future RRV burden. The current findings underscore the need to develop location-specific adaptation strategies for climate-sensitive disease control and prevention. Funding: Australian Research Council Discovery Program.

Improving the spatial and temporal resolution of burden of disease measures with Bayesian models.

This paper contributes to the field by addressing the critical issue of enhancing the spatial and temporal resolution of health data. Although Bayesian methods are frequently employed to address this challenge in various disciplines, the application of Bayesian spatio-temporal models to burden of disease (BOD) studies remains limited. Our novelty lies in the exploration of two existing Bayesian models that we show to be applicable to a wide range of BOD data, including mortality and prevalence, thereby providing evidence to support the adoption of Bayesian modeling in full BOD studies in the future. We illustrate the benefits of Bayesian modeling with an Australian case study involving asthma and coronary heart disease. Our results showcase the effectiveness of Bayesian approaches in increasing the number of small areas for which results are available and improving the reliability and stability of the results compared to using data directly from surveys or administrative sources.

The burden of occupational injury attributable to high temperatures in Australia, 2014-19: a retrospective observational study.

OBJECTIVES: To assess the population health impact of high temperatures on workplace health and safety by estimating the burden of heat-attributable occupational injury in Australia. STUDY DESIGN, SETTING: Retrospective observational study; estimation of burden of occupational injury in Australia attributable to high temperatures during 2014-19, based on Safe Work Australia (work-related traumatic injury fatalities and workers' compensation databases) and Australian Institute of Health and Welfare data (Australian Burden of Disease Study and National Hospital Morbidity databases), and a meta-analysis of climate zone-specific risk data. MAIN OUTCOME MEASURE: Burden of heat-attributable occupational injuries as disability-adjusted life years (DALYs), comprising the numbers of years of life lived with disability (YLDs) and years of life lost (YLLs), nationally, by Köppen-Geiger climate zone, and by state and territory. RESULTS: During 2014-19, an estimated 42 884 years of healthy life were lost to occupational injury, comprising 39 485 YLLs (92.1%) and 3399 YLDs (7.9%), at a rate of 0.80 DALYs per 1000 workers per year. A total of 967 occupational injury-related DALYs were attributable to heat (2.3% of occupational injury-related DALYs), comprising 890 YLLs (92%) and 77 YLDs (8%). By climate zone, the heat-attributable proportion was largest in the tropical Am (12 DALYs; 3.5%) and Aw zones (34 DALYs; 3.5%); by state and territory, the proportion was largest in New South Wales and Queensland (each 2.9%), which also included the largest numbers of heat-attributable occupational injury-related DALYs (NSW: 379 DALYs, 39% of national total; Queensland: 308 DALYs; 32%). CONCLUSION: An estimated 2.3% of the occupational injury burden in Australia is attributable to high ambient temperatures. To prevent this burden increasing with global warming, adaptive measures and industry-based policies are needed to safeguard workplace health and safety, particularly in heat-exposed industries, such as agriculture, transport, and construction.

Projection of high temperature-related burden of kidney disease in Australia under different climate change, population and adaptation scenarios: population-based study.

Background: The dual impacts of a warming climate and population ageing lead to an increasing kidney disease prevalence, highlighting the importance of quantifying the burden of kidney disease (BoKD) attributable to high temperature, yet studies on this subject are limited. The study aims to quantify the BoKD attributable to high temperatures in Australia across all states and territories, and project future BoKD under climatic, population and adaptation scenarios. Methods: Data on disability-adjusted-life-years (DALYs) due to kidney disease, including years of life lost (YLL), and years lived with disability (YLD), were collected during 2003-2018 (baseline) across all states and territories in Australia. The temperature-response association was estimated using a meta-regression model. Future temperature projections were calculated using eight downscaled climate models to estimate changes in attributable BoKD centred around 2030s and 2050s, under two greenhouse gas emissions scenarios (RCP4.5 and RCP8.5), while considering changes in population size and age structure, and human adaptation to climate change. Findings: Over the baseline (2003-2018), high-temperature contributed to 2.7% (Standard Deviation: 0.4%) of the observed BoKD in Australia. The future population attributable fraction and the attributable BoKD, projected using RCP4.5 and RCP8.5, showed a gradually increasing trend when assuming no human adaptation. Future projections were most strongly influenced by the population change, with the high temperature-related BoKD increasing by 18.4-67.4% compared to the baseline under constant population and by 100.2-291.2% when accounting for changes in population size and age structure. However, when human adaptation was adopted (from no to partial to full), the high temperature-related BoKD became smaller. Interpretation: It is expected that increasing high temperature exposure will substantially contribute to higher BoKD across Australia, underscoring the urgent need for public health interventions to mitigate the negative health impacts of a warming climate on BoKD. Funding: Australian Research Council Discovery Program.

The increasing significance of disease severity in a burden of disease framework.

Recent estimates have reiterated that non-fatal causes of disease, such as low back pain, headaches and depressive disorders, are amongst the leading causes of disability-adjusted life years (DALYs). For these causes, the contribution of years lived with disability (YLD) - put simply, ill-health - is what drives DALYs, not mortality. Being able to monitor trends in YLD closely is particularly relevant for countries that sit high on the socio-demographic spectrum of development, as it contributes more than half of all DALYs. There is a paucity of data on how the population-level occurrence of disease is distributed according to severity, and as such, the majority of global and national efforts in monitoring YLD lack the ability to differentiate changes in severity across time and location. This raises uncertainties in interpreting these findings without triangulation with other relevant data sources. Our commentary aims to bring this issue to the forefront for users of burden of disease estimates, as its impact is often easily overlooked as part of the fundamental process of generating DALY estimates. Moreover, the wider health harms of the COVID-19 pandemic have underlined the likelihood of latent and delayed demand in accessing vital health and care services that will ultimately lead to exacerbated disease severity and health outcomes. This places increased importance on attempts to be able to differentiate by both the occurrence and severity of disease.

Quantifying cause-related mortality in Australia, incorporating multiple causes: observed patterns, trends and practical considerations.

BACKGROUND: Mortality statistics using a single underlying cause of death (UC) are key health indicators. Rising multimorbidity and chronic disease mean that deaths increasingly involve multiple conditions. However, additional causes reported on death certificates are rarely integrated into mortality indicators, partly due to complexities in data and methods. This study aimed to assess trends and patterns in cause-related mortality in Australia, integrating multiple causes (MC) of death. METHODS: Deaths (n = 1 773 399) in Australia (2006-17) were mapped to 136 ICD-10-based groups and MC indicators applied. Age-standardized cause-related rates (deaths/100 000) based on the UC (ASRUC) were compared with rates based on any mention of the cause (ASRAM) using rate ratios (RR = ASRAM/ASRUC) and to rates based on weighting multiple contributing causes (ASRW). RESULTS: Deaths involved on average 3.4 causes in 2017; the percentage with >4 causes increased from 20.9 (2006) to 24.4 (2017). Ischaemic heart disease (ASRUC = 73.3, ASRAM = 135.8, ASRW = 63.5), dementia (ASRUC = 51.1, ASRAM = 98.1, ASRW = 52.1) and cerebrovascular diseases (ASRUC = 39.9, ASRAM = 76.7, ASRW = 33.5) ranked as leading causes by all methods. Causes with high RR included hypertension (ASRUC = 2.2, RR = 35.5), atrial fibrillation (ASRUC = 8.0, RR = 6.5) and diabetes (ASRUC = 18.5, RR = 3.5); the corresponding ASRW were 12.5, 12.6 and 24.0, respectively. Renal failure, atrial fibrillation and hypertension ranked among the 10 leading causes by ASRAM and ASRW but not by ASRUC. Practical considerations in working with MC data are discussed. CONCLUSIONS: Despite the similarities in leading causes under the three methods, with integration of MC several preventable diseases emerged as leading causes. MC analyses offer a richer additional perspective for population health monitoring and policy development.

Estimating the burden of disease attributable to high ambient temperature across climate zones: methodological framework with a case study.

BACKGROUND: With high temperature becoming an increasing health risk due to a changing climate, it is important to quantify the scale of the problem. However, estimating the burden of disease (BoD) attributable to high temperature can be challenging due to differences in risk patterns across geographical regions and data accessibility issues. METHODS: We present a methodological framework that uses Köppen-Geiger climate zones to refine exposure levels and quantifies the difference between the burden observed due to high temperatures and what would have been observed if the population had been exposed to the theoretical minimum risk exposure distribution (TMRED). Our proposed method aligned with the Australian Burden of Disease Study and included two parts: (i) estimation of the population attributable fractions (PAF); and then (ii) estimation of the BoD attributable to high temperature. We use suicide and self-inflicted injuries in Australia as an example, with most frequent temperatures (MFTs) as the minimum risk exposure threshold (TMRED). RESULTS: Our proposed framework to estimate the attributable BoD accounts for the importance of geographical variations of risk estimates between climate zones, and can be modified and adapted to other diseases and contexts that may be affected by high temperatures. CONCLUSIONS: As the heat-related BoD may continue to increase in the future, this method is useful in estimating burdens across climate zones. This work may have important implications for preventive health measures, by enhancing the reproducibility and transparency of BoD research.

Heat exposure and cardiovascular health outcomes: a systematic review and meta-analysis.

BACKGROUND: Heat exposure is an important but underappreciated risk factor contributing to cardiovascular disease. Warming temperatures might therefore pose substantial challenges to population health, especially in a rapidly aging population. To address a potential increase in the burden of cardiovascular disease, a better understanding of the effects of ambient heat on different types of cardiovascular disease and factors contributing to vulnerability is required, especially in the context of climate change. This study reviews the current epidemiological evidence linking heat exposures (both high temperatures and heatwaves) with cardiovascular disease outcomes, including mortality and morbidity. METHODS: In this systematic review and meta-analysis, we searched PubMed, Embase, and Scopus for literature published between Jan 1, 1990, and March 10, 2022, and evaluated the quality of the evidence following the Navigation Guide Criteria. We included original research on independent study populations in which the exposure metric was high temperatures or heatwaves, and observational studies using ecological time series, case crossover, or case series study designs comparing risks over different exposures or time periods. Reviews, commentaries, grey literature, and studies that examined only seasonal effects without explicitly considering temperature were excluded. The risk estimates were derived from included articles and if insufficient data were available we contacted the authors to provide clarification. We did a random-effects meta-analysis to pool the relative risk (RR) of the association between high temperatures and heatwaves and cardiovascular disease outcomes. The study protocol was registered with PROSPERO (CRD42021232601). FINDINGS: In total, 7360 results were returned from our search of which we included 282 articles in the systematic review, and of which 266 were eligible for the meta-analysis. There was substantial heterogeneity for both mortality (high temperatures: I2=93·6%, p<0·0001; heatwaves: I2=98·9%, p<0·0001) and morbidity (high temperatures: I2=98·8%, p<0·0001; heatwaves: I2=83·5%, p<0·0001). Despite the heterogeneity in environmental conditions and population dynamics among the reviewed studies, results showed that a 1°C increase in temperature was positively associated with cardiovascular disease-related mortality across all considered diagnoses. The overall risk of cardiovascular disease-related mortality increased by 2·1% (RR 1·021 [95%CI 1·020-1·023]), with the highest specific disease risk being for stroke and coronary heart disease. A 1°C temperature rise was also associated with a significant increase in morbidity due to arrhythmias and cardiac arrest and coronary heart disease. Our findings suggest heat exposure leads to elevated risk of morbidity and mortality for women, people 65 years and older, individuals living in tropical climates, and those in countries of lower-middle income. Heatwaves were also significantly associated with a 17% increase in risk of mortality (RR 1·117 [95% CI 1·093-1·141]), and increasing heatwave intensity with an increasing risk (RR 1·067 [95% CI 1·056-1·078] for low intensity, 1·088 [1·058-1·119] for middle intensity, and 1·189 [1·109-1·269] for high intensity settings). INTERPRETATION: This review strengthens the evidence on the increase in cardiovascular disease risk due to ambient heat exposures in different climate zones. The widespread prevalence of exposure to hot temperatures, in conjunction with an increase in the proportion of older people in the population, might result in a rise in poor cardiovascular disease health outcomes associated with a warming climate. Evidence-based prevention measures are needed to attenuate peaks in cardiovascular events during hot spells, thereby lowering the worldwide total heat-related burden of cardiovascular disease-related morbidity and death. FUNDING: Australian Research Council Discovery Program.

Education-related inequalities in cause-specific mortality: first estimates for Australia using individual-level linked census and mortality data.

BACKGROUND: Socioeconomic inequalities in mortality are evident in all high-income countries, and ongoing monitoring is recommended using linked census-mortality data. Using such data, we provide the first estimates of education-related inequalities in cause-specific mortality in Australia, suitable for international comparisons. METHODS: We used Australian Census (2016) linked to 13 months of Death Registrations (2016-17). We estimated relative rates (RR) and rate differences (RD, per 100 000 person-years), comparing rates in low (no qualifications) and intermediate (secondary school) with high (tertiary) education for individual causes of death (among those aged 25-84 years) and grouped according to preventability (25-74 years), separately by sex and age group, adjusting for age, using negative binomial regression. RESULTS: Among 13.9 M people contributing 14 452 732 person-years, 84 743 deaths occurred. All-cause mortality rates among men and women aged 25-84 years with low education were 2.76 [95% confidence interval (CI): 2.61-2.91] and 2.13 (2.01-2.26) times the rates of those with high education, respectively. We observed inequalities in most causes of death in each age-sex group. Among men aged 25-44 years, relative and absolute inequalities were largest for injuries, e.g. transport accidents [RR = 10.1 (5.4-18.7), RD = 21.2 (14.5-27.9)]). Among those aged 45-64 years, inequalities were greatest for chronic diseases, e.g. lung cancer [men RR = 6.6 (4.9-8.9), RD = 57.7 (49.7-65.8)] and ischaemic heart disease [women RR = 5.8 (3.7-9.1), RD = 20.2 (15.8-24.6)], with similar patterns for people aged 65-84 years. When grouped according to preventability, inequalities were large for causes amenable to behaviour change and medical intervention for all ages and causes amenable to injury prevention among young men. CONCLUSIONS: Australian education-related inequalities in mortality are substantial, generally higher than international estimates, and related to preventability. Findings highlight opportunities to reduce them and the potential to improve the health of the population.

Value of a national burden-of-disease study: a comparison of estimates between the Australian Burden of Disease Study 2015 and the Global Burden of Disease Study 2017.

BACKGROUND: Estimates of burden of disease are important for monitoring population health, informing policy and service planning. Burden estimates for the same population can be reported differently by national studies [e.g. the Australian Burden of Disease Study (ABDS) and the Global Burden of Disease Study (GBDS)]. METHODS: Australian ABDS 2015 and GBDS 2017 burden estimates and methods for 2015 were compared. Years of life lost (YLL), years lived with disability (YLD) and disability-adjusted life years (DALY) measures were compared for overall burden and 'top 50' causes. Disease-category definitions (based on ICD-10), redistribution algorithms, data sources, disability weights, modelling methods and assumptions were reviewed. RESULTS: GBDS 2017 estimated higher totals than ABDS 2015 for YLL, YLD and DALY for Australia. YLL differences were mainly driven by differences in the allocation of deaths to disease categories and the redistribution of implausible causes of death. For YLD, the main drivers were data sources, severity distributions and modelling strategies. Most top-50 diseases for DALY had a similar YLL:YLD composition reported. CONCLUSIONS: Differences in the ABDS and GBDS estimates reflect the different purposes of local and international studies and differences in data and modelling strategies. The GBDS uses all available evidence and is useful for international comparisons. National studies such as the ABDS have the flexibility to meet local needs and often the advantage of access to unpublished data. It is important that all data sources, inputs and models be assessed for quality and appropriateness. As studies evolve, differences should be accounted for through increased transparency of data and methods.

Hot weather as a risk factor for kidney disease outcomes: A systematic review and meta-analysis of epidemiological evidence.

BACKGROUND: The occurrence or exacerbation of kidney disease has been documented as a growing problem associated with hot weather. The implementation of effective prevention measures requires a better understanding of the risk factors that increase susceptibility. To fill gaps in knowledge, this study reviews the current literature on the effects of heat on kidney-disease outcomes (ICD-10 N00-N39), including morbidity and mortality. METHODS: Databases were systematically searched for relevant literature published between 1990 and 2020 and the quality of evidence evaluated. We performed random effects meta-analysis to calculate the pooled relative risks (RRs) of the association between high temperatures (and heatwaves) and kidney disease outcomes. We further evaluated vulnerability concerning contextual population characteristics. RESULTS: Of 2739 studies identified, 91 were reviewed and 82 of these studies met the criteria for inclusion in a meta-analysis. Findings showed that with a 1 °C increase in temperature, the risk of kidney-related morbidity increased by 1% (RR 1.010; 95% CI: 1.009-1.011), with the greatest risk for urolithiasis. Heatwaves were also associated with increased morbidity with a trend observed with heatwave intensity. During low-intensity heatwaves, there was an increase of 5.9% in morbidity, while during high-intensity heatwaves there was a 7.7% increase. There were greater RRs for males, people aged ≤64 years, and those living in temperate climate zones. Similarly, for every 1 °C temperature increase, there was a 3% (RR 1.031; 95% CI: 1.018-1.045) increase in the risk of kidney-related mortality, which also increased during heatwaves. CONCLUSIONS: High temperatures (and heatwaves) are associated with an elevated risk of kidney disease outcomes, particularly urolithiasis. Preventive measures that may minimize risks in vulnerable individuals during hot spells are discussed.

Is there an association between hot weather and poor mental health outcomes? A systematic review and meta-analysis.

BACKGROUND: Mental health is an important public health issue globally. A potential link between heat exposure and mental health outcomes has been recognised in the scientific literature; however, the associations between heat exposure (both high ambient temperatures and heatwaves) and mental health-related mortality and morbidity vary between studies and locations. OBJECTIVE: To fill gaps in knowledge, this systematic review aims to summarize the epidemiological evidence and investigate the quantitative effects of high ambient temperatures and heatwaves on mental health-related mortality and morbidity outcomes, while exploring sources of heterogeneity. METHODS: A systematic search of peer-reviewed epidemiological studies on heat exposure and mental health outcomes published between January 1990 and November 2020 was conducted using five databases (PubMed, Embase, Scopus, Web of Science and PsycINFO). We included studies that examined the association between high ambient temperatures and/or heatwaves and mental health-related mortality and morbidity (e.g. hospital admissions and emergency department visits) in the general population. A range of mental health conditions were defined using ICD-10 classifications. We performed random effects meta-analysis to summarize the relative risks (RRs) in mental health outcomes per 1 °C increase in temperature, and under different heatwaves definitions. We further evaluated whether variables such as age, sex, socioeconomic status, and climate zone may explain the observed heterogeneity. RESULTS: The keyword search yielded 4560 citations from which we identified 53 high temperatures/heatwaves studies that comprised over 1.7 million mental health-related mortality and 1.9 million morbidity cases in total. Our findings suggest associations between heat exposures and a range of mental health-related outcomes. Regarding high temperatures, our meta-analysis of study findings showed that for each 1 °C increase in temperature, the mental health-related mortality and morbidity increased with a RR of 1.022 (95%CI: 1.015-1.029) and 1.009 (95%CI: 1.007-1.015), respectively. The greatest mortality risk was attributed to substance-related mental disorders (RR, 1.046; 95%CI: 0.991-1.101), followed by organic mental disorders (RR, 1.033; 95%CI: 1.020-1.046). A 1 °C temperature rise was also associated with a significant increase in morbidity such as mood disorders, organic mental disorders, schizophrenia, neurotic and anxiety disorders. Findings suggest evidence of vulnerability for populations living in tropical and subtropical climate zones, and for people aged more than 65 years. There were significant moderate and high heterogeneities between effect estimates in overall mortality and morbidity categories, respectively. Lower heterogeneity was noted in some subgroups. The magnitude of the effect estimates for heatwaves varied depending on definitions used. The highest effect estimates for mental health-related morbidity was observed when heatwaves were defined as "mean temperature ≥90th percentile for ≥3 days" (RR, 1.753; 95%CI: 0.567-5.421), and a significant effect was also observed when the definition was "mean temperature ≥95th percentile for ≥3 days", with a RR of 1.064 (95%CI: 1.006-1.123). CONCLUSIONS: Our findings support the hypothesis of a positive association between elevated ambient temperatures and/or heatwaves and adverse mental health outcomes. This problem will likely increase with a warming climate, especially in the context of climate change. Further high-quality studies are needed to identify modifying factors of heat impacts.

History and development of national burden of disease assessment in Australia.

Australia's 1996 national burden of disease (BoD) study was one of the first in the world and updates have continued over the following two decades with the fifth study now underway. The studies adapt the global framework most recently implemented by the Global Burden of Disease Study and the World Health Organization to suit Australia's specific needs, producing estimates of fatal and non-fatal burden via the Disability Adjusted Life Year (DALY) metric, as well as attribution of the burden to many risk factors. Detailed Australian data are used with minimal reliance on modelling to fill data gaps. Comprehensive estimates are produced, including for the Indigenous population, for each of the eight states and territories, the five remoteness areas and five socioeconomic quintiles. A number of method developments have been made as part of these studies, including redistribution of deaths data and a detailed quality framework for describing the robustness of the underlying data and methods. Data and methods continue to be refined as part of the studies, and developments in global studies and other national studies are incorporated where appropriate.

Reflections on key methodological decisions in national burden of disease assessments.

BACKGROUND: Summary measures of population health are increasingly used in different public health reporting systems for setting priorities for health care and social service delivery and planning. Disability-adjusted life years (DALYs) are one of the most commonly used health gap summary measures in the field of public health and have become the key metric for quantifying burden of disease (BoD). BoD methodology is, however, complex and highly data demanding, requiring a substantial capacity to apply, which has led to major disparities across researchers and nations in their resources to perform themselves BoD studies and interpret the soundness of available estimates produced by the Global Burden of Disease Study. METHODS: BoD researchers from the COST Action European Burden of Disease network reflect on the most important methodological choices to be made when estimating DALYs. The paper provides an overview of eleven methodological decisions and challenges drawing on the experiences of countries working with BoD methodology in their own national studies. Each of these steps are briefly described and, where appropriate, some examples are provided from different BoD studies across the world. RESULTS: In this review article we have identified some of the key methodological choices and challenges that are important to understand when calculating BoD metrics. We have provided examples from different BoD studies that have developed their own strategies in data usage and implementation of statistical methods in the production of BoD estimates. CONCLUSIONS: With the increase in national BoD studies developing their own strategies in data usage and implementation of statistical methods in the production of BoD estimates, there is a pressing need for equitable capacity building on the one hand, and harmonization of methods on the other hand. In response to these issues, several BoD networks have emerged in the European region that bring together expertise across different domains and professional backgrounds. An intensive exchange in the experience of the researchers in the different countries will enable the understanding of the methods and the interpretation of the results from the local authorities who can effectively integrate the BoD estimates in public health policies, intervention and prevention programs.

Decoupling of host-symbiont-phage coadaptations following transfer between insect species.

Transferring endosymbiotic bacteria between different host species can perturb the coordinated regulation of the host and bacterial genomes. Here we use the most common maternally transmitted bacteria, Wolbachia pipientis, to test the consequences of host genetic background on infection densities and the processes underlying those changes in the parasitoid wasp genus Nasonia. Introgressing the genome of Nasonia giraulti into the infected cytoplasm of N. vitripennis causes a two-order-of-magnitude increase in bacterial loads in adults and a proliferation of the infection to somatic tissues. The host effect on W. pipientis distribution and densities is associated with a twofold decrease in densities of the temperate phage WO-B. Returning the bacteria from the new host species back to the resident host species restores the bacteria and phage to their native densities. To our knowledge, this is the first study to report a host-microbe genetic interaction that affects the densities of both W. pipientis and bacteriophage WO-B. The consequences of the increased bacterial density include a reduction in fecundity, an increase in levels of cytoplasmic incompatibility (CI), and unexpectedly, male-to-female transfer of the bacteria to uninfected females and an increased acceptance of densely infected females to interspecific mates. While paternal inheritance of the W. pipientis was not observed, the high incidence of male-to-female transfer in the introgressed background raises the possibility that paternal transmission could be more likely in hybrids where paternal leakage of other cytoplasmic elements is also known to occur. Taken together, these results establish a major change in W. pipientis densities and tissue tropism between closely related species and support a model in which phage WO, Wolbachia, and arthropods form a tripartite symbiotic association in which all three are integral to understanding the biology of this widespread endosymbiosis.