key: cord-0899300-vbp7qg4h
authors: Morrison, Robin E.; Mushimiyimana, Yvonne; Stoinski, Tara S.; Eckardt, Winnie
title: Rapid transmission of respiratory infections within but not between mountain gorilla groups
date: 2021-10-07
journal: Sci Rep
DOI: 10.1038/s41598-021-98969-8
sha: 48b959b7a732a58a56921e4ee41289b43a35fb07
doc_id: 899300
cord_uid: vbp7qg4h

Minimizing disease transmission between humans and wild apes and controlling outbreaks in ape populations is vital to both ape conservation and human health, but information on the transmission of real infections in wild populations is rare. We analyzed respiratory outbreaks in a subpopulation of wild mountain gorillas (Gorilla beringei beringei) between 2004 and 2020. We investigated transmission within groups during 7 outbreaks using social networks based on contact and proximity, and transmission between groups during 15 outbreaks using inter-group encounters, transfers and home range overlap. Patterns of contact and proximity within groups were highly predictable based on gorillas’ age and sex. Disease transmission within groups was rapid with a median estimated basic reproductive number (R0) of 4.18 (min = 1.74, max = 9.42), and transmission was not predicted by the social network. Between groups, encounters and transfers did not appear to have enabled disease transmission and the overlap of groups’ ranges did not predict concurrent outbreaks. Our findings suggest that gorilla social structure, with many strong connections within groups and weak ties between groups, may enable rapid transmission within a group once an infection is present, but limit the transmission of infections between groups.

Human exploitation of the world's natural resources has led to increased pressure on remaining habitats and increased contact between human and animal populations 1 . This has had a devastating effect on global biodiversity and also creates opportunities for cross-species disease transmission, posing considerable dangers to human and animal populations alike [2] [3] [4] . The One Health initiative has highlighted how deeply intertwined the health of human populations, animal populations and their environments are 5, 6 . This is particularly critical in relation to our close evolutionary relatives, the nonhuman great apes, who not only are vulnerable to many human pathogens 7 but have also been the origin of several novel human pathogens [8] [9] [10] [11] . The catastrophic effects of infectious diseases such as Ebola Virus Disease (EVD) and those caused by human respiratory pathogens in wild ape populations have been well documented over the past two decades 7, [12] [13] [14] [15] [16] [17] [18] . Understanding and limiting the transmission of these infectious diseases within ape populations is therefore of high priority when it comes to both the conservation of these critically endangered species 12, 19 and the health of human populations 20, 21 .

Across humans and other animals, the social structure of a population can strongly influence patterns of disease transmission 22, 23 . This is because social interactions can provide opportunities for pathogens to transmit, such as the transmission of respiratory pathogens through droplets or aerosols during close proximity interactions. By understanding the modes of disease transmission and quantifying these social interactions through networks we can estimate how likely an individual is to come into contact with a pathogen and to transmit it [24] [25] [26] . This can provide valuable knowledge to prepare for, prevent or respond to infectious disease outbreaks, and for reducing the prevalence of disease in a population. In chimpanzees, this type of analysis has highlighted the greater risk faced by highly connected community members 27, 28 and indicated that targeted vaccination of the most connected individuals could prevent outbreaks with fewer vaccinations 29 . Despite strong theoretical support linking social structure and disease transmission 24, 25, 30 , the modelling of real disease outbreaks transmitting through social interactions in wild populations is rarely possible as it requires the long-term collection of fine-scale social data in advance of a potential outbreak 26 .

The long-term data of the Dian Fossey Gorilla Fund includes fine-scale social data and health data on the daily presence or absence of grossly observable signs commonly associated with respiratory illness, for individuallyidentified mountain gorillas (Gorilla beringei beringei). This provides a rare opportunity for investigating the transmission of disease in an endangered, wild great ape, using the presence of these signs as a proxy for respiratory infection. Mountain gorilla numbers have been recovering due to intensive conservation strategies 31 .

