The Burden of Proof studies: assessing the evidence of risk

Exposure to risks throughout life results in a wide variety of outcomes. Objectively judging the relative impact of these risks on personal and population health is fundamental to individual survival and societal prosperity. Existing mechanisms to quantify and rank the magnitude of these myriad effects and the uncertainty in their estimation are largely subjective, leaving room for interpretation that can fuel academic controversy and add to confusion when communicating risk. We present a new suite of meta-analyses-termed the Burden of Proof studies-designed specifically to help evaluate these methodological issues objectively and quantitatively. Through this data-driven approach that complements existing systems, including GRADE and Cochrane Reviews, we aim to aggregate evidence across multiple studies and enable a quantitative comparison of risk-outcome pairs. We introduce the burden of proof risk function (BPRF), which estimates the level of risk closest to the null hypothesis that is consistent with available data. Here we illustrate the BPRF methodology for the evaluation of four exemplar risk-outcome pairs: smoking and lung cancer, systolic blood pressure and ischemic heart disease, vegetable consumption and ischemic heart disease, and unprocessed red meat consumption and ischemic heart disease. The strength of evidence for each relationship is assessed by computing and summarizing the BPRF, and then translating the summary to a simple star rating. The Burden of Proof methodology provides a consistent way to understand, evaluate and summarize evidence of risk across different risk-outcome pairs, and informs risk analysis conducted as part of the Global Burden of Diseases, Injuries, and Risk Factors Study.

SARS-CoV-2 Epidemiology on a Public University Campus in Washington State

BACKGROUND

We aimed to evaluate a testing program to facilitate control of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission at a large university and measure spread in the university community using viral genome sequencing.

METHODS

Our prospective longitudinal study used remote contactless enrollment, daily mobile symptom and exposure tracking, and self-swab sample collection. Individuals were tested if the participant was exposed to a known SARS-CoV-2-infected person, developed new symptoms, or reported high-risk behavior (such as attending an indoor gathering without masking or social distancing), if a member of a group experiencing an outbreak, or at enrollment. Study participants included students, staff, and faculty at an urban public university during the Autumn quarter of 2020.

RESULTS

We enrolled 16 476 individuals, performed 29 783 SARS-CoV-2 tests, and detected 236 infections. Seventy-five percent of positive cases reported at least 1 of the following: symptoms (60.8%), exposure (34.7%), or high-risk behaviors (21.5%). Greek community affiliation was the strongest risk factor for testing positive, and molecular epidemiology results suggest that specific large gatherings were responsible for several outbreaks.

CONCLUSIONS

A testing program focused on individuals with symptoms and unvaccinated persons who participate in large campus gatherings may be effective as part of a comprehensive university-wide mitigation strategy to control the spread of SARS-CoV-2.

Pathway-Based Approaches for Assessing Biological Hazards of Complex Mixtures of Contaminants: A Case Study in the Maumee River

Assessment of ecological risks of chemicals in the field usually involves complex mixtures of known and unknown compounds. We describe the use of pathway-based chemical and biological approaches to assess the risk of chemical mixtures in the Maumee River (OH, USA), which receives a variety of agricultural and urban inputs. Fathead minnows (Pimephales promelas) were deployed in cages for 4 d at a gradient of sites along the river and adjoining tributaries in 2012 and during 2 periods (April and June) in 2016, in conjunction with an automated system to collect composite water samples. More than 100 industrial chemicals, pharmaceuticals, and pesticides were detected in water at some of the study sites, with the greatest number typically found near domestic wastewater treatment plants. In 2016, there was an increase in concentrations of several herbicides from April to June at upstream agricultural sites. A comparison of chemical concentrations in site water with single chemical data from vitro high-throughput screening (HTS) assays suggested the potential for perturbation of multiple biological pathways, including several associated with induction or inhibition of different cytochrome P450 (CYP) isozymes. This was consistent with direct effects of water extracts in an HTS assay and induction of hepatic CYPs in caged fish. Targeted in vitro assays and measurements in the caged fish suggested minimal effects on endocrine function (e.g., estrogenicity). A nontargeted mass spectroscopy-based analysis suggested that hepatic endogenous metabolite profiles in caged fish covaried strongly with the occurrence of pesticides and pesticide degradates. These studies demonstrate the application of an integrated suite of measurements to help understand the effects of complex chemical mixtures in the field. Environ Toxicol Chem 2021;40:1098-1122. © 2020 SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 comparative risk assessment (CRA) is a comprehensive approach to risk factor quantification that offers a useful tool for synthesising evidence on risks and risk-outcome associations. With each annual GBD study, we update the GBD CRA to incorporate improved methods, new risks and risk-outcome pairs, and new data on risk exposure levels and risk-outcome associations.

