COVID-19 in the Brazilian Oil Capital: the successful experience of public and private engagement

Brazilian strategy to overcome the spread of COVID-19 has been particularly criticized due to the lack of a national coordinating effort and an appropriate testing program. Here, a successful approach to control the spread of COVID-19 transmission is described by the engagement of public (university and governance) and private sectors (hospitals and oil companies) in Macaé, state of Rio de Janeiro, Brazil, a city known as the National Oil Capital. and and were


Introduction
Coronaviridae is a family of RNA viruses well-known to infect a large variety of mammalian and avian hosts [1]. The genus betacoronavirus is particularly important for public health, since at least ve strains have already been reported to infect humans, leading to mild diseases, such as the common cold, pneumonia as well as other severe respiratory illness. In the past two decades, previous severe respiratory disease infections by coronavirus, SARS-CoV and MERS-CoV, have provided evidence that these viruses could lead to pandemic outbreaks. The World Health Organization (WHO) reported a new strain of a coronavirus arising in China at the end of December 2019, with a strong pandemic potential [2]. In Located in the interior of the state of Rio de Janeiro capital, the municipality of Macaé adopted a massive molecular (PCR) testing strategy followed by rapid isolation of positive individuals for COVID-19. This strategy was achieved by an alliance between municipality government and donations from private sectors which allowed a faster and cost-free molecular testing for the local population. Data analyses of 4,639 RT-qPCR tests, provide important insights of disease symptoms development, age and gender trends of positivity and deaths, as well as the neighborhood spatial distribution of COVID-19 in Macaé. Overall, Macaé displayed the lowest lethality rate in the state of Rio de Janeiro during the analyzed period (April-September 2020), indicating that the massive testing strategy was successful. These data suggest that a similar testing strategy between Brazilian public universities, research centers, government and private companies could have saved thousands of lives from COVID-19 deaths at other municipalities around the country.

Methods
Nasopharyngeal swab collection, RNA extraction and RT-qPCR From April 12th to September 12th, 2020 (17 th -38 th epidemiological weeks), a total of 4,639 inhabitants displaying symptoms indicative of COVID-19 were medically examined at the municipal Coronavirus Screening Center for COVID- 19 (CSC) or at one of the four hospitals involved in the study. At each place, healthcare professionals recorded the patient's vital signs, applied an admission questionnaire, which covered essential epidemiological data, including the presence-absence of 16 clinical symptoms, as well as relevant personal information (residential address and work activity). Written informed consent was obtained from all participants and/or their legal guardians before starting the study. All methods were performed in accordance with the relevant guidelines and regulations. The study protocols were reviewed and approved by the Research Ethics Committee of the Federal University of Rio de Janeiro -Macaé (Protocol number 32868720.4.0000.5699). Study inclusion criteria included subjects with a diagnosis of COVID-19 and positive PCR-based testing for SARS-CoV-2 and ongoing symptoms and/or clinical ndings consistent with acute COVID-19, who were separated by ages and gender. Study exclusion criteria included lack of willingness or ability to provide informed consent or lack of an appropriate legal guardian or representative to provide informed consent, or another medical contraindication to donate nasopharyngeal sample. Individual nasopharyngeal swabs were inserted into 15 mL conical tubes containing 2 mL Dulbecco's modi ed Eagle's culture medium (DMEM) and antibiotics (5,000 units per millimeter of penicillin and streptomycin). These conical tubes were delivered to the Institute NUPEM-UFRJ, where they were stored up to 72 hours at 4°C. 200 μL of each sample were employed for extraction of nucleic acids, which was performed using magnetic beads (Magmax Magnetic Kit -Thermo sher) following manufacturer's instructions. For the nal step, 50 μL of elution buffer was used to elute the puri ed nucleic acids. Alternatively, this kit was substituted by a fast direct a simple commercial extraction solution easy Extract TM , as recently described by similar methods [5,6]. 40uL of nasopharyngeal swab embedded in DMEM solution was added to 40uL of EasyExtract. The nal solution was mixed with an automatic pipette and was followed by a two-step heating process of 65°C for six minutes and 98°C for two minutes. Real-time reverse transcription polymerase reactions (RT-qPCR) for the identi cation of SARS-CoV-2 positive samples were performed using a TaqMan™ approach, as previously described in the Berlin [7] or the CDC [8] protocols. Reverse transcription and ampli cation reactions were performed using the Promega GoTaq® Probe 1-Step RT-qPCR System. Each 15 μL reaction contained 5 μL of RNA template and 10 μL of a solution containing reaction buffer, reverse transcriptase/Taq mixture and primer/probe mixture to a nal concentration of 0.66 μM. All oligonucleotides were synthesized by Integrated DNA Technologies (Iowa, USA) and ThermoScienti c. The thermal cycling conditions used were 50 °C for 10 min for the reverse transcription, followed by 95 °C for 2 min and then 45 cycles of 95°C for 10 s and 58 °C for 30 s. RT-qPCR tests were considered positive when two regions of SARS-CoV-2 genomes were ampli ed. Patients were considered uninfected when ampli cation was successful only for the human internal control (RNAse P). All RT-qPCR assays were performed on a StepOnePlus TM Real Time PCR System (Applied Biosystems).

