Diversity and Source of Airborne Microbial Communities at Differential Polluted Sites of Rome

Biogenic fraction of airborne PM10 which includes bacteria, viruses, fungi and pollens, has been proposed as one of the potential causes of the PM10 toxicity. The present study aimed to provide a comprehensive understanding of the microbial community variations associated to PM10, and their main local sources in the surrounding environment in three urban sites of Rome, characterized by differential pollution rate: green area, residential area and polluted area close to the traffic roads. We combined high-throughput amplicon sequencing of the bacterial 16S rRNA gene and the fungal internal transcribed spacer (ITS) region, with detailed chemical analysis of particulate matter sampled from air, paved road surfaces and leaf surfaces of Quercus ilex. Our results demonstrated that bacterial and fungal airborne communities were characterized by the highest alpha-diversity and grouped separately from epiphytic and road dust communities. The reconstruction of source-sink relationships revealed that the resuspension/deposition of road dust from traffic might contribute to the maximum magnitude of microbial exchanges. The relative abundance of extremotolerant microbes was found to be enhanced in epiphytic communities and was associated to a progressively increase of pollution levels as well as opportunistic human pathogenicity in fungal communities.


Introduction
Urban areas, being the heart of economic development, are characterized by enormous density of inhabitants. The majority of human's world population is already living in urban areas with the projection to reach 68% by 2050 [1]. As consequence, in the last decades, there has been an increase in anthropogenic pollution [2]. Non-homogenous distribution of air pollutants in terms of quality and quantity has been observed between different functional urban zones of European cities [3,4]. The highest concentrations of pollutants associated with typical industrial contaminants have been generally found in areas with heavy traffic loads or industrial ones. Residential areas were characterized by an increase of contaminants originating from human activities such as the household heating, whereas increase of typical soil elements were observed in urban green areas [3,4].
Particulate matter (PM) is recognized among the most harmful atmospheric pollutants for the human health in cities and mega cities with particles smaller than 10 µm being considered as the most dangerous. More than 92% of the world's human populations has been exposed to air pollution levels exceeding WHO guideline limits [5]. PM10 consists of a complex mixture of organic and inorganic components, with the organic matter generally being the major component. Microorganisms such as bacteria and fungi, constitute most of the organic matter of PM10 along with other biological materials including pollen, plant debris and viruses [6,7]. While considerable efforts have been taken for chemical characteristic of PM10 and the identification of its sources [3,[8][9][10], information on temporal and spatial variability/dispersal of pollution-associated microbiomes is still scarce. Bacterial and fungal communities associated with atmospheric aerosol are gaining increasing attention due to their ability to affect both, human health and atmospheric conditions [11][12][13]. In particular, airborne pathogenic microorganisms can cause important health issues in plants, wildlife and humans. They are responsible for skin and internal organs infections, cardiovascular diseases, asthma and other hyper immune responses [13,14].
Traditionally, most of the studies of airborne microorganisms have utilized culturing to highlight their abundance and identity, but capturing only a small fraction of the total microbial community. The recent development of next-generation sequencing techniques coupled to DNA metabarcoding analysis provided culture-independent robust tools for a fine-scale molecular characterization of microbial communities [15,16]. In this context, it has been observed that taxonomic composition of bacterial and fungal airborne communities varied with the air quality [17,18] as well as between rural and urban sites [19,20]. Although severe pollution levels differentially affected the diversity and composition of bacteria and fungi in Beijing, China [21], a growing number of studies revealed a general tendency of pathogenic microorganisms, both fungi and bacteria, to change and increase in relative abundance with the worsening of the city air quality [18,19,22,23]. These results implied a different susceptibility of pathogenic and non-pathogenic microbes to pollutants. However, the mechanisms underlying the responses of pathogenic microbiomes to air pollution as well as the role of different environmental and anthropogenic sources on shaping the diversity of airborne microbiomes in a site-specific manner, and contributing to their pathogenic components, are not fully elucidated [13]. The sources of airborne microbial communities can include soils, water bodies, and leaf surfaces but also animal feces and human activities such as crop cultivation, livestock operations, harvesting and biomass burning activities [20,[24][25][26]. Microorganisms could be also delivered with the dust storms and partially originate from long-distant sources. Indeed, wind was recognized as the principal driver for the microbial dispersal into the atmosphere [27]. However, the direct comparison of the airborne microbial structures with those detected in the surrounding habitats proved that the impact of the local environments on the airborne bacterial and fungal communities was often greater than the airflow transport, in both urban and rural sites [19,20,25]. Even for the indoor environment, where the city dwellers spend the majority of their lives, a consistent part of microorganisms derived from the outdoor air [28]. Hence, the identification of the factors that mostly influence the structural diversity of airborne microorganisms in the outdoor habits of densely populated urban sites is crucial for the development of effective air quality control strategies.
In this study, we present the fine-scale characterization of both bacterial and fungal communities for PM10 collected in three urban sites of Rome, different in pollution rate. We used a DNA metabarcoding approach to investigate the diversity and community compositions of PM10-associated microbes, which can provide insight into the following questions: (1) what are the microbial community structures associated with PM10?
(2) what are the predominant local sources of microbial communities in the urban sites? (3) to what scale vegetation and paved surface soil shape the diversity of urban microbiomes? (4) do the bacterial and fungal community compositions vary among different urban air quality levels? Our results provide baseline information on the presence of potentially human pathogenic species related to chemical pollutants in Rome.

Site Description
The sampling campaign was conducted in Rome (Italy) from 6 to 12 October 2020, covering the early autumn period, with instauration of the post vacation rates of traffic. Rome is the capital city of Italy, with 2,860,009 inhabitants in 1285 km 2 . It is the third most populous city in the European Union. Three functional zones different in terms of the anthropogenic load were identified in the city as areas suitable for sampling: South-Western Road Traffic zone (Enrico Fermi Square, Lat. 41.87184 Long. 12.46953), Central Residential zone (Arenula street, Lat. 41.864025 Long. 12.469712) and North-Eastern Urban Park zone (Villa Ada Park, Lat. 41.933001.12 Long. 12.507532) (Figure 1a). Overview of the three urban sites (Urban Park, Residential and Road Traffic) of Rome (a), PM10 mass concentration (b), three meteorological parameters (temperature, relative humidity and precipitation) (c-e) and wind conditions (speed and direction) (f), recorded during the sampling campaign of Early Autumn 2020 (yellow shaded area). For temperature and relative humidity, the solid lines represent the daily-averages (blue), the maxima (red) and the minima (dark blue). The standard deviations are reported as dashes lines.
Often referred to as a green lung of Rome due to its large scale and shape (~180 hectares), Villa Ada is the second largest urban Park in the city, rich of archeological, historical and recreational areas. Arenula is a 0.5 km long street, located in a residential zone with limited traffic, adjacent to a small public garden. Fermi Square is the heart of a densely populated district, crossed by high traffic roads. Differences in terms of anthropogenic load between the three sampling sites were partly confirmed by the trend of PM10 mass concentrations recorded at the air quality stations of Regional Agency for Environmental Protection (ARPA) Lazio, which were used for PM10 sampling procedures performed for the current work ( Figure 1b). Meteorological data including wind speed and direction, temperature, relative humidity and rainfall level were recorded and provided by ARPA Lazio (Figure 1c-f). During the period of sampling (6 days) the average minimum atmospheric temperature in center of Rome was 14.30 °C and average maximum 23.52 °C, with average minimum atmospheric humidity of 43.85% and average maximum of 80.48%. Low precipitation was registered in that period, with a mean value of 1.23 mm and intensity of 0.34 mm/h. During the whole sampling period in the urban sites, two main clusters of airflows were recorded: south (S) and northeast (NE) which accounted for 50.7% and 25.97% of records.

Sample Collection
The metabarcoding and the chemical analysis were performed on three biological matrices: (1) PM10 membrane filters (2) leaves of Quercus ilex L. (oak olm), one of the dominant tree species used for urban greening in Rome, and (3) dust from the paved road surface.
(1) Air sampling was performed at the 3 m above the ground level by using two active PM10 samplers installed within an environmental monitoring station of each urban site (property of ARPA Lazio). PM10 was sampled for 6 days, with a HSRS (High Spatial Resolution Sampler, Fai Instruments, Rome, Italy) at 2 L/min flow rate on PTFE membrane filters (47 mm diameter, 0.45 µm pore size). Due to the typically low environmental concentration of bio-aerosols and the challenging to collect sufficient genetic material, PM10 for molecular analyses was collected with an Impact Sampler (SKC, Dorset, UK) at 10 L/min flow rate on polycarbonate filters (47 mm diameter, 0.8 µm pore size) following the recommendations of Ferguson et al. [15]. Each polycarbonate filter was UV-sterilized (254 nm) for 15 min on each side within a laminar flow hood and individually stored in a UV-sterilized plastic filter holder. In order to improve the DNA recovering efficiency, polycarbonate filters were changed every 48 h with sterile gloves and transported in sterile conditions in a zipper cooling bag (4 °C) to the laboratory where immediately processed. Totally, 3 polycarbonate filters and 1 PTFE filter were available from each urban site. (2) Leaf collection was carried out at the end of the 6th day sampling period. Three mature healthy trees of Q. ilex were selected in each site, in the immediate vicinity of the environmental monitoring station. Thirty leaves were collected from the height of 3-4 m with telescopic pruning shears for each tree and stored at 4°C in sterile plastic bags until further analysis. (3) Similarly, at the end of sampling period, dust from the paved road surfaces close to the environmental monitoring station was collected from the area of 1 m 2 in three replicates in each site. Dust was placed with a sterile brush in 50 mL vials, transported in the cooling bag and conserved at 4 °C before processing.
All samples were processed not later than 3 days after the sampling procedure to prevent any microbial growth and change in community over the time. and Reverse: 805R5′-GACTACHVGGG-TATCTAATCC-3′ which target the hypervariable V3-V4 regions of the bacterial 16S rRNA gene [29] and Forward: ITS15′-TCCGTAGGTGAACCTGCGG-3′ and Reverse: ITS45′-TCCTCCGCTTATTGATATGC-3' targeting the ITS1-5.8S-ITS2 region of the fungal rRNA gene [30]. In order to integrate relevant flow-cell binding domains and unique indices, dual indexed libraries were generated using Nextera XT library preparation technology according to manufacturer's protocol (Illumina, San Diego, CA, USA). Eighteen equimolar single-stranded DNA libraries were pooled and then pair-ended sequenced (2 × 300 cycles) on Illumina MiSeq platform.

Bioinformatic Analysis
As a first step of the amplicon metagenomic bioinformatic analysis, a quality check was performed on the raw sequencing data, removing low-quality bases and adapters while preserving the longest high-quality part of the reads. In particular, demultiplexed 300 base pair (bp) paired-end reads were processed with the program Trimmomatic v0.36 [31] to remove adapter sequences, to discard low-quality reads (Phred quality score = 20) and to be trimmed to a minimum length of 50 bp. Secondly, the entire bioinformatic processing of the trimmed raw sequences was performed by GAIA 2.0 (http://gaia.sequentiabiotech.com, accessed on 28 December 2021), a metagenomics data analysis software developed and provided by Sequentia Biotech (Barcelona, Spain; https://metagenomics.sequentiabiotech.com/, accessed on 28 December 2021). As fully described in Paytuví et al. [32], trimmed paired-end reads were mapped against the custom-made databases obtained from the NCBI nr nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide/, accessed on 28 December 2021) by using BWA-MEM algorithm. Then the best alignments were extracted and classified into the most specific taxonomic level using an in-house Lowest Common Ancestor (LCA) algorithm implemented in GAIA 2.0. As described by Paytuví et al. [32], minimum identity thresholds were applied to classify reads into species, genus, family, order, class, and phylum levels and define the operational taxonomic units (OTUs) at each level.

Lifestyle, Extremotolerance and Human-Pathogenicity Classification of Bacteria and Fungi
Two-step classification scheme was introduced to describe the predominant lifestyles/habitats and extremotolerance of the microbial species identified in three sample types, Leaf Surface, Aerosol (PM10) and Road Dust, collected in three urban sites of Rome.
Two sub-datasets, including the overall bacterial and fungal OTUs detected at species level, were selected. As reported by Romano et al. [33], BacDive database (The Bacterial Diversity Metadatabase) [34], NCBI nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide/, accessed on 28 December 2021) and published sources were used to classify the identified bacterial species based on their main habitats: terrestrial and aquatic environments (TAE), soil (SOI), marine and fresh water (MFW), human, animal, and insect (HAI), plant-specific (PLANT) and other (OTHER). We further classified the same subset of bacterial species as extremotolerant when they were also recognized as microorganisms of hostile habitats (e.g., high/low temperature, dryness, UV irradiation, saline stress, soil contaminated with nitroaromatic compounds, soil from lead-zinc and uranium mine, arsenic-contaminated farmland soil and crude oilcontaminated soil). Similarly, MycoBank (https://www.mycobank.org/, accessed on 28 December 2021), the Encyclopedia of Life (EOL) (https://eol.org/, accessed on 28 December 2021) and NCBI nucleotide database were used to classified the identified fungal species as herbaceous (H), woody (W), Ubiquitous-plant-associated (U), saprophytic (sapro), parasitic (para) and surface inhabiting (inhab) fungi. As described by Gostinčar et al. [35], in order to identify the extremotolerant fungi, six approximative ecological categories were also considered: (1) ubiquitous fungi, (2) osmotolerant fungi, (3) rock-inhabiting fungi, mainly surface colonizers, which are subjected in their natural habitat to varying temperatures, dryness, and solar irradiation, (4) psychrotolerant fungi, represented mainly by species from deep ocean waters and polar regions, (5) xerotolerant fungi, including the dry fraction of the outdoor mycobiome, from desert soil and other similar habitats, and (6) other. Finally, the identified species of bacteria and fungi were classified as human opportunistic-pathogenic following the classification found in BacDive database [34] and de Hoog et al. [36], respectively, with support of NCBI nucleotide database. Pathogenic microorganisms were further divided into two groups, in accordance with their potential danger to human health: BSL1, and BSL2. Description of each extremotolerant and human-pathogenic bacterial and fungal species were available in Supplementary Materials.

Analyses of Chemical Components
PM10 membrane filters. PTFE membrane filters were used for the gravimetric assessment of PM10 mass concentration, using an automated microbalance (Sartorius R180D Analytical Balances, Sartorius AG, Göttingen, Germany). Separation and subsequent chemical analysis of both the water-soluble and insoluble fraction of PM10 were performed following the protocols of Canepari et al. [36,37]. This two-step chemical analytical method is highly efficient for increasing the selectivity of PM10 elemental components as specific source tracers, for both natural and anthropogenic emission sources [3,38]. Content of 38 elements in PM10 (Al, As, B, Ba, Be, Bi, Ca, Cd, Ce, Co, Cr, Cs, Cu, Fe, Ga, K, La, Li, Mg, Mn, Mo, Na, Nb, Ni, Pb, Rb, Sb, Se, Sn, Sr, Te, Ti, Tl, U, V, W, Zn, Zr) in the two solubility fractions was estimated by analyzing water-extracted and acid-digested solutions with inductively coupled plasma mass spectrometry (ICP-MS, Bruker 820-MS, Billerica, MA, USA), The limits of detection (LODs; Supplementary Materials)were set at 3 times the standard deviation (SD) of 5 replicate blank determinations. Further details on the preparation of external standards for the calibration and internal standards to control the nebulizer efficiency are reported in Massimi et al. [3].
Leaf surface. For the chemical characterization of leaf deposited PM10, freshly picked leaves (10 leaves from each tree, leaf area ca. 200 cm 2 ) were placed in sterile flask with 50 mL of deionized water (Milli-Q Reference Water Purification System, Merck Millipore, MA, USA) and shaken on a rotator at 300 rpm for 30 min. Solution was then filtered through 100 µm mesh to eliminate coarse impurities and placed at 65 °C in an oven (STF-F 52 FORCED VENTILATION, Falc Instruments, Bergamo, Italy) until water was fully evaporated. The obtained precipitate was weighted and then subjected to a two-step acid digestion. In the first step, 5 mL of HNO3 (67%, Promochem, Wesel, Germany) were added to the precipitate and kept for three hours at 95 °C. After cooling, 0.5 mL of H2O2 (30% Suprapur, Merck Millipore Ltd., Billerica, MA, USA) were added to the mixture and kept at 95 °C for an additional hour. Then, solutions were diluted in 50 mL of deionized water and filtrated with disposable syringes and NC filters (0.45 mm porosity, GVS Filter Technology, Indianapolis, IN, USA). These acid-digested solutions were chemically analyzed for determine the concentration of 34 elements (Al, As, B, Ba, Be, Bi, Ca, Cd, Ce, Co, Cr, Cs, Cu, Fe, La, Li, Mg, Mn, Mo, Na, Nb, Ni, Pb, Rb, Sb, Sn, Sr, Ti, Tl, U, V, W, Zn and Zr, by ICP-MS (Bruker 820-MS, Billerica, MA, USA), as previously reported for PM10 membrane filters.
Paved road surface. Dusts collected from the paved road surfaces were firstly sieved through a 100 µm mesh for homogenization. About 10 mg of each road dust sample were randomly selected and used for chemical analysis by adapting the two-step chemical fractionation procedure previously described for PM10 membrane filters. Both, water-extracted and acid-digested solutions, were analyzed by ICP-MS (Bruker 820, Bremen, Germany) for the quantification of 31 elemental components (As, B, Ba, Be, Bi, Cd, Ce, Co, Cr, Cs, Cu, Ga, La, Li, Mn, Mo, Nb, Ni, Pb, Rb, Sb, Sn, Sr, Te, Ti, Tl, U, V, W, Zn, Zr) in the two solubility fractions.

Statistical Analysis
All the exploratory statistical analyses were carried out using R v3.6.0 statistical programming language (R Development Core 2019) and the cited associated packages. In order to mitigate the impact of differences in the sequencing depth on the species richness evaluation, rarefaction and extrapolation curves were computed for each sequence library with iNEXT 2.0-12 R-package [39]. Hill numbers of order q [40], species richness (S*), and exponential of Shannon entropy ( 1 D*) were calculated on data rarefied to the smallest library size for all investigated type of samples (Leaf Surface, PM10 and Road Dust) in three urban sites of Rome (Urban Park, Residential and Road Traffic) for bacterial and fungal communities. All two Hill numbers were estimated as the mean of 100 replicate bootstrapping runs to estimate 95% confidence intervals. Significant differences in the alpha diversity estimators between sample types were tested by using one-way analysis of variance (ANOVA) with the aov function, and then post-hoc Tukey's HSD tests with R-package agricolae [41]. A linear regression analysis was performed to assess a possible correlation between the normalized PM10 mass concentration computed in the three urban sites for bacterial and fungal communities and the two normalized-alpha-diversity estimators, S*and 1 D*.
Beta-diversity measures were used to evaluate changes in bacterial and fungal community composition between sample types and urban sites. An UPGMA (unweighted pair group method with arithmetic mean) non-hierarchical clustering ordination of samples was performed using two pairwise distance matrices calculated using R-package vegan [42]: (1) based on the Jaccard index, a similarity coefficient which does not take into account the abundance of the OTUs but rather relies on presence/absence of OTUs, and (2) based on the Bray-Curtis index which consider the abundance of each OTU (i.e., structure of the microbiome) and estimates of the dissimilarity of the samples. Bootstrap support for each tree was determined by resampling loci 1000 times. To estimate the robustness of differences in microbial assemblies among sample types, we also used a non-parametric test of significant similarity differences between two or more groups (ANOSIM) implemented in vegan package (anosim function), with the number of permutations set to 1000. Subsequently, to reconstruct the source-sink relationships of aerosol (PM10)-associated microbial communities with microbes of leaf surface of Q. ilex, and surface-road dust in each urban site, we computed the normalized-differences in the proportional abundances (D) of each OTU at species-level between leaf and aerosol samples δ OTU was assigned to: (i) airborne PM10 when δ(leaf-aerosol) and δ(soil-aerosol) are both negative, (ii) road dust when δ(soil-aerosol) and δ(soil-leaf) are both positive, and (iii) to leaf surface when δ(leaf-aerosol) is positive and δ(soil-leaf) is negative. The normalizeddifferences in the proportional abundances of shared-OTUs were 3D-visualized for each sampling site by R-package plotly [43].
Finally, Principal Component Analysis was performed with R-package ggbiplot [44] to explore the urban site-specific patterns of selected chemical components in each sample type. The variables were selected depending on their ability to selectively trace PM10 emission sources [3]. Statistical differences in mean concentration of the elemental components between urban sites per sample type were investigated with the post-hoc Tukey's HSD tests after providing significance within the ANOVA analysis, or Mann-Whitney pairwise tests after providing significance using the non-parametric Kruskal Wallis test with Rpackage agricolae. In addition, Spearman's rank correlation analysis was applied to investigate the relationships between the chemical components used as pollution tracers and the relative abundance of extremotolerant microbes and human-pathogenic fungi in the three sampling types.

Richness and Diversity of Bacterial and Fungal Communities
A total of 4,613,980 bacterial 16S-amplicon sequences and 3,399,816 fungal ITS-amplicon sequences were detected from three aerosol, three road surface and three leaf surface samples. After quality filtering, our 16S and ITS datasets consisted of 4,611,848 and 3,364,968 sequences which clustered into 700 bacterial species-level OTUs and 683 fungal species-level OTUs, respectively, with 2,194,872 (47.6%) and 1,449,234 (43.10%) total reads assigned. The average number of reads (average ± SE) per sample was 243,874.7 ± 16,759.48 and 161,026 ± 12,210.04 for bacteria and fungi, respectively. In order to include all nine samples in each dataset, rarefaction was set to the lowest number of reads, 89,081 for bacteria and 59,622 for fungi (Supplementary Figure S1). All the rarefaction curves of OTUs abundance showed similar logarithmic shapes, tending to approach the saturation plateau with the increasing number of reads. This indicates that the sequencing depth was sufficient to recover the bacterial and fungal diversity in our different sample types.
The rarefied values of the alpha-diversity estimators for bacteria and fungi are summarized in the Table 1. For bacteria, mean species richness (S*) was higher for the PM10 and road dust samples (p < 0.05) than plant epiphytes, as well as the exponential of Shannon entropy ( 1 D*) d iversity index (p < 0.01) (Figure 2a). Table 1. Alpha-diversity parameters for all bacterial and fungal species-level OTUs found in the investigated type of samples (Leaf, Surface, aerosol PM10 and Road Dust) in three urban sites of Rome (Urban Park, Residential and Road Traffic): Hill numbers of order q [36], species richness (S), species richness (S*), and exponential of Shannon entropy ( 1 D*) calculated on data rarefied to the smallest library size were reported. Data sets are rarefied at the same sequencing depth for each library, 89,081 and 59,622 counts for bacteria and fungi, respectively. The 95% lower and upper confidence limits of diversity were reported in brackets.
The airborne PM10 harbored also higher fungal richness (p < 0.001) and diversity (p < 0.001) than paved road dust which in turn showed significant differences in terms of species richness and Shannon diversity (p < 0.001) in comparison to leaf samples (Figure 2b). Pairwise comparisons of least square means were performed using Tukey's HSD test, and the horizontal lines show significance levels between groups (* p < 0.05, ** p < 0.01, *** p < 0.001).

Fungal Communities
The statistical assignment of overall fungal OTUs at different taxonomic levels revealed 6 phyla, 23 classes and 78 orders (Supplementary Materials Table S2). The predominant fungal phyla were Ascomycota and Basidiomycota in all collected sample types. At phylum-level, UPGMA clustering analysis based on Bray-Curtis dissimilarity matrix revealed that the three fungal communities associated with PM10 fraction collected in Urban Park, Residential and Traffic Road were clearly differentiated from fungal microbes isolated from leaf surface and paved road in the same sites (Figure 3b). Indeed, the average relative abundance of Basidiomycota was 73.30 ± 0.84% in the airborne communities, about three time higher than in road dust (12.24 ± 5.32%) and phylloplane samples (20.70 ± 2.25%). PM10 communities held around three-time lower proportion of Ascomycota (15.98 ± 1.35%) than those in the epiphytic (61.94 ± 3.31%) and road dust (65.94 ± 9.62%). The predominant fungal classes of PM10 samples were Agaricomycetes (67.18 ± 1.12%) followed by Dothideomycetes (10.31 ± 0.62%), and Eurotiomycetes (2.24 ± 0.96%), in contrast with road dust and epiphytic samples mostly composed by Dothideomycetes (52.38-32.01%), Eurotiomycetes (15.82-4.17%), Tremellomycetes (6.17-5.54%), Exobasidiomycetes (5.92-1.04%) and Leotiomycetes (3.33-1.23%) (Supplementary Figure S3a). The fungal orders most responsible for these marked differences included the Polyporales, Agaricales, Corticiales, Hymenochaetales, Russulales, and Auriculariales, highly enriched in the airborne microbes (Supplementary Figure S3b). The composition on lifestyle-level also exhibited distinct differences between the particle fraction and leaf/road dust communities (Supplementary Figure S4). Applying a common lifestyle classification scheme to the 384 species taxa identified in the air communities, 239 (62.2%) were associated with dead wood or trees, so they predominantly either live as wood decomposers, tree parasites, or are found on the surface of trees, with a mean relative abundance of 71.93 ± 0.62% and 14.91 ± 0.85 for woodysaprophytic/parasitic and woody-surface inhabiting fungi, against the equivalent 46.23 ± 6.93% and 3.59 ± 4.38 for phylloplane and 15.99 ± 4.56% and 3.50 ± 3.04 for road dust. The remaining airborne taxa (37.8%) were associated to herbaceous material or non-woody plants, again being decomposers, parasites, or surface inhabitants (Supplementary Figure  S4).

Sink-Source Relationships between Microbial Communities
Three analytic approaches revealed the immigration extent of microbes from the potential sources (vegetation and road dust) to the sink (air). A non-hierarchical clustering ordination of sample types based on bacterial species-level compositions strongly differentiated the epiphytes from the two other biological matrix communities ( Figure 4). Aerosol and road dust samples were not readily distinguishable on the UPGMA cluster based on Jaccard index. However, taking into account the abundance of each specieslevel-OTU with Bray-Curtis computation, UPGMA cluster showed a clear divergence between the structure of the microbiome in the PM10 and road dust (Figure 4). A statistically significant separation of the three groups of sample types was further confirmed by the analysis of the similarities (ANOSIM, R = 0.95, p = 0.0035). However, Venn diagrams highlighted that the highest proportion of bacterial OTUs was shared among all the three sample types in each site, ranging from 313 (46.37%) for Urban Park to 472 (69.83%) for Traffic Road, where higher resuspension of road and soil dust is expected due to the influence of vehicular traffic. In addition, a total of 145 (21.48%), 218 (31.59%) and 77 (11.39%) bacterial OTUs overlapped between PM10 and road dust samples in Urban Park, Residential and Traffic Road site, respectively, whereas the epiphytic communities showed only 41 (6.07%), 13 (1.88%) and 25 (3.70%) of OTUs in common with airborne communities (Supplementary Figure S5a). The analysis of the two-way normalized-differences in the proportional abundances of OTUs revealed that several bacterial taxa detected in the aerosolassociated microbial communities were highly dominant on road dust (Figure 5a). Figure 5. The reconstruction of source-sink relationships of aerosol (PM10)-associated microbial communities with microbes of Q. ilex leaf surface and road dust detected in the Urban Park, Residential and Road Traffic sites of Rome. Prevalence of each bacterial (a) and fungal (b) OTU in the airborne PM10 (blue dot), road dust (red dot), and leaf surface (green dot) was displayed by plotting the normalized-differences in the proportional abundances between leaf and aerosol samples δ(leafaerosol), between road dust and aerosol samples δ(soil-aerosol), and between road dust and leaf samples δ(soil-leaf) on the x, y and z axes, respectively. All OTUs at species-level detected in the aerosol (PM10)-associated microbial communities were included in each analysis.
The proportion of airborne bacteria that were presumably of road dust origin ranged from 35.07% in Residential and 35.24% in Urban Park, with the highest peak, 37.48%, for Road Traffic site. Around half of the bacterial OTUs at species-level associated to PM10 microbial communities were still not prevalent in any of the two analyzed matrices. For fungal assemblies, statistical differences at OTUs-species-level were detected among three sample types (ANOSIM, R = 0.94, p = 0.005). Broad patterns of divergence in OTUs-species-level composition were observed between air samples and the remaining fungal microbes by using Jaccard index (Figure 4). In general, the beta diversities between road dust and leaf samples were more similar than air samples. UPGMA clustering analysis based and Bray-Curtis index divided the biological matrix types in three groups, with plant and road dust samples grouped separately. Similar to bacterial communities, mostly of the fungal OTUs at species-level were shared among all the three sample types, with the highest proportion in the Traffic Road site (37.69%) (Supplementary Figure S5b). Despite high incidence of unique OTUs in all air samples, from 18.05% for Urban Park to 24.14% for Residential area, an ample overlap of fungal OTUs between PM10 and road dust collections has been observed, with 21.35%, 23.55% and 14.11% for Urban Park, Residential and Road Traffic urban site, respectively. Venn diagrams also indicated that plant epiphytes were a weaker source of airborne fungi than road dust, with an overlapping of specieslevel OTUs from less than 3% to 6% (Supplementary Figure S5b). Taking into consideration the relative abundances of a given fungal taxa on the putative PM10 sink, we confirmed that road dust was indeed a substantial source of airborne fungi, with higher prevalence than epiphytes in all three urban areas, with the highest value at Road Traffic site (24.89%). It is noteworthy that 67.91%-68.79% of fungal taxa were dominant on-air-associated microbes, suggesting a wide contribution of various environmental unknown sources to the composition and abundance of outdoor air fungi (Figure 5b).

Effect of Air Pollutant on Bacterial and Fungal Communities
In agreement with ARPA-Lazio air quality parameters (Figure 1b), the PM10 mass concentration varied among the three sampling sites, progressively increasing from the Urban Park (20.76 ng/m 3 ) to Residential (35.47 ng/m 3 ) and Road Traffic (44.70 ng/m 3 ) area. This range was positively correlated with the rarefaction exponential of Shannon entropy index of the airborne bacterial (R 2 = 0.98, p < 0.001) and fungal (R 2 = 0.87, p < 0.001) communities of aerosol samples, with the highest alpha-diversity in the most polluted site, as well as with the fungal communities of road dust (R 2 = 0.85, p < 0.001). However, the epiphytic fungal communities showed an inverted relationship with the Shannon diversity estimator (R 2 = 0.99, p < 0.001) ( Figure S6).
The majority of the detected bacterial species were associated with soil (SOI) and the terrestrial and aquatic environment (TAE) habitats, with an overall mean relative abundance of 72.88 ± 6.37% for leaf, 71.73 ± 5.74% for PM10 and 80.16 ± 12.50% for road dust samples (Figure 6a). Out of the total 294 bacterial species identified in eight phyla, Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Cyanobacteria, Planctomycetes, Acidobacteria, and Abditibacteriota, 77 (26.1%) are classified as extremotolerant, adapted to grow in one or more extreme conditions such as low/high temperatures, high UV radiation, soil aridity and heavy metals-and hydrocarbons-contamination (Supplementary Materials Table S3). Their relative abundance was found to be enhanced in epiphytic communities (47.75%-55.55%), with Hymenobacter genus accounting for 35.42-39.93% of the reads. The extremotolerant bacteria, mostly SOI and TAE habitant, progressively increased in the epiphytic and airborne communities passing from Urban Park, to Residential and Traffic Road. The highest content of extremotolerant bacteria was also detected in the road dust collected in the polluted area close to the traffic roads (Figure 6a). Similarly, out of the total 424 fungal species identified in two phyla, Ascomycota and Basidiomycota, 83 (19.6%) were recognized as rock-inhabiting, osmotolerant, xerotolerant, psychotolerant and polytolerant fungi, subjected in their natural habitats to temperature fluctuations, variable dryness and solar irradiation (Supplementary Materials Table S4). Ubiquitous and osmotolerant fungi were prevalent in the epiphytic and airborne communities, whereas xerotolerance and rock-inhabiting lifestyles were preponderant on road dust microbes. All the biological matrices showed the same increasing trend of extremotolerant fungi in terms of relative abundance passing from Urban Park to most polluted site, with the highest mean value (67.70 ± 10.62%) on phylloplane communities (Figure 6b). In the latter category, phylloplane-associated osmotolerant Pseudomicrostroma phylloplanum and a black-yeast-like surface ubiquitous colonizer Aureobasidium pullulans, accounted for the majority of the reads (~44%). In addition, overall 32 (7.55%) species were classified as human-pathogenic fungi (Supplementary Materials Table S5). Among them, members of the genera Cladosporium (C. cladosporioides, C. herbarum, C. sphaerospermum, C. velox), Alternaria (A. alternata, A. infectoria, A. tenuissima), Aspergillus (A. chevalieri, A. niger), Didymella (D. glomerata, D. microchlamydospora) and Aureobasidium (A. pullulans, A. melanogenum) were the most abundant pathogens, together with Coprinellus domesticus for PM10 and Filobasidium magnum, Naganishia albida, and Phoma herbarum for road dust. Despite 14 bacteria (4.77%) were identified as human pathogens (Supplementary Materials Table S6), the polluted rate of urban sites seems to affect only the relative abundance of pathogenic fungi, which has been predominant in epiphytic and airborne communities of Road Traffic area (Figure 6c,d).
The Principal Component Analysis (PCA) applied to the elemental concentrations of PM10 samples revealed that the Road Traffic site was mostly characterized by a prevalence of the insoluble fractions of elements which are known as tracers of anthropogenic emission sources (Figure 7). Among these elemental components such as lead (Pb_Ins), anti-monium (Sb_Ins), zinc (Zn_Ins), nichel (Ni_Ins), chrome (Cr_Ins), copper (Cu_Ins), molybdenum (Mo_Ins), and tin (Sn_Ins), are tracers of the non-exhaust emission of vehicular traffic, being components of brake pads 'wear, disks and tires [3,45]. All these elements were abundant also in the road dust samples collected at Road Traffic area (Supplementary Materials, Tables S7-S9). In the same site, high concentrations of Cu, Sb and Mo were found on leaf surfaces whereas PM10 were also characterized by higher concentrations in the insoluble fractions of allumium (Al_Ins), rubidium (Rb_Ins), lithium (Li_Ins), cesium (Cs_Ins), tallium (Tl_Ins) and titanium (Ti_Ins). These elements are constituents of soils and often attributed to the removal and resuspension of deposited material from the ground to the atmosphere, both for natural (i.e., wind) and anthropogenic activities (i.e., vehicular traffic) [46][47][48]. In leaf surface and road dust samples, they were found at particularly high levels at the Urban Park. The soluble fraction of three main tracers of biomass burning emissions, cadmium (Cd_Sol), rubidium (Rb_Sol) and cesium (Cs_Sol) [44], were detectable in the PM10 of all urban sites, but mostly characterizing the Residential zone, as also suggested by the parallel analysis of road dusts. Two main tracers of sea salt source, magnesium (Mg_Sol) and sodium (Na_Sol), were also recorded in the overall airborne PM10 (Figure 7). Analyzing the correlation among the relative abundance of extremotolerant bacteria and fungi, as well as the human fungal pathogenics and chemical elements in the PM10 (Pb_Ins, Sb_Ins, Zn_Ins, Ni_Ins, Cr_Ins, Cu_Ins, Mo_Ins, and Sn_Ins,) a highly positively correlation (p < 0.05) was found. Similar positive correction (p < 0.05) was observed with the increasing concentration of Cu, Sb and Mo in the leaf surface (Supplementary Figure S7). The same eight chemical tracers of vehicular traffic also correlated (p < 0.05) with the fungal extremotolerance in road dust.

Discussion
In this study, we used a comprehensive multidisciplinary approach to provide airborne microbial fingerprints across a gradient of urbanization and traffic density in Rome. We reconstructed the source-sink relationships between the PM10 and the local sources, road dust and vegetation, and unveiled the effects of different urban environments (traffic, residential, and park) on bacterial and fungal community structure in terms of extremotolerance and human pathogenicity of microbes during Fall 2020.

Airborne Microbial Community Structure across a Gradient of Urbanization
Our results demonstrated that bacterial and fungal airborne communities were characterized by the highest alpha-diversity and grouped separately from epiphytic and road dust communities detected in different urban areas of Rome.
The major airborne bacterial phyla identified in this study were Proteobacteria, Actinobateria, Firmicutes, Bacteroidetes, and Cyanobacteria, similar to the taxa identified in recent bioaerosol studies performed in Southern and Northern Italy [33,49]. As dominant bacterial phyla, Proteobacteria, and Actinobateria, were ubiquitous in urban environments, related to anthropogenic activities such as urbanization and agriculture, with high atmospheric versatility [19][20][21][22]33,49]. The phyla Ascomycota and Basidiomycota dominated the PM10-associated fungal communities in Rome which is in accordance with the findings of several previous studies of fungi in the atmosphere using culture-independent methods [19][20][21][22]25]. Although Yamamoto et al. [50] demonstrated that the relative abundance of Basidiomycota increased with decreasing aerodynamic diameter in outdoor air, it is generally accepted that Ascomycota are characterized by single-celled or filamentous hyphal forms easily aerosolized, in contrast to the majority of Basidiomycota which are considered too large to be easily ejected and typical of rainforest and rural habits [11]. For these reasons, the airborne fungi of urban areas were often dominated by ascomycetes [19][20][21][22]25]. In this study, differences in the Ascomycota to Basidiomycota ratio were observed in all urban sites, with a strong prevalence of basidiomycetes. Plenty of studies reported that the composition of total fungal communities in the atmosphere varied seasonally, with peaks of basidiomycetes abundance in early autumn and winter [51][52][53]. This seasonal variation is driven by several factors including fruiting patterns of Basidiomycetes, and complex interactions between fungal dispersal abilities and meteorological factors such as moisture and wind speed [54]. Therefore, the observed airborne fungal patterns can reflect a snapshot in space and time of air composition in early fall of Rome. In particular, we can assume that during the sampling campaign, fruiting bodies developed with subsequent release of basidiospores favored by 60-70% of relative humidity.
In general, the large and temporal heterogeneity of air particulate might also explain the high species richness detected in our airborne microbial communities, almost equivalent to those computed for road dust samples which are considered one of the most biodiverse environments for microbes [25]. Instead, microbial diversity and richness observed in the surface soils are constantly higher than those on leaf surface. In particular, the Gram-negative bacteria, Proteobacteria and Bacteroidetes, and Gram-positive bacteria, Actinobateria, were the predominant phyla of the epiphytic communities associated to Q. ilex trees. Among the most abundant orders identified in all epiphytes, Sphingomonadales, Rhizobiales, Rhodospirillales, Enterobacterales, Burkholderiales, and Sphingobactereriales, are classified as plant-leaf associated bacteria and belonged to Proteobacteria and Bacteroidetes [55,56]. The remaining predominant orders of epiphytes, Cytophagales, Nakamurellales, Kineosporiales, Propionibacteriales, as well as Deinococcales belonging to Actinobateria and Deinococcus-Thermus phylum, respectively, have been recently recognized as tree phyllosphere or bark-associated residents [57,58]. Mostly of them, including Hymenobacter species, survive in unfavorable conditions, occurring in many extreme environments characterized by high temperatures, UV and gamma radiation, or by low levels of nutrients. Similarly, we also found that the epiphytic fungi of Q. ilex were dominated by five classes, ascomycetes Dothideomycetes, and Eurotiomycetes as well as basidiomycetes Tremellomycetes, Exobasidiomycetes and Cystobasidiomycetes, as detected in tropical mangrove species leaves [59] or in phyllosphere of Mediterranean plant species such as grapes [60] and olive trees [61]. The high incidence of the Dothideomycetes class was primarily determined by the abundant detection of the genera Aureobasidium (Order Dothideales), Cladosporium (Order Cladosporiales) and Alternaria (Order Pleosporales), characterized by melanized hyphae/spores. Melanin accumulation protects DNA from environmental stresses and confer a clear advantageous in areas with high levels of radiation and extreme temperatures, such as Mediterranean regions [61]. Although the taxonomic composition might be affected by the host-species identity [62,63], the current structure of phyllosphere bacteria and fungi of Q. ilex is not surprising given the hostile habitat of leaf surface, in terms of availability of nutrients and constantly exposition to environmental fluctuations. However, in the present study, the patterns of dissimilarity in bacterial composition observed between epiphytic, road dust and airborne communities allowed us to evaluate the putative local sources of microbial diversity of the atmosphere in our urban sites.