Predicting network position. From annual social networks constructed for 12 groups between 2004 and 2015, we extracted individual gorillas' eigenvector centrality as an estimate of their exposure to, and capacity to spread infectious disease through physical contact and proximity. Using generalized additive mixed models (GAMMs) with smoothed age terms for both sexes, we found that gorillas' social network positions were strongly predicted by their age and sex (Fig. 1 ). Linear mixed models using age/sex categories confirmed that infants were more central in both network types than all other age/sex categories (Fig. 1, Supplementary Table S1 ). Dominance also influenced adult male gorillas' centrality, with dominant adult males' having higher centrality than subordinate (t = − 13.951, Pnull < 0.001) and blackback males (t = − 7.005, Pnull < 0.001) in the proximity-based network, and higher centrality than subordinate males (t = − 5.675, Pnull < 0.001) in the contact-based network ( Fig. 1, Supplementary Table S1 ).

Within-group disease transmission. We examined within-group transmission in the 7 respiratory outbreaks that took place between 2004 and 2015 (Table 1 , Fig. 2 ). We built SIR (susceptible-infectious-removed) mathematical models to estimate epidemiologic parameters related to transmission. These suggested rapid transmission of respiratory pathogens was occurring, inferring a median basic reproductive number (R0) of 4.18 (min = 1.74, max = 9.42) across all outbreaks (Table 1) .

We used social networks constructed from contact and proximity data in the 3 months either side of the onset of 7 respiratory outbreaks (Fig. 2) to test whether occurrences of signs of infection were consistent with transmission through the social network. Path-based k-tests were used to calculate the mean path-length between individuals showing signs of respiratory infection in the social network at a given time point. This allowed us to determine whether individuals showing signs of infection were more closely connected to each other than expected by chance-evidence that transmission had occurred through the network. Across the six outbreaks in which not all group members showed signs of infection, those that showed signs were not more closely connected to each other in the social network than expected by chance ( Table 2 ). Those that showed signs were also no more central in their social network than expected by chance (linear mixed models, contact: est = 0.075 ± 0.717, z = 0.104, p = 0.917; proximity: est = − 0.243 ± 0.975, z = − 0.249, p = 0.803). Since it is possible that all group members were ultimately exposed to an infection, the presence or absence of signs of respiratory infection by the end of an outbreak might be determined entirely by an individual's immunity rather than the transmission pattern itself. We therefore examined the onset of signs more closely by running a path-based k-test at each outbreak time point-each day of the outbreak on which at least one group member showed signs of infection for the first time. If disease was transmitted through the social network but ultimately determined by individual immunity, we might expect a strong signal for transmission through the social network early in the outbreak, decreasing as all individuals were ultimately exposed. This pattern was only observed in the PAB Aug 14 outbreak (minimum p values on Day 1: p = 0.004 for proximity network and p = 0.116 for contact network) and the PAB Nov 06 outbreak (minimum p value on day 2: p = 0.022 for proximity network and p = 0.012 for contact network), out of the five www.nature.com/scientificreports/ outbreaks in which more than 2 infection time points occurred (Fig. 3) . Given the number of tests conducted, the support for transmission of respiratory pathogens through these social networks overall is fairly low. Examining across all 15 outbreaks (including those after 2015 for which detailed social data were not available), infants, juveniles and subadult females were the least likely to show any signs of respiratory infection (Fig. 4, Supplementary Table S2 ). Adults were the most likely to show signs, with the highest prevalence seen in dominant adult males. After including an individual's age/sex category to account for the potential for signs of infection to be more easily observed in some age/sex categories, network centrality for either network type in the smaller subset of outbreaks still did not predict whether a gorilla would show signs of infection during an outbreak (Supplementary Table S3 ).

Between-group transmission. Between 2004 and 2020, 15 respiratory outbreaks took place within the monitored groups (Fig. 5) . These fell into 9 distinct outbreaks at the population level, based on whether outbreaks occurred in multiple groups within a month of each other (referred to as concurrent outbreaks). Six outbreaks did not occur within a month of those in any other group. Outbreaks occurred once in two groups concurrently, once in three groups concurrently and once in four groups concurrently. This temporal overlap of www.nature.com/scientificreports/ outbreaks across groups could have occurred through inter-group transmission, through independent introductions from a common source or from independent introductions from multiple sources 17 .

