The Crystal Structure of a Biological Insulated Transmembrane Molecular Wire.

A growing number of bacteria are recognized to conduct electrons across their cell envelope, and yet molecular details of the mechanisms supporting this process remain unknown. Here, we report the atomic structure of an outer membrane spanning protein complex, MtrAB, that is representative of a protein family known to transport electrons between the interior and exterior environments of phylogenetically and metabolically diverse microorganisms. The structure is revealed as a naturally insulated biomolecular wire possessing a 10-heme cytochrome, MtrA, insulated from the membrane lipidic environment by embedding within a 26 strand ß-barrel formed by MtrB. MtrAB forms an intimate connection with an extracellular 10-heme cytochrome, MtrC, which presents its hemes across a large surface area for electrical contact with extracellular redox partners, including transition metals and electrodes.

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Effectiveness of Influenza Vaccination During Pregnancy Against Laboratory-Confirmed Seasonal Influenza Among Infants Under 6 Months of Age in Ontario, Canada.

BACKGROUND: Randomized trials conducted in low- and middle-income settings demonstrated efficacy of influenza vaccination during pregnancy against influenza infection among infants <6 months of age. However, vaccine effectiveness (VE) estimates from settings with different population characteristics and influenza seasonality remain limited. METHODS: We conducted a test-negative study in Ontario, Canada. All influenza virus tests among infants <6 months from 2010 to 2019 were identified and linked with health databases to ascertain information on maternal-infant dyads. VE was estimated from the odds ratio for influenza vaccination during pregnancy among cases versus controls, computed using logistic regression with adjustment for potential confounders. RESULTS: Among 23 806 infants tested for influenza, 1783 (7.5%) were positive and 1708 (7.2%) were born to mothers vaccinated against influenza during pregnancy. VE against laboratory-confirmed infant influenza infection was 64% (95% confidence interval [CI], 50%-74%). VE was similar by trimester of vaccination (first/second, 66% [95% CI, 40%-80%]; third, 63% [95% CI, 46%-74%]), infant age at testing (0 to <2 months, 63% [95% CI, 46%-75%]; 2 to <6 months, 64% [95% CI, 36%-79%]), and gestational age at birth (≥37 weeks, 64% [95% CI, 50%-75%]; < 37 weeks, 61% [95% CI, 4%-86%]). VE against influenza hospitalization was 67% (95% CI, 50%-78%). CONCLUSIONS: Influenza vaccination during pregnancy offers effective protection to infants <6 months, for whom vaccines are not currently available.

Site-specific cancer mortality after low level exposure to ionizing radiation: Findings from an update of the International Nuclear Workers Study (INWORKS).

A major update to the International Nuclear Workers Study was undertaken that allows us to report updated estimates of associations between radiation and site-specific solid cancer mortality. A cohort of 309,932 nuclear workers employed in France, the United Kingdom, and United States were monitored for external radiation exposure and associations with cancer mortality were quantified as the excess relative rate (ERR) per gray (Gy) using a maximum likelihood and a Markov chain Monte Carlo method (to stabilize estimates via a hierarchical regression). The analysis included 28,089 deaths due to solid cancer, the most common being lung, prostate, and colon cancer. Using maximum likelihood, positive estimates of ERR per Gy were obtained for stomach, colon, rectum, pancreas, peritoneum, larynx, lung, pleura/mesothelioma, bone and connective tissue, skin, prostate, testis, bladder, kidney, thyroid, and residual cancers; negative estimates of ERR per Gy were found cancers of oral cavity and pharynx, esophagus, and ovary. A hierarchical model stabilized site-specific estimates of association, including for lung (ERR per Gy=0.65; 95% credible interval [CrI]: 0.24, 1.07), prostate (ERR per Gy=0.44; 95% CrI: -0.06, 0.91), and colon cancer (ERR per Gy=0.53; 95% CrI: -0.07, 1.11). The results contribute evidence regarding associations between low dose radiation and cancer.

A Genomic Counseling Model for Population-Based Sequencing: A Pre-Post Intervention Study.

