In vivo model
Experiments were conducted on spf male BALB/c mice (Janvier Labs, Le Genest-Saint-Isle, France), 9 weeks old, 8 animals/group for each time points (3 and 6 month exposure) dispatched in 2 cages. The housing procedure respects the classical procedures with light and temperature control, free access to food and water and environmental enrichment. The number was determined according to our previous experiments with CS-exposed mice. The primary outcome measures were the alteration of the respiratory function and of the body weight. This mouse strain is described as sufficiently sensitive to oxidative stress and chemical induction of lung cancers [13].13 Animal procedures were in agreement with European directive 2010/63/EU for the protection of animals used for scientific purposes and obtained the Ethical Committee on Animal Experimentation (CEEA 75) approval (ref APAFIS #10363-2017062615002072v2). Animal randomisation was performed by the responsible of the animal facility at the delivery of the mice.
Animal body weights were recorded on Monday of each weak while clinical signs were monitored daily.
E-cigarettes and conventional cigarette
We chose the third generation “ModBox” model from NHOSS® (Innova, Bondues, France), used with the “Air Tank” clearomiser equipped with a 0.5 Ω kanthal coil and with a partially closed air flow. For experiments, we chose two power settings for the Modbox model: a “low” power of 18 W and a “high” power of 30 W. For the e-liquid, we chose the best-selling NHOSS® brand containing 65% propylene glycol, 35% glycerine, 16 mg/mL nicotine and the most common flavor, “blond tobacco”. Conventional 3R4F cigarettes were obtained from the University of Kentucky (Lexington, KY, USA)
To avoid chemical cross-contamination, two different pieces of equipment (dilution chamber, tubes, exposure towers and pipes) were used for e-cig and 3R4F exposures (Fig. 1). Aerosols from e-cigs and 3R4F cigarette were generated with an InExpose e-cigarette extension system on which we adapted the Modbox and a cigarette smoking robot (SCIREQ®, Emka technologies, Montreal, Quebec, Canada), respectively. Mice were exposed to aerosols by a nose-only tower (InExpose system, SCIREQ®, Emka technologies). In order to perform a comparative toxicological study of ECV and CS, we used the Health Canada Intense puff regime (55 mL puff volume, 2 s puff duration, 30 s puff period). Based on data from the literature and our preliminary study after a 4-day subacute exposure (data not shown), two exposure protocols were applied in this study for both e-cig and 3R4F emissions: a 3-month subchronic and a 6-month chronic exposure, 60 min/day and 5 days/week. For each exposure schedule, an additional group of mice was exposed to fresh conditioned air (negative control) and a period of acclimatisation of one week was performed before the exposure .
Sample Collection
Mice were euthanized 24 hours after the last exposure for both the protocols of 3 or 6 month exposures. Bronchoalveolar lavage (BAL) fluids, lungs, spleens and blood samples were collected and kept on ice until they were processed. In addition for lungs, one lobe was stored in RNA Later (Ambion,ThermoFisher Scientific, Illkirch, France) for transcriptomic analyses, posterior lobes were fixed with paraformaldehyde (PFA 4%, Labonord, Villeneuve d’Ascq, France) for histopathological analyses, and another lobe was frozen in liquid nitrogen and then stored at − 80°C until further analyses.
Broncho-alveolar lavage procedure, lung and spleen processing
BALs were performed by instilling 5 aliquots of 0,5 mL (final volume 2,5 mL) of sterile PBS. After centrifugation at 400g for 6 min at 4°C, supernatants (cell-free BAL fluid) from the first two aliquots of 0.5 ml were stored at − 80°C for cytokine analysis (ELISA), and cell pellets were used for flow-cytometry analysis.
The left lobe of the lung was mashed with a sterile blade and then digested with collagenase (Collagenase Type VI 17104–019 Gibco by Life technologies, Carlsbad, California United States) at 37°C. After 15 min of digestion, lungs were homogenized with an 18 G needle and further digested for 15 min. After centrifugation at 400g for 6 min at 4◦C, the pellets were resuspended in a 30% Percoll solution (Percoll TM GE Healthcare 17–0891-01, Chicago, IL, United States) and centrifuged at 500g for 15 min. Total spleen cells were also isolated from spleen and centrifuged at 400g for 6 min at 4°C. The lung and spleen pellets were resuspended in red blood cells (RBC) lysis buffer during 5 min at room temperature, to remove erythrocytes. The reaction of RBC lysis was stopped with PBS 2% FBS (Gibco by Life technologies, Carlsbad, California United States). After centrifugation at 400g for 6 min at 4 ◦C, pulmonary and spleen cells were resuspended in PBS 2% FBS, then enumerated and used for flow cytometry.
Flow cytometry
BAL, lung and spleen total cells were incubated with the appropriate panel of antibodies for 30 min in PBS 2% FCS. Conjugated antibodies were used against mouse CD5 (ref 130–102–574, FITC-conjugated), PBS57-loaded CD1 d Tetramer (NIH facility,PE-conjugated), NK1.1 (ref 130–103–963, PerCp-Cy5.5–conjugated), CD4 (ref 130–102– 411, PE-Cy7-conjugated), CD25 (ref 130–102–550, APC-conjugated), CD69 (ref 561–238, Alexa700-conjugated), TCRγδ (ref 130–104–016, APC-Vio770 conjugated), TCR-β (ref 130– 104–815, V450-conjugated), CD8 (ref 130–109–252, V500-conjugated), CD45 (ref BLE103140, BV605-conjugated), I-Ab (ref 130–102–168, FITC-conjugated), F4/80 (ref 130–102–422,PE conjugated), CD103 (ref 563–637, PerCP-Cy5.5-conjugated), CD11c (ref 558–079, PE Cy7-conjugated), CD86 (ref 560–581, Alexa-700 conjugated), Ly6G (ref 560–600, APC-H7conjugated), CD11b (ref 560–455, V45O conjugated), CD45 (ref 130–402–512, V500 conjugated), Ly6C (ref BLE128036, BV605-conjugated) (BD Biosciences, Franklin Lakes, United States; Biolegend, San Diego, United States and Myltenyi Biotech, Paris, France) and CCR2 (ref FAB 5538A, R&D systems, APC conjugated).
Data were acquired on a LSR Fortessa (BD Biosciences, Franklin Lakes, United States) and analyzed using FlowJo™ software v10.2 (Stanford, CA, USA). Gating strategy has been previously described [14]. Absolute cell numbers were calculated according to the total cell number and the frequency of CD45+ immune cells.
Cytokine Assay
Levels of IL-1β, IL-6, IL-22, IL-23, IFN-γ, CXCL1, CXCL2, CXCL5, CXCL17, TNF-α and Resistin were determined in BAL, lung and serum by enzyme-linked immunosorbent assay (ELISA) using the manufacturer’s recommendation (R&D systems, Biotechne, Minneapolis, MN, United States). In addition, concentrations of IL-17 and IL-23p19 were measured by an ELISA from Invitrogen (Waltham, MA, United States).
Lung Histology
Fixed lung lobes were paraffin-embedded and lung sections were stained with hematoxylin-eosin. Lung injury and inflammation were scored based on a scale evaluating bronchial damage, hyperplasia, inflammatory cell influx, alveolar exudates, damage, wall thickness, inflammation, and emphysema, all with a scale from 0 to 4. Moreover, hemorrhage, fibrinoid necrosis, leukocytoclasis and suppuration were noted as present or not (value of 1 and 0, respectively). This allows to generate a global histologic score from 0 to 36. The score was established in a blinded fashion by independent experts (OCV clinical research, Lille).
In order to assess the potential development of emphysema, we evaluated the mean linear intercept (MLI) on lung sections by using the Image J software (NIH).
Lung function
On the day of the sacrifice, each mouse was anaesthetized intraperitoneally with xylazine hydrochloride (15 mg/kg) / ketamine (100 mg/kg), tracheotomized, and cannulated with an 18 G metal cannula (resistance: 0.36–0.40 cmH2O.s/mL). Mice were then connected to a flexiVent FX system (SCIREQ Inc., Montreal, Qc, Canada) and operated by the flexiWare software v7.7 as previously described [15]. Immediately after connection to the ventilator, set at 150 breaths/min, two deep lung inflations were performed at least 6–12 seconds apart to recruit lung beyond any closed airway and to standardize lung volume history. This was done by inflating the lungs to 30 cm H2O over 3 seconds and holding that pressure for another 3 seconds to allow for the lungs to equilibrate after the inflation. Mice were then submitted to a 300 breaths/min hyperventilation in order to eliminate spontaneous breathing before a 150 breaths/min return. Next, the mechanical properties of the mouse respiratory system were assessed at baseline, i.e. before the construction of a full-range pressure-volume (PV) curve. This was done using a sequence of measurements integrated by default in the flexiVent operating software (flexiWare v7.7). The Area between the PV loop inflation and deflation limbs (Hysteresis (H)) A static compliance (Cst) was calculated. Cst, the parameters A (estimate of inspiratory capacity) and K (shape constant) can be extracted from the Salazar-Knowles equation [16]. Cst reflects the intrinsic elastic properties of the respiratory system (i.e. lung + chest wall).
Transcriptomic analyses
Total RNAs were extracted from RNA later treated tissue samples using the miRNeasy mini kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s protocol. The RNA concentration was measured with the Biospecnano spectrophotometer (Shimadzu, Marne-la-Vallée, France). Transcriptomic experiments were carried out using 8x60k OneColor microarrays (Agilent Technologies) coupled to 60-mer oligonucleotides covering the entire mouse genome. Labeling, hybridization and lncRNA detection were carried out according to the manufacturer's instructions (Agilent Technologies). For each microarray, Cyanine 3-coupled lncRNAs were synthesized by the QuickAmp Low Input Kit from 50 ng of total RNA. Spike-in RNAs were added to each tube and used as positive controls for the labeling and amplification steps. Labeled cDNAs were purified and 600 ng of each cDNA was then hybridized to the microarrays according to the manufacturer's instructions. After washing, the microarrays were scanned and the data exported using Agilent Feature Extraction Software© (FE version 10.7.3.1). Results were then interpreted by selecting those mRNAs whose expression was at least significantly 1.5-fold higher or 1.5-fold lower than that of unexposed control mice (p < 0.05). Statistical analyses were performed with the "linear models for microarray data" (limma) package for R, using moderated t statistics with standardized data. Functional analysis of selected deregulated lung RNAs was performed using Ingenuity Pathway Analysis software (Qiagen) by selecting the top 25 of canonical pathways for each exposure condition.
Statistical analyses
The data are expressed as mean ± Standard Error of the Mean (SEM). Statistical analyses were computed using GraphPad Prism software (5.00 version, GraphPad software, San Diego, USA) in accordance with the size of samples and the nature of the experiment. Values were compared to the controls with a bilateral and non-parametric Mann-Whitney test, or with a one-sample t-test, or a two-way ANOVA test when appropriate. Statistical significance was accepted for an error risk inferior to 5% and is represented as follows: *, p < 0.05; **, p < 0.01; ***, p < 0.001.