To investigate the potential for inter-group transmission, we examined transfers and encounters between groups. At least one transfer of a gorilla from an infected to an uninfected group took place during 3 out of 15 outbreaks (20%) (Supplementary Information). However, none of these gorillas were showing signs of respiratory infection at the time they transferred into these uninfected groups. Three transfer events took place from a group during an outbreak in which the dominant (and only) adult male died from the infection, resulting in group disintegration. At least one close proximity (< 50 m) encounter with another group occurred during 4 out of 15 group outbreaks (27%). In two of these cases, there was physical contact between members of different groups. In the first case, aggressive physical contact was made between two adult males from the infected group and an adult male of another group. One of these males had not shown any signs of infection during the outbreak and the other had not shown any in 4 days. In the second case, physical contact was made during play behavior between a subadult female from the infected group and a juvenile male from another group. The subadult female had not shown any signs of infection during the outbreak. None of the transfer events or inter-group encounters appear to have led to respiratory disease transmission between groups, as no consistent signs of infection were observed in any individuals in the previously uninfected groups within a month of the event.

As data on inter-group encounters are limited to when observers are present, range overlap during the outbreak period may provide an alternative measure of the probability of two groups encountering one another, and thus, opportunities for inter-group disease transmission. It also provides an estimate for potential environmental transmission, e.g. through feeding at the same location in close succession 46 . In cases of concurrent outbreaks, the ranges of pairs of groups which both had an outbreak did not overlap more than pairs of groups where only one of the pair had an outbreak (binomial generalized linear model: estimate = − 4.954 ± 4.583, z-value = − 1.081, Table 2 . Analysis of respiratory disease transmission during outbreaks (group code, month and year) between 2004 and 2015, excluding those where signs of infection were observed in every group member. Mean weighted path length (MWPL) calculated between all group members to show signs of infection during an outbreak. P null is the proportion of null models for which the mean weighted path length between individuals showing signs of infection was shorter than that observed in the real outbreak. For network diagrams of the distribution of signs across each network see Fig. 2 Figure S1 , n = 67 pairs of groups). Not only was the result non-significant but the estimated effect was in the opposite direction than expected if inter-group transmission was occurring, with greater range overlap between groups that didn't have concurrent outbreaks (Supplementary Figure S1 ). This lack of support for spatial overlap predicting transmission suggests that inter-group transmission of these respiratory pathogens through any of our proposed mechanisms is unlikely to be a cause of these concurrent outbreaks.

We identified strong, predictable variation in gorillas' network centrality based on their age and sex, and on adult males' dominance. This means that age, sex and adult male dominance could be used to estimate individuals' risk of exposure to infectious diseases and their risk of transmitting them. Selective approaches to minimizing disease transmission such as targeted vaccination strategies might therefore provide a valuable conservation tool in some circumstances. However, our analyses of respiratory outbreaks within these gorilla groups found only limited evidence that transmission was occurring through either the proximity or contact social networks www.nature.com/scientificreports/ ( Table 2) . Analysis of two outbreaks (out of a total of 5 outbreaks analysed) showed support for transmission through the network early in the outbreak (Fig. 3) . But individuals that ultimately became infected during an outbreak were not more closely connected with each other in either type of network than expected by chance and an individual's centrality in the social network did not predict whether they would become infected. This differs considerably from what has been observed in chimpanzees 27, 28 and many other species [22] [23] [24] [47] [48] [49] in which social network position can predict infectious disease risk. This could be for a number of reasons. A key caveat to this study, is that data on the true infection status of individual gorillas were not available. Infection status across this wild and endangered population could only be estimated using the presence or absence of observable signs commonly associated with respiratory illness-an imperfect measure of infection status 50 . Gorillas may vary in the onset of these signs relative to when they are exposed to the infection or in the detectability of these signs. For example, in these outbreaks, signs of infection were more commonly observed in adults, a pattern also observed in chimpanzee respiratory infections 15 . Adults may truly have higher rates of infection or could just be developing more detectable signs. However, even after accounting for age/sex categories and the potentially different detectability of signs within them, network centrality did not predict the presence of signs of respiratory infection. Variation in the immunity of group members could also play a part in obscuring any relationship between sociality and transmission 51, 52 . If all group members are rapidly exposed to the infection, the presence or absence of signs may have more to do with their susceptibility than their exposure.