Symptom data analyses
From April 12th to September 12th, 2020 (a total of 22 epidemiological weeks), detailed medical information of 3,495 individuals tested was recovered. The proportion of positive results were compared between gender among groups of individuals separated by age according to WHO, with the rst age class being from 0 to 14 years old, then from 15 to 80 years old the age class interval was every 5 years, totalizing 15 age classes. Additionally, comparisons were also made among age classes that encompass professional activities (between 15 and 74 years old) and among the days after the rst symptom onset, using the Test of Equal or Given Proportions, which assumes equal proportions in two or more groups.
The test of Equal or Given Proportions was also performed to compare the rate of positive outcomes between two temporal windows of 11 epidemiological weeks. The three possible types of cointegration (Type 1: no trend; Type 2: linear trend; Type 3: quadratic trend) between absolute mortality and positive results series of two earlier epidemiological weeks with the Engle-Granger Test were tested. The two weeks lag was chosen after checking the observation that 14 days (± 0.69 standard error) was the average time elapsed between RT-qPCR test and death of 75 fatal victims of COVID-19 in our cohort at Macaé.
Since the 16 symptoms included in the questionnaires presented a qualitative nature (1-0 or presenceabsence of a given symptom), we employed a multivariate analysis technique to explore the similarity of symptoms on CSC positive results between men and women from the 15 age classes. Multivariate analyses constitute the branch of mathematics that deals with the examination of numerous variables simultaneously and, therefore, they have been widely used in ecological studies with the aim of exploring the relationships of several characteristics measured on a number of individuals or communities, including microorganisms [9]; [10]. Herein, we dealt with a matrix of pairwise distances or dissimilarities between individuals based on their reported symptoms with the aim to ordinate them within a dissimilarity space limited to a low dimensional representation. Since multivariate analyses are generally sensitive to zero-in ated matrices, we excluded the three least reported symptoms (drowsiness, irritability and mental confusion), which were reported by less than 5% of the individuals. The Jaccard index was chosen to construct the distance-based matrix, since it is the most appropriate metric when variables are de ned as 0-1. The distance-based matrix was then submitted to a non-metric dimensional scaling (NMDS) with the function metaMDS in the 'vegan' package [11] in R 3.6.1, for visualization of the individuals on the two-dimensional space, according to their dissimilarity. Individuals presenting similar sets of symptoms are positioned close to one another, while those with dissimilar symptoms (e.g., not sharing any symptom) are positioned further apart. Compared to other ordination techniques, NMDS is recommended even when the relationship between dissimilarities and inter-object distances is nonlinear, as long as low stress values are obtained [12]. The stress value (from 0 to 1.0) measures how good the graphical representation is of the actual dissimilarities on the distance-based matrix. According to Quinn & Keough [12], stress values greater than 0.3 indicate that the con guration is no better than arbitrary and therefore should not be interpreted.
Signi cant effects of gender and/or age class and their interaction on the distance-based matrix were assessed with a two-way permutational analysis of variance (PERMANOVA) performed with the function adonis in the 'permute' package [13] in R 3.6.1. By expecting that the NMDS would be able to reveal structure on the most and less shared symptoms among individuals, we descriptively compare the scores between men and women among age classes along the rst and second NMDS axes (NMDS1 and NMDS2, respectively) and related them to the prevalence of the symptoms.

