Whole Exome Sequencing reveals new candidate genes in host genomic susceptibility to Respiratory Syncytial Virus Disease

A Genome-Wide Association Study of Respiratory Syncytial Virus Infection Severity in Infants

BACKGROUND

Respiratory syncytial virus (RSV) is a significant cause of infant morbidity and mortality worldwide. Most children experience at least one 1 RSV infection by the age of two 2 years, but not all develop severe disease. However, the understanding of genetic risk factors for severe RSV is incomplete. Consequently, we conducted a genome-wide association study of RSV severity.

METHODS

Disease severity was assessed by the ReSVinet scale, in a cohort of 251 infants aged 1 week to 1 year. Genotyping data were collected from multiple European study sites as part of the RESCEU Consortium. Linear regression models were used to assess the impact of genotype on RSV severity and gene expression as measured by microarray.

RESULTS

While no SNPs reached the genome-wide statistical significance threshold (P < 5 × 10-8), we identified 816 candidate SNPs with a P-value of <1 × 10-4. Functional annotation of candidate SNPs highlighted genes relevant to neutrophil trafficking and cytoskeletal functions, including LSP1 and RAB27A. Moreover, SNPs within the RAB27A locus significantly altered gene expression (false discovery rate, FDR P < .05).

CONCLUSIONS

These findings may provide insights into genetic mechanisms driving severe RSV infection, offering biologically relevant information for future investigations.

Risk Variants in the Exomes of Children With Critical Illness

IMPORTANCE

Diagnostic genetic testing can lead to changes in management in the pediatric intensive care unit. Genetic risk in children with critical illness but nondiagnostic exome sequencing (ES) has not been explored.

OBJECTIVE

To assess the association between loss-of-function (LOF) variants and pediatric critical illness.

DESIGN, SETTING, AND PARTICIPANTS

This genetic association study examined ES first screened for causative variants among 267 children at the Morgan Stanley Children's Hospital of NewYork-Presbyterian, of whom 22 were otherwise healthy with viral respiratory failure; 18 deceased children with bronchiolitis from the Office of the Chief Medical Examiner of New York City, of whom 14 were previously healthy; and 9990 controls from the Institute for Genomic Medicine at Columbia University Irving Medical Center. The ES data were generated between January 1, 2015, and December 31, 2020, and analyzed between January 1, 2017, and September 2, 2022.

EXPOSURE

Critical illness.

MAIN OUTCOMES AND MEASURES

Odds ratios and P values for genes and gene-sets enriched for rare LOF variants and the loss-of-function observed/expected upper bound fraction (LOEUF) score at which cases have a significant enrichment.

RESULTS

This study included 285 children with critical illness (median [range] age, 4.1 [0-18.9] years; 148 [52%] male) and 9990 controls. A total of 228 children (80%) did not receive a genetic diagnosis. After quality control (QC), 231 children harbored excess rare LOF variants in genes with a LOEUF score of 0.680 or less (intolerant genes) (P = 1.0 × 10-5). After QC, 176 children without a diagnosis harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 1.8; 95% CI, 1.3-2.5). After QC, 25 children with viral respiratory failure harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 2.8; 95% CI, 1.1-6.6). A total of 114 undiagnosed children were enriched for de novo LOF variants in genes without a known disease association (observed, 14; expected, 6.8; enrichment, 2.05).

CONCLUSIONS AND RELEVANCE

In this genetic association study, excess LOF variants were observed among critically ill children despite nondiagnostic ES. Variants lay in genes without a known disease association, suggesting future investigation may connect phenotypes to causative genes.

A multi-tissue study of immune gene expression profiling highlights the key role of the nasal epithelium in COVID-19 severity

