Gelatin filter capture-based high-throughput sequencing analysis of microbial diversity in haze particulate matter

Meiqing Sun, Zhanlin Ding, Hong Wang, Guangping Yu, Zhe Feng, Bingzhi Li, Penghui Li

1Wuqing District Center for Disease Control and Prevention, Tianjin 301700, China;

2Key Laboratory of Systems Bioengineering (Ministry of Education), School of Chemical Engineering and Technology,Tianjin University, Tianjin 300072, China;

3School of Environmental Science and Safety Engineering, Tianjin University of Technology, Tianjin 300384, China.

Abstract Airborne particulate matter (PM), especially PM2.5, can be easily adsorbed by human respiratory system. Their roles in carrying pathogens for spreading epidemic diseases has attracted great concern. Herein, we developed a novel gelatin filter-based and culture-independent method for investigation of the microbial diversity in PM samples during a haze episode in Tianjin, China. This method involves particle capture by gelatin filters, filter dissolution for DNA extraction, and high-throughput sequencing for analysis of the microbial diversity. A total of 584 operational taxonomic units (OTUs) of bacteria and 370 OTUs of fungi at the genus level were identified during hazy days. The results showed that both bacterial and fungal diversities could be evaluated by this method.This study provides a convenient strategy for investigation of microbial biodiversity in haze, facilitating accurate evaluation of airborne epidemic diseases.

Keywords: airborne particulate matter, gelatin filter, pathogen, high-throughput sequencing

Introduction

We breathe in, more or less, airborne particulate matter (PM),e.g., PM2.5and PM10. As compared to PM10that can be readily removed by mucosal clearance, PM2.5[1-3]can be adsorbed into human pulmonary parenchyma owing to their small sizes.Moreover, the higher specific surface area of PM2.5facilitates their role as carrier for numerous contaminants, such as water soluble inorganic ions,heavy metals, polycyclic aromatic hydrocarbons(PAHs), pathogens and bioaerosols[4-10]. After entering the human body, these particles can release antigenic compounds, microbial toxins, and viruses into the bloodstream, leading to serious cerebrovascular and cardiovascular diseases[11-12].

There is evidence suggesting that pathogenic microorganisms bound onto the particulate matter,such asStreptococcus pneumoniae,Streptococcus pyogenes,Mycoplasma pneumoniae,Haemophilus influenzae,Klebsiella pneumoniae,Pseudomonas aeruginosaandMycobacterium tuberculosis, are becoming a threat to public health[13]. Inhalation of these pathogens may cause allergic reactions and severe respiratory infections. High concentrations of these pathogens also cause skin diseases and even cancers[14]. Accordingly, their potential impacts on the environment and the human health have attracted great concerns.

Up to now, culture-based methods remain the prevalent strategy for the investigation of microbial diversity and evaluation of pathogen risks in airborne samples[15]. For example, Fanget alidentified 789 airborne bacteria, which were distributed across 55 genera and 184 species according to a culturable method[16]. However, most of the microorganisms cannot be cultured by common growth conditions,resulting in underestimation of microbial complexity in these samples. Therefore, it is urgent to develop culture-independent techniques for the investigation of airborne microorganisms. In terms of airborne sampling methods, traditional sampling processes involve air pump-mediated compression of air,followed by impact of microbes on the surface of quartz filters[17]. However, the impact process may lead to the escape of abundant microbes from the filters, also resulting in undervaluation of microbial integrity[18].

In this paper, we reported an approach by developing a novel gelatin filter-based and cultureindependent method for the investigation of the microbial diversity in PM samples. This method involves rapid gas compression by air pumps,thorough particle capture on gelatin filters, convenient dissolution of the filters for DNA extraction, and highthroughput sequencing for analysis of the microbial diversity. Hence, this study provides a new and convenient strategy for the investigation of biodiversity in haze, facilitating accurate evaluation of airborne epidemic diseases.

Methods

Sample collection

The PM samples were obtained each day by MD8 Air Scan sampling device (Air Scan Sartorius AG,Gottingen, Germany) with the average flow rate of 40 L h-1, through the sterile gelatin filters (diameter of 80 mm, pore size of 3 μm, 17528-80-ACD, Sartorius)during November 30 to December 21 of 2016 in different functional regions of Wuqing District,Tianjin (Table 1), including: the control groups sampled in clean days (Cle1 from hospitals, and Cle2 from traffic hubs), Group 1 sampled from hospitals(Hosp1, Hosp2, and Hosp3), Group 2 sampled from traffic hubs(Trans1, Trans2, and Trans3), Group 3 sampled from schools (Sch1, Sch2, and Sch3), and traffic hubs (Tran1, Tran2, and Tran3). Initially, we also sampled the materials in clean days from the school as the third control (Cle3). Unfortunately, the sequencing of this sample was failed. Therefore, we did not include this control sample to the manuscript.All the samplers were placed 1.5 m above the ground.The variations of PM2.5and PM10concentrations during the sampling periods were made available online (China national environmental monitoring center). The gelatin filters were stored at -20 °C prior to extraction of total DNA.