The density of connections within gorilla groups (Fig. 2 , Supplementary Videos) also hints at another reason why we may not be able to clearly predict disease transmission through these networks. When networks are so highly interconnected, fine-scale variation in rates of proximity and physical contact could have little real-world consequence, with all group members likely to be rapidly exposed 53 . These networks are also only based on interactions during the daytime and cannot account for potential transmission pathways during the night when group members nest in close proximity 54 . Finally, the relatively small group sizes of mountain gorillas may be limiting the power of such analyses compared to species living in larger social units such as chimpanzees 28 . Both outbreaks where transmission initially appeared to follow the social network were in groups of > 30 gorillas (54 and 31 group members), but no signal for transmission through the network was detected in groups formed of 8, 14 or 41 gorillas. Whilst there is no strong support for transmission following within-group social networks, this is by no means evidence against sociality influencing respiratory disease transmission in gorillas more broadly. At a population scale, sociality clearly remains an important predictor. Most outbreaks occurred in only one group at a time (Fig. 5) , with individuals in an infected group at far greater risk of becoming infected than individuals in neighbouring social groups.

Within-group transmission for the majority of outbreaks was extremely rapid. For example, in one outbreak, all 7 members of a group were showing signs of respiratory infection within 3 days (UGE Jul 14). In another, it took only 3 days for 45 out of 46 group members to show signs of infection (PAB Dec 10). Estimated basic reproductive numbers for outbreaks ranged between 1.74 and 9.42. For comparison, an outbreak of human metapneumovirus-one of the pathogens previously found to be responsible for respiratory outbreaks in the population, had a R0 of 1.27 in the Ngogo chimpanzee population 17 . Our findings suggest that the dense connections within gorilla groups may enable respiratory infections to spread so rapidly that we cannot effectively predict or mitigate their transmission within groups.

The highly modular social structure of gorilla populations, with many interactions within groups and few interactions between groups, may provide protection against wider disease transmission, as we found no evidence in support of transmission through the wider population. Transfers and inter-group encounters occurred during 20% and 27% of respiratory outbreaks respectively, suggesting that they could represent important pathways of inter-group transmission. However, neither mechanism appear to have caused any of the 15 outbreaks examined here, as signs of infection were not observed in previously uninfected groups following individuals transferring from, or encounters with an infected group. Where physical contact between members of different groups or transfers occurred during an outbreak, none of the individuals involved were displaying signs of infection at the time, suggesting that this pathway for inter-group transmission is rare, occurring in fewer than 1 in 15 outbreaks (< 6.6%). The ranges of infected groups during concurrent infections also did not overlap more with each other than with non-infected groups. This inability to predict the presence of signs of infection in a group from spatial overlap further supports the lack of inter-group transmission suggested by analysis of direct interactions (encounters and transfers). The minimal opportunities for transmission at this broader scale may prevent most respiratory infections from circulating more widely in the population, limiting their potential conservation impact. However, respiratory outbreaks remain common nonetheless and represent the second most likely cause of death 35, 36 . Since our findings suggest inter-group transmission is an unlikely source of these outbreaks, this implies that an alternative pathway is responsible for spreading infections to new groups.

Given the confirmed human origin of a number of respiratory infections in gorillas 37, 38 , humans represent the most plausible source of these outbreaks. The frequency of outbreaks and the lack of evidence for inter-group transmission suggests human-gorilla disease transmission may be a regular occurrence within this population. Repeated introduction of respiratory infections from humans, rather than through inter-group transmission has also been found in chimpanzees, where concurrent outbreaks in nearby communities were caused by two distinct respiratory viruses of human origin 17 . These outbreaks can result in considerable mortality, with a human respiratory pathogen outbreak resulting in 8.9% mortality in one chimpanzee population, a far higher rate than that caused by the same pathogen in humans 18 . This highlights the considerable risk from exposure to human diseases that habituated ape populations face, as well as the risk they pose as a potential source of zoonotic disease transmission to humans.

