Front. Immunol. Frontiers in Immunology Front. Immunol. 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2019.02714 Immunology Original Research Pulmonary Exposure to Magnéli Phase Titanium Suboxides Results in Significant Macrophage Abnormalities and Decreased Lung Function McDaniel Dylan K. 1 Ringel-Scaia Veronica M. 2 Morrison Holly A. 1 Coutermarsh-Ott Sheryl 1 Council-Troche McAlister 3 Angle Jonathan W. 4 Perry Justin B. 5 Davis Grace 5 Leng Weinan 6 Minarchick Valerie 7 Yang Yi 8 Chen Bo 9 Reece Sky W. 10 Brown David A. 5 Cecere Thomas E. 1 Brown Jared M. 7 Gowdy Kymberly M. 10 Hochella Michael F. Jr. 11 Allen Irving C. 1 2 * 1Department of Biomedical Sciences and Pathobiology, VA-MD College of Veterinary Medicine, Virginia Tech, Blacksburg, VA, United States 2Graduate Program in Translational Biology, Medicine and Health, Virginia Tech, Blacksburg, VA, United States 3Analytical Research Laboratory, VA-MD College of Veterinary Medicine, Virginia Tech, Blacksburg, VA, United States 4Department of Materials Science and Engineering, Virginia Tech, Blacksburg, VA, United States 5Department of Human Nutrition, Foods, and Exercise, Virginia Tech, Blacksburg, VA, United States 6National Center for Earth and Environmental Nanotechnology Infrastructure, Virginia Tech, Blacksburg, VA, United States 7Department of Pharmaceutical Sciences, University of Colorado, Anschutz Medical, Aurora, CO, United States 8Key Laboratory of Geographic Information Science of the Ministry of Education, School of Geographic Sciences, East China Normal University, Shanghai, China 9Department of Applied Physical Sciences, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States 10Department of Pharmacology & Toxicology, Brody School of Medicine, East Carolina University, Greenville, NC, United States 11Department of Geosciences, Virginia Tech, Blacksburg, VA, United States

Edited by: Janos G. Filep, Université de Montréal, Canada

Reviewed by: Nithyananda Thorenoor, College of Medicine, Pennsylvania State University, United States; Juli Bai, The University of Texas Health Science Center at San Antonio, United States

*Correspondence: Irving C. Allen icallen@vt.edu

This article was submitted to Molecular Innate Immunity, a section of the journal Frontiers in Immunology

28 11 2019 2019 10 2714 16 07 2019 05 11 2019 Copyright © 2019 McDaniel, Ringel-Scaia, Morrison, Coutermarsh-Ott, Council-Troche, Angle, Perry, Davis, Leng, Minarchick, Yang, Chen, Reece, Brown, Cecere, Brown, Gowdy, Hochella and Allen. 2019 McDaniel, Ringel-Scaia, Morrison, Coutermarsh-Ott, Council-Troche, Angle, Perry, Davis, Leng, Minarchick, Yang, Chen, Reece, Brown, Cecere, Brown, Gowdy, Hochella and Allen

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Coal is one of the most abundant and economic sources for global energy production. However, the burning of coal is widely recognized as a significant contributor to atmospheric particulate matter linked to deleterious respiratory impacts. Recently, we have discovered that burning coal generates large quantities of otherwise rare Magnéli phase titanium suboxides from TiO2 minerals naturally present in coal. These nanoscale Magnéli phases are biologically active without photostimulation and toxic to airway epithelial cells in vitro and to zebrafish in vivo. Here, we sought to determine the clinical and physiological impact of pulmonary exposure to Magnéli phases using mice as mammalian model organisms. Mice were exposed to the most frequently found Magnéli phases, Ti6O11, at 100 parts per million (ppm) via intratracheal administration. Local and systemic titanium concentrations, lung pathology, and changes in airway mechanics were assessed. Additional mechanistic studies were conducted with primary bone marrow derived macrophages. Our results indicate that macrophages are the cell type most impacted by exposure to these nanoscale particles. Following phagocytosis, macrophages fail to properly eliminate Magnéli phases, resulting in increased oxidative stress, mitochondrial dysfunction, and ultimately apoptosis. In the lungs, these nanoparticles become concentrated in macrophages, resulting in a feedback loop of reactive oxygen species production, cell death, and the initiation of gene expression profiles consistent with lung injury within 6 weeks of exposure. Chronic exposure and accumulation of Magnéli phases ultimately results in significantly reduced lung function impacting airway resistance, compliance, and elastance. Together, these studies demonstrate that Magnéli phases are toxic in the mammalian airway and are likely a significant nanoscale environmental pollutant, especially in geographic regions where coal combustion is a major contributor to atmospheric particulate matter.

cytotoxicity air pollution nanoparticle TixO2x−1 in vivo environmental exposure