Table 1 Sampling groups of the different functional areas

Electron microscopy

The morphology of PM samples on the gelatin membranes was characterized by a Phenom scanning electron microscope (SEM) supported with image software (SEM, ProX, Holland).

DNA extraction

The gelatin filters were dissolved in 5 mL sterile water (preheated to 37 °C). The pathogenic microorganisms captured by gelatin filter were lysed by the alkaline lysing liquid (NaOH 50 mmol/L, SDS 1%,Protease K 10 mg/L, RNase 20 mg/L), the total DNA was enriched using magnetic nanoparticles modified by polyquaternary amino salt polymers, and the total DNA was eluted by elution (Jinping Biotech., China)according to the instruction. The extracted DNA was stored at -80 °C for further use.

Amplification of 16S and 18S rDNA gene sequences

The 16S rDNA V3+V4 region was amplified using the forward primer 338F (5′-ACTCCTACGGGAGGC AGCA-3′) and the reverse primer 806R (5′-GGACTAC HVGGGTWTCTAAT-3′), the 18S rDNA ITS region was amplified using primers ITS1 (CTTGGTCAT TTAGAGGAAGTAA) and ITS2 (GCTGCGTTCTTC ATCGATGC). The primers were synthesized by Sangon Ltd., Beijing. The PCR reaction procedure is shown inSupplementary Table 1, available online.

High-throughput sequencing of 16S rDNA

The obtained total 16S rDNA and 18S rDNA were sequenced by double end sequencing (Paired-End)method through FLASH v1.2.7 using the Illumina HiSeq 2500 sequencing platform (BioMarker Technologies Co., Ltd., Beijing, China)[19].Bioinformatics analysis of sequences was conducted using the QIIME (V1.7.0) software package.Sequences with similarities greater than or equal to 97% were grouped into operational taxonomic units(OTUs). The Shannon index was used to estimate the biodiversity of bacteria in a single sample[20]. The principal component analysis (PCA) and the cluster heatmap analysis were also performed to assess the bacterial composition of the samples[21].

Results

Concentration and diameter of fine particles during the haze events

Fig. 1shows the variations of PM2.5and PM10concentrations during the sampling periods. The concentrations of PM2.5were higher than 200 μg·m-3in 8 days, and the concentrations of PM10were higher than 200 μg·m-3in 10 days, with the highest PM10concentration being almost 350 μg·m-3, which far exceeded the PRC National Standard PM standards:50 μg·m-3for PM10, and 35 μg·m-3for PM2.5,respectively (PRC National Standard, 2012). As shown inFig. 2, the statistical mean diameters of the sampled particles were in the ranges of 100-900 nm and 1 -2.5 μm, confirming that the particles were in the scope of fine particulate matter.

Fig. 1 Daily average PM2.5 and PM10 concentrations estimated from the samples collected from November 30 to December 21, 2016.

Fig. 2 Scanning electron microscope characterization of the gelatin filter after PM sampling in different functional areas and the control group, where the green marks indicated the diameter of fine particles.

Microbial species richness and diversity

Rarefaction curves derived from the observed OTU number and shannon index were flattened (Fig. 3),showing that our sequencing depth was sufficient to cover the vast majority of bacteria and fungi in the samples. The statistics of OTU species in different ranks of the bacteria and the fungi are shown inTable 2and3, and a total of 584 OTUs of bacteria and 370 OTUs of fungi at the genus level were identified during hazy days.

Fig. 3 Rarefaction curves and Shannon index derived from the observed OTU numbers, (A) and (B) for bacteria, (C) and (D) for fungi.

Table 2 Statistics of OTU species of bacteria in different ranks

As shown inFig. 4, in all samples,Proteobacteriawas the most abundant phylum, and four other dominant phyla wereFirmicutes,Bacteroidetes,Actinobacteria, andCyanobacteria(Fig. 4A). Compared to the samples from the clear days, the samples from hospitals and traffic hubs during hazy days contain much higher levels of human health-related bacteria, such asAcinetobacter,Staphylococcus,Corynebacterium,Lactobacillus, andDuganella(Fig. 4B). Especially,StaphylococcusandCorynebacteriumwere the predominant pathogens around the high-density-population functional areas, such as hospital and transportation areas. For example,Staphylococcuson PM collected around the hospital areas was the dominant bacteria (18% in Hosp2 samples), which may be derived from the high density of patients and health care personnel[22]. Nevertheless,the average abundance ofStaphylococcusin Sch2 was only 1.2% because of the regular population. The dominant bacteria wereCorynebacterium(34.3%) in Tran2 andLachnospiraceae_Nk4A136(10.5%) in Tran3. In Sch3 areas the dominant bacteria wereLactobacillus(13.5%) andKlebsiella(11.8%).Interestingly, the genusSphingomonasis detected on both the clear days and hazy days, which are the environment-associated bacteria widely found in surface water, the rhizosphere, sediments, and even soils.Fig. 4Cshows that the dominant fungi wereMalassezia(58.0%) in Hosp1 and (54.9%) in Hosp2,Alternaria(55.5%) andCladosporium(9.17%) in Tran3. The genusMalasseziahas been associated with a number of diseases affecting the human skin, such as pityriasis versicolor,Malassezia(Pityrosporum)folliculitis, seborrheic dermatitis and dandruff, atopic dermatitis, and psoriasis.Malasseziayeasts are a part of the normal micro-flora, but under certain conditions they can cause superficial skin infection[23].Cladosporium,AspergillusandAlternariacould enter the deep lung and cause respiratory diseases[24].