Tourism, research, daily protection and monitoring have been crucial to the intensive conservation strategy responsible for the recovery of the mountain gorilla population 31 www.nature.com/scientificreports/ great ape populations. However, these methods also pose some risks of disease transmission between humans and nonhuman great apes 33, 55 . Continuing the ongoing efforts to protect wild great apes through extensive testing, vaccination and mask-wearing for all those that come into close proximity with them will be vital, and priority vaccination of human populations close to great ape habitats could further reduce these risks 56 . Given the risk now posed to great apes by the coronavirus pandemic [41] [42] [43] , it is more vital than ever to minimize pathways of human-ape disease transmission, which pose a risk to wild great apes and humans alike. In addition to the respiratory outbreaks examined here, there are a number of other diseases that could pose a threat to gorilla conservation in the coming years including COVID-19 and ebolavirus disease (EVD). Our analyses provide information on disease transmission within wild gorillas which could help in preparation for and in response to outbreaks of such diseases. However, it is important to consider that the mode of transmission and the species between which transmission can occur may differ greatly between pathogens. Given these limitations, our findings can provide only rough guidance for conservation approaches for limiting disease transmission in gorillas more broadly. The predictable variation in contact and proximity based on age, sex and dominance could enable strategies such as targeted vaccination to be effective for pathogens with low transmissibility. But the density of connections within gorilla groups suggests that with most infectious pathogens, all group members are likely to be exposed in a relatively short period of time. Opportunities for disease transmission between groups through transfers, encounters or range overlap appear to be considerably rarer-a pattern also seen in chimpanzees 17 . However, contact with carcasses during outbreaks of high mortality 57 or contact with other infected species might still enable wider transmission of certain diseases through a population. Measures to prevent or control disease in great apes may therefore be best targeted at preventing disease reaching a group where possible e.g. by preventing initial transmission from humans in the case of anthroponotic disease. This may be considerably more effective than attempting to limit transmission within a group once it is there, where transmission is rapid and unpredictable, or preventing transmission between groups, which is likely to be rare if close inter-group contact is required for transmission.

Ethics statement. This study was non-invasive and strictly observational. All data collection followed the ethical guidelines for gorilla research put in place by the Dian Fossey Gorilla Fund and the Rwanda Development Board.

Rwanda were monitored daily by trained field assistants and researchers at the Dian Fossey Gorilla Fund. These groups were also regularly visited by tourists and park rangers. Groups were formed of one dominant adult male, multiple adult females and their offspring, and often also contained additional subordinate adult males. Researchers visited each group for up to 4 h per day and individually identified group members by physical characteristics. Long-term monitoring has enabled demographic data, such as dates of birth, sex and maternal kinship to be recorded. All gorillas could therefore be classified into the following age/sex categories: infants 0-3.5 years, juveniles 3.5-6 years, subadult males or females from 6 to 8 years of age, adult females from 8 years, blackback males from 8 to 12 years and adult males from 12 years. Adult males were further categorized into either subordinate adult males or dominant adult males based on their dominance hierarchy. Male dominance was calculated using the Elo-rating method 58,59 based on observed displacement and avoidance events using the R package 'EloRating' , version 0.43 60 as described by Wright et al. 61 .

Within-group social data. Between 2004 and 2015 data on the social relationships within habituated gorilla groups were recorded through focal sampling. During daily monitoring, researchers chose a focal individual for a 50-min follow by systematically working their way through a randomly ordered list of all group members above the age of one year. If an individual could not be observed (e.g. obscured by dense vegetation), the researcher moved on to the next individual on the list and returned to them the subsequent day. During these focal follows, a scan was completed every 10 min recording all gorillas in physical contact with the focal individual and within 5 m of the focal individual.

Predicting network position. We constructed weighted, undirected social networks based on proximity within 5 m and physical contact for each group in each year using the focal sampling data. The simple ratio index was used to estimate network edges (pairwise relationships between gorillas), such that edges represented the proportion of scans of either individual in which the pair were observed associating (either in physical contact or within 5 m). Networks included only individuals that were present in the group for at least 11 months of the year, were > 1 year old at the start of the year, and for whom more than 24 focal scans had been recorded that year. We extracted eigenvector centrality (hereafter, centrality) for individuals in both the proximity and contact networks using the 'igraph' package in R 62 .