METHODS

We used the CRA framework developed for previous iterations of GBD to estimate levels and trends in exposure, attributable deaths, and attributable disability-adjusted life-years (DALYs), by age group, sex, year, and location for 84 behavioural, environmental and occupational, and metabolic risks or groups of risks from 1990 to 2017. This study included 476 risk-outcome pairs that met the GBD study criteria for convincing or probable evidence of causation. We extracted relative risk and exposure estimates from 46 749 randomised controlled trials, cohort studies, household surveys, census data, satellite data, and other sources. We used statistical models to pool data, adjust for bias, and incorporate covariates. Using the counterfactual scenario of theoretical minimum risk exposure level (TMREL), we estimated the portion of deaths and DALYs that could be attributed to a given risk. We explored the relationship between development and risk exposure by modelling the relationship between the Socio-demographic Index (SDI) and risk-weighted exposure prevalence and estimated expected levels of exposure and risk-attributable burden by SDI. Finally, we explored temporal changes in risk-attributable DALYs by decomposing those changes into six main component drivers of change as follows: (1) population growth; (2) changes in population age structures; (3) changes in exposure to environmental and occupational risks; (4) changes in exposure to behavioural risks; (5) changes in exposure to metabolic risks; and (6) changes due to all other factors, approximated as the risk-deleted death and DALY rates, where the risk-deleted rate is the rate that would be observed had we reduced the exposure levels to the TMREL for all risk factors included in GBD 2017.

FINDINGS

In 2017, 34·1 million (95% uncertainty interval [UI] 33·3-35·0) deaths and 1·21 billion (1·14-1·28) DALYs were attributable to GBD risk factors. Globally, 61·0% (59·6-62·4) of deaths and 48·3% (46·3-50·2) of DALYs were attributed to the GBD 2017 risk factors. When ranked by risk-attributable DALYs, high systolic blood pressure (SBP) was the leading risk factor, accounting for 10·4 million (9·39-11·5) deaths and 218 million (198-237) DALYs, followed by smoking (7·10 million [6·83-7·37] deaths and 182 million [173-193] DALYs), high fasting plasma glucose (6·53 million [5·23-8·23] deaths and 171 million [144-201] DALYs), high body-mass index (BMI; 4·72 million [2·99-6·70] deaths and 148 million [98·6-202] DALYs), and short gestation for birthweight (1·43 million [1·36-1·51] deaths and 139 million [131-147] DALYs). In total, risk-attributable DALYs declined by 4·9% (3·3-6·5) between 2007 and 2017. In the absence of demographic changes (ie, population growth and ageing), changes in risk exposure and risk-deleted DALYs would have led to a 23·5% decline in DALYs during that period. Conversely, in the absence of changes in risk exposure and risk-deleted DALYs, demographic changes would have led to an 18·6% increase in DALYs during that period. The ratios of observed risk exposure levels to exposure levels expected based on SDI (O/E ratios) increased globally for unsafe drinking water and household air pollution between 1990 and 2017. This result suggests that development is occurring more rapidly than are changes in the underlying risk structure in a population. Conversely, nearly universal declines in O/E ratios for smoking and alcohol use indicate that, for a given SDI, exposure to these risks is declining. In 2017, the leading Level 4 risk factor for age-standardised DALY rates was high SBP in four super-regions: central Europe, eastern Europe, and central Asia; north Africa and Middle East; south Asia; and southeast Asia, east Asia, and Oceania. The leading risk factor in the high-income super-region was smoking, in Latin America and Caribbean was high BMI, and in sub-Saharan Africa was unsafe sex. O/E ratios for unsafe sex in sub-Saharan Africa were notably high, and those for alcohol use in north Africa and the Middle East were notably low.