Comparison of viral cycle thresholds (Cts)
The ∆Ct approach was chosen to compare the approximate viral load of nasopharyngeal swabs, where Cts values of RNAse P (RP) endogenous control were subtracted from the arithmetic mean of N1 and N2 targets of SARS-CoV-2 (∆Ct = (CtN1 + CtN2)/2 -CtRP), for each positive patient [14,15]. The arithmetic mean of Ct values is usually employed when using multiple reference genes in relative analyses of gene expression studies [16], and was used here to the targeted genes under the assumption that the expression of N1 and N2 viral genes are not differentially regulated between infected individuals. When compared to the presentation of crude Ct values, the ∆Ct approach has the advantage of providing a normalized Ct score, discounting for initial variations on swab collected biological material and also for the presence of possible polymerase inhibitors in the RT-qPCR reactions.
RT-qPCR positive patients were classi ed as mild or severe based on criteria previously de ned by Liu et al. [14], which include any of these conditions: respiratory distress (respiratory frequency ≥ 30 breaths/min), oxygen saturation at rest ≤ 93%, severe disease complications leading to hospitalization and death in consequence of the infection. We also attempted to classify these patients in early and late infected cases based on the day of appearance of the rst COVID-19 symptoms, as reported in the admission questionnaires. The boundary up to and after the 5th day from manifestation of the rst symptoms was chosen to establish early and late classes (respectively <5 and >5 days), as the proportions of mild and severe cases were comparatively balanced between these two intervals.
The ∆Ct of individuals were compared between mild and severe cases and between < 5 and > 5 days after the rst symptom. Additionally, interactions among them were compared by analyses of variance (ANOVA), after checking for the residuals homoscedasticity with the function leveneTest, available in the package carData [17].
Geocodi cation, spatial distribution and spatial cluster analysis The noti cation forms of the Brazilian Ministry of Health's Noti cation and Surveillance System were extracted from the Municipal Health Department of Macaé. The address data (street address, number, neighborhood and zip code) were used for the geocoding process. To describe the spatial distribution of COVID-19 occurrences in the municipality, each case location was plotted together with the case density level curves, estimated using the Kernel smoothing method. The smoothed relative risk (SRR) was used to assess the spatial distribution of COVID-19 cases (deaths) [18]. This analysis allowed the comparison of deaths caused by COVID-19 among neighborhoods. To estimate the SRR, the geocoded cases (deaths) were grouped by neighborhood, and indirect standardization [35] was used to compute the expected number of cases for each neighborhood. The SRR then follows as the ratio of the observed number of events (cases and deaths) over the expected number: where Oi represents the observed number of COVID-19 cases (deaths) in the area (neighborhood) i, and Ei represents the expected number of COVID-19 cases (deaths) for the area i.
To assess the spatial dependence of the distribution of COVID-19 cases (deaths), the Moran's I coe cient of autocorrelation was calculated to obtain the SRR. Autocorrelation statistics for aggregated data provide an estimate of the degree of spatial similarity observed among neighboring values of an attribute over a study area [19].
The spatial analyses of the data and generation of the maps were performed using R 3.6.1. The Google Maps API was adopted from the R software package ggmap [20] to geocoding, the tmap package [21] was used to plot the maps, the smooth relative risks were estimated using the Dcluster package [22] and the package spdep [23] were used for the estimates and tests the Moran's I coe cient.