The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH, NSF, EPA, or any other funding agency.

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Coal combustion is a significant component of the global energy portfolio and accounts for approximately 30% of the world's overall energy needs (1). The global use of coal has continued to increase over the last half century, especially in developing countries, such as China and India, due in part to abundant supplies and favorable economic factors (1). The burning of coal on such a massive global scale has also resulted in a range of detrimental environmental and human health concerns. While the long-term effects of increased coal burning, such as rising greenhouse gas emissions and climate change are clearly a concern, the much more visible short-term consequences associated with air pollution must also be addressed. Coal-burning is a significant contributor to atmospheric particulate matter with aerodynamic diameter <2.5 μM (PM2.5), which is the most concerning fraction for human health (2).

      The smallest fraction of PM2.5 is ultrafine nanoscale PM associated with air pollution. These nanoparticles have highly significant biological impacts and can have detrimental effects on human health (3). These mostly incidental nanomaterials are unintentionally produced as byproducts of anthropogenic processes and typically enter the environment shortly after generation (4). Until recently, the majority of incidental nanomaterials either went unnoticed or were not evaluated for potential detrimental impacts on environmental health. However, we now know that many of these nanomaterials can contribute disproportionately to environmental chemistry impacting multiple earth systems (4). Likewise, due to their nano-scale size, these particles can readily penetrate many of the physiological and cellular barriers used to protect biological systems. For example, depending on size, shape, and composition, many nanoscale particles can readily enter the airway and gastrointestinal tract of mammals, translocate epithelial cell barriers, access the bloodstream, and affect the function of vital body systems (3). The convergence of nanoscience, environmental health, and biomedical research are beginning to provide significant insight into the impacts of incidental nanomaterials across earth systems and in human health (4).

      As evidence of the progress being made in identifying and characterizing new incidental nanomaterials, recent work identified novel titanium suboxides, defined as Magnéli phases (TixO2x−1), which are generated during coal combustion (1). TixO2x−1 is produced during coal combustion at specific temperatures and oxygen fugasities from stoichiometric TiO2 titanium oxide that is a common accessory mineral in nearly all coals worldwide (1). Magnéli phases were originally discovered during the investigation of a coal ash spill into the Dan River (North Carolina, USA), where these novel titanium suboxides were found downstream from the spill site (5). Magnéli phases were observed in the size ranges of a few tens to hundreds of nanometers and exhibit a unique electron diffraction pattern compared to anatase and rutile, the two most common titanium oxide minerals found in coal (5). Since the discovery of Magnéli phases in this riverine environment, subsequent analyses of coal ash samples from power plants throughout the United States and China, each utilizing various types and compositions of coal, revealed that Magnéli phases are widespread and were present at every site tested (5). Magnéli phases were found with compositional ranges between Ti4O7 and Ti9O17, with the most frequent being Ti6O11 (5). Subsequent analyses of power plant stack emissions, sludge from waste water treatment plants, and road dust from Chinese urban areas were all found to contain Magnéli phases (5).

      Together, these data illustrate that Magnéli phases generated as incidental nanoparticles through the combustion of coal are widespread in the environment. However, there is currently a paucity of data related to the in vivo relevance and physiological effects of Magnéli phases in the mammalian respiratory system. In the current manuscript, we demonstrate that Magnéli phases are concentrated and ultimately sequestered in lung associated macrophages. Magnéli phase phagocytosis significantly impairs mitochondrial function and stimulates reactive oxygen species (ROS) production by the macrophages. Ultimately, these trigger pathways associated with apoptosis and lung injury. Consistent with these findings, mice chronically exposed to Magnéli phases demonstrate significantly decreased lung function. Together, these data reveal the significant impact of these incidental nanoparticles on overall respiratory function and provide further evidence of the need for improved environmental monitoring to screen for these and similar materials.

      Materials and Methods Magnéli Phase Fabrication and Characterization

      Magnéli phases were synthesized using a tube furnace (diameter = 8.9 cm) with a heating and cooling rate of 5°C min−1 and an N2 atmosphere (flow rate = 0.28 m3 min−1) as previously described (1). Heating and cooling processes were isothermal at the target temperature for 2 hours (h). The process includes heating pulverized coal with TiO2 nanoparticles. Magnéli phases were produced using commercial P25 nanoparticles, which is a mixture of the 80% anatase and 20% rutile forms of TiO2. Magnéli phase samples were characterized using a scanning transmission electron microscope operating at 200 kV and equipped with a silicon drift detector-based Energy Dispersive X-ray Spectroscopy (EDS) system as previously described (1).

      Experimental Animals

      All mouse studies were approved by the respective Institutional Animal Care and Use Committees (IACUC) at Virginia Tech and East Carolina University, and were conducted in accordance with the Federal NIH Guide for the Care and Use of Laboratory Animals. All studies utilized age and gender matched wild type C57Bl/6J mice.

      Cytokine and Lactate Acid Dehydrogenase (LDH) Assessments in Bone Marrow-Derived Macrophages (BMDMs)

      BMDMs were harvested from mice using standard protocols (6). Briefly, harvested bone marrow cells were cultured for 6 days in Dulbecco's Modified Eagle Media (DMEM) supplemented with 10% fetal bovine serum (FBS), L-Glutamine, 1 × penicillin/streptomycin, and 20% L929-conditioned medium. After this 6-day incubation, BMDMs were re-plated at 250,000 cells/well and were allowed to adhere to the plate for 24 h. Cells were then treated with different concentrations of Magnéli phase particles (0, 1, 10, 100, or 1,000 ppm) overnight. Following these treatments, cell-free supernatant was collected, and cytokine levels were assessed via enzyme-linked immunosorbent assay (ELISA). Lipopolysaccharide (LPS) was used at 1 μg/mL as a positive control for cytokine release. LDH activity was measured in the supernatant to assess cytotoxicity.

      Transmission Electron Microscopy to Assess Magnéli Phase Phagocytosis

      Following treatment with Magnéli phases, treatment media was discarded and BMDMs (1 × 106 cells/well) were fixed with 2.5% Gluteraldehyde in 0.1 M Sodium Cacodylate directly in the cell culture plate for at least 1 h. Following fixation in 2.5% glutaraldehyde in 0.1 M Na cacodylate, samples were washed two times in 0.1 M Na cacodylate for 15 min each, post-fixed in 1% OsO4 in 0.1 M Na cacodylate for 1 h and washed two times in 0.1 M Na cacodylate for 10 min each. Fixative was removed and cells were dehydrated with increasing concentrations of ethanol for 15 min each as follows (15, 30, 50, 70, 95, and 100%). The dehydration was completed with propylene oxide for 15 min. The samples were then infiltrated with a 50:50 solution of propylene oxide:Poly/Bed 812 (Polysciences Inc., Warrington, Pennsylvania, USA) for 6–24 h followed by complete infiltration with a 100% mixture of Poly/Bed 812 for 6–12 h. The samples were embedded in fresh 100% Poly/Bed 812 using Beem embedding capsules (Ted Pella Inc., Redding, California, USA) overnight. These were then placed into a 60°C oven for at least 48 h to cure and harden. The embedded samples were trimmed and thin (60–90 nm) sections were cut and collected on copper grids (Electron Miscroscopy Sciences). These sections were stained with aqueous uranyl acetate and lead citrate and examined and photographed using a JOEL JEM 1400 transmission electron microscope.

      Mitochondrial Stress Test

      Oxygen Consumption Rate (OCR) of BMDMs after exposure to nanoparticles was determined using a Seahorse XF96 Extracellular Flux Analyzer (Agilent Technologies). Seahorse 96-well cell culture microplates were seeded at 100,000 cells per well and subjected to a 24-h incubation with Magnéli phases. Following the incubation, cells were washed with OCR assay media (1 mM pyruvate, 2 mM glutamine, and 10 mM glucose) twice and then immersed in a total volume of 180 μL of media in each well immediately before the assay. Plates were then placed in the XF96 to establish a basal level of respiration. Afterwards, three separate injections were performed: the ATPase inhibitor oligomycin (1 μM), the mitochondrial uncoupler carbonilcyanide p-triflouromethoxyphenylhydrazone (FCCP) (3 μM) and the mitochondrial electron transport chain complex III inhibitor antimycin A (2 μM), respectively. Each condition was measured three times before moving to the next injection. Data is represented as pmol of oxygen per minute.

      Mitochondrial Membrane Potential and ROS Assessments

      Cells were grown following the conditions described above and plated at 100,000 cells per well in black walled, clear, flat bottom Corning Costar 96-well microplates. Images were collected in three different channels with a 60 × objective using Hoechst 33258 (100 nM) to counterstain nuclei, dichlorofluorescein (DCF) (100 μM) to monitor reaction to H2O2 injury, and tetramethylrhodamine, methyl ester (TMRM) (10 nM) to assess mitochondrial membrane potential. The following excitation/emission filters were selected to collect the fluorescence signal for each channel on the instrument; DAPI (Excitation 390.0 nm/18.0 nm bandpass, Emission 435.0/48.0 nm bandpass), FITC_511 (Excitation 475.0 nm/28.0 nm bandpass, 511.0 nm/23.0 nm bandpass), Cy3 (Excitation 542.0 nm/27.0 bandpass, Emission 587.0 nm/45.0 nm bandpass), respectively. Prior to the start of the imaging protocol, wells were treated with TMRM and DCF for 30 min, protected from light, in a 37°C/5% CO2 incubator and washed twice with phosphate buffered saline (PBS). Sequential images were taken in 10 fields of view for each channel in each well containing cells. After a basal level of fluorescence was obtained for each field, all wells were treated with 500 μM H2O2 to induce injury and a second round of images were obtained following the same procedure. Images were analyzed in an automated fashion using GE's InCarta software version 1.6. Multiple parameters were collected from each plate by creating custom “masks” that incorporated the fluorescent target of interest. Data are expressed as the mean pixel value under the mask minus mean pixel value for the local background or, more simply, Intensity-Background.

      Quantification of Titanium Concentration in Tissues Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

      ICP-MS analysis was performed on Agilent 7900 ICP-MS operating in helium collision mode to reduce/remove polyatomic interferences. Data was acquired in spectrum mode and utilized a 3-point peak pattern for quantitative analysis. Each data point acquisition had three replicates and utilized 250 sweeps per replicate. Titanium levels were monitored at both 47 and 49 m/z to additionally check for potential interferences at either value. The instrumental parameters for ICP-MS are shown in Supplemental Figure S2A.

      Scandium (45 m/z) was used as an on-line addition internal standard. A calibration curve that ranged from 0.04 to 10 ppb titanium in solution was used to quantify titanium in the samples according to their response ratio with the scandium internal standard (Ti/Sc). An Anton-Parr Multiwave GO microwave system that utilized modified Teflon (PTFE-TFM) microwave vessels was used in the microwave digestions. This microwave oven provides feedback from the digestion conditions based on temperatures and whether excessive venting was detected.

      Ultrapure water was purified with a Millipore MilliQ water purification system set at 18.2 MOhm. Single element stock solutions of 1,000 ppm titanium (Ti) and scandium (Sc) for ICP-MS analysis were obtained from Inorganic Ventures. Nitric acid (67–70 w/w%) and hydrofluoric acid (47–51 w/w%) were both trace metals grade and obtained from Fisher Scientific. Suprapur hydrogen peroxide (30 w/w%) was also trace metals grade and obtained from Merck. Lung tissue samples were stored in a freezer at −20°C until analysis. When ready, the samples were removed from the freezer and allowed to come to room temperature. The entirety of the provided samples was weighed directly into PTFE-TFM microwave vessels and ranged from approximately 200–300 mg.

      One milliliter of ultrapure water (18.2 MΩ), 1.5 mL of 70% (w/w) nitric acid (HNO3), 0.3 mL of 10% (w/w) hydrofluoric acid (HF), and 0.2 mL of 30% w/w hydrogen peroxide (H2O2) were added to each vessel. The vessels were then capped appropriately and loaded into the carousel of the microwave digestion system. The samples were digested using a temperature gradient that ramped from room temperature to 180°C over 20 min, then held at 180°C for 10 min before cooling back to 60°C over 10 min. All samples digested clear and colorless, with no visible particulates in solution.

      Each of the digested samples was decanted into fresh, separate 15 mL polypropylene (PP) tubes, and their corresponding digestion vessels were rinsed with 2 mL of 2% (w/v) HNO3 and decanted into their respective tubes for a 5 mL total volume. The concentrated samples were then vortexed and shaken to mix immediately before being diluted 1:100 in 2% (w/v) HNO3 to ensure homogeneity of the dilute samples. The 1:100 dilution was necessary to reduce the total concentration of HF within the samples themselves so that it did not attack the quartz portions of the ICP-MS instrumentation.

      Intratracheal (i.t.) Administration and Imaging of Magnéli Phases in the Lung

      Wild-type C57Bl/6 mice were treated with one dose of 100 ppm of Magnéli phases using i.t. administration, following anaesthetization with isoflurane. Mice were euthanized using CO2 and lungs were harvested 7 days following i.t. administration. To assess chronic effects of Magnéli phases in the lung, mice were treated i.t. with 100 ppm of Magnéli phases three times per week. Lungs were harvested at 6 weeks after the initial treatment. Hematoxylin and eosin (H&E) stained lung slides were evaluated using darkfield microscopy (Cytoviva, Auburn, AL, USA) and images were collected at a magnification of 100X. Magnéli phase localization and lung histopathology were also evaluated by a board-certified veterinary pathologist using H&E stained sections prepared from formalin fixed, paraffin embedded tissues. Immunohistochemistry was conducted on the formalin fixed, paraffin embedded tissues using an APAF1 polyclonal antibody (Invitrogen; PA5-85121) and rabbit specific HRP/DAB (ABC) Detection IHC kit (Abcam).

      Gene Expression Profiling and Pathway Analysis

      Total RNA was harvested from whole lungs using FastRNA Pro Green Kit following manufacturer's protocols (MP Biomedicals). Total RNA was pooled from 3-8 individual mice per group for RT2 Profiler PCR Array Platform (QIAgen) cDNA reaction. Gene expression was evaluated using PAMM-212Z and PAMM-052Z arrays (Qiagen) following manufacturer's protocols. The gene list for these arrays and functional categories are available through the manufacturer (Qiagen). Ingenuity Pathways Analysis (IPA) and the manufacturer's array software (Qiagen) was used to analyze gene expression data and conduct pathway analysis. IPA data were ranked and evaluated based on z-score as previously described by the research team (79).

      Pulmonary Function Testing

      Male C57Bl/6J mice were exposed to PBS or Magnéli phases three times/week for 6 weeks (18 total doses). Briefly, mice were anesthetized with isoflurane and received Magnéli phases (100 ppm) in 50 μl by i.t. administration as described above for the chronic exposure assessments. Pulmonary function testing was performed 24 h following the last exposure. All mice were anesthetized with tribromoethanol (TBE) (400 mg/kg), tracheostomized, and baseline pulmonary function was recorded using the FlexiVent system (SCIREQ, Montreal, QC, Canada) (10, 11).

      Mice were ventilated with a tide volume of 10 ml/kg at a frequency of 150 breaths/min and a positive end expiratory pressure (PEEP) of 3 cm H2O to prevent alveolar collapse. In addition, mice were paralyzed with pancuronium bromide (0.8 mg/kg) to prevent spontaneous breathing. Electrocardiogram (EKG) was monitored for all mice to determine anesthetic depth and potential complications that could have arisen during testing. A deep inflation perturbation was used to maximally inflate the lungs to a standard pressure of 30 cm H2O followed by a breath hold of a few seconds to establish a consistent volume history. A snapshot perturbation maneuver consisting of a three-cycle sinusoidal wave of inspiration and expiration to measure total respiratory system resistance (R), dynamic compliance (C), and elastance (E). A Quickprime-3 perturbation, which produced a broadband frequency (0.5–19.75 Hz) over 3 s, was used to measure Newtonian resistance (Rn), which is a measure of central airway resistance. All perturbations were performed until three acceptable measurements with coefficient of determination (COD) ≥0.9 were recorded in each individual subject.

      After three baseline measurements of lung function using “snapshot” and “Quickprime-3” perturbations, mice were challenged with aerosolized saline or cumulative doses of 1.5, 3, 6, 12, and 24 mg/ml methacholine (MCh) (Sigma-Aldrich, St. Louis, MO) generated by an ultrasonic nebulizer for 10 s without altering the ventilation pattern (1012). These challenges were followed by assessments with “snapshot” and “Quickprime-3” perturbations every 30 s. Between each dose of methacholine, a deep inflation perturbation was initiated to reset lung hysteresis. For each animal, all perturbations were performed until six acceptable measurements with COD ≥0.9 were recorded.

      Replicates

      All studies were repeated at least three independent times unless noted.

      Statistical Analysis

      Single data point comparisons were evaluated by Student's two-tailed t-test. Multiple comparisons were evaluated for significance using Analysis of Variance (ANOVA) followed by either Tukeys or Newman–Keuls post-test, as appropriate. Graphs and statistical analyses were performed and generated using GraphPad Prism software 8 (GraphPad, San Diego, CA). All data are presented as mean ± the standard error of the mean (SEM) with p ≤ 0.05 considered significant. Data shown are representative of at least three independent studies.

      Results Magnéli Phases Are Cytotoxic in Murine Bone Marrow-Derived Macrophages

      In initial toxicity studies, Magnéli phases (Ti6O11) were found to be toxic to dechorionated zebrafish embryos at 100 ppm (1). Thus, we chose this formulation and dose as the focus of our subsequent studies. Utilizing previously described methods, we generated Magnéli phases that were predominately composed of Ti6O11 (1) (Figure 1A). The resultant nanoparticles were confirmed by electron microscopy analysis and X-ray diffraction patterns, as previously described (1) (Figures 1B,C). These nanoparticles have excellent light absorption in the near-infrared, UV, and visible light range (1). Likewise, the Magnéli phases also display a low amount of photocatalytic activity, as previously reported (1). The resultant nanoparticles used in our studies ranged in size from a few tens of nm to hundreds of nm (1).

      Magnéli phase phagocytosis results in increased cell death in bone marrow-derived macrophages. (A–C) Characterization of Magnéli phases used in this study. (A) Schematic illustrating Magnéli phase generation. (B,C) TEM images of Magnéli phases formed by annealing P25 TiO2 nanoparticles with coal in a pure N2 atmosphere for 2 h at 900°C. Electron diffraction patterns were characteristic of Magnéli phases and confirmed these as predominately Ti6O11 particles. Particles were between 10 and 200nm in size. (D) Un-treated bone marrow-derived macrophages (1 × 106 cells/well) and (E) macrophages treated with Ti6O11 (1, 10, 100, or 1,000 ppm) were visualized using TEM (Scale bar: 5 μm). Magnéli phases appear as punctate dark dots in the macrophages and are concentrated in phagolysosomes. (F) Macrophages containing Magnéli phases demonstrate morphological features consistent with apoptosis, including cell shrinking and membrane blebbing (red arrows) (Scale bar: 1μm). (G,H) Cytotoxicity was evaluated using trypan blue exclusion (G) across the Ti6O11 dose range and (H) at 100 ppm over a 24 h time course. (I,J) Inflammation was evaluated by assessing the production of pro-inflammatory cytokines, such as (I) IL-1β, (J) IL-6, and (K) TNF in the cell free supernatant following exposure to different doses of Magnéli phases. (L) Macrophage function was evaluated by assessing the ability of macrophages to phagocytose fluorescent Escherichia coli 24 h post-exposure to 100 ppm Magnéli phases. Data are expressed as mean ± SEM (n = 3 independent experiments). **p < 0.01.

      Once inhaled, nanoparticles are rapidly phagocytosed by resident macrophages in the lungs, which represent the predominate cell type responsible for neutralization and clearance (13). To evaluate the effect of Magnéli phases on these cells, we generated primary bone marrow derived macrophages (BMDMs) using established protocols (13). The BMDMs were treated with varying doses of Ti6O11 Magnéli phases across a 0–1,000 ppm range and over a 24-h time course (Figures 1D–H). Transmission electron microscopy (TEM) assessments revealed Magnéli phases, concentrated in phagolysosomes within treated macrophages at all concentrations and timepoints evaluated (Figures 1D,E). Morphology assessments did not appear to show high levels of macrophage activation following phagocytosis. However, at 100 ppm, a significant number of macrophages demonstrated morphology consistent with increased cell death, specifically apoptosis (Figure 1F). Several features were noted, including cell shrinkage, loss of membrane integrity, and membrane blebbing (Figure 1F). To quantify cell death, we utilized trypan blue exclusion and manual counting using a hemacytometer. Here, we observed a dose dependent increase in cytotoxicity between 10 and 1,000 ppm, ranging from 19.75 ± 2.63% to 44.25 ± 4.787%, respectively (Figure 1G). The average cell death over the same 24-h time range in the mock treated macrophages was 16.25 ± 2.062% (Figure 1G). Macrophage cell death did not initiate immediately, requiring 24 h post-exposure to peak (Figure 1H). Using 100 ppm, cell death was assessed over a time course, with almost double the number of macrophages undergoing cell death at 24 h post-exposure (15.25 ± 0.957% mock vs. 31.25 ± 3.403% Ti6O11 Magnéli phase treatment) (Figure 1G). Together, these data show that Ti6O11 Magnéli phases, at levels observed in the environment as incidental nanoparticles generated through the industrial burning of coal, are cytotoxic to mammalian macrophages.

      Nanoparticle exposure typically drives macrophage activation and the production of inflammatory mediators (13). The cells generated using the methods detailed here are skewed toward M1 polarized, pro-inflammatory macrophages (6). Likewise, the cell death induced by the Magnéli phases were originally predicted to generate damage associated molecular patterns, such as aberrant ATP production or Calcium efflux similar to that described for A549 cells (14), that should drive innate immune system signaling in the macrophages. Here, we monitored the generation of several potent pro-inflammatory mediators, including IL-1β, IL-6, and TNF (Figures 1I–K). Consistent with electron microscopy findings indicating apoptosis as the method of cell death and the lack of morphological features associated with macrophage activation, we did not observe any statistically significant increases in either gene transcription by rtPCR or at the protein level by ELISA for any pro-inflammatory mediators evaluated (Figures 1I–K). Indeed, at the gene transcription level, we observed several mediators and transcription factors actively down-regulated (Figure 2). These data indicate that cell death and inflammation induced by Magnéli phases are not linked. Despite higher levels of cytotoxicity, cell death and Magnéli phase exposure are not sufficient to significantly activate elements of the innate immune system associated with inflammation.

      Macrophage exposure to Magnéli phases activates gene expression profiles associated with apoptosis and mitochondria dysfunction. (A) Gene expression was evaluated using real time PCR based Superarrays (Qiagen). For each gene target on the array, fold change was calculated based on ΔΔCt values. Genes found to be ±2-fold change in expression from untreated specimens were defined as significant. The resultant data was evaluated using Ingenuity Pathway Analysis (IPA) to define pathways and global correlations between gene expression profiles and biological functions. Cell death signaling, specifically apoptosis, was the top pathway up-regulated in the macrophages following Magnéli phase exposure. IPA further identified gene expression profiles consistent with mitochondria dysfunction as a potential factor associated with the increased apoptosis signaling. (B) Pathways associated with inflammation and specifically inflammation associated cell death (i.e., pyroptosis and necroptosis) were significantly down-regulated at the level of gene transcription. Either fold change in gene expression or z-score values are displayed for each node as appropriate.

      Previous studies have shown that nanoparticle phagocytosis by macrophages can impair their function, resulting in greater susceptibility to a myriad of pathogens (15, 16). To evaluate macrophage functionality 24 h following either mock or Magnéli phase exposure, live cells containing nanoparticles were sorted and co-exposed to fluorescently labeled E. coli. In macrophages lacking the Magnéli phases, we observed a significant increase in mean fluorescence intensity (MFI: 10,630 ± 290) compared to those containing the nanoparticles (7,826 ± 114.2) (Figure 1L). Thus, cells containing Magnéli phases were significantly less efficient at bacterial phagocytosis, implying that nanoparticle sequestration has a deleterious impact on macrophage function.

      Magnéli Phase Phagocytosis Activates Apoptosis Signaling Pathways Modulated Through the Mitochondria

      To better define mechanisms underlying Magnéli phase cytotoxicity, we profiled gene expression and utilized Ingenuity Pathway Analysis to decipher changes in gene networks and pathways associated with a range of biological functions in macrophages. BMDMs were collected 24 h post-exposure to 100 ppm Ti6O11. Total RNA was collected and gene expression was profiled using rtPCR based gene expression arrays (Qiagen). Changes in gene expression were defined as significant if the change was ±2-fold different between the mock and treated cells. Based on the pattern of gene expression changes, apoptosis signaling was the top canonical pathway that was significantly increased following Magnéli phase exposure (Figure 2; Supplemental Figure S1A). Individually, there were no clear networks or relationships identified between the top 10 dysregulated genes (Supplemental Figure S1B). However, analysis of the congregate gene expression patterns and networks revealed a strong relationship between genes associated with apoptosis signaling with a focus on mitochondria (Figure 2A). Here, we observed significant increases in gene expression associated with both pro-apoptotic BCL-2 family members, specifically Bax (Figure 2A). Up-stream analysis revealed a significant increase in p53 signaling and down-stream analysis revealed a significant increase in Apaf1 (Figure 2A). Consistent with these findings, Caspases-3,−6, and−9 were all found up-regulated following Magnéli phase exposure, while Caspase-7 was significantly downregulated (Figure 2A). We also observed a significant up-regulation in the gene encoding ICAD, which is associated with DNA fragmentation during apoptosis (Figure 2A). In addition to apoptosis, pathway analysis further revealed that inflammation and inflammatory forms of cell death, such as necroptosis, were significantly down-regulated following Magnéli phase exposure (Figures 2A,B). While we realize the general limitations of utilizing gene expression data to define cell death pathways, the data indicates significant changes in genes associated with apoptosis pathways and is consistent with the electron microscopy and inflammation data (Figure 1). Beyond cell death, very few additional pathways were found to be significantly impacted following Magnéli phase phagocytosis at the time points and conditions evaluated here, with a notable exception being a general downregulation of genes associated with NF-κB signaling (Figure 2A). This finding is consistent with the lack of inflammation and cytokine responses observed following exposure (Figures 1I,J). To better evaluate our findings compared to those previously reported for TiO2 exposure, we compared our findings to gene expression patterns associated for TiO2 in the IPA database (Supplemental Figure S1C). In general, many of the same pathways are altered for both Magnéli phases and TiO2, including the down-regulation of Tnf and the up-regulation of Ifng (Supplemental Figure S1C). However, the major difference between Magnéli phases and TiO2 is associated with Caspase-3/-7 regulation. Caspase-7 is down regulated in both cases; but, Caspase-3 is significantly up-regulated following Magnéli phase exposure and is central in the apoptosis pathways identified (Figure 2A; Supplemental Figure S1C). This is in contrast to the down-regulation of Caspase-3 and cell death mechanisms previously described following TiO2 exposure (17). Together, these data suggest that macrophage cell death induced by Magnéli phase exposure occurs though the mitochondrial pathway of apoptosis and is associated with the activation of pro-apoptotic BCL-2 family members, including a potential p53-BAX-APAF1 axis.

      Magnéli Phase Exposure Results in Significant Mitochondrial Dysfunction

      Previous studies have revealed that nanoparticles can significantly impair mitochondria function and the convergence of signaling pathways on the mitochondria identified in the gene expression profiling and pathway analysis together warrant higher resolution functional assessments. We first utilized an Agilent Seahorse XF cell mitochondrial stress test to evaluate the oxygen consumption rate (OCR) in BMDMs treated for 24 h with either mock, 100 ppm, or 1,000 ppm Magnéli phases (Figure 3A). At both doses, the basal respiration and the maximal respiratory capacity (MRC) were significantly increased in BMDMs treated with Magnéli phases compared to mock (Figure 3A). Cells treated with 1,000 ppm of Magnéli phases did not show the typical decline in OCR after the addition of oligomycin (Oligo), a known ATP synthase inhibitor, indicating an increase in proton leak in the mitochondria (Figure 3A). Proton leak is typically associated with mitochondrial damage and stress (18, 19). Likewise, these cells only showed a minimal change in MRC following the addition of the proton ionophore FCCP (Figure 3A), which likely reflects the significant increase in cells undergoing cell death. The BMDMs treated with 100 ppm Magnéli phases showed the most dramatic shifts, with the highest level of basal respiration and the largest MRC increase, following FCCP treatment (Figure 3A). The impressive increase in MRC suggests that Magnéli phase exposure resulted in a substantially higher reserve capacity in BMDMs, which is also likely associated with the increased energy demands of the cells following increased damage, stress, and cell death. Complementing the OCR data, we also evaluated the extracellular acidification rate (ECAR), which is an assessment of glycolysis. While we did not observe significant differences in ECAR following exposure to Magnéli phases at 1,000 ppm, we did observe a significant increase following exposure to 100 ppm (Figure 3B). These data indicate a shift in the substrate utilization, likely due to increased cellular metabolism in cells treated with 100 ppm Magnéli phases. Together, these data indicate that Magnéli phase exposure results in altered mitochondrial bioenergetics that can significantly contribute to damage associated processes, such as reactive oxygen species (ROS) generation and the activation of cell death cascades.

      Alterations in cellular energetics and mitochondrial membrane potential in macrophages treated with Ti6O11. (A) Oxygen consumption rate (OCR) and (B) extracellular acidification rate (ECAR) were measured using Agilent Seahorse XF96 Analyzer in mock and Ti6O11 treated BMDMs at either 100 ppm or 1,000 ppm. (C) Mitochondrial membrane potential was evaluated using tetramethylrhodamine (TMRM) in macrophages treated with vehicle (blue) or 100 ppm Magnéli phases (green), with H2O2 used to induce maximum loss of membrane potential in each group. (D) Cellular oxidative stress was evaluated using dihydrodichlorofluorescein (DCF) fluorescence. All data are expressed as mean ± SEM (n = 3/group). AU, Arbitrary Units. *p < 0.05, #p < 0.05, **p < 0.01.

      Because proton leak is also a sign of mitochondrial damage and stress, we next investigated mitochondrial membrane potential (MMP), as decreases in MMP are associated with mitochondrial membrane damage and changes in mitochondrial function. The cell-permeant dye, tetramethylrhodamine (TMRM) accumulates in active mitochondria with healthy MMPs. A decrease in TMRM fluorescence indicates a loss of mitochondrial membrane potential. BMDMs treated with 100 ppm Magnéli phases had significantly lower TMRM fluorescence (86.26 ± 13.06 AU) compared to mock (613.60 ± 51.89 AU) (Figure 3C), indicating treatment significantly altered mitochondrial membrane potential. The subsequent addition of H2O2 to each sample was used as a positive control and resulted in significant decreased TMRM fluorescence compared to basal fluorescence in the mock treated cells, with a much smaller decrease in the cells treated with Magnéli phases (Figure 3C). Together, these data suggest that Magnéli phase exposure results in a significant decrease in mitochondrial membrane potential, and confirms severe mitochondrial membrane damage.

      The production and modulation of ROS is a common biological response in macrophages to nanoparticle exposure, is associated with cytotoxicity, and is commonly involved in mitochondrial stress (2022). Increased ROS production is a significant driver of cell death and stress following TiO2 exposure (20). A previous study assessed ROS levels in A549 alveolar epithelial cell lines following TiO2 and undefined Magnéli phase exposure (14). Under the conditions evaluated, the TiOx nanoparticles increased intracellular Ca2+ that was associated with low levels of cytotoxicity, but no significant changes on ROS production were observed (14). Due to the inherent differences in nanoparticle responses between different cell types and the potent cytotoxic effects of dysregulated ROS production, we next evaluated this response following exposure to Ti6O11 nanoparticles in macrophages. For this we utilized another cell-permeant dye, 2′,7′-dichlorofluorescin diacetate (DCFDA), which is oxidized by ROS into the fluorescent DCF. BMDMs challenged with 100 ppm Magnéli phases for 24 h had a significantly higher level of ROS production (1,105 ± 94.17) compared to mock exposed cells (76.39 ± 7.22) (Figure 3D). Taken together, these data suggest that macrophages treated with Magnéli phases show increased ROS production and loss of mitochondrial membrane potential after 24 h. All of which are indications that these particles induce mitochondrial stress and are consistent with the induction of mitochondria-mediated apoptosis following Magnéli phase exposure.

      Magnéli Phases Are Sequestered in Lung-Associated Macrophages and Concentrate in the Tissue After a Single Exposure

      In the only in vivo toxicity study conducted to date, Magnéli phases were shown to be bioactive and toxic to zebrafish at 100 ppm (1). These studies provided the first indication that these incidental nanoparticles could be biologically active and potentially toxic to a variety of animal species, including humans. However, the impact of Magnéli phases in mammals and the mechanism of toxicity were not explored (1). Thus, we next sought to extend these findings using mice and a respiratory exposure model. Mice were exposed to 100 ppm Ti6O11 Magnéli phases in 1 × PBS via intra-tracheal administration using methods previously described to evaluate nanoparticle toxicity and trafficking in the respiratory tract (13, 23). To quantify the concentration of Magnéli phases in the lungs over time, we utilized inductively coupled plasma mass spectrometry analysis to detect titanium levels (Figure 4A; Supplemental Figure S2A). Immediately following i.t. administration, tissues including the lungs, liver, spleen, and kidneys were removed, digested, and processed for ICP-MS. These data represent the maximum deposition control (Figure 4). As expected, no titanium was detected in tissues outside of the lungs or in any tissues that were not exposed to the Magnéli phases. However, in exposed lungs, we observed 14.46 ± 4.32 μg/g of lung tissue immediately after deposition (Figure 4A). Within 24 h post-exposure, titanium levels detected in the lungs dropped to 3.40 ± 1.19 μg/g of lung tissue (Figure 4A). However, surprisingly, after this single exposure to the Ti6O11 Magnéli phases, the levels of titanium in the lungs were maintained over the course of a week with levels holding steady at 5.30 ± 1.23 μg/g (Day 3) and 5.91 ± 1.30 μg/g (Day 7) (Figure 4A). Titanium levels were below the level of detection in all of the other tissues evaluated at the time points assessed (Supplemental Figure S2B). Together, these data emphasize the potential consequences of even a single exposure to Magnéli phases, which are retained in the lungs at significant levels for an extended duration following inhalation.

      Magnéli phases concentrate in pulmonary macrophages following a single exposure and are retained in the lung. (A) Following a single airway exposure to 100 ppm, titanium associated with the Magnéli phases were retained in the lungs and detected using ICP-MS for over 7 days post-expsure. Control tissues were collected from animals exposed to 100 ppm Ti6O11 and immediately euthanized. (B–D) Representative H&E stained tissue sections used for histopathology evaluation. (B) Vehicle control showing no nanoparticles or pathology. (C) Significant airway inflammation was observed in LPS treated animals, but no evidence of inflammation was found in the mice treated with Magnéli phases at any timepoint evaluated. (D) Larges areas of macrophages containing Magnéli phases were found in all treated animals (yellow box). (E–G) The majority of pulmonary macrophages contain Magnéli phases, with no other cell types appearing to contain or associate with the nanoparticles determined using dark field microscopy. Darkfield images were taken of (E) Ti6O11 nanoparticles alone or the lungs of mice 7 days post-exposure to either (F) PBS or (G) 100 ppm Magnéli phases. (G) Particles in the tissues were identified by bright punctate dots, almost exclusively localized in macrophages. Representative H&E stained lung sections (Scale bar = 50 μm). All data are expressed as mean ± SEM (n = 7/group). **p < 0.01.

      We next sought to identify where in the lungs the Magnéli phases were being sequestered. We initially utilized histopathology assessments of H&E stained lung sections to evaluate both nanoparticle deposition and pathology. Lungs were fixed by gravity inflation with formalin and paraffin embedded for staining as previously described (13, 23). Lung pathology was evaluated by a board-certified veterinary pathologist (S.C.O.). Control mice received either vehicle or were challenged with a 1 mg/ml dose of LPS (Figures 4B,C). As anticipated, lung histopathology was normal in the vehicle treated animals (Figure 4B). Conversely, we observed significant airway inflammation in the LPS treated animals, characterized by predominately macrophage and neutrophil infiltration concentrated around the airways and vasculature (Figure 4C). In the experimental animal groups, we did not observe significant evidence of inflammation in the mice exposed to Magnéli phases at 1, 3, or 7 days post-exposure (Figure 4D). However, Magnéli phases were readily observed as punctate foci within cells in the lungs under light microscopy conditions (Figure 4D). Further assessments at higher magnification revealed that the vast majority of Magnéli phases were sequestered in macrophages, with no nanoparticles observed under light microscopy conditions associated with or within any other cell type in the lung or observed in the vehicle control specimens (Figures 4B,D). To confirm the association of Magnéli phases with the macrophages in the lungs, we conducted enhanced darkfield microscopy (Cytoviva) (Figures 4E–G). A reference image of Magnéli phases was generated and used to identify the nanoparticles in vivo (Figure 4E). As anticipated, we did not detect any Magnéli phases in the vehicle control specimens (PBS) (Figure 4F). However, consistent with our pathology assessments, Magnéli phases were confirmed in macrophages throughout the lungs and were readily identified by their distinctive punctate staining (Figure 4G; green arrow). Together, these data show that Magnéli phases are effectively phagocytosed by lung-associated macrophages; however, the nanoparticles are not effectively cleared from the tissue, remaining at steady levels over the course of at least 7 days following a single exposure.

      Based on our ex vivo findings that Magnéli phases are cytotoxic, we next evaluated the in vivo impact of nanoparticle sequestration on macrophages and lung tissue. Bronchoalveolar lavage fluid (BALF) was collected following euthanasia and cellularity was evaluated following differential staining. In both the PBS and Magnéli phase treated specimens, macrophages were the dominate population of cells recovered (Figure 5A). This is in contrast to the highly inflammatory conditions observed following LPS treatment (data not shown) or following exposure to nanoparticles with immunostimulatory properties, where neutrophils are the dominate cell type observed (2427). Magnéli phases were readily identified in macrophages collected from the lungs of exposed mice at all time points evaluated, with 71 ± 10.16% of the macrophages containing at least 1 readily observable nanoparticle under 40 × light microscopy (Figure 5A). Consistent with our in vitro cytotoxicity data, many of the macrophages recovered from the airways demonstrated morphological features consistent with cellular stress and death (Figure 5A). The cell populations recovered in the BALF were further quantified using a hemocytometer (Figure 5B). Macrophages, neutrophils, and lymphocytes were recovered in the BALF from mice exposed to either PBS or Magnéli phases (Figure 5B). The quantity and populations of cells recovered in BALF are routinely utilized as surrogates to characterize lung inflammation (28). Compared to findings commonly reported following exposure to immunostimulatory agents, such as LPS, and consistent with the pathology evaluations shown in Figures 4B–D, the number and populations of cells recovered following exposure to Magnéli phases indicate that these nanoparticles do not induce significant lung inflammation (Figure 5B). Indeed, assessments of inflammatory mediator levels, such as IL-6, in the BALF following either Magnéli phase exposure or LPS exposure, further confirm these findings (Supplemental Figure S3A). Consistent with the increased cytotoxicity, we observed a significant decrease in the number of macrophages recovered in the BALF from mice exposed to Magnéli phases (Figure 5B). No other cell populations appeared to be impacted (Figure 5B). Together, these data indicate that airway exposure to a single dose of Magnéli phases at a range previously found in environmental samples following incidental nanoparticle release and cytotoxic to zebrafish embryos indeed facilitates macrophage cell death in the mammalian lung. Consistent with the apoptosis findings (shown in Figures 1, 2), neither in vivo macrophage cell death nor exposure to the Magnéli phases result in significant inflammation.

      Magnéli phases concentrate in pulmonary macrophages, resulting in significant dysfunction. (A) Representative images of BALF cytology from PBS and Ti6O11 (100 ppm) treated mice. (B) Differential cell counts in the BALF after treatment with PBS, Ti6O11, or LPS. (C) Gene expression profiling and analysis using Ingenuity Pathway Analysis software identified gene transcription patterns associated with increases in apoptosis, formation of reactive oxygen species, and wound healing/fibrosis. Data are expressed as mean ± SEM (n = 3 mice per PBS treated group, n = 7 mice per Ti6O11 treated group, and n = 4 per LPS treated group). *p < 0.05.

      To better define the mechanisms underlying the in vivo response to Magnéli phase exposure, lung tissue was collected 7 days post-exposure. Total RNA was isolated and processed for gene expression profiling using the SuperArray platform and IPA, as previously described (79). The gene expression profiling analysis revealed significant changes in several pathways associated with cell death and lung function. The three top pathways predicted by IPA analysis to be significantly activated 7 days post-exposure to the Magnéli phases were (1) apoptosis—specifically apoptosis of monocytes; (2) formation of ROS; and (3) fibrosis. In all three cases the combination of genes up-regulated and down-regulated were predicted to have an activating effect on the pathways identified (Figure 5C). Also, of note, no pathways were predicted to be significantly inhibited based on the gene expression patterns evaluated and no pathways associated with inflammation were significantly modulated following Magnéli phase exposure. The predicted activation of apoptosis of monocytes and formation of ROS in vivo are highly consistent with the findings from our in vitro studies and further emphasize the detrimental impact that the Magnéli phases are having on the macrophage populations in the lungs following exposure. The fibrosis gene expression signature is intriguing. The 7-day timeframe of the exposure model is far too soon for lung fibrosis to present in the mice and no pathologic features associated with fibrosis were detected in our specimens by histopathology, collagen, or fibronectin assessments (Figure 4). It is possible that we are identifying gene signatures associated with fibrosis at a very early time in the lungs, prior to pathological disease onset. However, this gene expression pattern would also be consistent with increased regeneration and repair processes in the lungs, perhaps associated with the increased cell death.

      Repeated Exposure to Magnéli Phases Results in Significant Nanoparticle Accumulation in the Lungs

      In most cases of incidental nanoparticle release into the environment, biological exposure to Magnéli phases are expected to occur repeatedly and for an extended duration. Thus, we next sought to determine the impact of chronic airway exposure to the Ti6O11 nanoparticles over a 30-day time period. Under these conditions, we observed a significant increase in titanium concentrations in the lungs, which accumulated to levels that were ~6-fold higher (35.82 μg/g of lung) than levels observed after a single exposure (5.92 μg/g of lung) (Figure 6A). Similar to the findings from the single exposures, pathology assessments of the lungs again revealed that the majority of the Magnéli phases detectable under bright-field microscopy were sequestered in lung-associated macrophages (Figure 6B). Despite the increased concentrations of nanoparticles in the lungs, we did not observe any pathologic features associated with increased inflammation or differences in pro-inflammatory mediators, such as IL-6, IL-1β, or TNF levels, in the lungs (Figure 6B; Supplemental Figures S3A–F). Pathology assessments revealed no evidence of interactions with other cell types beyond the macrophages under bright-field microscopy. To better define nanoparticle-cell interactions, we evaluated specimens using enhanced dark field microscopy (Cytoviva) (Figure 6C). Consistent with our bright field data and similar to what we observed following the single exposure, the Magnéli phases are overwhelmingly concentrated in macrophages in the lung (Figure 6C). It should be noted that we did identify punctate staining associated with alveolar epithelial cells, indicating that under chronic exposure conditions, there are some interactions between other cell types in the lungs and the Ti6O11 nanoparticles (Figure 6C). However, compared to the macrophages, Magnéli phase interactions with these cells were more sporadic and involved far fewer nanoparticles per cell (Figure 6C). Using cells recovered following BALF and evaluated following differential staining assessments, the number of macrophages containing Magnéli phases per 100 cells was determined and used to generate a phagocytic index (29). Following the single exposure to Magnéli phases, we calculated the phagocytic index to be 71.14 ± 3.84 on day 7, with the majority of macrophages recovered containing >1 nanoparticle per cell (Figure 6D). However, following the multiple exposures, the phagocytic index significantly increased to 91.00 ± 2.17 on day 30, again with almost all macrophages recovered containing >1 nanoparticle per cell (Figure 6D). Together, these data are consistent with the Magnéli phases being sequestered in the macrophages. Following cell death, the nanoparticles likely continue to concentrate in the lungs as cells undergoing apoptosis and associated debris (including the Magnéli phases) are further phagocytosed by additional macrophages, creating a positive feedback loop.

      Repeated exposures concentrates Magnéli phases in pulmonary macrophages. (A) Using ICP-MS, we quantified the amount of titanium in the lungs following multiple airway exposures to 100 ppm Ti6O11 over a 30 day period. (B) Histopathology using H&E stained tissue sections reveled significant airway inflammation LPS treated animals, but no evidence of inflammation in the mice treated with Magnéli phases. Larges areas of macrophages containing Magnéli phases were found in all treated animals (Scale bar = 20 μm). (C) Darkfield images were taken of mouse lungs post-repeated exposures, on day 30. Particles in the tissues were identified by bright punctate dots. Almost all particles were concentrated in macrophages (green arrows). However, small numbers of Ti6O11 nanoparticles were found associated with alveolar epithelial cells (green arrows). Representative H&E stained lung sections (Scale bar = 50 μm). (D) Phagocytic Index (number of BALF macrophages containing Magnéli phases per 100 cells). (E–G) APAF1 immunohistochemistry staining from (E) saline, (F) LPS, and (G) Magnéli Phase exposed lungs. (G) Following Magnéli Phase exposure, only endothelial cells and macrophages containing particles were broadly positive for APAF1. (E) In saline exposed mice, only endothelial cells were positive for APAF1; (F) whereas, LPS exposure resulted in a range of cell types staining positive for APAF1. Macrophages in each image are identified by red arrows. All data are expressed as mean ± SEM (n = 7/group). **p < 0.01.

      Our gene expression profiling data indicated that this sequestration resulted in apoptosis. To further evaluate this, we conducted immunohistochemistry targeting APAF1, which is an important mediator in mitochondria mediated apoptosis and was significantly up-regulated at the gene expression level (Figure 2A). Immunohistochemistry revealed that Ti6O11 treatment resulted in significant increases in APAF1+ cells in the lungs. Consistent with previous data, endothelial cells have high levels of APAF1+ staining, which was observed in the Saline treated animals (Figure 6E). However, we did not observe any significant staining in other cell types, including macrophages (Figure 6E). In LPS treated lungs, we observed robust staining in multiple cell types in the lungs, including endothelial cells, epithelial cells, macrophages, and neutrophils (Figure 6G). However, in the Magnéli phase treated animals, we did not observe the same level of widespread staining (Figure 6F). Rather, the APAF+ cells were predominately limited to the endothelial cells and macrophages containing Magnéli phases (Figure 6F). Also, of note, only the macrophages containing Magnélli phases were APAF+. Together, these data are consistent with our other findings localizing the in vivo effects of the Magnélli phases to the macrophage compartment in the lungs.

      Magnéli Phase Exposure Attenuates Lung Function and Impacts Airway Pathophysiology

      Based on the data above, we predicted that exposure to Magnéli phases would negatively impact mammalian lung function. To directly test this hypothesis, we determined the impact of multiple exposures of Ti6O11 nanoparticles over a 30 day period of time on basal airflow and airway mechanics utilizing a computer-controlled small-animal ventilator with highly sensitive pressure transducers (Flexivent) to quantify airway opening pressures, volume, and airflow (10, 11). In addition to basal measurements, we also determined whether exposure to these Magnéli phases modulates lung pathophysiology and airway mechanics following exposure to the bronchoconstricting agent methacholine. To evaluate airway mechanics, a snapshot perturbation maneuver that consisted of a three-cycle sinusoidal wave of inspiration and expiration was used to measure total respiratory system resistance (R), dynamic compliance (C), and elastance (E). A Quickprime-3 perturbation, which produced a broadband frequency (0.5–19.75 Hz) over 3 s, was used to measure Newtonian resistance (Rn), which is a measure of central airway resistance. For the methacholine studies, after three baseline measurements of lung function, mice were challenged with aerosolized saline or increasing doses of 1.5, 3, 6, 12, and 24 mg/ml of methacholine, generated by an ultrasonic nebulizer for 10 s without altering the ventilation pattern, followed by assessments every 30 s (30, 31). Between each dose of methacholine, a deep inflation perturbation was initiated to reset lung hysteresis.

      Consistent with our hypothesis, exposure to Magnéli phases resulted in significant increases in basal resistance and Newtonian resistance (Figures 7A,B). These data indicate that the accumulation of Ti6O11 nanoparticles significantly attenuates baseline lung function. The increased basal resistance and Newtonian resistance suggests that the airway diameter and lung volumes are reduced under baseline conditions. We also observed an increase in basal elastance and the converse decrease in compliance (Figures 7C,D). These findings reflect defects in alveolar and elastic tissue in the lungs and indicate that exposure to Magnéli phases results in increased lung stiffness and reduced lung expansion abilities. Beyond baseline assessments, we also sought to evaluate lung function following methacholine challenge, which is commonly used to induce smooth muscle constriction in the lungs. Despite the increased baseline resistance and Newtonian resistance, the lungs of mice exposed to Magnéli phases showed minimal increases in airway constriction following methacholine challenge (Figures 7E,F). This is in contrast to the unexposed animals, which demonstrated the anticipated dose dependent increases (Figures 7E,F). These data indicate that the lungs of these mice lack the capacity to properly function and are likely at maximum for resistance and Newtonian resistance. Methacholine challenge also results in increased elastance and decreased compliance in untreated mice (Figures 7G,H). However, similar to the resistance measures, Magnéli phase exposure resulted in minimal changes in elastance and dynamic compliance, indicating a stiffening of the lungs (Figures 7G,H). Together, these data are consistent with Magnéli phase exposure resulting in significant detrimental changes in lung mechanics and function.

      Chronic exposure to Magnéli phases significantly attenuates lung function. (A–D) Baseline pulmonary function parameters including (A) total respiratory system resistance (R), (B) Newtonian resistance (Rn), (C) elastance (E), and (D) dynamic compliance (C) were measured by FlexiVent in control and Ti6O11 exposed mice. (E–H) Assessment of AHR to methacholine (MCh) in control and Ti6O11 exposed mice. (E) Total respiratory system resistance (R), (F) Newtonian resistance (Rn), (G) elastance (E), and (H) dynamic compliance (C) dose response curves following challenge with increasing concentrations of MCh in aerosolized saline (0, 1.5, 3, 6, 12, and 24 mg/ml) for control and Ti6O11 exposed mice were examined on FlexiVent. All data are expressed as mean ± SEM (n = 7/group). *p < 0.05, ****p < 0.0001, compared to the PBS treated control group.

      Discussion

      In initial toxicity studies, Magnéli phases (Ti6O11) were found to be toxic to dechorionated zebrafish embryos at 100 ppm, without concurrent solar radiation (1). This is in contrast to TiO2 nanoparticles, which are photocatalytically active and toxic at the same concentrations in conditions where solar radiation is present (1). The mechanism of cytotoxicity associated with TiO2 has previously been associated with the increase in ROS generation (32, 33). Prior to data presented in the current study, the physiological impact in mammalian cells and cytotoxicity mechanisms associated with Magnéli phases were unclear. To date, only one additional study has explored the cytotoxicity of these nanoparticles in mammalian cells (14). The cytotoxicity of three different Magnéli phases and two different TiO2 formulations were evaluated utilizing A549 cells, which are a human lung alveolar epithelial cell line (14). In these studies, rather than using specific formulations of Magnéli phases, the nanoparticles were evaluated based on size, ranging from 192 to 795 nm in average diameter (14). In the A549 cells, the TiOx Magnéli phases demonstrated low levels of cytotoxicity in some of the size ranges evaluated and increased intracellular Ca2+, but no significant changes in ROS levels were observed (14). Thus, the authors concluded that in A549 cells Magnéli phases had either similar or reduced cytotoxicity compared to TiO2 nanoparticles with cell death associated with different mechanisms related to Ca2+ flux (14).

      The prior studies using zebrafish embryos also provides critical insight into general in vivo toxicity and immediate relevance to aquatic organism exposure (1). However, the impact of Magnéli phases on other terrestrial biological systems and the extension of these findings to human physiology is somewhat limited. Similarly, the in vitro studies using the A549 model provides important insight related to the interaction of Magnéli phases with a human alveolar epithelial cell line and a comparison with TiO2 nanoparticles (14). However, while the A549 cells are certainly relevant, when nanoparticles are internalized through the airway, several other cell types beyond the alveolar epithelial cells are arguably more critical to the host response (13). Based on data presented here, interactions between Magnéli phases and alveolar epithelial cells appear quite limited in vivo. Our data indicates that the vast majority of nanoparticles are sequestered in lung-associated macrophages. These cells are the primary phagocytic cell in the airway and these findings are similar to other studies evaluating nanoparticle deposition in the lungs (13, 23, 24). Likewise, our data indicates that the macrophages are not able to effectively breakdown and clear the concentrated Magnéli phases from the lungs. Even more concerning, this sequestration occurs even after a single challenge and suggests that even short term, low dose exposure to Magnéli phases may pose a potential health risk.

      Once phagocytosed by macrophages, our data reveals that Ti6O11 Magnéli phases are cytotoxic and induce apoptosis, which is considered a non-inflammatory form of cell death. This finding is consistent with the lack of inflammatory signaling in the ex vivo primary macrophages and in vivo findings. However, the lack of inflammation, both in vitro and in vivo, following nanoparticle exposure was quite unexpected. As discussed above, we observed active suppression of gene expression pathways associated with necroptosis (Figure 2B), which is a proinflammatory form of cell death. While this was also unexpected, similar findings have been reported for other metal-based nanoparticles. Indeed, the ability to attenuate inflammation while inducing cell death is a favorable characteristic for engineered nanomaterials generated for biomedical and therapeutic applications. For example, zinc oxide nanoparticles induce high levels of apoptosis associated with increased ROS production and significant depletion of glutathione in human pulmonary adenocarcinoma cell lines and have well defined anti-inflammatory activities (34, 35). Thus, it is tempting to speculate a similar model applies to the Magnéli phases evaluated here. Similarly, silver nanoparticles have shown effects in neutrophils that are consistent to those described here for the Magnéli phases in macrophages (36). Silver nanoparticle exposure triggers activation and maturation in specific subsets of neutrophils, followed by increased levels of IL-8 (36). However, classical pro-inflammatory pathway activation, overall inflammation and necrotic cell death are significantly attenuated following exposure (36). The implication of these data indicate that silver nanoparticles stimulate neutrophil up-take as the cell attempts to minimize inflammation and stimulate nanoparticle clearance (36). It is possible that a similar mechanism is also occurring in macrophages, where Magnéli phase nanoparticle phagocytosis is occurring, but the cells are not able to effectively degrade or clear the nanoparticles.

      Mitochondria are critical in the activation of apoptosis in mammalian cells and the Bcl-2 family is central in the regulation of these processes (37). Our data reveals that this pathway is significantly dysregulated following Magnéli phase exposure and provides insight into the mechanism underlying the increased apoptosis observed following nanoparticle exposure. These data are consistent with findings associated with a diverse range of metal-based nanoparticles. For example, gadolinium oxide (Gd2O3) nanoparticles have been shown to induce apoptosis through dysregulation of BCL2/BAX signaling and ROS production in response to DNA damage in human neuronal cells (38). Likewise, similar to the gene expression pattern observed following Magnéli phase exposure, silver nanoparticle exposure also increases caspase 3, Bax, and P53 expression, increased ROS signaling, and apoptosis in the MCF-7 human adenocarcinoma cell line (39). But, as discussed above and also similar to Magnéli phases, silver nanoparticle exposure results in an attenuated immune system response (36). Interestingly, lower doses of Magneli phases (100 ppm) increased maximal respiration, consistent with the concept of mitochondrial hormesis, whereby ROS upregulate mitochondria as a compensatory mechanism to fight cellular stressors (40). This is consistent with our observation that higher doses (1,000 ppm) were associated with functional decrements, including uncoupled respiration and a loss in respiratory reserve in macrophages. Likewise, the dramatic shifts observed in the OCR and ECAR at 100 ppm further emphasize that the mitochondria are dysfunctional and indicates that Magnéli phase exposure results in significantly altered mitochondrial bioenergetics. This likely contributes to damage associated processes, such as ROS generation and the activation of cell death cascades. Together, these data support a model whereby Magnéli phases directly impair mitochondrial function following exposure and induce cell death in macrophages through the mitochondrial-mediated intrinsic apoptosis pathway.

      Despite the lack of inflammation in the lungs following Magnéli phase exposure, it is clear that the sequestration and concentration of these nanoparticles and resultant macrophage cell death significantly affects lung function. We observed significant changes in resistance, compliance, and elastance at baseline and following methacholine stimulation (Figure 7). At baseline, the increase in resistance following repeated Magnéli phase exposure reflects a narrowing of the central and conducting airways, including the trachea, bronchi, bronchioles, and terminal bronchioles. Similarly, the reduction in compliance and corresponding increase in elastance at baseline reflect the need for greater pressures to increase lung volume in animals exposed to Magnéli phases. Following stimulation with methacholine, we observed predicted changes in R, Rn, E, and C in animals unexposed to Magnéli phases, reflecting normal lung function. However, in mice exposed to the nanoparticles, we observed relatively minimal changes in lung mechanics from baseline (Figure 7). This suggests that the lungs have already reached their near maximum potential for the parameters assessed in this study. Together, these consequential alterations in airway mechanics indicate that the lungs of mice exposed to Ti6O11 Magnéli phases are significantly stiffer and more constricted compared to unexposed animals. This influences both baseline lung function and normal airway contractility following stimulation. These findings are consistent with restrictive lung disease and are similar to findings reported for silica associated models of fibrosis, which are characterized by low levels of inflammation and the production of pro-fibrotic growth factors that ultimately result in fibrotic lesions due to diminished lung clearance of the silica particles (4144). These fibrotic lesions are directly associated with diminished macrophage function. Macrophages that concentrate the silica particles are highly effective in directly stimulating fibrogenesis associated growth factors in the lungs that can impair lung function. These findings are consistent with those presented here following exposure to Magnélli phases and we speculate that a similar mechanism underlies the observed decrease in lung function.

      Due to the widespread environmental distribution of incidental Magnéli phases in the air and water associated with global coal combustion, it is critical that we have a better understanding of their biological impacts. Improved understanding of the dynamics associated with nanoparticle physicochemistry, cellular, and organismal responses following exposure, environmental transport, and global distribution patterns are all critical for providing a basis for establishing future frameworks to determine potential risks to human health and inform environmental policy decisions. Our findings here are highly cautionary, as a single exposure to Magnéli phases can result in potentially harmful long-term health effects. Unfortunately, there is currently no practical way to limit the formation of Magnéli phases and other nanoparticle formation during the coal burning process. However, in countries with strong environmental regulations, most of these nanoparticles can be captured by particle traps prior to final emission of exhaust gas. While we have provided the first evidence of the detrimental impacts of Magnéli phase exposure in mammals, it is clear that further toxicology studies are necessary and future assessments of the impact of these nanoparticles on human health is warranted.

      Data Availability Statement

      The datasets generated for this study are available on request to the corresponding author.

      Ethics Statement

      The animal study was reviewed and approved by Virginia Tech IACUC East Carolina University IACUC.

      Author Contributions

      DM and IA designed, performed the in vitro and in vivo experiments, analyzed and interpreted the data, prepared the figures, and wrote the paper. VR-S, HM, WL, VM, YY, BC, JB, and MH were responsible for designing and executing specific studies associated with Magnéli phase synthesis and characterization, localization, and macrophage studies. MC-T and JA were responsible for the ICP-MS analysis. TC conducted the TEM analysis. JP, GD, and DB were responsible for the mitochondria studies. SR and KG conducted the lung function studies. Pathology assessments were conducted by a board certified veterinary pathologist, SC-O. All authors have read and approved the manuscript.

      Conflict of Interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      We would like to thank Dr. Kristin Eden, Dr. Daniel Rothschild, Melissa Makris, and Dr. Mitsu Murayama for insightful comments, critical reagents, and instrumentation assistance associated with this publication. We would also like to acknowledge the staff and personnel related to our Teaching and Research Animal Care Support Service (TRACSS). Center support was provided by the Virginia Tech Institute for Critical Technology and Applied Science, which maintains critical support for the Virginia Tech Center for Engineered Health, the Virginia Tech Center for Sustainable Nanotechnology, and the National Center for Earth and Environmental Nanotechnology Infrastructure.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fimmu.2019.02714/full#supplementary-material

      Apoptosis signaling was significantly increased in macrophages following exposure to Magnéli phases. (A) IPA revealed that cell death and apoptosis signaling pathways were significantly up-regulated following Ti6O11 exposure. The top molecular and cellular functions predicted to be impacted by exposure were associated with cell death, cellular function and maintenance, DNA repair, and cellular organization. (B) The 10 genes most impacted (the top 5 up-regulated and down-regulated) by Magnéli phase exposure in macrophages have broad and diverse biological functions. (C) Macrophage responses to Magnéli phases produce a gene expression pattern similar to that observed following exposure and activation following titanium dioxide. However, Magnéli phase exposure is significantly less inflammatory (lower levels of IFNG, IGF1, and FASLG), while inducing higher levels of apoptosis and increased mitochondria dysfunction. The activation of gene networks associated with Caspase 3 is a significant difference observed here for Magnéli phases, which is down-regulated following titanium dioxide exposure.

      Titanium was not found in tissues and organs outside of the lungs. (A) Table showing ICP-MS instrument parameters used for the detection of titanium in tissue. (B) Blood, liver, and kidneys were evaluated for titanium levels over the course of 7 days post-exposure to a single dose of Magnéli phases at 100 ppm. All levels in all tissues were below the level of detection, with the exception of a single liver sample, shown here 3 days post-exposure.

      Magnéli phases are non-inflammatory. (A,B) Pro-inflammatory cytokine levels were evaluated by gene expression using rtPCR and protein evaluations using ELISA in all in vitro and in vivo studies. However, no significant increases in the majority of potent pro-inflammatory mediators were detected. (A–F) IL-6, IL-1β, and TNF levels were evaluated in bronchoalveolar lavage fluid (BALF) by ELISA following either (A–C) a single exposure of 100 ppm; (D–F) chronic exposures of 100 ppm Magnéli phases; or (A–F) airway exposure to 1 mg/ml of lipopolysaccharide (LPS). LPS induced a robust IL-6, IL-1β, and TNF response; however, no exposures to the Magnéli phases significantly induced these proinflammatory mediators. *p < 0.05.

      References Yang Y Chen B Hower J Schindler M Winkler C Brandt J . Discovery and ramifications of incidental Magneli phase generation and release from industrial coal-burning. Nat Commun. (2017) 8:194. 10.1038/s41467-017-00276-228790379 Rohde RA Muller RA. Air pollution in China: mapping of concentrations and sources. PLoS ONE. (2015) 10:e0135749. 10.1371/journal.pone.013574926291610 Oberdorster G Oberdorster E Oberdorster J. Nanotoxicology: an emerging discipline evolving from studies of ultrafine particles. Environ Health Perspect. (2005) 113:82339. 10.1289/ehp.733916002369 Hochella MF Jr Mogk DW Ranville J Allen IC Luther GW Marr LC . Natural, incidental, and engineered nanomaterials and their impacts on the Earth system. Science. (2019) 363:eaau8299. 10.1126/science.aau829930923195 Yang Y Colman BP Bernhardt ES Hochella MF. Importance of a nanoscience approach in the understanding of major aqueous contamination scenarios: case study from a recent coal ash spill. Environ Sci Technol. (2015) 49:337582. 10.1021/es505662q25688977 Davis BK. Derivation of macrophages from mouse bone marrow. Methods Mol Biol. (2019) 1960:4155. 10.1007/978-1-4939-9167-9_330798519 Allen IC Moore CB Schneider M Lei Y Davis BK Scull MA . NLRX1 protein attenuates inflammatory responses to infection by interfering with the RIG-I-MAVS and TRAF6-NF-kappaB signaling pathways. Immunity. (2011) 34:85465. 10.1016/j.immuni.2011.03.02621703540 Coutermarsh-Ott S Simmons A Capria V LeRoith T Wilson JE Heid B . NLRX1 suppresses tumorigenesis and attenuates histiocytic sarcoma through the negative regulation of NF-kappaB signaling. Oncotarget. (2016) 7:33096110. 10.18632/oncotarget.886127105514 Eden K Rothschild DE McDaniel DK Heid B Allen IC. Noncanonical NF-kappaB signaling and the essential kinase NIK modulate crucial features associated with eosinophilic esophagitis pathogenesis. Dis Model Mech. (2017) 10:151727. 10.1242/dmm.03076729259025 Allen IC Pace AJ Jania LA Ledford JG Latour AM Snouwaert JN . Expression and function of NPSR1/GPRA in the lung before and after induction of asthma-like disease. Am J Physiol Lung Cell Mol Physiol. (2006) 291:L10051017. 10.1152/ajplung.00174.200616829631 Cyphert JM Kovarova M Allen IC Hartney JM Murphy DL Wess J . Cooperation between mast cells and neurons is essential for antigen-mediated bronchoconstriction. J Immunol. (2009) 182:74309. 10.4049/jimmunol.090003919494266 Allen IC Hartney JM Coffman TM Penn RB Wess J Koller BH. Thromboxane A2 induces airway constriction through an M3 muscarinic acetylcholine receptor-dependent mechanism. Am J Physiol Lung Cell Mol Physiol. (2006) 290:L52633. 10.1152/ajplung.00340.200516243899 McDaniel DK Jo A Ringel-Scaia VM Coutermarsh-Ott S Rothschild DE Powell MD . TIPS pentacene loaded PEO-PDLLA core-shell nanoparticles have similar cellular uptake dynamics in M1 and M2 macrophages and in corresponding in vivo microenvironments. Nanomedicine. (2017) 13:125566. 10.1016/j.nano.2016.12.01528040495 Kononenko V Drobne D. In vitro cytotoxicity evaluation of the magneli phase titanium suboxides (TixO2x−1) on A549 human lung cells. Int J Mol Sci. (2019) 20:E196. 10.3390/ijms2001019630625978 Tsai CY Lu SL Hu CW Yeh CS Lee GB Lei HY. Size-dependent attenuation of TLR9 signaling by gold nanoparticles in macrophages. J Immunol. (2012) 188:6876. 10.4049/jimmunol.110034422156340 Kodali V Littke MH Tilton SC Teeguarden JG Shi L Frevert CW . Dysregulation of macrophage activation profiles by engineered nanoparticles. ACS Nano. (2013) 7:69977010. 10.1021/nn402145t23808590 Rahmani Kukia N Rasmi Y Abbasi A Koshoridze N Shirpoor A Burjanadze G . Bio-effects of TiO2 nanoparticles on human colorectal cancer and umbilical vein endothelial cell lines. Asian Pac J Cancer Prev. (2018) 19:28219. 10.22034/APJCP.2018.19.10.282130361551 Echtay KS Roussel D St.-Pierre J Jekabsons MB Cadenas S Stuart JA . Superoxide activates mitochondrial uncoupling proteins. Nature. (2002) 415:969. 10.1038/415096a11780125 Shabalina IG Nedergaard J. Mitochondrial ('mild') uncoupling and ROS production: physiologically relevant or not? Biochem Soc Trans. (2011) 39:13059. 10.1042/BST0391305 Manke A Wang L Rojanasakul Y. Mechanisms of nanoparticle-induced oxidative stress and toxicity. Biomed Res Int. (2013) 2013:942916. 10.1155/2013/94291624027766 Sonkusre P Cameotra SS. Biogenic selenium nanoparticles induce ROS-mediated necroptosis in PC-3 cancer cells through TNF activation. J Nanobiotechnol. (2017) 15:43. 10.1186/s12951-017-0276-328592284 Sun Y Miao H Ma S Zhang L You C Tang F . FePt-Cys nanoparticles induce ROS-dependent cell toxicity, and enhance chemo-radiation sensitivity of NSCLC cells in vivo and in vitro. Cancer Lett. (2018) 418:2740. 10.1016/j.canlet.2018.01.02429331422 McDaniel DK Ringel-Scaia VM Coutermarsh-Ott SL Allen IC. Utilizing the lung as a model to study nanoparticle-based drug delivery systems. Methods Mol Biol. (2018) 1831:17990. 10.1007/978-1-4939-8661-3_1330051432 Roberts RA Shen T Allen IC Hasan W DeSimone JM Ting JP. Analysis of the murine immune response to pulmonary delivery of precisely fabricated nano- and microscale particles. PLoS ONE. (2013) 8:e62115. 10.1371/journal.pone.006211523593509 Allen IC. The utilization of oropharyngeal intratracheal PAMP administration and bronchoalveolar lavage to evaluate the host immune response in mice. J Vis Exp. (2014). 10.3791/5139124747836 Han Y Lee DK Kim SH Lee S Jeon S Cho WS. High inflammogenic potential of rare earth oxide nanoparticles: the New Hazardous Entity. Nanotoxicology. (2018) 12:71228. 10.1080/17435390.2018.147231129848123 Mo Y Jiang M Zhang Y Wan R Li J Zhong CJ . Comparative mouse lung injury by nickel nanoparticles with differential surface modification. J Nanobiotechnol. (2019) 17:2. 10.1186/s12951-018-0436-030616599 Yoshii C Nagata N Tao Y Suematsu R Nikaido Y Kido M. Relationship between inflammatory cells in bronchoalveolar lavage fluid and pathologic changes in the lung interstitium. Respiration. (1998) 65:38692. 10.1159/0000293009782222 Akilesh HM Buechler MB Duggan JM Hahn WO Matta B Sun X . Chronic TLR7 and TLR9 signaling drives anemia via differentiation of specialized hemophagocytes. Science. (2019) 363:eaao5213. 10.1126/science.aao521330630901 Wang X Katwa P Podila R Chen P Ke PC Rao AM . Multi-walled carbon nanotube instillation impairs pulmonary function in C57BL/6 mice. Part Fibre Toxicol. (2011) 8:24. 10.1186/1743-8977-8-2421851604 Beamer CA Girtsman TA Seaver BP Finsaas KJ Migliaccio CT Perry VK . IL-33 mediates multi-walled carbon nanotube (MWCNT)-induced airway hyper-reactivity via the mobilization of innate helper cells in the lung. Nanotoxicology. (2013) 7:107081. 10.3109/17435390.2012.70223022686327 Olmedo DG Tasat DR Guglielmotti MB Cabrini RL. Effect of titanium dioxide on the oxidative metabolism of alveolar macrophages: an experimental study in rats. J Biomed Mater Res A. (2005) 73:1429. 10.1002/jbm.a.3023015742364 Long TC Tajuba J Sama P Saleh N Swartz C Parker J . Nanosize titanium dioxide stimulates reactive oxygen species in brain microglia and damages neurons in vitro. Environ Health Perspect. (2007) 115:16317. 10.1289/ehp.1021618007996 Wang C Hu X Gao Y Ji Y. ZnO nanoparticles treatment induces apoptosis by increasing intracellular ROS levels in LTEP-a-2 cells. Biomed Res Int. (2015) 2015:423287. 10.1155/2015/42328726339612 Jiang J Pi J Cai J. The advancing of Zinc Oxide nanoparticles for biomedical applications. Bioinorg Chem Appl. (2018) 2018:1062562. 10.1155/2018/106256230073019 Fraser JA Kemp S Young L Ross M Prach M Hutchison GR . Silver nanoparticles promote the emergence of heterogeneic human neutrophil sub-populations. Sci Rep. (2018) 8:7506. 10.1038/s41598-018-25854-229760395 Wang C Youle RJ. The role of mitochondria in apoptosis*. Annu Rev Genet. (2009) 43:95118. 10.1146/annurev-genet-102108-13485019659442 Alarifi S Ali H Alkahtani S Alessia MS. Regulation of apoptosis through bcl-2/bax proteins expression and DNA damage by nano-sized gadolinium oxide. Int J Nanomed. (2017) 12:454151. 10.2147/IJN.S13932628684914 Plackal Adimuriyil George B Kumar N Abrahamse H Ray SS. Apoptotic efficacy of multifaceted biosynthesized silver nanoparticles on human adenocarcinoma cells. Sci Rep. (2018) 8:14368. 10.1038/s41598-018-32480-530254325 Alleman RJ Katunga LA Nelson MA Brown DA Anderson EJ. The Goldilocks Zone from a redox perspective-Adaptive vs. deleterious responses to oxidative stress in striated muscle. Front Physiol. (2014) 5:358. 10.3389/fphys.2014.0035825278906 Davis GS Holmes CE Pfeiffer LM Hemenway DR. Lymphocytes, lymphokines, and silicosis. J Environ Pathol Toxicol Oncol. (2001) 20(Suppl 1):5365. 10.1615/JEnvironPatholToxicolOncol.v20.iSuppl.1.5011570674 Barbarin V Arras M Misson P Delos M McGarry B Phan SH . Characterization of the effect of interleukin-10 on silica-induced lung fibrosis in mice. Am J Respir Cell Mol Biol. (2004) 31:7885. 10.1165/rcmb.2003-0299OC14975940 Misson P van den Brule S Barbarin V Lison D Huaux F. Markers of macrophage differentiation in experimental silicosis. J Leukoc Biol. (2004) 76:92632. 10.1189/jlb.010401915292275 Barbarin V Xing Z Delos M Lison D Huaux F. Pulmonary overexpression of IL-10 augments lung fibrosis and Th2 responses induced by silica particles. Am J Physiol Lung Cell Mol Physiol. (2005) 288:L8418. 10.1152/ajplung.00329.200415608148

      Funding. This work was financially supported by the Virginia-Maryland College of Veterinary Medicine, the National Natural Science Foundation of China (41522111 and 4176114462), The Center for the Environmental Implications of Nanotechnology (NSF Cooperative Agreement EF-0830093), the Virginia Tech Institute for Critical Technology and Applied Science (ICTAS), the Open Foundation of East China Normal University, the Virginia Tech National Center for Earth and Environmental Nanotechnology Infrastructure (NSF Cooperative Agreement ECCS-1542100), the STAR Fellowship Assistance Agreement no. FP-91780101-1 awarded by the U.S. Environmental Protection Agency (EPA), The National Institutes of Health (R01 ES019311). Student work on this publication was supported by the National Institute of Allergy and Infectious Diseases Animal Model Research for Veterinarians (AMRV) training grant (T32-OD010430) and the American Association of Immunologist Careers in Immunology Fellowship Program (VR-S).

      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016fg500.com.cn
      hldmslc.com.cn
      www.mbservice.com.cn
      fbklqm.com.cn
      hyhid.com.cn
      fwdbcb.com.cn
      shengyu123.com.cn
      www.rlxeyt.com.cn
      wsxti.com.cn
      www.wzpetq.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p