Fig. 4 The bacterial richness and diversity of the 16S rDNA samples at the phylum level (A) and (B) at the genus level. C: The fungal richness and diversity at the genus level. Only nodes with the first-tenth relative abundance are labeled. D: The PCA analysis of different functional areas at the genus level.

The variation of bacteria in collected samples were analyzed by two methods: the PCA analysis and the clustering heatmap analysis. PCA analysis showed that the high-population-density functional areas(Hosp1, Hosp2, Tran2, and Tran3 but not Tran1) had a similar bacterial diversity (Fig. 4D). However, the samples from Trans1 and Sch1 displayed little similarity in the bacterial composition. The clustering heat map analysis also demonstrated the similarity of the bacterial composition between the hospital areas and the transportation areas (Fig. 5). One cluster was composed of the hospital areas and the transportation areas; another cluster was composed of the samples from control and the school areas (samples Cle1, Cle2,Sch1, and Sch3). The bacterial diversity around school areas (Sch1 and Sch2) is similar to the control group,but not always the case: the sample Sch3 showed partially similar bacterial diversity to the samples from the hubs and hospitals, which may be attributed to transportation of PM from hubs and hospitals.

Fig. 5 The cluster heatmap analysis of different functional areas.

Table 3 Statistics of OTU species of fungi in different ranks

Discussion

Numerous epidemiologic studies have documented that PM is associated with inflammation-related diseases. For example, PM may result in the alteration of immune functions, such as IgA, IgG, IgM, IgE and lymphocyte profiles (CD4 T cells, CD8 T cells,CD4/CD8 T cells) in blood[25], which might link to various adverse health effects, including asthma,pulmonary infectious diseases, diabetes mellitus and obstructive bronchitis[26]. In China, numerous studies have been conducted to monitor regional air pollution[27], and the exposure-response relationships between airborne PM and cardiovascular diseases. All these studies have focused on the particle concentration effects, the airborne gas and other chemical substances[28]. However, little information has been available on the association of particle-bound pathogenic bacteria with the prevalence of infectious diseases and epidemics. In this paper, we found that the amounts ofStaphylococcusandCorynebacteriumare high on the PM collected around hospital areas and traffic hubs.Staphylococcusmay cause both endemics and epidemic infections, such as neumonia,pseudomembranous colitis, pericarditis, and sepsis,especially outbreak of severe acute pneumonia[29].Corynebacteriumgenus is also able to cause various types of healthcare-associated infections in immunocompromised hosts, such as pneumonia,pharyngitis and epidemic disease[30]. ExceptStaphylococcusandCorynebacterium, other prevailing pathogenic bacteria areAcinetobacter,Pseudomonas, which often cause large multifacility,nosocomial outbreaks and are frequently evolved into drug-resistant bacteria[31]. Our data indicate that densely populated traffic areas show strong evidence of bacteria and the presence of some potentially pathogenic organisms. In addition, PM2.5exposure may alter and impair the normal immune responses of the lung, and decrease the phagocytosis of alveolar macrophages through disrupting the normal physical and immunological function of the lung, induce disorder of inflammatory cytokine[32]. Plus these more active pathogens, all of these effects would lead to the decline in immunity and facilitate infectious diseases.4-7

In conclusions, the pathogenic microorganisms in PM are closely related to their surrounding environmental conditions and people's density,mobility and activities. In this study, we developed a convenient culture-independent method for the investigation of the microbial diversity in PM samples, which is composed of gas compression,particle capture, DNA extraction, and high-throughput sequencing. The results showed that both bacterial and fungal diversities could be accurately evaluated by this method. The results demonstrated the presence of some pathogenic bacteria, such asStaphylococcus,Corynebacterium,AcinetobacterandPseudomonas,which may affect immunocompromised populations(e.g., the elderly, children and postoperative convalescence patients). This study provides a new strategy for the evaluation and surveillance of airborne epidemic diseases that increasingly threaten urbanized human societies.

Acknowledgments

This work is supported by Project of Science and Technology Development in Wuqing District, Tianjin(No. WQKJ201614), Tianjin 131 innovative talent training project, Postdoctoral Science Foundation.