Coronavirus Disease-19 (COVID-19) symptoms range from mild to severe illness; the cause for this differential response to infection remains unknown. Unravelling the immune mechanisms acting at different levels of the colonization process might be key to understand these differences. We carried out a multi-tissue (nasal, buccal and blood; n = 156) gene expression analysis of immune-related genes from patients affected by different COVID-19 severities, and healthy controls through the nCounter technology. Mild and asymptomatic cases showed a powerful innate antiviral response in nasal epithelium, characterized by activation of interferon (IFN) pathway and downstream cascades, successfully controlling the infection at local level. In contrast, weak macrophage/monocyte driven innate antiviral response and lack of IFN signalling activity were present in severe cases. Consequently, oral mucosa from severe patients showed signals of viral activity, cell arresting and viral dissemination to the lower respiratory tract, which ultimately could explain the exacerbated innate immune response and impaired adaptative immune responses observed at systemic level. Results from saliva transcriptome suggest that the buccal cavity might play a key role in SARS-CoV-2 infection and dissemination in patients with worse prognosis. Co-expression network analysis adds further support to these findings, by detecting modules specifically correlated with severity involved in the abovementioned biological routes; this analysis also provides new candidate genes that might be tested as biomarkers in future studies. We also found tissue specific severity-related signatures mainly represented by genes involved in the innate immune system and cytokine/chemokine signalling. Local immune response could be key to determine the course of the systemic response and thus COVID-19 severity. Our findings provide a framework to investigate severity host gene biomarkers and pathways that might be relevant to diagnosis, prognosis, and therapy.

Respiratory Syncytial Virus-Associated Neurologic Complications in Children: A Systematic Review and Aggregated Case Series

OBJECTIVES

To describe the features and frequency of respiratory syncytial virus (RSV)-associated severe acute neurologic disease in children.

STUDY DESIGN

We performed a systematic review of the literature to identify reports of severe acute neurologic complications associated with acute RSV infection in children aged <15 years (PROSPERO Registration CRD42019125722). Main outcomes included neurologic, clinical, and demographic features of cases and the frequency of disease. We aggregated available case data from the published literature and from the Australian Acute Childhood Encephalitis (ACE) study.

RESULTS

We identified 87 unique studies from 26 countries describing a spectrum of RSV-associated severe acute neurologic syndromes including proven encephalitis, acute encephalopathy, complex seizures, hyponatremic seizures, and immune-mediated disorders. The frequency of RSV infection in acute childhood encephalitis/encephalopathy was 1.2%-6.5%. We aggregated data from 155 individual cases with RSV-associated severe acute neurologic complications; median age was 11.0 months (IQR 2.0-21.5), most were previously healthy (71/104, 68%). Seizure was the most frequently reported neurologic feature (127/150, 85%). RSV was detected in the central nervous system of 12 cases. Most children recovered (81/122, 66%); however, some reports described partial recovery (33/122, 27%) and death (8/122, 7%).

CONCLUSIONS

RSV-associated neurologic complications have been widely reported, but there is substantial heterogeneity in the design and quality of existing studies. The findings from our study have implications for the investigation, management, and prevention of RSV-associated neurologic complications. Further, this systematic review can inform the design of future studies aiming to quantify the burden of childhood RSV-associated neurologic disease.

Clinical and molecular characterization of COVID-19 hospitalized patients

Clinical and molecular characterization by Whole Exome Sequencing (WES) is reported in 35 COVID-19 patients attending the University Hospital in Siena, Italy, from April 7 to May 7, 2020. Eighty percent of patients required respiratory assistance, half of them being on mechanical ventilation. Fiftyone percent had hepatic involvement and hyposmia was ascertained in 3 patients. Searching for common genes by collapsing methods against 150 WES of controls of the Italian population failed to give straightforward statistically significant results with the exception of two genes. This result is not unexpected since we are facing the most challenging common disorder triggered by environmental factors with a strong underlying heritability (50%). The lesson learned from Autism-Spectrum-Disorders prompted us to re-analyse the cohort treating each patient as an independent case, following a Mendelian-like model. We identified for each patient an average of 2.5 pathogenic mutations involved in virus infection susceptibility and pinpointing to one or more rare disorder(s). To our knowledge, this is the first report on WES and COVID-19. Our results suggest a combined model for COVID-19 susceptibility with a number of common susceptibility genes which represent the favorite background in which additional host private mutations may determine disease progression.

Decoding Susceptibility to Respiratory Viral Infections and Asthma Inception in Children