We used generalized additive mixed models (GAMMs) in the 'gamm4' R package 63 to model how centrality in both network types changed with age for both sexes. The response variable was the centrality of each individual per group per year (combining measures across all annual networks). A smoothing term for age by each sex was included as a predictor as well as a smoothing term for the number of focal scans per individual per year to account for any potential differences in centrality driven by sampling effort. Smoothing parameters were fit using REML (restricted maximum likelihood) to avoid overfitting and checked using the gam.check function in the 'mgcv' R package 64 . The network (for each group in each year) and the individual identity were included as random effects to account for the non-independence of centrality metrics within the same network and for the same individual in different networks. Age was taken as the individuals' age at the midpoint of the year. www.nature.com/scientificreports/ We ran linear mixed models in the 'lme4' R package 65 to assess whether centrality (and the disease risk with which it is linked 24 ) could be predicted by age/sex categories (described above). The response variable was the centrality of each individual in each group in each year. Predictors included the age/sex category and both a linear and quadratic term for the number of focal scans to account for any potential differences in centrality driven by sampling effort. The network and the individual were again included as random effects. Age/sex categories were determined based on the individuals' age at the midpoint of the year.

As network centrality measures are non-independent, we calculated p-values using node-based permutations. The age/sex category of individuals was permuted within each network to create 10,000 sets of permutations in which individuals' age/sex categories had been randomized. The same model was then run on these permuted data sets. P-values represented the proportion of these null models with an effect size greater than or equal to that observed in the model run on the real data. Models were run firstly with infants as the reference age/sex category, identifying whether infants had significantly higher centrality than other age/sex categories; and then with dominant males as the reference age/sex category to verify whether dominant adult males had higher centrality than subordinate adult males.

Between January 2004 and August 2020, the presence of any signs of respiratory disease in the monitored gorillas, including coughing, sneezing and nose discharge were recorded each day. From these data, we extracted all cases of respiratory outbreaks within a gorilla group. These were defined as time periods during which at least 25% of group members showed signs of infection and no more than 2 days passed without at least one group member showing signs. For groups of fewer than 12, at least 3 group members had to show signs of respiratory infection to be classified as an outbreak, to avoid chronic cases in a small number of group members being incorrectly classified as a respiratory outbreak. The start of the outbreak was classified as the first day a gorilla showed signs of respiratory infection. The end of the outbreak was classified as the last day on which signs were observed, with no signs in any individuals for at least 3 days subsequently. In total, 15 outbreaks were identified across the 16 groups monitored between 2004 and 2020 (Fig. 5 ).

Fine-scale data on social contact and proximity within groups were available for the 7 outbreaks that took place between 2004 and 2015. We constructed susceptible-infectious-removed (SIR) models following Scully et al. 18 for each group in each outbreak using the following set of differential equations:

where S is the number of susceptible gorillas, I is the number of infectious gorillas, R is the number of recovered gorillas, β is the daily transmission probability, γ is the recovery rate, and N is the number of gorillas in the group. β and γ were estimated by fitting SIR models to the observed incidence data for the presence of signs of respiratory infection during each outbreak. The basic reproductive number, R0, was calculated as R0 = β γ . For graphs of SIR model predictions plotted over data from observations of infection signs see Supplementary Figure S2 .

Social networks based on 5 m proximity and physical contact were constructed for each outbreak using focal scans collected in the 3 months either side of the date on which the outbreak started (6 months total). The simple ratio index was again used to estimate network edges. Networks included only individuals that were present during the outbreak, were > 1 year old during the outbreak, and for whom more than 24 focal scans were available during the 6-month period. We ran path-based k-tests 50 using the test.path function from the 'k-test' github repository 66 with 10,000 iterations. These analyses assessed whether the mean weighted path length (MWPL) between individuals showing signs of respiratory infection was smaller than expected by chance, and therefore the extent to which the social network predicted the occurrence of these signs. These were run on both the proximity and contact-based networks for each group outbreak. These models were first used to determine whether individuals that had shown any signs during the outbreak were more closely connected than expected by chance. This was not done for the UGE Jul 14 outbreak since all group members had shown signs by the end of the outbreak. As all group members may ultimately be exposed to an infection due to the dense connections within group networks, we then examined the onset of signs of respiratory infection in more detail. This was done by running path-based k-tests for each day of the outbreak on which at least one group member showed signs for the first time, assessing whether the group members that had shown signs so far in the outbreak were more closely connected than expected by chance. This enabled us to detect whether there was any evidence for transmission through the social network at any time point during the outbreak.