INTERPRETATION

By quantifying levels and trends in exposures to risk factors and the resulting disease burden, this assessment offers insight into where past policy and programme efforts might have been successful and highlights current priorities for public health action. Decreases in behavioural, environmental, and occupational risks have largely offset the effects of population growth and ageing, in relation to trends in absolute burden. Conversely, the combination of increasing metabolic risks and population ageing will probably continue to drive the increasing trends in non-communicable diseases at the global level, which presents both a public health challenge and opportunity. We see considerable spatiotemporal heterogeneity in levels of risk exposure and risk-attributable burden. Although levels of development underlie some of this heterogeneity, O/E ratios show risks for which countries are overperforming or underperforming relative to their level of development. As such, these ratios provide a benchmarking tool to help to focus local decision making. Our findings reinforce the importance of both risk exposure monitoring and epidemiological research to assess causal connections between risks and health outcomes, and they highlight the usefulness of the GBD study in synthesising data to draw comprehensive and robust conclusions that help to inform good policy and strategic health planning.

FUNDING

Bill & Melinda Gates Foundation.

Population and fertility by age and sex for 195 countries and territories, 1950-2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND

Population estimates underpin demographic and epidemiological research and are used to track progress on numerous international indicators of health and development. To date, internationally available estimates of population and fertility, although useful, have not been produced with transparent and replicable methods and do not use standardised estimates of mortality. We present single-calendar year and single-year of age estimates of fertility and population by sex with standardised and replicable methods.

METHODS

We estimated population in 195 locations by single year of age and single calendar year from 1950 to 2017 with standardised and replicable methods. We based the estimates on the demographic balancing equation, with inputs of fertility, mortality, population, and migration data. Fertility data came from 7817 location-years of vital registration data, 429 surveys reporting complete birth histories, and 977 surveys and censuses reporting summary birth histories. We estimated age-specific fertility rates (ASFRs; the annual number of livebirths to women of a specified age group per 1000 women in that age group) by use of spatiotemporal Gaussian process regression and used the ASFRs to estimate total fertility rates (TFRs; the average number of children a woman would bear if she survived through the end of the reproductive age span [age 10-54 years] and experienced at each age a particular set of ASFRs observed in the year of interest). Because of sparse data, fertility at ages 10-14 years and 50-54 years was estimated from data on fertility in women aged 15-19 years and 45-49 years, through use of linear regression. Age-specific mortality data came from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 estimates. Data on population came from 1257 censuses and 761 population registry location-years and were adjusted for underenumeration and age misreporting with standard demographic methods. Migration was estimated with the GBD Bayesian demographic balancing model, after incorporating information about refugee migration into the model prior. Final population estimates used the cohort-component method of population projection, with inputs of fertility, mortality, and migration data. Population uncertainty was estimated by use of out-of-sample predictive validity testing. With these data, we estimated the trends in population by age and sex and in fertility by age between 1950 and 2017 in 195 countries and territories.

FINDINGS

