Identification of Novel Genes and Biological Pathways That Overlap in Infectious and Nonallergic Diseases of the Upper and Lower Airways Using Network Analyses

Previous genetic studies on susceptibility to otitis media and airway infections have focused on immune pathways acting within the local mucosal epithelium, and outside of allergic rhinitis and asthma, limited studies exist on the overlaps at the gene, pathway or network level between the upper and lower airways. In this report, we compared [1] pathways identified from network analysis using genes derived from published genome-wide family-based and association studies for otitis media, sinusitis, and lung phenotypes, to [2] pathways identified using differentially expressed genes from RNA-sequence data from lower airway, sinus, and middle ear tissues, in particular cholesteatoma tissue compared to middle ear mucosa. For otitis media, a large number of genes (n = 1,806) were identified as differentially expressed between cholesteatoma and middle ear mucosa, which in turn led to the identification of 68 pathways that are enriched in cholesteatoma. Two differentially expressed genes CR1 and SAA1 overlap in middle ear, sinus, and lower airway samples and are potentially novel genes for otitis media susceptibility. In addition, 56 genes were differentially expressed in both tissues from the middle ear and either sinus or lower airways. Pathways that are common in upper and lower airway diseases, whether from published DNA studies or from our RNA-sequencing analyses, include chromatin organization/remodeling, endocytosis, immune system process, protein folding, and viral process. Taken together, our findings from genetic susceptibility and differential tissue expression studies support the hypothesis that the unified airway theory wherein the upper and lower respiratory tracts act as an integrated unit also applies to infectious and nonallergic airway epithelial disease. Our results may be used as reference for identification of genes or pathways that are relevant to upper and lower airways, whether common across sites, or unique to each disease.

[1]  Wenbin Liu,et al.  Long noncoding RNAs show differential expression profiles and display ceRNA potential in cholesteatoma pathogenesis , 2018, Oncology reports.

[2]  Levi Waldron,et al.  Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples , 2012, Genome Biology.

[3]  Nadezhda T. Doncheva,et al.  Topological analysis and interactive visualization of biological networks and protein structures , 2012, Nature Protocols.

[4]  T. Matsunaga,et al.  Immunohistochemical studies of complement receptor (CR1) in cases with normal sinus mucosa and chronic sinusitis , 2004, Archives of oto-rhino-laryngology.

[5]  Ivana V. Yang,et al.  FUT2 Variants Confer Susceptibility to Familial Otitis Media. , 2018, American journal of human genetics.

[6]  Brian Bushnell,et al.  BBMerge – Accurate paired shotgun read merging via overlap , 2017, PloS one.

[7]  W. Ahmad,et al.  Nasal polyps. , 2018, British medical journal.

[8]  A. Matsuda,et al.  Distinct gene expression profiles and regulation networks of nasal polyps in eosinophilic and non‐eosinophilic chronic rhinosinusitis , 2018, International forum of allergy & rhinology.

[9]  Bjarni V. Halldórsson,et al.  HLA class II sequence variants influence tuberculosis risk in populations of European ancestry , 2016, Nature Genetics.

[10]  Bjarni V. Halldórsson,et al.  A loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis , 2019, Nature Genetics.

[11]  E. Silverman,et al.  RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD , 2019, Respiratory Research.

[12]  Othman Soufan,et al.  NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis , 2019, Nucleic Acids Res..

[13]  Rob Patro,et al.  Salmon provides fast and bias-aware quantification of transcript expression , 2017, Nature Methods.

[14]  Gao T. Wang,et al.  Rare A2ML1 variants confer susceptibility to otitis media , 2015, Nature Genetics.

[15]  Kristin E. Burnum-Johnson,et al.  Application of multiplexed ion mobility spectrometry towards the identification of host protein signatures of treatment effect in pulmonary tuberculosis , 2018, Tuberculosis.

[16]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[17]  Y. Bossé,et al.  Novel genes for airway wall thickness identified with combined genome-wide association and expression analyses. , 2015, American journal of respiratory and critical care medicine.

[18]  G. Cheng,et al.  Transcriptome Analysis of Bronchoalveolar Lavage Fluid From Children With Mycoplasma pneumoniae Pneumonia Reveals Natural Killer and T Cell-Proliferation Responses , 2018, Front. Immunol..

[19]  S. Fujieda,et al.  Expression and Functional Analysis of CST1 in Intractable Nasal Polyps , 2018, American journal of respiratory cell and molecular biology.

[20]  S. Kinnamon,et al.  Rna Sequencing and Pathway Analysis Identify Tumor Necrosis Factor Alpha Driven Small Proline-Rich Protein Dysregulation in Chronic Rhinosinusitis , 2017, American journal of rhinology & allergy.

[21]  COPDGene and Eclipse Investigators,et al.  Genetic susceptibility for chronic bronchitis in chronic obstructive pulmonary disease , 2014 .