BACKGROUND: Novel uses of genome sequencing (GS) present an opportunity for return of results to healthy individuals, prompting the need for scalable genetic counseling strategies. We evaluate the effectiveness of a genomic counseling model (GCM) and explore preferences for GS findings in the general population. METHODS: Participants (N=466) completed GS and our GCM (digital genomics platform and group-based webinar), and indicated results preferences. Surveys were administered pre- (T0) and post- (T1) GCM. Change in knowledge and decisional conflict (DC) were evaluated using paired-sample T and Wilcoxon tests. Factors influencing knowledge and results preferences were evaluated using linear and logistic regression models. RESULTS: Participants were 56% female, 58% white, and 53% ≥40 years of age. Mean knowledge scores increased (Limitations: 3.73 to 5.63; benefits: 3.73 to 5.48, p<0.0001) and DC decreased (-21.9, p<0.0001) at T1 versus T0. Eighty-six percent of participants wished to learn all GS findings at T1 vs 78% at T0 (p<0.0001). Older age, negative/mixed attitudes toward genetics, and greater DC were associated with change in preferences post-intervention. CONCLUSION: In a population-based cohort undergoing GS interested in learning GS findings, our GCM increased knowledge and reduced DC, illustrating the GCM's potential effectiveness for GS counseling in the general population.

Universal Difference-in-Differences for Causal Inference in Epidemiology.

Difference-in-differences is undoubtedly one of the most widely used methods for evaluating the causal effect of an intervention in observational (i.e., nonrandomized) settings. The approach is typically used when pre- and postexposure outcome measurements are available, and one can reasonably assume that the association of the unobserved confounder with the outcome has the same absolute magnitude in the two exposure arms and is constant over time; a so-called parallel trends assumption. The parallel trends assumption may not be credible in many practical settings, for example, if the outcome is binary, a count, or polytomous, as well as when an uncontrolled confounder exhibits nonadditive effects on the distribution of the outcome, even if such effects are constant over time. We introduce an alternative approach that replaces the parallel trends assumption with an odds ratio equi-confounding assumption under which an association between treatment and the potential outcome under no treatment is identified with a well-specified generalized linear model relating the pre-exposure outcome and the exposure. Because the proposed method identifies any causal effect that is conceivably identified in the absence of confounding bias, including nonlinear effects such as quantile treatment effects, the approach is aptly called universal difference-in-differences. We describe and illustrate both fully parametric and more robust semiparametric universal difference-in-differences estimators in a real-world application concerning the causal effects of a Zika virus outbreak on birth rate in Brazil. A supplementary digital video is available at: http://links.lww.com/EDE/C90.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird.

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Longitudinal gut microbiome dynamics in relation to age and senescence in a wild animal population.

In humans, gut microbiome (GM) differences are often correlated with, and sometimes causally implicated in, ageing. However, it is unclear how these findings translate in wild animal populations. Studies that investigate how GM dynamics change within individuals, and with declines in physiological condition, are needed to fully understand links between chronological age, senescence and the GM, but have rarely been done. Here, we use longitudinal data collected from a closed population of Seychelles warblers (Acrocephalus sechellensis) to investigate how bacterial GM alpha diversity, composition and stability are associated with host senescence. We hypothesised that GM diversity and composition will differ, and become more variable, in older adults, particularly in the terminal year prior to death, as the GM becomes increasingly dysregulated due to senescence. However, GM alpha diversity and composition remained largely invariable with respect to adult age and did not differ in an individual's terminal year. Furthermore, there was no evidence that the GM became more heterogenous in senescent age groups (individuals older than 6 years), or in the terminal year. Instead, environmental variables such as season, territory quality and time of day, were the strongest predictors of GM variation in adult Seychelles warblers. These results contrast with studies on humans, captive animal populations and some (but not all) studies on non-human primates, suggesting that GM deterioration may not be a universal hallmark of senescence in wild animal species. Further work is needed to disentangle the factors driving variation in GM-senescence relationships across different host taxa.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

A Feasibility Randomized Controlled Trial of Prehabilitation During Neoadjuvant Chemotherapy for Women with Breast Cancer: A Mixed Methods Study.