Results
Macaé fatality rate was lower compared with other populous municipalities Macaé is located in the State of Rio de Janeiro, southeast of Brazil ( Figure 1). Comparison of COVID-19 death rates showed that in six months of epidemic (up to the 38 th epidemiological week), Macaé had the lowest mortality rate among the 23 most populous municipalities (>125,000 inhabitants) in the state of Rio de Janeiro ( Figure 1).
Men actively working contracted COVID-19 at higher rates than women Essential activities of offshore companies and transport are still generally men-based worldwide, and these sectors did not stop during the quarantine in Macaé. On the other hand, public services, educational system and street commerce, which usually employs genders more equally were strongly restrained by municipal decrees. Only after ve months, some of those sectors (i.e., non-essential commerce) were restarting traditional activities, whereas by October 18 th 2020, presential activities in schools, universities and public departments had not yet restarted.  Fig. 2).
Taken together, these results suggest that men actively working were more exposed to infection and contracted COVID-19 at higher rates than women.
Two distinct stages of COVID-19 spreading were identi ed in Macaé The increase in the number of tests over time at least partially re ects the carrying capacity of sampling processing at the NUPEM-UFRJ laboratory, which increased from an average of 38 to 430 tests per week, until reaching a stable rate of 400 tests, from the 31 th epidemiological week onwards (Fig. 3A). The ratio of positive RT-qPCR decreases throughout the epidemiological weeks, even though high uctuations were recorded along the rst 11 weeks (Fig. 3A). There were signi cant differences in the positive rates between two temporal windows, the rst characterized by a low number of tests (Stage 1) and the other by a high number of tests (Stage 2) (Chi2=126.62; d.f.=1; p<0.001) (Fig. 3BC). This fact may partially re ect the selection criteria of patients for RT-qPCR that arrived more severely ill at CSC during the rst weeks of the pandemic and the effectiveness of the governance policies against the spread of the virus.
A congruence between the number of positive tested and absolute deaths started to be observed only after the 32 th epidemiological week (Fig. 5D) and this trend may explain the lack of cointegration between the number of positive results and the number of deaths for Types 1(EG=-2.07; p=0.10), 2 (EG=0.01; p=0.10) and 3 (EH=0.95; p=0.10). On Stage 2, a constant ratio of 1:10 between fatality and positive outcomes was observed (Fig. 5E), even though both variables diminished with time.
RT-qPCR tests were performed during ideal window for virus detection Irrespective of gender, most positive RT-qPCR results were from sample collections between the fourth and sixth day after the rst symptoms (55% of the women and 53% of the men) and there was no signi cant difference between gender on the overall representativeness among the days after the rst symptoms (Chi2=18.077; d.f.=19; p=0.517) (Fig. 4). This nding reinforces that the ideal window for doing the RT-qPCR is between 3 and 8 days after symptom onset, regardless of the gender tested.

Respiratory symptoms are associated with hospitalization
Among the 3,751 tested for COVID-19 at the health center CSC and the hospital settings, 878 tested positive (51% from CSC and 49% from hospital). Symptoms related to respiratory disorders (i.e., cough and shortness of breath) were the most frequently reported by those requiring mostly intensive care in hospital settings, whereas varied and less frequent symptoms were more frequently reported by those with milder symptoms (Fig. 5). Among the latter, typical symptoms of u, such as headache, myalgia, runny nose and sore throat were similarly reported by women and men, whereas loss of taste, anosmia and nausea and vomit were more typically reported by women ( Fig. 5; Supl.. Fig. 3). But irrespective of gender, individuals younger than 40 years old reported more frequently headache, loss of taste, anosmia and sore throat, while individuals older than 40 reported more frequently myalgia and fever (Supl. Fig. 3). Independently of gender deaths were concentrated in individuals older than 50 years (Sup. Fig 4).

Identity and frequency of symptoms vary depending on age and gender
An analysis using the multidimensional space of the NMDS (see methods for details) showed a widely spread ordination of the positive tested men and women, reinforcing a general pattern of heterogeneity of their symptoms (Fig. 6A-B). The stress value (from 0 to 1.0) measures how good the graphical representation is of the actual dissimilarities on the distance-based matrix. According to Quinn & Keough [12], stress values greater than 0.3 indicate that con guration is no better than arbitrary and therefore should not be interpreted. Since a stress value of 0.23 was observed in our dataset, the ordination was considered adequate. The centroids of the ordinated individuals signi cantly differ between gender and age classes (Fig. 6C). Since PERMANOVA analysis detected signi cant differences within both effects, age (<40 and >40) and gender classes, it is possible to conclude that in general the identity and the frequency of the reported symptoms differ between men and women, but it is also dependent on age ( Table 2; Supl. Fig. 3).
Detailed analysis showed that those signi cant differences among individuals are associated with the milder symptoms, such as loss of taste, anosmia and nausea-vomit, that prevailed among women (Supl. Fig. 3). In addition, there was also an effect of age on the identity and frequency of the reported symptoms. Interestingly, irrespective of gender, individuals younger than 40 years old reported more frequently headache, loss of taste, anosmia and sore throat, while individuals older than 40 reported more frequently myalgia and fever (Table 2; Supl. Fig. 3). These trends are reinforced by the comparison of the scores of positive tested individuals along the rst two NMDS axes (Fig. 6C). Age has a stronger effect along NMDS1 (younger individuals with positive and older with negative scores, respectively) and gender along NMDS2, where men presented positive scores whereas women presented negative scores.