[22]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[23]  Fidel Ramírez,et al.  Computing topological parameters of biological networks , 2008, Bioinform..

[24]  Alissa S. Hanshew,et al.  Integrating the microbiota of the respiratory tract with the unified airway model. , 2017, Respiratory medicine.

[25]  D. Briles,et al.  Complement Receptor 1 Expression on Mouse Erythrocytes Mediates Clearance of Streptococcus pneumoniae by Immune Adherence , 2010, Infection and Immunity.

[26]  M. Robinson,et al.  Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences , 2015, F1000Research.

[27]  Kyle Bittinger,et al.  Topographical continuity of bacterial populations in the healthy human respiratory tract. , 2011, American journal of respiratory and critical care medicine.

[28]  Mimi P Macias,et al.  Differential Gene Expression in Cholesteatoma by DNA Chip Analysis , 2013, The Laryngoscope.

[29]  A. Rosli,et al.  Bioinformatics analysis of differentially expressed genes in liver cancer for identification of key genes and pathways , 2018 .

[30]  Raúl San José Estépar,et al.  A Genome-Wide Association Study of Emphysema and Airway Quantitative Imaging Phenotypes. , 2015, American journal of respiratory and critical care medicine.

[31]  Anne W. Ndungu,et al.  Polymorphism in a lincRNA Associates with a Doubled Risk of Pneumococcal Bacteremia in Kenyan Children , 2016, American journal of human genetics.

[32]  C. Carenfelt,et al.  CR1-expression and C3b-mediated phagocytosis of granulocytes in purulent maxillary secretion and peripheral blood from patients with sinusitis. , 1989, Acta Oto-Laryngologica.

[33]  H. Tettelin,et al.  Whole‐Exome Sequencing Identifies the 6q12‐q16 Linkage Region and a Candidate Gene, TTK, for Pulmonary Nontuberculous Mycobacterial Disease , 2017, American journal of respiratory and critical care medicine.

[34]  The Virulence Function of Streptococcus pneumoniae Surface Protein A Involves Inhibition of Complement Activation and Impairment of Complement Receptor-Mediated Protection1 , 2004, The Journal of Immunology.

[35]  M. Fraczek,et al.  A microarray study of gene expression profiles in nasal polyps. , 2011, Auris, nasus, larynx.

[36]  H. Jung,et al.  Analysis of gene expression profiles of normal human nasal mucosa and nasal polyp tissues by SAGE. , 2006, The Journal of allergy and clinical immunology.

[37]  G. Cheng,et al.  Transcriptome analysis of bronchoalveolar lavage fluid from children with severe Mycoplasma pneumoniae pneumonia reveals novel gene expression and immunodeficiency , 2017, Human Genomics.

[38]  H. Jung,et al.  Analysis of gene expression profiles in cholesteatoma using oligonucleotide microarray , 2006, Acta oto-laryngologica.

[39]  A. Salas,et al.  Whole Exome Sequencing reveals new candidate genes in host genomic susceptibility to Respiratory Syncytial Virus Disease , 2017, Scientific Reports.

[40]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[41]  Astrid Gall,et al.  Ensembl 2018 , 2017, Nucleic Acids Res..

[42]  J. Negri,et al.  Microarray Analysis of Distinct Gene Transcription Profiles in Non-eosinophilic Chronic Sinusitis with Nasal Polyps , 2008, American journal of rhinology.

[43]  Anushya Muruganujan,et al.  PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools , 2018, Nucleic Acids Res..

[44]  S. Fujieda,et al.  Periostin as a novel biomarker for postoperative recurrence of chronic rhinosinitis with nasal polyps , 2018, Scientific Reports.

[45]  V. Ramakrishnan,et al.  Upper and lower airways associations in patients with chronic rhinosinusitis and bronchiectasis , 2013, International forum of allergy & rhinology.

[46]  Q. Hamid,et al.  Similar allergic inflammation in the middle ear and the upper airway: evidence linking otitis media with effusion to the united airways concept. , 2004, The Journal of allergy and clinical immunology.

[47]  Weiqing Wang,et al.  Transcriptome Analysis Reveals Distinct Gene Expression Profiles in Eosinophilic and Noneosinophilic Chronic Rhinosinusitis with Nasal Polyps , 2016, Scientific Reports.

[48]  E. Burchard,et al.  Functional genomics of CDHR3 confirms its role in HRV-C infection and childhood asthma exacerbations. , 2019, The Journal of allergy and clinical immunology.

[49]  David Haussler,et al.  UCSC Data Integrator and Variant Annotation Integrator , 2016, Bioinform..

[50]  R. Wunderink,et al.  Multidimensional Assessment of the Host Response in Mechanically Ventilated Patients with Suspected Pneumonia. , 2019, American journal of respiratory and critical care medicine.

[51]  I. Wang,et al.  Gene expression profiles in human nasal polyp tissues studied by means of DNA microarray. , 2004, The Journal of allergy and clinical immunology.