Human Respiratory Syncytial Virus and Human Rhinovirus are the most frequent cause of respiratory tract infections in infants and children and are major triggers of acute viral bronchiolitis, wheezing and asthma exacerbations. Here, we will discuss the application of the powerful tools of systems biology to decode the molecular mechanisms that determine risk for infection and subsequent asthma. An important conceptual advance is the understanding that the innate immune system is governed by a Bow-tie architecture, where diverse input signals converge onto a few core pathways (e.g., IRF7), which in turn generate diverse outputs that orchestrate effector and regulatory functions. Molecular profiling studies in children with severe exacerbations of asthma/wheeze have identified two major immunological phenotypes. The IRF7hi phenotype is characterised by robust upregulation of antiviral response networks, and the IRF7lo phenotype is characterised by upregulation of markers of TGFβ signalling and type 2 inflammation. Similar phenotypes have been identified in infants and children with severe viral bronchiolitis. Notably, genome-wide association studies supported by experimental validation have identified key pathways that increase susceptibility to HRV infection (ORMDL3 and CHDR3) and modulate TGFβ signalling (GSDMB, TGFBR1, and SMAD3). Moreover, functional deficiencies in the activation of type I and III interferon responses are already evident at birth in children at risk of developing febrile lower respiratory tract infections and persistent asthma/wheeze, suggesting that the trajectory to asthma begins at birth or in utero. Finally, exposure to microbes and their products reprograms innate immunity and provides protection from the development of allergies and asthma in children, and therefore microbial products are logical candidates for the primary prevention of asthma.

The MUC6/AP2A2 Locus and Its Relevance to Alzheimer's Disease: A Review

We recently reported evidence of Alzheimer's disease (AD)-linked genetic variation within the mucin 6 (MUC6) gene on chromosome 11p, nearby the adaptor-related protein complex 2 subunit alpha 2 (AP2A2) gene. This locus has interesting features related to human genomics and clinical research. MUC6 gene variants have been reported to potentially influence viral-including herpesvirus-immunity and the gut microbiome. Within the MUC6 gene is a unique variable number of tandem repeat (VNTR) region. We discovered an association between MUC6 VNTR repeat expansion and AD pathologic severity, particularly tau proteinopathy. Here, we review the relevant literature. The AD-linked VNTR polymorphism may also influence AP2A2 gene expression. AP2A2 encodes a polypeptide component of the adaptor protein complex, AP-2, which is involved in clathrin-coated vesicle function and was previously implicated in AD pathogenesis. To provide background information, we describe some key knowledge gaps in AD genetics research. The "missing/hidden heritability problem" of AD is highlighted. Extensive portions of the human genome, including the MUC6 VNTR, have not been thoroughly evaluated due to limitations of existing high-throughput sequencing technology. We present and discuss additional data, along with cautionary considerations, relevant to the hypothesis that MUC6 repeat expansion influences AD pathogenesis.

RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans

There is a growing interest in unraveling gene expression mechanisms leading to viral host invasion and infection progression. Current findings reveal that long non-coding RNAs (lncRNAs) are implicated in the regulation of the immune system by influencing gene expression through a wide range of mechanisms. By mining whole-transcriptome shotgun sequencing (RNA-seq) data using machine learning approaches, we detected two lncRNAs (ENSG00000254680 and ENSG00000273149) that are downregulated in a wide range of viral infections and different cell types, including blood monocluclear cells, umbilical vein endothelial cells, and dermal fibroblasts. The efficiency of these two lncRNAs was positively validated in different viral phenotypic scenarios. These two lncRNAs showed a strong downregulation in virus-infected patients when compared to healthy control transcriptomes, indicating that these biomarkers are promising targets for infection diagnosis. To the best of our knowledge, this is the very first study using host lncRNAs biomarkers for the diagnosis of human viral infections.

A Meta-Analysis of Multiple Whole Blood Gene Expression Data Unveils a Diagnostic Host-Response Transcript Signature for Respiratory Syncytial Virus