BACKGROUND: Limited data exist regarding the role of multimodal prehabilitation during neoadjuvant chemotherapy (NACT) for breast cancer. Determining large trial feasibility and identifying signals of prehabilitation benefit are needed. PATIENTS AND METHODS: We conducted a randomized controlled feasibility trial of multimodal prehabilitation versus usual care during NACT among women diagnosed with non-metastatic breast cancer. Intervention participants received an individualized exercise program, dietetic support, and stress management counseling during NACT. The trial assessed feasibility via rates of recruitment, attrition, adherence, and study-related adverse events. Physical fitness (Six Minute Walk Test, grip strength, anthropometrics) and patient-reported outcomes were assessed at baseline, after NACT completion, and 6 months after surgery as exploratory outcomes, and analyzed using linear mixed effects models. Qualitative data were collected from a subsample to understand feasibility and acceptability of prehabilitation. RESULTS: A total of 72 participants were enrolled from the 123 eligible patients (recruitment rate of 53%). There was a 13% attrition rate and no intervention-related adverse events. Participants in the prehabilitation group had better 6-min walk distance at the post-chemotherapy timepoint [between group difference of 49.43 m, 95% confidence interval (CI) - 118.1, 19.2] and at the post-surgery timepoint (27.3, 95% CI -96.8, 42.2) compared with the control group. Prehabilitation participants reported better quality of life, less fatigue, and improved physical activity levels compared with usual care participants. Interviews revealed that the intervention had a positive impact on the treatment experience. CONCLUSIONS: This study demonstrated feasibility and improvement in physical and psychosocial outcomes. Larger trials assessing intervention efficacy to confirm indications of prehabilitation benefit are warranted.

Coherent O-band WDM transmission of DP-16QAM over a 50-km BDFA-amplified link.

We present wavelength-division multiplexed coherent transmission in an O-band amplified link enabled by bismuth-doped fiber amplifiers (BDFAs). Transmission of 4 × 25 GBd DP-16QAM (4 × 200 Gb/s) is demonstrated over a single span of 50-km length, occupying a bandwidth of 4.7 THz across the wavelengths 1323 nm to 1351 nm.

Transient gas-induced differential refractive index effects in as-drawn hollow core optical fibers.

When a hollow core fiber is drawn, the core and cladding holes within the internal cane geometry are pressurized with an inert gas to enable precise control over the internal microstructure of the fiber and counteract surface tension forces. Primarily by considering the temperature drop as the fiber passes through the furnace and the geometrical transformation of the internal microstructure from preform-to-fiber, we recently established that the gas pressure within the final 'as-drawn' fiber is substantially below atmospheric pressure. We have also established that slight changes in the gas refractive index within the core and surrounding cladding holes induced by changes in gas pressure are sufficient to significantly affect both the modality and loss of the fiber. Here we demonstrate, through both simulations and experimental measurements, that the combination of these effects leads to transient changes in the fiber's attenuation when the fibers are opened to atmosphere post-fabrication. It is important to account for this phenomenon for accurate fiber characterization, particularly when long lengths of fiber are drawn where it could take many weeks for every part of the internal microstructure to reach atmospheric pressure.

Considerations for developing and implementing a safe list for alien taxa.

Many species have been intentionally introduced to new regions for their benefits. Some of these alien species cause damage, others do not (or at least have not yet). There are several approaches to address this problem: prohibit taxa that will cause damage, try to limit damages while preserving benefits, or promote taxa that are safe. In the present article, we unpack the safe list approach, which we define as "a list of taxa alien to the region of interest that are considered of sufficiently low risk of invasion and impact that the taxa can be widely used without concerns of negative impacts." We discuss the potential use of safe lists in the management of biological invasions; disentangle aspects related to the purpose, development, implementation, and impact of safe lists; and provide guidance for those considering to develop and implement such lists.

Life history correlations and trade-offs resulting from selection for dispersal in Tribolium castaneum.