We used linear mixed models with a binomial distribution to predict whether an individual would show signs of infection during an outbreak. In the first set of models, we predicted whether signs would be observed based only on an individual's centrality including the specific outbreak as a random effect. These assessed whether more central individuals were more likely to become infected and used only outbreaks for which social data was available. In the second model, we predicted whether signs of infection would be observed in a gorilla based on their age/sex category, this time expanding the analyses to include all 15 outbreaks (2004-2020), as fine-scale Annual centrality model 1:

gamm4 Centrality ∼ s Age, by = Sex + Sex + s(Focal scans), random = ∼ (1|Network) + (1|Individual) .

Annual centrality model 2:

lmer Centrality ∼ age/sex category + poly(Focal scans, 2) + (1|Network) + (1|Individual) . www.nature.com/scientificreports/ social data was not required. Again, the specific outbreak was included as a random effect. As age/sex categorybased differences in the likelihood of observing signs in an infected individual could potentially be obscuring an effect of centrality on the likelihood of becoming infected, we then ran a third model for the 7 outbreaks with social data available, predicting whether signs of respiratory infection would be observed based on both the individual's network centrality and their age/sex category, again with the outbreak included as a random effect.

Between-group transmission. Across the entire study period (2004-2020) up to 4 GPS locations (location of nest site, location when researchers found the group, location at noon and location when researchers left the group), any encounters with other groups or solitary males and any transfer of individuals into or out of the group were recorded for each group each day. All transfers that took place from a group during an outbreak were extracted from this dataset. Health records of the group that they transferred into were checked to determine whether the transfer event could have enabled disease transmission between the two groups. Health records of the individual transferring were also checked to determine whether they were showing signs of respiratory infection leading up to and after the transfer. All encounters with another group or solitary at proximities of less than 50 m during an outbreak were also extracted. The health records of the groups involved were then checked to verify whether the presence of any signs of infection in these groups indicated that respiratory disease transmission could have occurred during the encounter. Cases where outbreaks were detected in more than one group within a month of each other were defined as concurrent cases 35 (Fig. 5) . To investigate whether there was a spatial correlation between the ranges of concurrently infected groups, GPS locations of all groups during concurrent outbreaks were extracted. These were extracted from the start of the first group outbreak (date on which a member of that group first showed signs of infection) until the start of the last group outbreak within the overall set of concurrent outbreaks. This covered the total time period during which groups could potentially have infected each other. When this period was less than one month in total, the period was extended to one month from the start of the first outbreak to enable adequate location data points for range estimation. Home range utilization distributions for all groups during these concurrent outbreak periods were estimated using the biased random bridge approach in the 'adehabi-tatHR' R package 67 . This approach explicitly modelled successive relocations following a biased random walk. We restricted this to relocations of within 72 h and of more than 50 m. The smoothing parameter was set to 100 following Caillaud et al. 44 , with the diffusion coefficient estimated using maximum likelihood. The overlap of the home range utilization distributions of each pair of groups was calculated using the utilization distribution overlap index in the 'adehabitatHR' R package 67 . A binomial generalized linear model was then used to predict whether the ranges of infected groups overlapped more with other infected groups than they did with uninfected groups. This enabled us to determine whether concurrent outbreaks were occurring in groups that had greater home range overlap during the period of the outbreaks, suggestive of inter-group transmission through either direct encounters or spatial overlap.

Pervasive human-driven decline of life on Earth points to the need for transformative change

Anthropogenic environmental change and the emergence of infectious diseases in wildlife

Wildlife disease prevalence in human-modified landscapes

Emerging zoonotic diseases: Should we rethink the animal-human interface

One Health: A New Professional Imperative. One Health Initiative Task Force: Final Report

Parallel pandemics illustrate the need for One Health solutions

Pandemic human viruses cause decline of endangered great apes

Ape origins of human malaria

Origin of the human malaria parasite Plasmodium falciparum in gorillas

Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science (80-)

Wild great apes as sentinels and sources of infectious disease

Consequences of non-intervention for infectious disease in African great apes

Ebola outbreak killed 5000 gorillas

Catastrophic ape decline in western equatorial Africa

Risk factors for respiratory illness in a community of wild chimpanzees

Causes of death in the Kasekela chimpanzees of Gombe National Park, Tanzania

Simultaneous outbreaks of respiratory disease in wild chimpanzees caused by distinct viruses of human origin

Lethal respiratory disease associated with human rhinovirus C in wild Chimpanzees