Respiratory syncytial virus (RSV) is one of the major causes of acute lower respiratory tract infection worldwide. The absence of a commercial vaccine and the limited success of current therapeutic strategies against RSV make further research necessary. We used a multi-cohort analysis approach to investigate host transcriptomic biomarkers and shed further light on the molecular mechanism underlying RSV-host interactions. We meta-analyzed seven transcriptome microarray studies from the public Gene Expression Omnibus (GEO) repository containing a total of 922 samples, including RSV, healthy controls, coronaviruses, enteroviruses, influenzas, rhinoviruses, and coinfections, from both adult and pediatric patients. We identified > 1500 genes differentially expressed when comparing the transcriptomes of RSV-infected patients against healthy controls. Functional enrichment analysis showed several pathways significantly altered, including immunologic response mediated by RSV infection, pattern recognition receptors, cell cycle, and olfactory signaling. In addition, we identified a minimal 17-transcript host signature specific for RSV infection by comparing transcriptomic profiles against other respiratory viruses. These multi-genic signatures might help to investigate future drug targets against RSV infection.

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.

Ancestry patterns inferred from massive RNA-seq data

There is a growing body of evidence suggesting that patterns of gene expression vary within and between human populations. However, the impact of this variation in human diseases has been poorly explored, in part owing to the lack of a standardized protocol to estimate biogeographical ancestry from gene expression studies. Here we examine several studies that provide new solid evidence indicating that the ancestral background of individuals impacts gene expression patterns. Next, we test a procedure to infer genetic ancestry from RNA-seq data in 25 data sets where information on ethnicity was reported. Genome data of reference continental populations retrieved from The 1000 Genomes Project were used for comparisons. Remarkably, only eight out of 25 data sets passed FastQC default filters. We demonstrate that, for these eight population sets, the ancestral background of donors could be inferred very efficiently, even in data sets including samples with complex patterns of admixture (e.g., American-admixed populations). For most of the gene expression data sets of suboptimal quality, ancestral inference yielded odd patterns. The present study thus brings a cautionary note for gene expression studies highlighting the importance to control for the potential confounding effect of ancestral genetic background.

Personalizing the Management of Pneumonia

Pneumonia is a highly prevalent disease with considerable morbidity and mortality. However, diagnosis and therapy still rely on antiquated methods, leading to the vast overuse of antimicrobials, which carries risks for both society and the individual. Furthermore, outcomes in severe pneumonia remain poor. Genomic techniques have the potential to transform the management of pneumonia through deep characterization of pathogens as well as the host response to infection. This characterization will enable the delivery of selective antimicrobials and immunomodulatory therapy that will help to offset the disorder associated with overexuberant immune responses.

TIPICO IX: report of the 9th interactive infectious disease workshop on infectious diseases and vaccines

The Ninth Interactive Infectious Disease workshop TIPICO was held on November 22–23, 2018, in Santiago de Compostela, Spain. This 2-day academic experience addressed current and topical issues in the field of infectious diseases and vaccination. Summary findings of the meeting include: cervical cancer elimination will be possible in the future, thanks to the implementation of global vaccination action plans in combination with appropriate screening interventions. The introduction of appropriate immunization programs is key to maintain the success of current effective vaccines such as those against meningococcal disease or rotavirus infection. Additionally, reduced dose schedules might improve the efficiency of some vaccines (i.e., PCV13). New vaccines to improve current preventive alternatives are under development (e.g., against tuberculosis or influenza virus), while others to protect against infectious diseases with no current available vaccines (e.g., enterovirus, parechovirus and flaviviruses) need to be developed. Vaccinomics will be fundamental in this process, while infectomics will allow the application of precision medicine. Further research is also required to understand the impact of heterologous vaccine effects. Finally, vaccination requires education at all levels (individuals, community, healthcare professionals) to ensure its success by helping to overcome major barriers such as vaccine hesitancy and false contraindications.

Single Nucleotide Polymorphisms in Chemosensory Pathway Genes GNB3, TAS2R19, and TAS2R38 Are Associated with Chronic Rhinosinusitis

BACKGROUND

Chronic rhinosinusitis (CRS) is a multifaceted disease with a significant genetic component. The importance of taste receptor signaling has recently been highlighted in CRS; single nucleotide polymorphisms (SNPs) of bitter tastant-responsive G-protein-coupled receptors have been linked with CRS and with altered innate immune responses to multiple bacterially derived signals.

OBJECTIVE

To determine in CRS the frequency of six SNPs in genes with known bitter tastant signaling function.

METHODS