Dispersal is an important facet of the life history of many organisms and is, therefore, subject to selective pressure but does not evolve in isolation. Across nature, there are examples of dispersal syndromes and life history strategies in which suites of traits coevolve and covary with dispersal in combinations that serve to maximize fitness in a given ecological context. The red rust flour beetle, Tribolium castaneum, is a model organism and globally significant post-harvest pest that relies on dispersal to reach new patches of ephemeral habitat. Dispersal behaviour in Tribolium has a strong genetic basis. However, a robust understanding of the relationship between dispersal and other life-history components, which could elucidate evolutionary processes and allow pest managers to control their spread and reduce the impact of infestation, is currently lacking. Here, we use highly replicated lines of T. castaneum previously artificially selected for divergent small-scale dispersal propensity to robustly test several important life history components: reproductive strategy, development time, and longevity. As predicted, we find that a suite of important changes as a result of our selection on dispersal: high dispersal propensity is associated with a lower number of longer mating attempts by males, lower investment in early life reproduction by females, slower development of later-laid offspring, and longer female life span. These findings indicate that correlated intraspecific variation in dispersal and related traits may represent alternative life history strategies in T. castaneum. We therefore suggest that pest management efforts to mitigate the species' agro-economic impact should consider the eco-evolutionary dynamics within multiple life histories. The benefits of doing so could be felt both through improved targeting of efforts to reduce spread and also in forecasting how the selection pressures applied through pest management are likely to affect pest evolution.

Single proxy control.

Negative control variables are sometimes used in nonexperimental studies to detect the presence of confounding by hidden factors. A negative control outcome (NCO) is an outcome that is influenced by unobserved confounders of the exposure effects on the outcome in view, but is not causally impacted by the exposure. Tchetgen Tchetgen (2013) introduced the Control Outcome Calibration Approach (COCA) as a formal NCO counterfactual method to detect and correct for residual confounding bias. For identification, COCA treats the NCO as an error-prone proxy of the treatment-free counterfactual outcome of interest, and involves regressing the NCO on the treatment-free counterfactual, together with a rank-preserving structural model, which assumes a constant individual-level causal effect. In this work, we establish nonparametric COCA identification for the average causal effect for the treated, without requiring rank-preservation, therefore accommodating unrestricted effect heterogeneity across units. This nonparametric identification result has important practical implications, as it provides single-proxy confounding control, in contrast to recently proposed proximal causal inference, which relies for identification on a pair of confounding proxies. For COCA estimation we propose 3 separate strategies: (i) an extended propensity score approach, (ii) an outcome bridge function approach, and (iii) a doubly-robust approach. Finally, we illustrate the proposed methods in an application evaluating the causal impact of a Zika virus outbreak on birth rate in Brazil.

Higher-order evidence.

Higher-order evidence is evidence about evidence. Epidemiologic examples of higher-order evidence include the settings where the study data constitute first-order evidence and estimates of misclassification comprise the second-order evidence (e.g., sensitivity, specificity) of a binary exposure or outcome collected in the main study. While sampling variability in higher-order evidence is typically acknowledged, higher-order evidence is often assumed to be free of measurement error (e.g., gold standard measures). Here we provide two examples, each with multiple scenarios where second-order evidence is imperfectly measured, and this measurement error can either amplify or attenuate standard corrections to first-order evidence. We propose a way to account for such imperfections that requires third-order evidence. Further illustrations and exploration of how higher-order evidence impacts results of epidemiologic studies is warranted.

Tobacco retailer density and its association with birth outcomes in the USA: 2000-2016.

INTRODUCTION: Significant progress has been made in reducing maternal exposure to tobacco smoke and subsequent adverse birth outcomes, however, reductions may require strategies that reduce the availability of tobacco retailers. In this study, we investigated the relationship between tobacco retailer density and birth outcomes across the USA and predicted the potential impact of a tobacco retailer density cap on these outcomes. METHODS: Annual US county (n=3105), rates of preterm birth, low birth weight, small-for-gestational age, all-cause infant mortality and sudden infant death syndrome (SIDS) were calculated using National Vital Statistics System data. Tobacco retailers were identified from the National Establishment Time-Series Database. We used Poisson regression to estimate the effect of capping retailer density at 1.4 retailers per 1000 population, controlling for county demographics and air pollution, using propensity score weighting. RESULTS: Tobacco retailer density was positively associated with most adverse birth outcomes. We estimate that a nationwide cap on tobacco retailer density, implemented in 2016, would have resulted in a reduction of 4275 (95% CI 2210 to 6392) preterm births, 6096 (95% CI 4421 to 7806) small-for-gestational-age births, 3483 (95% CI 2615 to 4378) low birthweight births, 538 (95% CI 345 to 733) all-cause infant deaths and 107 (95% CI 55 to 158) SIDS deaths in that year. CONCLUSION: Higher rates of adverse birth outcomes were seen in counties with high tobacco retailer density compared with those with low density. These results provide further support for regulating tobacco retail density to reduce adverse health outcomes associated with tobacco use.

