Airway inflammation is highly prevalent in horses, with the majority of non-infectious cases being defined as equine asthma. Currently, cytological analysis of airway derived samples is the principal method of assessing lower airway inflammation. Samples can be obtained by tracheal wash (TW) or by lavage of the lower respiratory tract (bronchoalveolar lavage fluid; BALF). Although BALF cytology carries significant diagnostic advantages over TW cytology, sample acquisition is invasive, making it prohibitive for routine and sequential-screening of airway health. The aim of this study was to establish a robust protocol to isolate macrophages, protein and RNA for molecular characterisation of TW samples and demonstrate the applicability of sample handling to rodent and human pediatric bronchoalveolar lavage fluid isolates. TW samples provided a good quality and yield of both RNA and protein for downstream transcriptomic/proteomic analyses. The sample handling methodologies were successfully applicable to BALF for rodent and human research. TW samples represent a rich source of airway cells, and molecular analysis to facilitate and study airway inflammation, based on both transcriptomic and proteomic analysis. This study provides a necessary methodological platform for future transcriptomic and/or proteomic studies on equine lower respiratory tract secretions and BALF samples from humans and mice.
Research Article
A One Health Approach to Simple Liquid Sample Handling for Respiratory Based -Omics Analysis; Tracheal Wash and Bronchoalveolar Lavage Fluid
https://doi.org/10.21203/rs.3.rs-279312/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 11 Jul, 2021
You are reading this latest preprint version
The European Thoroughbred horseracing sector has an annual economic impact of €12 billion, including approximately 155,000 employees in a variety of roles (British Horseracing Authority). Racehorses commonly develop airway inflammation and /or exercise induced pulmonary haemorrhage during training, with prevalence rates as high as 70-80% and 100% respectively, resulting in a significant impact on animal welfare and the racehorse industry 1-3. Cytological analysis of tracheal wash (TW) and bronchoalveolar lavage fluid (BALF) is the most commonly applied method of assessing lower respiratory tract health of the horse and helps to guide therapy 4,5. The diagnostic application of this technique includes the identification of common equine airway disorders of the horse, including equine asthma and exercise induced pulmonary haemorrhage 6,7. Tracheal wash, BALF, peritoneal fluid and blood samples are commonly obtained by equine veterinarians and analysed for diagnostic purposes 8-11. Despite that, a limited number of equine airway derived gene expression studies have been reported in the literature, largely aimed at defining the pathogenesis of equine asthma 12-14. In each of these studies, the airway cells were derived from BALF samples, which are more difficult to obtain due to the invasive collection procedure. In comparison, the analyses of TW samples rarely extend beyond cytology and subjective scoring of mucus content. However, recent technological advances in sample collection and processing procedures have made it possible to determine whether a wider range of analyses might be applied to TW samples. This would be advantageous as equine tracheal secretions represent a rich source of: (i) cells, (ii) host-derived microbiome, (iii) transcriptome and (iv) local environment-derived proteins; and (v) are relatively simple to collect with minimal discomfort to the horse.
Similar approaches have already been applied to a comparable sample set in humans15,16; namely, induced sputum. Induced sputum in humans provides a non-invasive method to sample airway biofluids that are enriched in proteins and may reveal novel mechanisms in the pathogenesis of respiratory diseases, such as asthma and chronic obstructive pulmonary disease 16,17.For equine derived samples, the proteomic approach is novel and follows a recent upsurge of interest in proteomic analysis of human airway secretions as a highly informative non-invasive method of assessing the mechanisms associated with lung disease 15. Application of both a transcriptomic and proteomic approach to equine TW samples has the potential to reveal novel biomarkers of early disease development or disease associated tissue remodelling. Furthermore, such investigations could provide additional criteria on which to base disease subcategorization (e.g., equine asthma) and/or identify novel therapeutic targets.
Moreover, the horse is widely acknowledged as an attractive model for human asthma 3,18,19, addressing needs that cannot be addressed through small animal models: mice are widely recognised as poor models of the human immune response, particularly in relation to the lung 20,21. In addition to the important similarities between horses and humans with regard to both pathophysiology and macrophage/monocyte biology, the horse also represents an ideal source of large volumes of various biological sample types 18,19,22-25. For example, the volume and the associated cell retrieval rates derived from horses represent at least 2‐3 orders of magnitude greater than those obtained from rodent models. This is of particular relevance in light of the extremely limited availability of human derived samples.
Thus, in this study we investigated the most appropriate means of (a) isolating airway macrophages from TW and (b) maximising the yield and quality of RNA and protein from equine TW samples. Furthermore, we compared the relevance of our results to murine and paediatric samples, determining whether a similar approach might be adopted.
Samples used in the study
Equine samples
Tracheal secretions were collected from 39 (38 male and 1 female) Thoroughbred horses [7+0.2 years (mean+SEM); range 4-12 years]. All horses were housed at the same racing yard and samples were obtained as part of a routine assessment of respiratory health (differential cytology). Residual sample was retained for an ongoing proteomic and transcriptomic study. Sample collection procedures were approved by the Veterinary Ethical Review Committee of the School of Veterinary Medicine, University of Edinburgh, following the relevant guidelines and regulations relating to the provisions of the Animals (Scientific Procedures) Act 1986. All methods reported are in accordance with ARRIVE guidelines (https://arriveguidelines.org) 40. Informed consent was obtained from the trainer for the use of residual sample material. Four BALF samples were collected at the Royal (Dick) School of Veterinary Studies, University of Edinburgh. All horses were clinically healthy and treated according to standard welfare procedures.
Murine samples
Bronchoalveolar lavage fluid were generously provided by Prof Jurgen Schwarze (Medical Research Council Centre for Inflammation Research, University of Edinburgh, UK). Samples were collected from four wild type BALB/c female mice, aged 8-12 weeks, at four days following intranasal administration of UV treated Respiratory syncytial Virus (RSV) as described by Jurgen Schwarze’s et al (2107) 41. These mice represented controls to compare against models of RSV infection for investigating the influence of RSV infection on tight junction integrity. All experimental procedures were carried out in accordance with ARRIVE guidelines (https://arriveguidelines.org)40,41. All experimental protocols were approved by the Ethical Review Committee of the Medical Research Council Centre for Inflammation Research, University of Edinburgh, UK, as previously described 41. All methods were carried out in accordance with relevant guidelines and regulations of the Home Office regulations 41.
Mice were housed
at the Medical Research Council Centre for Inflammation
Research, University of Edinburgh, UK in individually ven-
tilated cages for the duration of the study, and all experi-
mental procedures were carried out in accordance with
Home Office regulations
Mice were housed
at the Medical Research Council Centre for Inflammation
Research, University of Edinburgh, UK in individually ven-
tilated cages for the duration of the study, and all experi-
mental procedures were carried out in accordance with
Home Office regulations
Paediatric samples
Bronchoalveolar lavage fluid was collected from three children, with no evidence of lung disease, who were undergoing planned procedures under general anaesthetic. Children were attending for planned surgical procedures (inguinal hernia and hypospadias repair and developmental dysplasia of the hip repair) and were aged between six months and five years (2.5+1 years old, 2 female, 1 male). Children were sampled at the Royal Hospital for Sick Children Edinburgh as part of the Lower respiratory tract infection Kids Study approved by the South East Scotland Research Ethics Committee 02. Informed consent was obtained from appropriate adult with parental responsibility. The clinical study and all associated procedures were conducted within Good Clinical Practice guidelines and the regulations stipulated by the MHRA (Medicines Healthcare Regulation Agency).
Sample collection
Equine samples
All TW samples were collected one hour following morning exercise (interval training) using a trans-endoscopic technique. Horses were restrained and a nose twitch applied when necessary. A 1500mm working length, 9.2 mm outer diameter video endoscope (2.8mm biopsy channel; Aohua, China), was passed via the ventral nasal meatus to the pharynx and then advanced into the trachea via the rima glottidis. Following assessment of the amount and nature of mucus deposits within the tracheal lumen, approximately 20 mL of sterile 0.9% saline at room temperature was instilled via a catheter passed via the biopsy channel of the endoscope at the proximal aspect of the trachea. The endoscope was then further advanced to the level of the thoracic inlet where the pool of instilled fluid had gravitated. As much fluid as possible was subsequently aspirated via the transendoscopic catheter. Samples were stored on ice and submitted for laboratory analysis at the Roslin Institute and the Royal (Dick) School of Veterinary Studies and processed within 4-5 hours of collection. Three BALF samples were collected at the Royal (Dick) School of Veterinary studies from horses admitted to the hospital, as previously described 42; these procedures were performed for diagnostic purposes, with residual BALF made available for this study with owner consent.
Murine samples
Bronchoalveolar lavage fluid was collected by injecting 1ml of PBS containing protease inhibitor via tracheal cannula – 650-900µl of fluid being retrieved. BALF supernatant was then collected and stored at -80oC till further use.
Paediatric samples
While anaesthetised a suction catheter was inserted blindly down the endotracheal tube and 2ml/kg of 0.9% sodium chloride instilled and around 1ml/kg of fluid being retrieved. Following a short period of time, and without removing the suction catheter, sample was collected by applying 10-15kpa of suction to said catheter. Resultant BALF was spun at 400g at 4oC for 10 minutes, the supernatant collected and then frozen at -80°C for storage and further analysis.
Total and differential cell count and total cell and macrophage isolation for equine tracheal wash samples
Differential cell count
An aliquot of 0.5 ml was submitted to the pathology lab at the Royal (Dick) School of Veterinary Studies for DCC analysis. Differential leucocyte count (minimum of 200 cells) was performed and expressed as a percentage of total non-squamous and non-epithelial nucleated cells. Horses were considered free of inflammatory airway disorders based on the following DCC cut off values: neutrophils <20%.
Total cell isolation
First, a 1ml aliquot of sample was immediately stored at -80°C for future proteomic analyses and the remaining sample incubated for 15 minutes at room temperature in 0.1% dithiothreitol (DTT) to depolymerize secreted mucin. Dithiothreitol has been demonstrated to cause no deleterious effects on human sputum derived cells or interfere with surface marker measurements using flow cytometry. Following gravity filtration through a 70µm pore mesh filter, the sample was centrifuged at 400g for 10min at 4oC. Supernatant was carefully removed, and the cell pellet re-suspended in Dubelcco’s PBS, from which a total cell count (excluding squamous epithelial cells) and cell viability (Trypan Blue exclusion staining) was performed using a haemocytometer. Afterwards, the cell pellet was resuspended in 1ml of Trizol and stored at -80oC for future RNA analysis.
Macrophage isolation
Following total cell isolation, immunomagnetic separation was used to isolate macrophages from DTT treated TW samples from 16 of the 39 animals, using magnetic beads coated with mouse anti-human CD163 antibody (GHI/61, BioLegend, cat no 333605), according to the manufacturer’s instructions. We have previously demonstrated cross-reactivity of this antibody with equine macrophages/monocytes. Briefly, the cells were first stained with a R-Phycoerythrin (PE)-conjugated CD163 primary antibody; subsequently, the cells were magnetically labelled with anti-PE MicroBeads (Miltenyi Biotec, cat no 130-105-639) before the cell suspension was loaded on a MACS® LS Column (Miltenyi Biotec Ltd., cat n0 130-042-401), placed in the magnetic field of a MACS Separator. The magnetically labeled cells were retained within the column, while the un-labeled cells passed through. Following removal of the column from the magnetic field, magnetically retained cells were eluted as the positively selected cell fraction.
Flow cytometry analysis of tracheal wash derived cells
Total TW cells were stained with anti-CD163 antibody against macrophages (BioLegend, Catalogue No. 333605). Cells were also stained with Zombie aqua dye (Biolegend, Catalogue No. 423101) to test for viability. A negative/isotype control was used to control for non-specific binding (PE Mouse IgG1, κ Isotype Ctrl, Biolegend, Catalogue No. 400111). Data were acquired on 10,000 live cell events and cells gated according to size (FSC-A: forward scatter) and granularity (SSC-A: side scatter) after removing artifacts, debris, doublet discrimination and dead cells. Data analysis was performed in FlowJo® v10.5.3 (Tree Star).
RNA analysis of equine tracheal wash derived cells
Total RNA extraction
Total RNA was extracted using a combination of Trizol reagent (Thermo Scientific™, 15596026) and an RNAeasy plus micro kit (Qiagen, cat no 74034), according to manufacturer’s instructions. Briefly, following removal of the supernatant, the remaining cell pellet was lysed by the addition of 1mL Trizol Reagent. Subsequently, 0.2 mL 1-Bromo-3-chloropropane (BCP) (Sigma Aldrich, cat no B9673-200ML) was added per 1 mL of Trizol. The sample was shaken vigorously for approximately 30 sec and left at room temperature for 5 min to allow complete dissociation of nuclear-protein complexes. The homogenate was then centrifuged at 18,000g for 15 min at 40C resulting in the formation of a lower red phenol-chloroform phase, an interphase, and an upper colourless aqueous phase. The aqueous phase contained the RNA and had almost 50% of the volume of the Trizol used, plus the volume of the sample. Following transfer to a clean tube for the precipitation step, 0.5 mL of 70% ethanol was added, the sample stored for 2 hours at -200C and then transferred to an RNeasy spin column and centrifuged at 18,000g for 5 min at 40C. Following centrifugation, the flow through was removed, the RNA washed once with RW1 buffer and DNA treatment was performed using the RNase-Free DNase Set (Qiagen, cat no 79254) according to the manufacturer’s instruction. Due to potential DNA contamination, this step was performed twice and samples were run through gDNA Eliminator Spin Columns twice after the elution step. Afterwards, the RNA membrane was washed with RW1, RPE and 80% ethanol. Finally, RNA was eluted in 20 µl RNase-free water and RNA samples were stored at -80oC until further use.
RNA quality assessment
RNA concentration and purity were measured using ND-1000 Nanodrop spectrophotometer (Thermo Scientific, Wilmington, USA) by measuring absorbance at 260 and 280 nm (A260, and A280 respectively). Purity of RNA was determined using the A260/A280 ratio. A ratio close to two was considered to be indicative of pure RNA. RNA integrity was confirmed with the High Sensitivity RNA ScreenTape system (Agilent Technologies). A RIN number greater than seven was considered appropriate for RT-qPCR and RNAseq analysis.
Protein Analysis
Briefly, an aliquot of 500µl of untreated TW or BALF samples was homogenized in protein extraction buffer (100mM Tris, pH 7.6 and 4%w/v SDS) +1% Halt Protease Inhibitor Cocktail, EDTA-Free (Thermo Scientific, cat no 87785). Following homogenization, samples were centrifuged at 20,000 g for 20 min at 100C. The supernatant containing the solubilized protein was removed and stored at -80oC. Protein concentration of samples was determined using a Micro BCA Protein Assay Kit (Thermo Scientific™, cat no 23235) according to the manufacturer’s instructions. Finally, total protein analysis was carried out for quality control purposes and to determine equivocal protein load between samples. Samples were separated by electrophoresis on gradient gels (NuPAGE 4-12% Bis-Tris Protein Gels, 1.0 mm, 12-well, Fisher Scientific, cat no NP0322BOX) and stained with InstantBlue™ Protein Stain (Expedeon Ltd, cat no ISB1L) as previously described 43. The stained gel was then imaged using the LICOR Odyssey imager to visualise and quantify the total protein load within each lane of the gel using associated Image Studio Software.
Trachea wash sample preparation for LC-MS
Sample preparation for LC-MS was carried out using S-Trap micro spin column digestion protocol (Protifi, Huntington, NY). Samples containing 10 μg of protein, were diluted to 50 μL with lysis buffer (5% SDS in 50-mM TEAB) and reduced with 20mM Dithiothreitol at 95 °C for 10 min. After cooling to room temperature, samples were alkylated in the dark for 30 min with 40 mM of iodoacetamide. This mix was then acidified to 1% phosphoric acid and proteins were precipitated with 350 μL of binding buffer [90% methanol, 100-mM Triethyl ammonium bicarbonate (TEAB)]. This protein suspension was loaded onto S-Trap column columns and centrifuged at 4000×g for 30 s. After four washes with 125 μL of 50-mM TEAB, trypsin (0.5 μg in 50mM TEAB) was added to the trap and incubated overnight. Peptides were then eluted with 80 μL of each of the following: 50-mM TEAB, 0.2% formic acid (FA), and 50% acetonitrile in 0.2% FA with centrifugation at 1000g after each elution step. All elution fractions were pooled and dried under vacuum and cleaned up by stage-tipping 44
LC-MS analysis on tracheal wash samples
NanoflowLC-MS/MS was performed on a micrOTOF-II mass spectrometer (Bruker, Germany) coupled to a RSLCnano LC system (Thermo) following earlier method 45 with minor modifications. Raw spectral data were processed with DataAnalysis (Bruker) software and the resulting peak lists were searched using Mascot 2.4 server (Matrix Science) against Uniprot horse sequence database (Uniprot ID UP000002281), containing 44,485 entries. Precursor and fragment ion mass tolerance were set at 25 ppm and 0.06 Da respectively and false discovery rate was set at < 1% for peptide IDs after searching decoy databases.
Protein function annotation and pathway analysis
As previously described 46, identification of enriched KEGG pathways and biological processes in the protein list was performed with Database for Annotation, Visualization and Integrated Discovery (DAVID) software (v6.8) 27,47,48. KEGG pathway database contains regulatory and metabolic pathways, representing global knowledge on molecular interactions and reaction networks. Pathway analysis was performed using Ingenuity Pathway Analysis (IPA Spring 2021 release, https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis) 28 to infer the functional roles and relationships of the detected proteins.
The aim of this project was to highlight the usefulness of TW samples for several biological applications (Figure 1) and to establish the conditions for optimal RNA and protein isolation, in order to use for transcriptomic and proteomic analysis. Thus, TW samples from 39 Thoroughbred racehorses were collected. Furthermore, the proteomic analysis protocol was successfully applied to equine, murine and human BALF samples.
Establishing a protocol for tracheal wash cell isolation
One aspect of the current project was the isolation of macrophages from TW samples. Analysis had confirmed cross-reactivity of an anti-human CD163 antibody, a classical marker for mature macrophages, and successful application of magnetic bead separation to dithiothreitol (DTT)-treated equine tracheal secretion samples. Flow cytometry analysis revealed that almost half of the cell population of TW samples were CD163 positive (Figure 2a-c). Isolation of CD163 positive cells was also confirmed with light microscopy on cytospin slide preparations stained with Leishman stain (Figure 2d).
The present study used a series of 39 horse-derived samples and an average of 5.9x106 (+1.7) cells were isolated from TW samples. RNA average yield extracted from the total cell population of TW samples was 244 + 43ng/µl. RNA yield concentration of CD163 positive cells was 43 + 19ng/µl. RNA integrity number (RIN) greater than 7 is recommended for RNAseq and qPCR analysis. RNA samples derived from the total population of TW-derived cells had an average RIN number of 7.92 + 0.14 (SEM), thus rendering samples suitable for sequencing or RNA analysis. Differential cell counts (DCC) of TW samples were as follows (mean ± SEM): 38.5 ± 3 % macrophages, 17.8 + 2.8 haemosiderophages, 30.4 ± 1.7% lymphocytes, 13.3 + 2.4 % neutrophils, and 0.01 ± 0.02 eosinophils. Results of the DCC are shown in Figure 3. In agreement with the data derived from flow cytometric analysis (Figure 2b) almost half the number of TW cells were macrophages (including haemosiderophages). Since all samples were collected from racehorses, the presence of haemosiderophages in these samples was unavoidable as almost all racing Thoroughbred horses will bleed into the airways; this condition is termed Exercise-Induced Pulmonary Haemorrhage 6.
Protein extraction of tracheal wash samples, has been successfully performed
In order to expand on the RNA studies, we planned to define the airway protein profiles (total proteome), thus revealing the mechanisms which may underpin any alterations in immune function. Protein extraction was successfully performed on TW samples from 39 horses. An average of 1.3 ± 0.2 (+SEM) mg protein has been isolated from 500µl of TW per animal. Protein extraction of equine BALF samples resulted in 0.8+ 0.07 (+SEM) mg protein isolated from 500µl. The same protocol was successfully applied on the same volume of human (n=3) and murine (n=4) BALF samples, from which a protein yield of 0.6± 0.1 (+SEM) and 0.4± 0.1 (+SEM) mg total protein were isolated, respectively.
In order to visualise the total protein load, all samples were run on gradient gels and stained with instant blue protein stain. Figure 4 is representative of a gel stain of three equine TW, two BALF samples, as well as two human and four murine samples; all sample types showed a similar pattern in each specie. Protein samples were also submitted for proteomic analysis.
Lane 1 shows the protein ladder. Total protein stain of equine (Eq) TW samples (Lanes 2-4) and BALF samples (Lanes 5-6). Lanes 7-8 show total protein stain from human (Hum) BALF samples and lanes 9-12 contain total protein stain of four murine (Mur) BALF samples. Selected prominent bands isolated from the gel and analysed in the Mass Spectometry Facility are shown in white boxes (Lane 2, B1-B7). Note the increased diversity of proteins derived from equine TW samples compared to BALF.
In order to investigate the variety of proteins present in the samples, selected prominent bands were isolated from the gel and analysed in the Mass Spectometry Facility at the Roslin Institute. Uniprot, a comprehensive resource for protein sequence and annotation data was used (https://www.uniprot.org/proteomes/UP000002281) for horse protein annotation.
Initial quality control analysis based on the staining patterns on gel clearly shows protein samples were intact and showed no signs of degradation, in spite of the sample preparation and addition of extraction buffer components directly into liquid biological samples. Also, a range of protein identities were revealed from a few bands analysed by LC-MS (Supplementary Data 1).
Furthermore, in an attempt to investigate whether DTT treatment has an effect on TW derived protein; a comparison was also made between DTT treated and untreated TW samples. The addition of DTT during sputum processing is widely used to release cells from mucus. However, it is reported that DTT treatment may result in decreased detection of inflammatory cytokines in human sputum compared to untreated samples 26. In line with this, our results suggested that the yield of unique protein identifications from the untreated samples was double that derived from DDT treated samples (Supplementary Data 1). Therefore, we decided to proceed with proteomic analysis of TW samples that had not been treated with DTT. In total, 436 unique proteins were detected from all samples. Some example proteins detected in the TW samples are shown in Table 1. These are related to immune response and lung biology, including CD14, complement factors, immunoglobulins, surfactant proteins and mucins, suggesting that samples isolated and processed in this way are likely to be representative of the in vivo composition. The complete lists of mapped proteins are presented in the Supplementary Data 1.
KEGG pathway enrichment analysis of the detected proteins was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) annotation software (https://david.ncifcrf.gov/) 27. The detected protein list included proteins involved in metabolic pathways and immune response to infection (Supplementary Data 2). Ingenuity Pathway Analysis (IPA) was also used to functionally analyse the list of the detected proteins 28. The main canonical pathways and diseases/biofunctions related to the protein list are shown in Figure 5. The IPA-generated network (Figure 6) contains most of the molecules related to “Cellular Function” and “Maintenance”, “Humoral Immune response” and “Inflammatory Response”.
Table 1. List of proteins related to immune response and lung biology derived from tracheal wash samples.
To the best of our knowledge, this study has shown for the first time that it is feasible to carry out both transcriptomic and proteomic analysis on a single TW collected from equines, in this case UK national hunt Thoroughbreds. Such samples are routinely collected in veterinary practice, in contrast to the less readily available and more technically invasive BALF-derived samples, which have previously been favoured for transcriptomic and proteomic analyses14,29,30. Moreover, macrophages were successfully isolated from TW samples based on magnetic bead separation and isolated cells were deemed suitable for downstream transcriptomic analysis. Similar experiments have been performed on human isolated macrophages from induced sputum 31 but this is the first time it has been carried out in the equine species.
Proteomic analysis revealed a wide range of proteins in the equine TW samples, many of which have already been detected and reported in previously published studies on human and equine BALF and human sputum 30,32-34. Some of the detected proteins have previously been proposed as potential biomarkers for human lung cancer; these include carbonic anhydrase 1 (CA1), apolipoprotein A2 (APOA2), transthyretin (TTR), clusterin CLU, ceruloplasmin (CP) and gelsolin (GSN) 33. In a similar context, TW proteomics can be envisioned in terms of biomarker screening for the development of ‘real’ biomarkers with potential future clinical utility in relation to two common respiratory diseases in horses; namely, equine asthma and exercise induced pulmonary haemorrhage 6,35. Moreover, results may have translational application to human respiratory disorders, such as equine asthma or respiratory infections.
Several proteins were related to important biological processes of the respiratory system, such as cellular function/differentiation [annexin (ANXA4, ANXA7), cathepsin B (CTSB), transgelin 2 (TAGLN2)], immune defense/response [CD14 molecule (CD14), lysozyme (LYZ), complement components (C1QB, C1QC, C3, C4A, CFB)] and antigen processing and presentation [MHC class II DR alpha chain (DRA), beta-2-microglobulin (B2M), heat shock proteins (HSPA1A, HSP90AB1)]. Interestingly, MUC5AC, also detected in our analyses, has previously been shown to be one of the predominant mucins in airway secretions derived from healthy horses and humans. Both MUC5AC and MUC5B are products of epithelial goblet cells and submucosal glands 36-38. Finally, consistent with the cytological profile of the TW samples and the high presence of haemosiderophages, several proteins related to haematopoietic homeostasis were also present in the equine TW samples; examples include haemoglobin, zeta (HBZ1) and glucose-6-phosphate dehydrogenase (G6PD).
The non-bronchoscopic method of collecting human BALF is well recognised in paediatrics, where a formal bronchoscopy is both technically challenging and difficult to justify without clear clinical need. As it is relatively less invasive, this alternative technique has several advantages over formal bronchoscopy in the paediatric setting. Furthermore, it requires relatively minimal training, avoids the requirement for expensive equipment and is quick and well tolerated by both patients with and without chest infection. Its adoption is therefore suited to medical research in resource limited settings. Although previously used to successfully obtain dendritic cells from healthy children and those with RSV bronchiolitis 39, this is the first report of its use to obtain samples for proteomic analysis.
In conclusion, we have developed a robust protocol for obtaining a high quantity and quality of RNA and protein for transcriptomic and proteomic analysis from a single equine TW sample. Our “proof of principle” results should facilitate future gene-expression and extensive proteomic studies offering deep proteome coverage utilizing TW samples. These technological advances have the real capacity to improve our understanding of the function of the pulmonary immune system of the horse. They may also permit a more definitive and comprehensive characterization of common diseases such as equine asthma, potentially revealing not only disease susceptibility-associated biomarkers but also novel therapeutic targets. Finally, comparative global analysis of gene/protein expression across airway derived samples between horse and human would allow further investigation of the horse as an appropriate candidate model for human disease.
Acknowledgments
The authors would like to thank all the horse owners that gave their consent to use their animals for the purpose of the current study and the clinicians from the Royal (Dick) School of Veterinary Studies for suppling us with the tracheal wash samples; the pathology lab at the Royal (Dick) School of Veterinary Studies for differential cell count analysis and the Mass Spectometry Facility at the Roslin Institute for conducting the proteomic experiments; Professor Jurgen Schwarze (Medical Research Council Centre for Inflammation Research, University of Edinburgh, Edinburgh, UK) for kindly providing the murine samples for this study. Finally, we would like to thank the Horserace Betting Levy Board for funding this study.
Author contributions
A.E.K., S.R.P. and T.M.W. designed the study. A.E.K. led the laboratory work, the analyses and the manuscript writing. E.C.G., S.L.E., J.T.P. and A.R. contributed to the laboratory methods and sample collection. D.K. performed the proteomic analysis. All authors reviewed and approved the final manuscript.
Additional information
Supplementary Information:
Supplementary Data 1 - Complete lists of mapped proteins
Supplementary Data 2 - KEGG pathway and Gene Ontology enrichment analysis of the detected proteins
Competing financial interests:
One of the authors of this manuscript is an academic editor for Scientific Reports.
Dr T M Wishart , Reader in Molecular Anatomy, The Roslin Institute, Royal (Dick) School of Veterinary Studies, College of Medicine and Veterinary Medicine, University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, SCOTLAND, UK, T.M.Wishart@Roslin.ed.ac.uk
Research explorer: https://www.research.ed.ac.uk/portal/en/persons/thomas-wishart(3d17d4e1-df57-4c19-8d26-2cd1e5a5aaac)/publications.html
Academic Lead - Proteomics and Metabolomics Facility; Roslin Institute
https://www.ed.ac.uk/roslin/facilities-resources/proteomics-and-metabolomics-facility
Co-Head - Translational Biomarker Development; Centre for Dementia Prevention; University of Edinburgh http://centrefordementiaprevention.com/research/translational-research-groups/
None of the co-authors have any conflict of interest
- Wood, J. L. N., Newton, J. R., Chanter, N. & Mumford, J. A. Inflammatory airway disease, nasal discharge and respiratory infections in young British racehorses. Equine Vet. J. 37, 236–242 (2005).
- Morley, P. S., Bromberek, J. L., Saulez, M. N., Hinchcliff, K. W. & Guthrie, A. J. Exercise-induced pulmonary haemorrhage impairs racing performance in Thoroughbred racehorses. Equine Vet J. 47, 358–365 https://doi.org/10.1111/evj.12368 (2015).
- Couetil, L. et al. Equine Asthma: Current Understanding and Future Directions. Front Vet Sci. 7, 450 https://doi.org/10.3389/fvets.2020.00450 (2020).
- Rossi, H. et al. Comparison of Tracheal Wash and Bronchoalveolar Lavage Cytology in 154 Horses With and Without Respiratory Signs in a Referral Hospital Over 2009 – 2015. Frontiers in Veterinary Science. 5, https://doi.org/10.3389/fvets.2018.00061 (2018).
- Couetil, L. L. & Thompson, C. A. Airway Diagnostics: Bronchoalveolar Lavage, Tracheal Wash, and Pleural Fluid. Vet Clin North Am Equine Pract. 36, 87–103 https://doi.org/10.1016/j.cveq.2019.12.006 (2020).
- Hinchcliff, K. W. et al. Exercise induced pulmonary hemorrhage in horses: American College of Veterinary Internal Medicine consensus statement. Journal of veterinary internal medicine. 29, 743–758 https://doi.org/10.1111/jvim.12593 (2015).
- Couetil, L. L. et al. Inflammatory Airway Disease of Horses–Revised Consensus Statement. 30,503–515, doi:10.1111/jvim.13824 (2016).
- Hoffman, A. M. Bronchoalveolar lavage technique and cytological diagnosis of small airway inflammatory disease. Equine Veterinary Education. 11, 330–336 (1999).
- Harrold, S. M. et al. Tissue sites of persistent infection and active replication of equine infectious anemia virus during acute disease and asymptomatic infection in experimentally infected equids. Journal of Virology. 74, 3112–3121 https://doi.org/10.1128/jvi.74.7.3112-3121.2000 (2000).
- Richard, E. A. et al. Sub-clinical diseases affecting performance in Standardbred trotters: Diagnostic methods and predictive parameters. Vet. J. 184, 282–289 https://doi.org/10.1016/j.tvjl.2009.04.016 (2010).
- Christley, R. M. et al. Coughing in thoroughbred racehorses: risk factors and tracheal endoscopic and cytological findings. Vet. Rec. 148, 99–104 (2001).
- Lavoie, J. P. et al. Bronchoalveolar Lavage Fluid Cytology and Cytokine Messenger Ribonucleic Acid Expression of Racehorses with Exercise Intolerance and Lower Airway Inflammation. Journal of veterinary internal medicine. 25, 322–329 https://doi.org/10.1111/j.1939-1676.2010.0664.x (2011).
- Hughes, K. J. et al. Evaluation of cytokine mRNA expression in bronchoalveolar lavage cells from horses with inflammatory airway disease. Veterinary Immunology and Immunopathology. 140, 82–89 https://doi.org/10.1016/j.vetimm.2010.11.018 (2011).
- Beekman, L., Tohver, T. & Léguillette, R. Comparison of cytokine mRNA expression in the bronchoalveolar lavage fluid of horses with inflammatory airway disease and bronchoalveolar lavage mastocytosis or neutrophilia using REST software analysis. J. Vet. Intern. Med./Am. Coll. Vet. Intern. Med. 26, 153–161 https://doi.org/10.1111/j.1939-1676.2011.00847.x (2012).
- Iadarola, P. & Viglio, S. Spit it out! How could the sputum proteome aid clinical research into pulmonary diseases? Expert review of proteomics. 14, 391–393 https://doi.org/10.1080/14789450.2017.1317246 (2017).
- Zhu, T. et al. Induced sputum metabolomic profiles and oxidative stress are associated with chronic obstructive pulmonary disease (COPD) severity: potential use for predictive, preventive, and personalized medicine. Epma j. 11, 645–659 https://doi.org/10.1007/s13167-020-00227-w (2020).
- Gharib, S. A. et al. Induced sputum proteome in healthy subjects and asthmatic patients. Journal of Allergy and Clinical Immunology. 128, 1176–11841176 https://doi.org/10.1016/j.jaci.2011.07.053 (2011).
- Bullone, M. & Lavoie, J. P. Asthma "of horses and men"--how can equine heaves help us better understand human asthma immunopathology and its functional consequences? Mol Immunol. 66, 97–105 https://doi.org/10.1016/j.molimm.2014.12.005 (2015).
- Hulliger, M. F. et al. An Integrative miRNA-mRNA Expression Analysis Reveals Striking Transcriptomic Similarities between Severe Equine Asthma and Specific Asthma Endotypes in Humans. Genes (Basel). 11, https://doi.org/10.3390/genes11101143 (2020).
- Persson, C. G., Erjefält, J. S., Korsgren, M. & Sundler, F. The mouse trap. Trends Pharmacol Sci. 18, 465–467 https://doi.org/10.1016/s0165-6147(97)01142-5 (1997).
- Shin, Y. S., Takeda, K. & Gelfand, E. W. Understanding asthma using animal models. Allergy Asthma Immunol Res. 1, 10–18 https://doi.org/10.4168/aair.2009.1.1.10 (2009).
- Karagianni, A. E. et al. Comparative transcriptome analysis of equine alveolar macrophages. Equine veterinary journal. https://doi.org/10.1111/evj.12584 (2016).
- Karagianni, A. E., Kapetanovic, R., McGorum, B. C., Hume, D. A. & Pirie, S. R. The equine alveolar macrophage: Functional and phenotypic comparisons with peritoneal macrophages. Veterinary immunology and immunopathology. 155, 219–228 https://doi.org/10.1016/j.vetimm.2013.07.003 (2013).
- Young, R. et al. Species-Specific Transcriptional Regulation of Genes Involved in Nitric Oxide Production and Arginine Metabolism in Macrophages. ImmunoHorizons 2, 27–37, doi:10.4049/immunohorizons.1700073 (2018).
- Karagianni, A. E., Lisowski, Z. M. & Hume, D. A. & Scott Pirie, R. The equine mononuclear phagocyte system: The relevance of the horse as a model for understanding human innate immunity. Equine Vet J. https://doi.org/10.1111/evj.13341 (2020).
- Woolhouse, I. S., Bayley, D. L. & Stockley, R. A. Effect of sputum processing with dithiothreitol on the detection of inflammatory mediators in chronic bronchitis and bronchiectasis. Thorax. 57, 667–671 https://doi.org/10.1136/thorax.57.8.667 (2002).
- Huang da, W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 4, 44–57 https://doi.org/10.1038/nprot.2008.211 (2009).
- Kramer, A., Green, J., Pollard, J. Jr. & Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 30, 523–530 https://doi.org/10.1093/bioinformatics/btt703 (2014).
- Beekman, L., Tohver, T., Dardari, R. & Léguillette, R. Evaluation of suitable reference genes for gene expression studies in bronchoalveolar lavage cells from horses with inflammatory airway disease. BMC Mol Biol. 12, 5 https://doi.org/10.1186/1471-2199-12-5 (2011).
- Bright, L. A. et al. Modeling the pasture-associated severe equine asthma bronchoalveolar lavage fluid proteome identifies molecular events mediating neutrophilic airway inflammation. Vet Med (Auckl). 10, 43–63 https://doi.org/10.2147/vmrr.s194427 (2019).
- Sikkeland, L. I., Kongerud, J., Stangeland, A. M., Haug, T. & Alexis, N. E. Macrophage enrichment from induced sputum. Thorax. 62, 558–559 https://doi.org/10.1136/thx.2006.073544 (2007).
- Kumar, Y. et al. Molecular analysis of serum and bronchoalveolar lavage in a mouse model of influenza reveals markers of disease severity that can be clinically useful in humans. PLoS One. 9, e86912 https://doi.org/10.1371/journal.pone.0086912 (2014).
- Zamay, T. N. et al. Current and Prospective Protein Biomarkers of Lung Cancer. Cancers (Basel). 9, 155 https://doi.org/10.3390/cancers9110155 (2017).
- Bishwal, S. C. et al. Sputum Proteomics Reveals a Shift in Vitamin D-binding Protein and Antimicrobial Protein Axis in Tuberculosis Patients. Sci Rep. 9, 1036–1036 https://doi.org/10.1038/s41598-018-37662-9 (2019).
- Allen, K. J., Tremaine, W. H. & Franklin, S. H. Prevalence of inflammatory airway disease in national hunt horses referred for investigation of poor athletic performance. Equine Vet J Suppl. 529–534 https://doi.org/10.1111/j.2042-3306.2006.tb05599.x (2006).
- Rousseau, K. et al. Muc5b and Muc5ac are the major oligomeric mucins in equine airway mucus. Am J Physiol Lung Cell Mol Physiol. 292, L1396–1404 https://doi.org/10.1152/ajplung.00444.2006 (2007).
- Wickström, C., Davies, J. R., Eriksen, G. V., Veerman, E. C. & Carlstedt, I. MUC5B is a major gel-forming, oligomeric mucin from human salivary gland, respiratory tract and endocervix: identification of glycoforms and C-terminal cleavage. Biochem J. 334 (Pt 3), 685–693 https://doi.org/10.1042/bj3340685 (1998).
- Hovenberg, H. W., Davies, J. R. & Carlstedt, I. Different mucins are produced by the surface epithelium and the submucosa in human trachea: identification of MUC5AC as a major mucin from the goblet cells. Biochem J. 318 (Pt 1), 319–324 https://doi.org/10.1042/bj3180319 (1996).
- Kerrin, A. et al. Differential lower airway dendritic cell patterns may reveal distinct endotypes of RSV bronchiolitis. Thorax. 72, 620–627 https://doi.org/10.1136/thoraxjnl-2015-207358 (2017).
- Kilkenny, C., Browne, W. J., Cuthill, I. C., Emerson, M. & Altman, D. G. Improving Bioscience Research Reporting: The ARRIVE Guidelines for Reporting Animal Research. PLOS Biology. 8, e1000412 https://doi.org/10.1371/journal.pbio.1000412 (2010).
- Kast, J. I. et al. Respiratory syncytial virus infection influences tight junction integrity. Clin Exp Immunol. 190, 351–359 https://doi.org/10.1111/cei.13042 (2017).
- Depecker, M. et al. Bronchoalveolar lavage fluid in Standardbred racehorses: Influence of unilateral/bilateral profiles and cut-off values on lower airway disease diagnosis. Vet. J. 199, 150–156 https://doi.org/10.1016/j.tvjl.2013.10.013 (2014).
- Eaton, S. L. et al. Total protein analysis as a reliable loading control for quantitative fluorescent Western blotting. PloS one. 8, e72457 https://doi.org/10.1371/journal.pone.0072457 (2013).
- Rappsilber, J., Mann, M. & Ishihama, Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nature protocols. 2, 1896–1906 https://doi.org/10.1038/nprot.2007.261 (2007).
- Dwyer, B. J. et al. TWEAK/Fn14 signalling promotes cholangiocarcinoma niche formation and progression. Journal of Hepatology. https://doi.org/10.1016/j.jhep.2020.11.018 (2020).
- Karagianni, A. E. et al. Transcriptional Response of Ovine Lung to Infection with Jaagsiekte Sheep Retrovirus. J Virol. 93, e00876–00819 https://doi.org/10.1128/JVI.00876-19 (2019).
- Ashburner, M. et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 25, 25–29 https://doi.org/10.1038/75556 (2000).
- Ogata, H. et al. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic acids research. 27, 29–34 https://doi.org/10.1093/nar/27.1.29 (1999).
Competing interest reported. One of the authors of this manuscript is an academic editor for Scientific Reports.
Dr T M Wishart , Reader in Molecular Anatomy, The Roslin Institute, Royal (Dick) School of Veterinary Studies, College of Medicine and Veterinary Medicine, University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, SCOTLAND, UK, T.M.Wishart@Roslin.ed.ac.uk
Research explorer: https://www.research.ed.ac.uk/portal/en/persons/thomas-wishart(3d17d4e1-df57-4c19-8d26-2cd1e5a5aaac)/publications.html
Academic Lead - Proteomics and Metabolomics Facility; Roslin Institute
https://www.ed.ac.uk/roslin/facilities-resources/proteomics-and-metabolomics-facility
Co-Head - Translational Biomarker Development; Centre for Dementia Prevention; University of Edinburgh http://centrefordementiaprevention.com/research/translational-research-groups/ None of the co-authors have any conflict of interest