Genomic DNA was isolated from 74 CRS volunteers in West Virginia, and allele frequency was determined and compared with demographically matched data from the 1,000 Genomes database.

RESULTS

For two SNPs in a gene recently associated with bitterant signaling regulation, RGS21, there were no associations with CRS (although the frequency of the minor allele of RGS21, rs7528947, was seen to increase with increasing Lund-Mackay CT staging score). Two TAS2R bitter taste receptor gene variants (TAS2R19 rs10772420 and TAS2R38 rs713598), identified in prior CRS genetics studies, were found to have similar associations in this study.

CONCLUSION

Unique to our study is the establishment of an association between CRS in this patient population and GNB3 SNP rs5443, a variation in an established G protein component downstream of bitterant receptor signal transduction.

A genome-wide association study of tramadol metabolism from post-mortem samples

Phase I tramadol metabolism requires cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) to form O-desmethyltramadol (M1). CYP2D6 genetic variants may infer metabolizer phenotype; however, drug ADME (absorption, distribution, metabolism, and excretion) and response depend on protein pathway(s), not CYP2D6 alone. There is a paucity of data regarding the contribution of trans-acting proteins to idiosyncratic phenotypes following drug exposure. A genome-wide association study identified five markers (rs79983226/kgp11274252, rs9384825, rs62435418/kgp10370907, rs72732317/kgp3743668, and rs184199168/exm1592932) associated with the conversion of tramadol to M1 (M1:T). These SNPs reside within five genes previously implicated with adverse reactions. Analysis of accompanying toxicological meta-data revealed a significant positive linear relationship between M1:T and degree of sample polypharmacy. Taken together, these data identify candidate loci for potential clinical inferences of phenotype following exposure to tramadol and highlight sample polypharmacy as a possible diagnostic covariate in post-mortem genetic studies.

The natural selection that shapes our genomes

Most of the variation in the human genome (∼95%) is constrained, directly or indirectly, by purifying selection and GC-biased gene conversion, according to a recent article by Pouyet et al. (2018). The use of 'non-neutral' variation to infer human demographies can lead to undesirable biases; for example, in estimation of the time of the most recent common ancestor. Further examination of 'neutral' variation in entire human genomes from The 1000 Genomes Project reveals that ∼99% of this variation lacks exonic function, but ∼35% of it falls in introns. In addition, estimates of biogeographical ancestry using 'non-neutral' SNPs differ very marginally from inferences obtained from 'neutral' variation. Additional investigations should be carried out before establishing the roadmap for future human population and forensic genetic studies.

RSV infection and respiratory sequelae

INTRODUCTION

The association between respiratory syncytial virus (RSV) infections and long-term respiratory sequelae has long been recognized. It is estimated that individuals with a history of RSV bronchiolitis have 2- to 12-fold higher risk of developing asthma. Although this risk tends to decrease with age, persistent airway obstruction and hyperresponsiveness are observed even 30 years after RSV infection.

EVIDENCE ACQUISITION

Our data search strategy was designed to address the following questions: What is the epidemiological evidence available on the association between RSV infection and long-term respiratory morbidity? What are the potential pathogenic pathways linking RSV infection to long-term respiratory morbidity? Are there any host genetic backgrounds that can predispose to both severe RSV lower respiratory tract infection and asthma? Are antiviral therapies and RSV prevention measures effective in reducing respiratory morbidities?

EVIDENCE SYNTHESIS

This article reviews the recent scientific literature on the epidemiological association and pathogenic links between early RSV infection and long-term respiratory morbidities.

CONCLUSIONS

Nowadays, asthma is increasingly considered a heterogeneous disease, caused by interactions between several host and environmental factors. Understanding the specific causative role of respiratory viruses, and the pathogenic mechanisms through which bronchiolitis predisposes to asthma, is a challenging, but essential starting point for the development of prevention and treatment strategies potentially capable of preserving lung function.

Respiratory syncytial virus