Lack of correlation between ∆Ct values from nasopharyngeal RT-qPCR swab tests and disease severity
To understand if the differences of patients with mild and severe COVID-19 symptoms previously described (see Figure 5) can also be correlated with changes in the viral load from nasopharyngeal swabs, a RT-qPCR comparative analysis was performed using the ∆Ct values as a response variable (see [24] and references therein). Lower ∆Ct values directly correspond to higher viral load in nasopharyngeal swabs. The results showed that the viral load clearly reduced with time after the appearance of the rst symptoms, but this is perceptible only for mild COVID-19 patients (Figure 7). In contrast, patients presenting severe symptoms showed a lower viral load (higher ∆Ct values), independent of the day of the RT-qPCR test was performed (<5 or >5) ( Table 3).

Highly populated neighborhoods contain COVID-19 RTq-PCR positive hotspots and deaths
To understand if the RT-qPCR positive results and deaths were evenly distributed over Macaé city or were concentrated in speci c neighborhoods, the smoothed relative risk (SRR) values (see methodology for details) were calculated. The analysis allowed to compare positivity and deaths rates among neighborhoods. SRR quintiles of positive RT-qPCR cases are displayed in Figure 8A with dark and light regions representing high and low SRR values, respectively. Similarly, Figure 8B shows the SRR quintiles of con rmed urban deaths. These maps indicate the distribution of positive RT-qPCR and death SRRs values along the city, and measures of spatial autocorrelation could be analyzed. In general, if the SRR values in the dataset are clustered spatially (high values cluster near other high values; low values cluster near other low values), the autocorrelation method of Moran's Index will be positive. When high values repel other high values, and tend to be near low values, the Index will be negative. The Moran Index statistic results were -0.035 (p-value = 0.5264) for SRR positive RT-PCR and -0.019 (p-value = 0.3555) for SRR deaths, con rming that COVID-19 infections and deaths were evenly distributed along the city and not spatially concentrated. Maps of weekly geocoded RT-qPCR positivity and deaths were developed (Sup. Fig 5 and 6) and aggregated into single maps ( Figure 8C and 8D). Our weekly comparative analysis shows that at the 17 th epidemiological week positive RT-PCR cases were already spread along the city.
Interestingly, there is a large correspondence between Figures 8 C and 8 D (positive RT-PCRs and death) showing that regions with high RT-PCR positivity region were also the ones corresponding to high death occurrences. Comparison of these data with population density (Sup. Fig 7) also highlights that positive RT-PCRs and deaths were concentrated in the neighborhoods in which the greatest number of inhabitants are found.