[52]  Ralf Hofestädt,et al.  Identification of Novel Cholesteatoma-Related Gene Expression Signatures Using Full-Genome Microarrays , 2012, PloS one.

[53]  Xinchun Chen,et al.  Genome-wide association study identifies two risk loci for tuberculosis in Han Chinese , 2018, Nature Communications.

[54]  Karin Breuer,et al.  InnateDB: systems biology of innate immunity and beyond—recent updates and continuing curation , 2012, Nucleic Acids Res..

[55]  K. Bromek,et al.  Structure of the C3b Binding Site of CR1 (CD35), the Immune Adherence Receptor , 2002, Cell.

[56]  P. Burgel United airway diseases. Should we add upper airway inflammatory disorders to the list of chronic obstructive pulmonary disease comorbidities? , 2015, Annals of the American Thoracic Society.

[57]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[58]  J. Krouse The unified airway--conceptual framework. , 2008, Otolaryngologic clinics of North America.

[59]  Robert E. W. Hancock,et al.  NetworkAnalyst - integrative approaches for protein–protein interaction network analysis and visual exploration , 2014, Nucleic Acids Res..

[60]  A. Uitterlinden,et al.  Genome-wide association study for acute otitis media in children identifies FNDC1 as disease contributing gene , 2016, Nature Communications.

[61]  Jianguo Xia,et al.  NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data , 2015, Nature Protocols.

[62]  Katherine S. Pollard,et al.  Integrating host response and unbiased microbe detection for lower respiratory tract infection diagnosis in critically ill adults , 2018, Proceedings of the National Academy of Sciences.

[63]  E. Silverman,et al.  Susceptibility to Childhood Pneumonia: A Genome‐Wide Analysis , 2017, American journal of respiratory cell and molecular biology.

[64]  M. Decramer,et al.  Original article: Sinonasal pathology in nonallergic asthma and COPD: ‘united airway disease’ beyond the scope of allergy , 2007, Allergy. European Journal of Allergy and Clinical Immunology.

[65]  F. Gudat,et al.  Glomerulonephritis in a patient with complement factor I deficiency. , 1999, American journal of kidney diseases : the official journal of the National Kidney Foundation.

[66]  J. Barrett,et al.  Susceptibility to tuberculosis is associated with variants in the ASAP1 gene encoding a regulator of dendritic cell migration , 2015, Nature Genetics.

[67]  D. Weeks,et al.  A Genome-Wide Association Study of Chronic Otitis Media with Effusion and Recurrent Otitis Media Identifies a Novel Susceptibility Locus on Chromosome 2 , 2013, Journal of the Association for Research in Otolaryngology.

[68]  Joseph K. Pickrell,et al.  Detection and interpretation of shared genetic influences on 42 human traits , 2015, Nature Genetics.

[69]  Feng-jun Wang,et al.  Identification of key genes and pathways in chronic rhinosinusitis with nasal polyps using bioinformatics analysis. , 2019, American journal of otolaryngology.

[70]  S. Jesic,et al.  Perimatrix of middle ear cholesteatoma: A granulation tissue with a specific transcriptomic signature , 2020, The Laryngoscope.

[71]  R. Simon,et al.  The united allergic airway: Connections between allergic rhinitis, asthma, and chronic sinusitis , 2012, American journal of rhinology & allergy.

[72]  E. Arellano,et al.  C-Reactive Protein and Serum Amyloid A Overexpression in Lung Tissues of Chronic Obstructive Pulmonary Disease Patients: A Case-Control Study , 2013, International journal of medical sciences.

[73]  R. Metson,et al.  Gene Expression Profiling of Nasal Polyps Associated With Chronic Sinusitis and Aspirin‐Sensitive Asthma , 2008, The Laryngoscope.

[74]  Thomas Blicher,et al.  A genome-wide association study identifies CDHR3 as a susceptibility locus for early childhood asthma with severe exacerbations , 2013, Nature Genetics.

[75]  P. Giavina-Bianchi,et al.  United airway disease: current perspectives , 2016, Journal of asthma and allergy.

[76]  M. Segal,et al.  Biomarkers of Tuberculosis Severity and Treatment Effect: A Directed Screen of 70 Host Markers in a Randomized Clinical Trial , 2017, EBioMedicine.

[77]  A. Salas,et al.  Whole Exome Sequencing Identifies New Host Genomic Susceptibility Factors in Empyema Caused by Streptococcus pneumoniae in Children: A Pilot Study , 2018, Genes.

[78]  Stuart M. Brown,et al.  Gene expression profiles of bronchoalveolar cells in pulmonary TB. , 2008, Tuberculosis.

[79]  N. Eriksson,et al.  Genome-wide association and HLA region fine-mapping studies identify susceptibility loci for multiple common infections , 2016, Nature Communications.

[80]  J. Kere,et al.  Genome-wide association analysis reveals variants on chromosome 19 that contribute to childhood risk of chronic otitis media with effusion , 2016, Scientific Reports.