Respiratory syncytial virus (RSV) is the most common cause of infant hospitalization and causes a high burden of disease in the elderly, too. This enveloped negative-stranded RNA virus has been recently reclassified in the Pneumoviridae family. Infections of the respiratory cells happens when the two major surface glycoproteins, G and F, take contact with the cell receptor CX3CR1 and mediate entry by fusion, respectively. Viral mRNA transcription, genomic RNA synthesis and nucleocapsid formation occur in large cytoplasmic inclusion bodies to avoid recognition by the host innate immune response. Most progeny virions remain associated to the infected cell surface; fusion of infected with adjacent cells results in the formation of large multinucleated syncytia that eventually undergo apoptosis. Desquamated epithelial cells form the plugs that with mucus and fibrin may cause lower airway obstructions. Pathogenetic mechanism of severe RSV disease likely involve both the extent of viral replication and the host immune response. Regarding the latter, single nucleotide polymorphism analysis and genome-wide association studies showed that genetic susceptibility to severe RSV infection is likely a complex trait, in which many different host genetic variants contribute. Recent studies pointed to the fact that bronchiolitis severity depends more on the specific infecting RSV genotypes than on the amount of viral loads. A population-based surveillance system to better define RSV burden of disease would be of valuable help for implementing future vaccination programs.

A 2-transcript host cell signature distinguishes viral from bacterial diarrhea and it is influenced by the severity of symptoms

Recently, a biomarker signature consisting of 2-transcript host RNAs was proposed for discriminating bacterial from viral infections in febrile children. We evaluated the performance of this signature in a different disease scenario, namely a cohort of Mexican children (n = 174) suffering from acute diarrhea of different infectious etiologies. We first examined the admixed background of the patients, indicating that most of them have a predominantly Native American genetic ancestry with a variable amount of European background (ranging from 0% to 57%). The results confirm that the RNA test can discriminate between viral and bacterial causes of infection (t-test; P-value = 6.94×10⁻¹¹; AUC = 80%; sensitivity: 68% [95% CI: 55%–79%]; specificity: 84% [95% CI: 78%–90%]), but the strength of the signal differs substantially depending on the causal pathogen, with the stronger signal being that of Shigella (P-value = 3.14 × 10⁻¹²; AUC = 89; sensitivity: 70% [95% CI: 57%–83%]; specificity: 100% [95% CI: 100%–100%]). The accuracy of this test improves significantly when excluding mild cases (P-value = 2.13 × 10⁻⁶; AUC = 85%; sensitivity: 79% [95% CI: 58%–95%]; specificity: 78% [95% CI: 65%–88%]). The results broaden the scope of previous studies by incorporating different pathogens, variable levels of disease severity, and different ancestral background of patients, and add confirmatory support to the clinical utility of these 2-transcript biomarkers.

Whole Exome Sequencing Identifies New Host Genomic Susceptibility Factors in Empyema Caused by Streptococcus pneumoniae in Children: A Pilot Study

Pneumonia is the leading cause of death amongst infectious diseases. Streptococcus pneumoniae is responsible for about 25% of pneumonia cases worldwide, and it is a major cause of childhood mortality. We carried out a whole exome sequencing (WES) study in eight patients with complicated cases of pneumococcal pneumonia (empyema). An initial assessment of statistical association of WES variation with pneumonia was carried out using data from the 1000 Genomes Project (1000G) for the Iberian Peninsula (IBS) as reference controls. Pseudo-replication statistical analyses were carried out using different European control groups. Association tests pointed to single nucleotide polymorphism (SNP) rs201967957 (gene MEIS1; chromosome 2; p-value_IBS = 3.71 × 10^-13) and rs576099063 (gene TSPAN15; chromosome 10; p-value_IBS = 2.36 × 10^-8) as the best candidate variants associated to pneumococcal pneumonia. A burden gene test of pathogenicity signaled four genes, namely, OR9G9, MUC6, MUC3A and APOB, which carry significantly increased pathogenic variation when compared to controls. By analyzing various transcriptomic data repositories, we found strong supportive evidence for the role of MEIS1, TSPAN15 and APOBR (encoding the receptor of the APOB protein) in pneumonia in mouse and human models. Furthermore, the association of the olfactory receptor gene OR9G9 has recently been related to some viral infectious diseases, while the role of mucin genes (MUC6 and MUC3A), encoding mucin glycoproteins, are well-known factors related to chronic obstructive airway disease. WES emerges as a promising technique to disentangle the genetic basis of host genome susceptibility to infectious respiratory diseases.