From 1950 to 2017, TFRs decreased by 49·4% (95% uncertainty interval [UI] 46·4-52·0). The TFR decreased from 4·7 livebirths (4·5-4·9) to 2·4 livebirths (2·2-2·5), and the ASFR of mothers aged 10-19 years decreased from 37 livebirths (34-40) to 22 livebirths (19-24) per 1000 women. Despite reductions in the TFR, the global population has been increasing by an average of 83·8 million people per year since 1985. The global population increased by 197·2% (193·3-200·8) since 1950, from 2·6 billion (2·5-2·6) to 7·6 billion (7·4-7·9) people in 2017; much of this increase was in the proportion of the global population in south Asia and sub-Saharan Africa. The global annual rate of population growth increased between 1950 and 1964, when it peaked at 2·0%; this rate then remained nearly constant until 1970 and then decreased to 1·1% in 2017. Population growth rates in the southeast Asia, east Asia, and Oceania GBD super-region decreased from 2·5% in 1963 to 0·7% in 2017, whereas in sub-Saharan Africa, population growth rates were almost at the highest reported levels ever in 2017, when they were at 2·7%. The global average age increased from 26·6 years in 1950 to 32·1 years in 2017, and the proportion of the population that is of working age (age 15-64 years) increased from 59·9% to 65·3%. At the national level, the TFR decreased in all countries and territories between 1950 and 2017; in 2017, TFRs ranged from a low of 1·0 livebirths (95% UI 0·9-1·2) in Cyprus to a high of 7·1 livebirths (6·8-7·4) in Niger. The TFR under age 25 years (TFU25; number of livebirths expected by age 25 years for a hypothetical woman who survived the age group and was exposed to current ASFRs) in 2017 ranged from 0·08 livebirths (0·07-0·09) in South Korea to 2·4 livebirths (2·2-2·6) in Niger, and the TFR over age 30 years (TFO30; number of livebirths expected for a hypothetical woman ageing from 30 to 54 years who survived the age group and was exposed to current ASFRs) ranged from a low of 0·3 livebirths (0·3-0·4) in Puerto Rico to a high of 3·1 livebirths (3·0-3·2) in Niger. TFO30 was higher than TFU25 in 145 countries and territories in 2017. 33 countries had a negative population growth rate from 2010 to 2017, most of which were located in central, eastern, and western Europe, whereas population growth rates of more than 2·0% were seen in 33 of 46 countries in sub-Saharan Africa. In 2017, less than 65% of the national population was of working age in 12 of 34 high-income countries, and less than 50% of the national population was of working age in Mali, Chad, and Niger.

INTERPRETATION

Population trends create demographic dividends and headwinds (ie, economic benefits and detriments) that affect national economies and determine national planning needs. Although TFRs are decreasing, the global population continues to grow as mortality declines, with diverse patterns at the national level and across age groups. To our knowledge, this is the first study to provide transparent and replicable estimates of population and fertility, which can be used to inform decision making and to monitor progress.

FUNDING

Bill & Melinda Gates Foundation.

Genetic and Genomic Response to Selection for Food Consumption in Drosophila melanogaster

Food consumption is an essential component of animal fitness; however, excessive food intake in humans increases risk for many diseases. The roles of neuroendocrine feedback loops, food sensing modalities, and physiological state in regulating food intake are well understood, but not the genetic basis underlying variation in food consumption. Here, we applied ten generations of artificial selection for high and low food consumption in replicate populations of Drosophila melanogaster. The phenotypic response to selection was highly asymmetric, with significant responses only for increased food consumption and minimal correlated responses in body mass and composition. We assessed the molecular correlates of selection responses by DNA and RNA sequencing of the selection lines. The high and low selection lines had variants with significantly divergent allele frequencies within or near 2081 genes and 3526 differentially expressed genes in one or both sexes. A total of 519 genes were both genetically divergent and differentially expressed between the divergent selection lines. We performed functional analyses of the effects of RNAi suppression of gene expression and induced mutations for 27 of these candidate genes that have human orthologs and the strongest statistical support, and confirmed that 25 (93 %) affected the mean and/or variance of food consumption.