Discussion
Rio de Janeiro is one of the most populous Brazilian states and even though it has a diversi ed economy, the state relies mainly on the extraction of natural resources, such as oil and gas [25]. Located 200 km north of the state capital, Macaé has attracted dozens of oil and hundreds of outsourced companies in the last 50 years, which renders the municipality the title of the National Oil Capital [26]. Since the beginning of March 2020, when Macaé recorded the rst COVID-19 cases, the municipal governance began publishing a series of decrees (available at [27]) that included the closure of schools and temples, restriction of all non-essential activities, rigorous control of the arrivals by installation of sanitary barriers on all major municipality entries and the creation of the municipal Coronavirus Screening Center for COVID-19 patients (CSC), in the city center. Additionally, an alliance between the municipal governance, the Institute of Biodiversity and Sustainability of the Federal University of Rio de Janeiro in Macaé (NUPEM-UFRJ) and several oil industry related companies and hospitals allowed the implementation of a new laboratory that was urgently adapted for such a sanitary emergency, with the aim to set up the gold standard technique for SARS-Cov2 identi cation, the RT-qPCR. Consequently, from the 17 th epidemiological week onwards (April 12 th , 2020), samples from nasopharyngeal swabs from CSC and four associated hospitals were collected and delivered to a recently created laboratory at NUPEM-UFRJ for RT-qPCR analyses.
In the present study it was observed a higher percentage of positive COVID-19 men at working ages (between 25-54) than women in Macaé ( Figure 2). Research has shown that women are highly underrepresented in most extractive industries. The U.S. Department of Labor de nes a male-dominated sector as one where women constitute less than a quarter of the total workforce [28]. Macaé is one of those oilproducing cities, thus, presenting a highly men biased workforce. Therefore, this result can be directly correlated with the presence of male-dominated workforce in the city.
Another interesting observation from our data regards the distribution of symptoms among the positive cases ( Figure 5). COVID-19 symptoms have been described over six months by now [29] and our data support previous observations that symptoms related to respiratory disorders (i.e cough and shortness of breath) were most frequently reported by those requiring most intensive care in the hospital settings ( Figure 5). As previously observed [30], absolute deaths increased within age classes independently of gender in Macaé, Rio de Janeiro, Brazil, between April 12 th and September 12 th 2020. This increase in deaths with ages has been reported by several studies [30,31] highlighting that age is one of the most important death risk factors. Altogether, our data provides a framework for clinical doctors to access the most frequent symptoms which might lead to hospitalization.
The correlation between ∆Ct values and COVID-19 severity can also be discussed from our data. Brie y, there was no correlation between viral load, evaluated by ∆Ct values, and disease severity (Figure 7).
Although in principle counterintuitive, our data suggests that patients with mild symptoms may display higher viral load at the time of the exam than the severely affected patients. Interestingly, [32] have reported a similar nding when comparing viral loads between hospitalized and non-hospitalized New York patients. These authors argued that higher viral loads observed in patients presenting mild symptoms may re ect the elapsed time from infection onset, as viral loads usually peak during the presymptomatic stage or shortly after the manifestation of rst symptoms of COVID-19 and then show a slow decline during the following two-weeks [33]. Then, while many patients recover from the disease in this two-week period, a small proportion of them will suffer from progressive health deterioration, even with lower viral load using nasopharyngeal swabs. Thus, most patients identi ed as severe cases may have been infected and replicating SARS-CoV-2 several days before presenting the rst symptoms.
Together with the report from [32] our observations corroborate previous studies showing that long lasting COVID-19 syndrome is not necessarily correlated with higher viral loads, but rather with lower ones [33]. The larger sample number presented here and the more robust statistical tests than most of the previous studies reinforces this conclusion. This is of special concern for health authorities, as asymptomatic or mildly symptomatic individuals may offer more risk of spreading SARS-CoV-2 than severely compromised patients.
Another noteworthy nding from the current study comprised the COVID-19 positivity distribution along Macaé city (Figure 8, Sup. Fig 5-7). Highly populated areas from the city containing banks, supermarkets and pharmacies concentrate both RT-qPCR positive and death cases, suggesting that these essential trade activities might foster disease spreading along the city. Importantly, population density and SSR RT-qPCR positivity were directly correlated with spatial death occurrence along the city over the whole analyzed period (17 th and 38 th epidemiologic week), thus, authorities should focus on containing the spread of the disease mainly at these speci c locations. Previous analysis at two states of Brazil, Ceará and São Paulo, also provided evidence that metropolitan highly populated areas showed the greatest number of COVID-19 cases [34,35]. Spread of the virus can be highly associated with the mobility of people by means of transportation, such as planes and buses [35]. Importantly, one of the hotspots of Macaé COVID-19 RT-qPCR data overlaps with the largest bus station in the city, providing further evidence that people´s mobility and population concentration are essential factors for disease spreading. To sum up, a well-designed strategy of testing and isolation of positive individuals was successful to mitigate the effects of COVID-19 in an important oil production city of Brazil until September, 2020.
Symptom analysis showed that respiratory effects are associated with a worse disease prognosis. Geocoded spatial analysis indicated that highly populated areas display the largest number of RT-qPCR positive cases and deaths. Altogether, based on death comparison with cities lacking extensive molecular tests suggest that more than 50,000 lives could have been spared if similar strategies of massive testing would have been applied country wise. Figure 1 Macaé, the National Oil Capital, occupying the lowest fatality rate due to COVID-19 among the 23 most populous municipalities (>125,000 inhabitants) of the state of Rio de Janeiro, Brazil. Source: Rio de Janeiro Health State Department up to the 38th epidemiological week. Note: The designations employed and the presentation of the material on this map do not imply the expression of any opinion whatsoever on the part of Research Square concerning the legal status of any country, territory, city or area or of its authorities, or concerning the delimitation of its frontiers or boundaries. This map has been provided by the authors.  Table 1 to identify age classes with signi cant differences between gender (in bold).   Supplementary Files