The role of infectious diseases in biological conservation

One health, vaccines and ebola: The opportunities for shared benefits

Ebola in great apes-current knowledge, possibilities for vaccination, and implications for conservation and human health

Social networks and the spread of Salmonella in a sleepy lizard population

Linking social and pathogen transmission networks using microbial genetics in giraffe (Giraffa camelopardalis)

Using social network measures in wildlife disease ecology, epidemiology, and management

Infectious disease transmission and contact networks in wildlife and livestock

Network models: An underutilized tool in wildlife epidemiology?

Social network analysis of wild chimpanzees provides insights for predicting infectious disease risk

Social network predicts exposure to respiratory infection in a wild chimpanzee group

Network-based vaccination improves prospects for disease control in wild chimpanzees

Unraveling the disease consequences and mechanisms of modular structure in animal social networks

Extreme conservation leads to recovery of the virunga mountain gorillas

Estimating abundance and growth rates in a wild mountain gorilla population

Lack of rule-adherence during mountain gorilla tourism encounters in Bwindi Impenetrable National Park, Uganda, places gorillas at risk from human disease. Front. Public Health

Habituating the great apes: The disease risks

Respiratory disease in mountain gorillas (gorilla beringei beringei) in rwanda

Examining causes of death for mountain gorillas (Gorilla beringei beringei and G.b. undecided) from 1968-2004: An aid to conservation programs

Human metapneumovirus infection in wild mountain gorillas

Human respiratory syncytial virus detected in Mountain Gorilla respiratory outbreaks

Lethal pneumonia in a captive juvenile chimpanzee (Pan troglodytes) due to human-transmitted human respiratory syncytial virus (HRSV) and infection with Streptococcus pneumoniae

Codetection of respiratory syncytial virus in habituated wild western lowland gorillas and humans during a respiratory disease outbreak

Risk of human-to-wildlife transmission of SARS-CoV-2

Comparative ACE2 variation and primate COVID-19 risk

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates

Violent encounters between social units hinder the growth of a high-density mountain gorilla population

Gorilla susceptibility to Ebola virus: The cost of sociality

Persistence of human rhinovirus infectivity under diverse environmental conditions

Networks and the epidemiology of infectious disease

Infectious disease and grouping patterns in mule deer

Badger social networks correlate with tuberculosis infection

Evaluating empirical contact networks as potential transmission pathways for infectious diseases

Innate immunity to respiratory infection in early life

Susceptibility and immunity to common upper respiratory viral infections-the common cold

The evolutionary and ecological consequences of animal social networks: Emerging issues

An analysis of gorilla nesting sites of the Mt. Kahuzi Region (Zaire)

Keep your distance: Using social media to evaluate the risk of disease transmission in gorilla ecotourism

Vaccinate in biodiversity hotspots to protect people and wildlife from each other

Behavioral responses around conspecific corpses in adult eastern gorillas (Gorilla beringei spp

Elo-rating as a tool in the sequential estimation of dominance strengths

Assessing dominance hierarchies: Validation and advantages of progressive evaluation with Elo-rating

EloRating: Animal dominance hierarchies by Elo rating. R Package Version 0

Male body size, dominance rank and strategic use of aggression in a group-living mammal

The igraph software package for complex network research

Generalized additive mixed models using 'mgcv' and 'lme4' . R Package Version 0.2-6

Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models

Fitting linear mixed-effects models using lme4

The package 'adehabitat' for the R software: A tool for the analysis of space and habitat use by animals

We are grateful to the Rwanda Development Board (RDB) for their long-term support of the Dian Fossey Gorilla Fund and thank the many field staff involved in the collection of the data. We also thank Holly Tibble for her feedback on a draft of this manuscript and the editors and reviewers at Scientific Reports for their valuable comments. This research was made possible by the many donors who have supported the Dian Fossey Gorilla Fund and in particular the long-term support provided through U.S. Fish & Wildlife Service's Great Ape Conservation Fund.

This research was conceived by R.E.M., T.S.S. and W.E. Y.M. and W.E. were involved in collecting and managing data, and data analyses were conducted by R.E.M. The manuscript was written by R.E.M. with review and edits from Y.M., T.S.S. and W.E.

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-98969-8.Correspondence and requests for materials should be addressed to R.E.M.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.