Local Sources of Airborne Microbial Taxa in Urban Environment
As recently demonstrated in US [64] and China [19], urban bacterial and fungal communities tend to be less variable than those in the rural areas. The authors suggested that the urbanization and the invariability of land use might promoted the homogenization of the airborne microbiota in urban areas as well as led to a decrease of species richness. In this study, the identified microbial taxa varied by sampling type. Our data also supports the hypothesis that habitat heterogeneity, presence of urban green infrastructures and anthropogenic disturbances can affect bacterial and fungal structure in the three selected urban environments. At the beginning of autumn 2020, the airborne microbial community of the three urban sites of Rome were in fact not panmictic, showing a positive correlation between the alpha-diversity richness and the concentration of PM10, that increased passing from Urban Park to Road Traffic sampling site. The source-sink analysis indicated that road dust was a primary local source of outdoor air bacteria and fungi in urban environment of Rome. This is in agreement with previous studies [19,20,25] in which soil-derived microbes reached the peak in relative abundance during late summer and fall. The highest levels of dust loads were detected under conditions of prolonged drought. In our case, the most wind-dispersible soil particles were very rich of microorganisms and might have served as one of the most important sources of airborne bacteria and fungi when soil moisture was still low and temperatures started declining. It is worth noticing that the most polluted urban site showed the predominance of the insoluble fractions of elements which are known tracers of anthropogenic emission sources [45]. They are considered reliable tracers for soil resuspension, which can be favored, as in this specific case, by the passage of vehicles, particularly during warm and dry periods [46,47]. Thus, we can argue that that the resuspension/deposition of road dust from traffic might have contributed to the maximum magnitude of microbial exchanges observed at Road Traffic site and its high observed airborne diversity. However, in agreement with Lymperopoulou et al. [56], road dust appeared in turn to be a consistent but not the absolute dominant source of microbes to the near-surface atmosphere. The detection of two main tracers of sea salt source [48] and three main tracers of biomass burning emissions [3], suggested a putative roles of bubble bursting associated with sea spray and biomass burning as putative sources of airborne bacteria and fungi in Rome. The long-distance transportation of dust-associated microorganisms by wind from the surrounding Southern rural and coastal geographic zones might has affected the bacterial and fungal community composition in the atmosphere. In particular, Camacho et al. [65] revealed that the concentration of fungal spores increased in the atmosphere in around two weeks after an urban or forest fire occurrence. Convective clouds associated to the intense heat from wildfires promoted the emission of fungal spores from their original habitats and the subsequent germination of fungal colonizers of burnt areas. The European Forest Fire Information System (EFFIS) registered > 70 of actively burning fires from 1 August to 30 September 2020, mostly located 80-170 Km South of Rome (https://effis.jrc.ec.europa.eu/, accessed on 28 December 2021). We cannot rule out that the contribution of combustion-aerosolized and viable microbial organisms influenced the observed airborne communities.

Functional Profiles: Distribution of Extremotolerant and Human-Pathogenic Microbes across Urban Sites
Change of pathogenicity of airborne microorganisms with worsening of the air quality have been highlighted indirectly in epidemiologic studies and directly by analyzing airborne microbiome characteristics. An increment of total respiratory outpatient visits in China [66] and US [67] for both, children and adults, have been associated to prolonged exposure to elevated PM levels. Increased risk of influenza hospitalizations was linked to source-specific particulate matter originated from traffic and other combustion sources [68]. In agreement with these data, Korneykova et al. [23] found the largest proportion of airborne opportunistic fungi in traffic and residential areas and lowest ones in green urban zones. Most dangerous fungi were detected only in the air of traffic areas of Moscow (temperate climate) and Murmansk (subarctic climate) [23]. In the present study, the relative abundance of the pathogenic fungi in total airborne and phylloplane mycobiome increased with the air quality worsening, following the increasing trend of pollution from Urban Park to Traffic Road. In parallel, PM10 was significantly correlated with fungal and bacterial composition, in both cases affecting the relative abundance of extremotolerant microbes, especially for epiphytic communities. Indeed, urban environment is characterized by multiple anthropogenic stressors (pollution, heat islands, limited water, etc.) hence, providing conditions for selection and proliferation of specialized microorganisms such as poly-extremotolerant microbes [35]. In particular, the contribution of the phyllosphere microbiome to degrade air pollutants (e.g., aromatic compounds) via oxidation mechanisms has been recently recognized [58]. Significant co-occurrence of extremotolerance (e.g., osmotolerance and psychrotolerance) and mammal opportunism was shown at the level of fungal orders in a kingdom-wide phylogenetic analysis [35]. To infect the human body, opportunistic and pathogenic fungi should be capable to overcome multiple barriers such as elevated body temperature, unfavorable pH and humidity, epithelial and mucus obstacles [69,70]. In case of true pathogens, having mammals hosts as a primary habitat, resistance mechanisms have been developed for allowing successful infections. The opportunist fungi are not specialized for infection of mammals but their acquired adaptability and stress tolerance allow the establishment in the host [71]. For instance, Aureobasidium pullulans, the most abundant opportunistic-pathogenic fungus in our phylloplane samples, is capable to colonize almost all environmental niches including soil, water, air, and limestone due to its thermotolerance and osmotolerance [72]. It has been reported as one on the most abundant fungal colonizers of phyllosphere in several plant species as olive trees [61]. Aureobasidium pullulans had also emerged as an opportunistic human pathogen, especially among immunocompromised patients [72]. Our data confirmed that such characteristics can contribute to microbial survival under conditions of heat, dryness and solar irradiation in harsh-polluted urban environments, but also ensure the possibility of opportunism. However, this interconnection between extremotolerance and opportunism has not been clearly detected in the bacterial microbiome in either substrate. In accordance with data recorded in Beijing during severe smog event [73], part of our inhalable microorganisms was soil-associated and nonpathogenic. We did not unveil a gradual enhance of the airborne opportunistic pathogenic bacterial fraction in response to increased levels of PM as reported in previous studies [19][20][21][22][23]. We can speculate that lower variation of the air quality between different sites of Rome (PM mass concentration doubled from Urban Park to Traffic Road) in respect to a wide range of pollution (5 to 20 times increase) covered by seasonal sampling in these studies may explain the observed differences in pathogenic patterns.

Supplementary Materials:
The following are available online at www.mdpi.com/article/10.3390/at-mos13020224/s1, Figure S1: Rarefaction curves of OTUs abundance grouped by types of samples (leaf surface of Q. ilex trees, aerosol and road dust) for bacteria (a) and fungal (b) communities sampled in three urban sites of Rome. Solid curves represent the observations while the dashed ones show the interpolation. The dashed vertical line indicates the number of sequences subsampled from each sample, 89081 for bacteria and 59622 for fungi; Figure S2: The relative abundance (%) of the most pronounced different classes (a) and orders (b) of bacteria recovered in three urban sites of Rome (Urban Park, Residential and Road Traffic) and in three sample types (Leaf Surface, Aerosol and Road Dust). The classes/orders with relative abundance ≥ 1% of the reads per given sample were included, accounting for ≥ 40% of all sequences recovered; Figure S3: The relative abundance (%) of the most pronounced different classes (a) and orders (b) of fungi recovered in three urban sites of Rome (Urban Park, Residential and Road Traffic) and in three sample types (Leaf Surface, Aerosol and Road Dust). The classes/orders with relative abundance ≥ 1% of the reads per given sample were included, accounting for ≥ 35% of all sequences recovered; Figure S4: Percentage contribution of the life-style associated with the fungal species identified in three sample types, Leaf Surface, Aerosol (PM10) and Road Dust, collected in the Urban Park, Residential and Road Traffic sites of Rome. Ubiquitous = U, H = herbaceous-plant-associated; W = woody; sapro = saprophytic; para = parasitic; inhab = surface inhabiting; Figure S5: Venn diagrams illustrating the number of unique and shared bacterial (a) and fungal (b) OTUs among three sample types, leaf surface samples of Q. ilex, aerosol (PM10), and road dust detected in three urban sites of Rome: green area (Urban Park), residential area (Residential) and polluted area close to the traffic roads (Road Traffic); Figure  S6: Linear regression analysis between normalized PM10 mass concentration and two normalizedalpha-diversity estimators, Species Richness (S*) (a) and exponential of Shannon entropy ( 1 D*) (b) computed in the three urban sites for bacterial and fungal communities of leaf surface, aerosol (PM10) and road dust. Statistical significance of R-squared coefficient was reported (ns = Not significant, * p < 0.05, ** p < 0.01, *** p < 0.001); Figure S7: Heatmap of Spearman's rank correlation between the chemical species used as pollution tracers and the relative abundance of extremotolerant microbes and human-pathogenic fungi in three sample types: Leaf Surface, Aerosol (PM10) and Road Dust. Only significant Spearman's rank correlation values for p < 0.05 were reported. Table S1: Frequency distribution of bacterial classes/orders detected in three urban sites of Rome (Urban Park, Residential and Road Traffic) and in three sample types (Leaf Surface, Aerosol and Road Dust); Table S2: Frequency distribution of fungal classes/orders detected in three urban sites of Rome (Urban Park, Residential and Road Traffic) and in three sample types (Leaf Surface, Aerosol and Road Dust). Table S3: List of 77 extremotolerant bacteria detected in three sample types, Leaf Surface, Aerosol (PM10) and Road Dust, collected in the Urban Park (Park), Residential (Res) and Road Traffic (Traf) sites of Rome. Isolation_sources were collected from published sources, "The Bacterial Diversity Metadatabase BacDive" (Reimer et al. 2019) and NCBI nucleotide database. Table S4: List of 83 extremotolerant fungi detected in three sample types, Leaf Surface, Aerosol (PM10) and Road Dust, collected in the Urban Park (Park), Residential (Res) and Road Traffic (Traf) sites of Rome. "Type-classification" was assigned from MycoBank and NCBI nucleotide databases. Table S5: The opportunistic human pathogenic fungi isolated in the investigated type of samples (Leaf, Surface, aerosol PM10 and Road Dust) in three urban sites of Rome (Urban Park, Residential and Road Traffic). The identification of human pathogenic fungi, Biosafety Levels (BSL) classification and the description of the human pathogenicity were performed using MycoBank Database and the Atlas of Clinical Fungi [10]. Table S6: Mean concentration ± standard deviation (ng/cm 2 ) of 34 elemental components of PM deposited on the Quercus ilex leaf surfaces collected at the three sampling sites. Table S7: Concentration (ng/m 3 ) in the water-soluble and insoluble fractions of 38 elemental components of atmospheric PM10 samples collected at the three sampling sites. The limit of detection (LOD) for each element was also reported (-= under the limit of detection). Table S8: Mean concentration ± standard deviation (ng/mg) in the water-soluble and insoluble fractions of 31 elemental components of road dust samples collected in three sampling sites. Table S9: Mean concentration ± standard deviation (ng/mg) in the water-soluble and insoluble fractions of 31 elemental components of road dust samples collected in three sampling sites.