Long-term ambient sulfur dioxide exposure during gestation and preterm birth in North Carolina, 2003-2015.

BACKGROUND: Coal-fired power plants are major contributors of ambient sulfur dioxide (SO2) air pollution. Epidemiological literature suggests an adverse association between SO2 exposure during gestation and preterm birth (PTB; <37 weeks completed gestation). PTB is strongly associated with infant mortality and increased risk for later life morbidities. OBJECTIVE: We investigated associations between SO2 and PTB in North Carolina and evaluated whether the associations were modified by race/ethnicity. METHODS: We assembled a retrospective, administrative cohort of singleton births in North Carolina from 2003-2015. We used US EPA EQUATES data to assign long-term SO2 gestational exposures to eligible births for the entire pregnancy and by trimester. We used multivariable generalized linear regression to estimate risk differences (RD (95%CI)) per 1-ppb increase in SO2, adjusted for gestational parent education, Medicaid status, marital status, and season of conception. Multi-pollutant models were additionally adjusted for other criteria air co-pollutants (O3, PM2.5, NO2). RESULTS: The median SO2 (24-hour average) across exposure windows was ~1.5 (IQR: 1.8) ppb. The overall baseline risk for PTB was 8,756 per 100,000 live births. When stratified by race/ethnicity, the baseline risk for PTB was 12215, 7824, and 7187 per 100,000 live births among non-Hispanic Black, non-Hispanic white, and Hispanic births, respectively. RDs per 1-ppb increase in SO2 averaged across the entire pregnancy were 317.0 (95%CI: 279.4, 354.5) and 568.2 (95%CI: 500.3, 636.1) per 100,000 live births for single- and multi-pollutant models, respectively. For the PTB multi-pollutant models, we observed similar RDs for non-Hispanic Black participants (669.6 [95%CI: 573.9, 765.2]) and non-Hispanic white participants (635.4 [95%CI: 557.2, 713.6]) with smaller RDs for Hispanic participants (336.8 [95%CI: 241.3, 432.2]). SIGNIFICANCE: The results for our adjusted single- and multi-pollutant models showed adverse associations between SO2 and PTB, with some evidence of effect measure modification by race/ethnicity within subcategories of PTB.

Twenty-five year occupational homicide mortality trends in North Carolina: 1992-2017.

INTRODUCTION: Determining industry of decedents and victim-perpetrator relationships is crucial to inform and evaluate occupational homicide prevention strategies. In this study, we examine occupational homicide rates in North Carolina (NC) by victim characteristics, industry and victim-perpetrator relationship from 1992 to 2017. METHODS: Occupational homicides were identified from records of the NC Office of the Chief Medical Examiner system and the NC death certificates. Sex, age, race, ethnicity, class of worker, manner of death, victim-perpetrator relationship and industry were abstracted. Crude and age-standardised homicide rates were calculated as the number of homicides that occurred at work divided by an estimate of worker-years (w-y). Rate ratios and 95% CIs were calculated, and trends over calendar time in occupational homicide rates were examined overall and by industry. RESULTS: 456 homicides over 111 573 049 w-y were observed. Occupational homicide rates decreased from 0.82 per 100 000 w-y for the period 1992-1995 to 0.21 per 100 000 w-y for the period 2011-2015, but increased to 0.32 per 100 000 w-y in the period 2016-2017. Fifty-five per cent (252) of homicides were perpetrated by strangers. Taxi drivers experienced an occupational homicide rate that was 110 times (95% CI 76.52 to 160.19) the overall occupational homicide rate in NC; however, this rate declined by 76.5% between 1992 and 2017. Disparities were observed among workers 65+ years old, racially and ethnically minoritised workers and self-employed workers. CONCLUSION: Our findings identify industries and worker demographics that experienced high occupational homicide fatality rates. Targeted and tailored mitigation strategies among vulnerable industries and workers are recommended.