Bulk RNA sequencing of human pediatric lung cell populations reveals unique transcriptomic signature associated with postnatal pulmonary development

Postnatal lung development results in an increasingly functional organ prepared for gas exchange and pathogenic challenges. It is achieved through cellular differentiation and migration. Changes in the tissue architecture during this development process are well-documented and increasing cellular diversity associated with it are reported in recent years. Despite recent progress, transcriptomic and molecular pathways associated with human postnatal lung development are yet to be fully understood. In this study, we investigated gene expression patterns associated with healthy pediatric lung development in four major enriched cell populations (epithelial, endothelial, and nonendothelial mesenchymal cells, along with lung leukocytes) from 1-day-old to 8-yr-old organ donors with no known lung disease. For analysis, we considered the donors in four age groups [less than 30 days old neonates, 30 days to < 1 yr old infants, toddlers (1 to < 2 yr), and children 2 yr and older] and assessed differentially expressed genes (DEG). We found increasing age-associated transcriptional changes in all four major cell types in pediatric lung. Transition from neonate to infant stage showed highest number of DEG compared with the number of DEG found during infant to toddler- or toddler to older children-transitions. Profiles of differential gene expression and further pathway enrichment analyses indicate functional epithelial cell maturation and increased capability of antigen presentation and chemokine-mediated communication. Our study provides a comprehensive reference of gene expression patterns during healthy pediatric lung development that will be useful in identifying and understanding aberrant gene expression patterns associated with early life respiratory diseases.NEW & NOTEWORTHY This study presents postnatal transcriptomic changes in major cell populations in human lung, namely endothelial, epithelial, mesenchymal cells, and leukocytes. Although human postnatal lung development continues through early adulthood, our results demonstrate that greatest transcriptional changes occur in first few months of life during neonate to infant transition. These early transcriptional changes in lung parenchyma are particularly notable for functional maturation and activation of alveolar type II cell genes.

Single-Cell Transcriptomic Profiling Identifies Molecular Phenotypes of Newborn Human Lung Cells

While animal model studies have extensively defined the mechanisms controlling cell diversity in the developing mammalian lung, there exists a significant knowledge gap with regards to late-stage human lung development. The NHLBI Molecular Atlas of Lung Development Program (LungMAP) seeks to fill this gap by creating a structural, cellular and molecular atlas of the human and mouse lung. Transcriptomic profiling at the single-cell level created a cellular atlas of newborn human lungs. Frozen single-cell isolates obtained from two newborn human lungs from the LungMAP Human Tissue Core Biorepository, were captured, and library preparation was completed on the Chromium 10X system. Data was analyzed in Seurat, and cellular annotation was performed using the ToppGene functional analysis tool. Transcriptional interrogation of 5500 newborn human lung cells identified distinct clusters representing multiple populations of epithelial, endothelial, fibroblasts, pericytes, smooth muscle, immune cells and their gene signatures. Computational integration of data from newborn human cells and with 32,000 cells from postnatal days 1 through 10 mouse lungs generated by the LungMAP Cincinnati Research Center facilitated the identification of distinct cellular lineages among all the major cell types. Integration of the newborn human and mouse cellular transcriptomes also demonstrated cell type-specific differences in maturation states of newborn human lung cells. Specifically, newborn human lung matrix fibroblasts could be separated into those representative of younger cells (n = 393), or older cells (n = 158). Cells with each molecular profile were spatially resolved within newborn human lung tissue. This is the first comprehensive molecular map of the cellular landscape of neonatal human lung, including biomarkers for cells at distinct states of maturity.

Single-cell resolution of human airway epithelial cells exposed to bronchiolitis obliterans-associated chemicals

Bronchiolitis obliterans (BO) is a fibrotic lung disease characterized by progressive luminal narrowing and obliteration of the small airways. In the nontransplant population, inhalation exposure to certain chemicals is associated with BO; however, the mechanisms contributing to disease induction remain poorly understood. This study's objective was to use single-cell RNA sequencing for the identification of transcriptomic signatures common to primary human airway epithelial cells after chemical exposure to BO-associated chemicals-diacetyl or nitrogen mustard-to help explain BO induction. Primary airway epithelial cells were cultured at air-liquid interface and exposed to diacetyl, nitrogen mustard, or control vapors. Cultures were dissociated and sequenced for single-cell RNA. Differential gene expression and functional pathway analyses were compared across exposures. In total, 75,663 single cells were captured and sequenced from all exposure conditions. Unbiased clustering identified 11 discrete phenotypes, including 5 basal, 2 ciliated, and 2 secretory cell clusters. With chemical exposure, the proportion of cells assigned to keratin 5+ basal cells decreased, whereas the proportion of cells aligned to secretory cell clusters increased compared with control exposures. Functional pathway analysis identified interferon signaling and antigen processing/presentation as pathways commonly upregulated after diacetyl or nitrogen mustard exposure in a ciliated cell cluster. Conversely, the response of airway basal cells differed significantly with upregulation of the unfolded protein response in diacetyl-exposed basal cells, not seen in nitrogen mustard-exposed cultures. These new insights provide early identification of airway epithelial signatures common to BO-associated chemical exposures.NEW & NOTEWORTHY Bronchiolitis obliterans (BO) is a devastating fibrotic lung disease of the small airways, or bronchioles. This original manuscript uses single-cell RNA sequencing for identifying common signatures of chemically exposed airway epithelial cells in BO induction. Chemical exposure reduced the proportion of keratin 5+ basal cells while increasing the proportion of keratin 4+ suprabasal cells. Functional pathways contributory to these shifts differed significantly across exposures. These new results highlight similarities and differences in BO induction across exposures.

New insights into the natural history of bronchopulmonary dysplasia from proteomics and multiplexed immunohistochemistry

Bronchopulmonary dysplasia (BPD) is a disease of prematurity related to the arrest of normal lung development. The objective of this study was to better understand how proteome modulation and cell-type shifts are noted in BPD pathology. Pediatric human donors aged 1-3 yr were classified based on history of prematurity and histopathology consistent with "healed" BPD (hBPD, n = 3) and "established" BPD (eBPD, n = 3) compared with respective full-term born (n = 6) age-matched term controls. Proteins were quantified by tandem mass spectroscopy with selected Western blot validations. Multiplexed immunofluorescence (MxIF) microscopy was performed on lung sections to enumerate cell types. Protein abundances and MxIF cell frequencies were compared among groups using ANOVA. Cell type and ontology enrichment were performed using an in-house tool and/or EnrichR. Proteomics detected 5,746 unique proteins, 186 upregulated and 534 downregulated, in eBPD versus control with fewer proteins differentially abundant in hBPD as compared with age-matched term controls. Cell-type enrichment suggested a loss of alveolar type I, alveolar type II, endothelial/capillary, and lymphatics, and an increase in smooth muscle and fibroblasts consistent with MxIF. Histochemistry and Western analysis also supported predictions of upregulated ferroptosis in eBPD versus control. Finally, several extracellular matrix components mapping to angiogenesis signaling pathways were altered in eBPD. Despite clear parsing by protein abundance, comparative MxIF analysis confirms phenotypic variability in BPD. This work provides the first demonstration of tandem mass spectrometry and multiplexed molecular analysis of human lung tissue for critical elucidation of BPD trajectory-defining factors into early childhood.NEW & NOTEWORTHY We provide new insights into the natural history of bronchopulmonary dysplasia in donor human lungs after the neonatal intensive care unit hospitalization. This study provides new insights into how the proteome and histopathology of BPD changes in early childhood, uncovering novel pathways for future study.

Building Career Paths for Ph.D., Basic and Translational Scientists in Clinical Departments in the United States: An Official American Thoracic Society Workshop Report

Rationale: To identify barriers and opportunities for Ph.D., basic and translational scientists to be fully integrated into clinical units. Objectives: In 2022, an ad hoc committee of the American Thoracic Society developed a project proposal and workshop to identify opportunities and barriers for scientists who do not practice medicine to develop successful careers and achieve tenure-track faculty positions in clinical departments and divisions within academic medical centers (AMCs) in the United States. Methods: This document focuses on results from a survey of adult and pediatric pulmonary, critical care, and sleep medicine division chiefs as well as a survey of workshop participants, including faculty in departmental and school leadership roles in both basic science and clinical units within U.S. AMCs. Results: We conclude that full integration of non-clinically practicing basic and translational scientists into the clinical units, in addition to their traditional placements in basic science units, best serves the tripartite mission of AMCs to provide care, perform research, and educate the next generation. Evidence suggests clinical units do employ Ph.D. scientists in large numbers, but these faculty are often hired into non-tenure track positions, which do not provide the salary support, start-up funds, research independence, or space often associated with hiring in basic science units within the same institution. These barriers to success of Ph.D. faculty in clinical units are largely financial. Conclusions: Our recommendation is for AMCs to consider and explore some of our proposed strategies to accomplish the goal of integrating basic and translational scientists into clinical units in a meaningful way.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

A Trisomy 21 Lung Cell Atlas

Trisomy 21 (T21), resulting in Down Syndrome (DS), is the most prevalent chromosomal abnormality worldwide. While pulmonary disease is a major cause of morbidity and mortality in DS, the ontogeny of pulmonary complications remains poorly understood. We recently demonstrated that T21 lung anomalies, including airway branching and vascular lymphatic abnormalities, are initiated in utero. Here, we aimed to describe molecular changes at the single cell level in prenatal T21 lungs. Our results demonstrate differences in the proportion of cell populations and detail changes in gene expression at the time of initiation of histopathological abnormalities. Notably, we identify shifts in the distribution of alveolar epithelial progenitors, widespread induction of key extracellular matrix molecules in mesenchymal cells and hyper-activation of IFN signaling in endothelial cells. This single cell atlas of T21 lungs greatly expands our understanding of antecedents to pulmonary complications and should facilitate efforts to mitigate respiratory disease in DS.

Preparation of noninfectious scRNAseq samples from SARS-CoV-2-infected epithelial cells

Coronavirus disease (COVID-19) is an infectious disease caused by the SARS coronavirus 2 (SARS-CoV-2) virus. Direct assessment, detection, and quantitative analysis using high throughput methods like single-cell RNA sequencing (scRNAseq) is imperative to understanding the host response to SARS-CoV-2. One barrier to studying SARS-CoV-2 in the laboratory setting is the requirement to process virus-infected cell cultures, and potentially infectious materials derived therefrom, under Biosafety Level 3 (BSL-3) containment. However, there are only 190 BSL3 laboratory facilities registered with the U.S. Federal Select Agent Program, as of 2020, and only a subset of these are outfitted with the equipment needed to perform high-throughput molecular assays. Here, we describe a method for preparing non-hazardous RNA samples from SARS-CoV-2 infected cells, that enables scRNAseq analyses to be conducted safely in a BSL2 facility-thereby making molecular assays of SARS-CoV-2 cells accessible to a much larger community of researchers. Briefly, we infected African green monkey kidney epithelial cells (Vero-E6) with SARS-CoV-2 for 96 hours, trypsin-dissociated the cells, and inactivated them with methanol-acetone in a single-cell suspension. Fixed cells were tested for the presence of infectious SARS-CoV-2 virions using the Tissue Culture Infectious Dose Assay (TCID50), and also tested for viability using flow cytometry. We then tested the dissociation and methanol-acetone inactivation method on primary human lung epithelial cells that had been differentiated on an air-liquid interface. Finally, we performed scRNAseq quality control analysis on the resulting cell populations to evaluate the effects of our virus inactivation and sample preparation protocol on the quality of the cDNA produced. We found that methanol-acetone inactivated SARS-CoV-2, fixed the lung epithelial cells, and could be used to obtain noninfectious, high-quality cDNA libraries. This methodology makes investigating SARS-CoV-2, and related high-containment RNA viruses at a single-cell level more accessible to an expanded community of researchers.

Use of Host Response to Refine the Diagnosis of Group A Streptococcal Pharyngitis

BACKGROUND

Current diagnostic tests for pharyngitis do not distinguish between symptomatic group A Streptococcus (GAS) infection and asymptomatic colonization, resulting in over-diagnosis and unnecessary use of antibiotics. We assessed whether measures of host response could make this distinction.

METHODS

We enrolled 18 children with pharyngitis having Centor scores of 4 or 5 and 21 controls without pharyngitis or other acute infections. Both groups had throat cultures, molecular tests for GAS and respiratory viruses and IgM serology for Epstein-Barr virus. Host response was evaluated with white blood cell count (WBC), C-reactive protein (CRP), procalcitonin (PCT), and sequencing of RNA from peripheral blood leukocytes.

RESULTS

Of 18 cases, 11 had GAS pharyngitis, 3 had adenovirus pharyngitis and 4 had other pharyngitis. Among asymptomatic controls, 5 were positive for GAS. WBC, CRP, and PCT were higher in subjects with pharyngitis compared to asymptomatic controls including those with GAS. Transcriptional profiles from children with symptomatic GAS were clearly distinct from those of children in all other groups. The levels of two genes, CD177 and TLR5 each individually accurately distinguished between symptomatic and asymptomatic GAS. Optimal diagnostic sensitivity and specificity were achieved by the combination of CRP and PCT, and by each of the two gene markers.

CONCLUSION

In this exploratory study, we showed that traditional measures of inflammation and markers of host gene expression distinguish between symptomatic and asymptomatic GAS. These results point to future rapid molecular approaches for improving the diagnosis of GAS pharyngitis, that may help reduce unnecessary antibiotic use.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection.. [PMID: 36172120 PMCID: PMC9516852 DOI: 10.1101/2022.09.22.509123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots were identified for DVG recombination and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single cell RNA-seq analysis indicated the IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection.

Proteomic Analysis of Human Lung Development

Rationale: The current understanding of human lung development derives mostly from animal studies. Although transcript-level studies have analyzed human donor tissue to identify genes expressed during normal human lung development, protein-level analysis that would enable the generation of new hypotheses on the processes involved in pulmonary development are lacking. Objectives: To define the temporal dynamic of protein expression during human lung development. Methods: We performed proteomics analysis of human lungs at 10 distinct times from birth to 8 years to identify the molecular networks mediating postnatal lung maturation. Measurements and Main Results: We identified 8,938 proteins providing a comprehensive view of the developing human lung proteome. The analysis of the data supports the existence of distinct molecular substages of alveolar development and predicted the age of independent human lung samples, and extensive remodeling of the lung proteome occurred during postnatal development. Evidence of post-transcriptional control was identified in early postnatal development. An extensive extracellular matrix remodeling was supported by changes in the proteome during alveologenesis. The concept of maturation of the immune system as an inherent part of normal lung development was substantiated by flow cytometry and transcriptomics. Conclusions: This study provides the first in-depth characterization of the human lung proteome during development, providing a unique proteomic resource freely accessible at Lungmap.net. The data support the extensive remodeling of the lung proteome during development, the existence of molecular substages of alveologenesis, and evidence of post-transcriptional control in early postnatal development.

Reduced Notch1 Cleavage Promotes the Development of Pulmonary Hypertension

Clinical trials of Dll4 (Delta-like 4) neutralizing antibodies (Dll4nAbs) in cancer patients are ongoing. Surprisingly, pulmonary hypertension (PH) occurs in 14% to 18% of patients treated with Dll4nAbs, but the mechanisms have not been studied. Here, PH progression was measured in mice treated with Dll4nAbs. We detected Notch signaling in lung tissues and analyzed pulmonary vascular permeability and inflammation. Notch target gene array was performed on adult human pulmonary microvascular endothelial cells (ECs) after inhibiting Notch cleavage. Similar mechanisms were studied in PH mouse models and pulmonary arterial hypertension patients. The rescue effects of constitutively activated Notch1 in vivo were also measured. We observed that Dll4nAbs induced PH in mice as indicated by significantly increased right ventricular systolic pressure, as well as pulmonary vascular and right ventricular remodeling. Mechanistically, Dll4nAbs inhibited Notch1 cleavage and subsequently impaired lung endothelial barrier function and increased immune cell infiltration in vessel walls. In vitro, Notch targeted genes' expression related to cell growth and inflammation was decreased in human pulmonary microvascular ECs after the Notch1 inactivation. In lungs of PH mouse models and pulmonary arterial hypertension patients, Notch1 cleavage was inhibited. Consistently, EC cell-cell junction was leaky, and immune cell infiltration increased in PH mouse models. Overexpression activated Notch1-attenuated progression of PH in mice. In conclusion, Dll4nAbs led to PH development in mice by impaired EC barrier function and increased immune cell infiltration through inhibition of Notch1 cleavage in lung ECs. Reduced Notch1 cleavage in lung ECs could be an underlying mechanism of PH pathogenesis.

A systems genomics approach uncovers molecular associates of RSV severity

Respiratory syncytial virus (RSV) infection results in millions of hospitalizations and thousands of deaths each year. Variations in the adaptive and innate immune response appear to be associated with RSV severity. To investigate the host response to RSV infection in infants, we performed a systems-level study of RSV pathophysiology, incorporating high-throughput measurements of the peripheral innate and adaptive immune systems and the airway epithelium and microbiota. We implemented a novel multi-omic data integration method based on multilayered principal component analysis, penalized regression, and feature weight back-propagation, which enabled us to identify cellular pathways associated with RSV severity. In both airway and immune cells, we found an association between RSV severity and activation of pathways controlling Th17 and acute phase response signaling, as well as inhibition of B cell receptor signaling. Dysregulation of both the humoral and mucosal response to RSV may play a critical role in determining illness severity.

This paper presents a novel approach to understanding the localized molecular responses to respiratory syncytial virus (RSV) and the system-level correlates of clinical outcomes. To do this, we developed a novel statistical method able to integrate high dimensional molecular data characterizing the host airway microbota and immune and nasal gene expression. We show that this integrative approach facilitates superior performance in estimating clinical outcome as opposed to any single data type. Using this approach, we identified both cell type-specific and shared biomarkers and regulatory pathways associated with RSV severity. Specifically, we identified an association between RSV severity, activation of pathways controlling Th17, and inhibition of B cell receptor signaling, which were present in both the site of infection airway and in peripheral immune cells. These results can guide future efforts to identify biomarkers for identifying or predicting illness severity following infant RSV infection. They may also be useful as biomarkers to inform the efficacy of future interventions (e.g., therapies) or preventative measures to suppress the rate of severe disease (e.g., vaccines).

Gene Expression Risk Scores for COVID-19 Illness Severity

BACKGROUND

The correlates of coronavirus disease 2019 (COVID-19) illness severity following infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are incompletely understood.

METHODS

We assessed peripheral blood gene expression in 53 adults with confirmed SARS-CoV-2 infection clinically adjudicated as having mild, moderate, or severe disease. Supervised principal components analysis was used to build a weighted gene expression risk score (WGERS) to discriminate between severe and nonsevere COVID-19.

RESULTS

Gene expression patterns in participants with mild and moderate illness were similar, but significantly different from severe illness. When comparing severe versus nonsevere illness, we identified >4000 genes differentially expressed (false discovery rate < 0.05). Biological pathways increased in severe COVID-19 were associated with platelet activation and coagulation, and those significantly decreased with T-cell signaling and differentiation. A WGERS based on 18 genes distinguished severe illness in our training cohort (cross-validated receiver operating characteristic-area under the curve [ROC-AUC] = 0.98), and need for intensive care in an independent cohort (ROC-AUC = 0.85). Dichotomizing the WGERS yielded 100% sensitivity and 85% specificity for classifying severe illness in our training cohort, and 84% sensitivity and 74% specificity for defining the need for intensive care in the validation cohort.

CONCLUSIONS

These data suggest that gene expression classifiers may provide clinical utility as predictors of COVID-19 illness severity.

Lung disease manifestations in Down syndrome

Down syndrome (DS) is one of the most prevalent chromosomal abnormalities worldwide, affecting 1 in 700 live births. Although multiple organ systems are affected by the chromosomal defects, respiratory failure and lung disease are the leading causes of morbidity and mortality observed in DS. Manifestations of DS in the respiratory system encompass the entire lung starting from the nasopharynx to the trachea/upper airways to the lower airways and alveolar spaces, as well as vascular and lymphatic defects. Most of our knowledge on respiratory illness in persons with DS arises from pediatric studies; however, many of these disorders present early in infancy, supporting developmental mechanisms. In this review, we will focus on the different lung phenotypes in DS, as well as the genetic and molecular pathways that may be contributing to these complications during development.

The effect of air pollution on the transcriptomics of the immune response to respiratory infection

Combustion related particulate matter air pollution (PM) is associated with an increased risk of respiratory infections in adults. The exact mechanism underlying this association has not been determined. We hypothesized that increased concentrations of combustion related PM would result in dysregulation of the innate immune system. This epidemiological study includes 111 adult patients hospitalized with respiratory infections who underwent transcriptional analysis of their peripheral blood. We examined the association between gene expression at the time of hospitalization and ambient measurements of particulate air pollutants in the 28 days prior to hospitalization. For each pollutant and time lag, gene-specific linear models adjusting for infection type were fit using LIMMA (Linear Models For Microarray Data), and pathway/gene set analyses were performed using the CAMERA (Correlation Adjusted Mean Rank) program. Comparing patients with viral and/or bacterial infection, the expression patterns associated with air pollution exposure differed. Adjusting for the type of infection, increased concentrations of Delta-C (a marker of biomass smoke) and other PM were associated with upregulation of iron homeostasis and protein folding. Increased concentrations of black carbon (BC) were associated with upregulation of viral related gene pathways and downregulation of pathways related to antigen presentation. The pollutant/pathway associations differed by lag time and by type of infection. This study suggests that the effect of air pollution on the pathogenesis of respiratory infection may be pollutant, timing, and infection specific.

Gene Expression Risk Scores for COVID-19 Illness Severity

BACKGROUND

The correlates of COVID-19 illness severity following infection with SARS-Coronavirus 2 (SARS-CoV-2) are incompletely understood.

METHODS

We assessed peripheral blood gene expression in 53 adults with confirmed SARS-CoV-2-infection clinically adjudicated as having mild, moderate or severe disease. Supervised principal components analysis was used to build a weighted gene expression risk score (WGERS) to discriminate between severe and non-severe COVID.

RESULTS

Gene expression patterns in participants with mild and moderate illness were similar, but significantly different from severe illness. When comparing severe versus non-severe illness, we identified >4000 genes differentially expressed (FDR<0.05). Biological pathways increased in severe COVID-19 were associated with platelet activation and coagulation, and those significantly decreased with T cell signaling and differentiation. A WGERS based on 18 genes distinguished severe illness in our training cohort (cross-validated ROC-AUC=0.98), and need for intensive care in an independent cohort (ROC-AUC=0.85). Dichotomizing the WGERS yielded 100% sensitivity and 85% specificity for classifying severe illness in our training cohort, and 84% sensitivity and 74% specificity for defining the need for intensive care in the validation cohort.

CONCLUSION

These data suggest that gene expression classifiers may provide clinical utility as predictors of COVID-19 illness severity.

Serotonin system genes contribute to the susceptibility to obesity in Black adolescents

OBJECTIVE

The importance of the central and peripheral serotonin systems in regulating energy balance and obesity development has been highlighted in animal models. Yet, the role of both serotonin systems has not been systematically assessed in humans. The purpose of this study was to investigate the association of genes within both serotonin systems with obesity outcomes in black adolescents.

METHODS

African-American adolescents (n = 1052) whose mothers participated the Memphis New Mother's Study were assessed. In total, 110 polymorphisms mapped to 10 serotonin genes were examined for their associations with standardized body mass index (BMI-z) scores and waist circumferences using generalized estimating equation models.

RESULTS

Over 39% of adolescents were overweight or had obesity. Three single nucleotide polymorphisms (SNPs) within TPH2, HTR3B, and SLC6A4, were significantly associated with BMI-z scores (p < 1.7 × 10-3). Two SNPs in TPH2 were nominally associated with waist circumferences. One SNP in HTR2C was associated with BMI-z scores (p = 0.001) and waist circumferences (p = 0.005) only in girls. Tissue-specific expression indicates that three identified genes are predominantly expressed in the brain.

CONCLUSION

The central serotonin system may play a key role in obesity development in black adolescents. Future studies are warranted to explore additional serotonin system genes and their potential obesogenic mechanisms in humans.

Prenatal histological, cellular, and molecular anomalies in trisomy 21 lung

Down syndrome (DS), also known as trisomy 21 (T21), is the most common human chromosomal anomaly. Although DS can affect many organ systems, lung and heart disease are the leading causes of death. An abundance of existing data suggests that lung abnormalities originate postnatally in DS. However, a single report of branching insufficiency in DS has inferred a potential prenatal origin. The histology of T21 fetal lungs (n = 15) was assessed by an experienced pathologist. Spatial differences in cellular phenotypes were examined using immunohistochemistry (IHC). Comprehensive gene expression in prenatal T21 lungs (n = 19), and age-matched controls (n = 19), was performed using high-throughput RNA sequencing (RNAseq) and validated by RT-qPCR. Histopathological abnormalities were observed in approximately half of T21 prenatal lung samples analyzed, which included dilated terminal airways/acinar tubules, dilated lymphatics, and arterial wall thickening. IHC for Ki67 revealed significant reductions in epithelial and mesenchymal cell proliferation, predominantly in tissues displaying pathology. IHC demonstrated that airway smooth muscle was reduced and discontinuous in the proximal airway in conjunction with reduced SOX2. RNAseq identified 118 genes significantly dysregulated (FDR < 0.05) in T21 lung when unadjusted and 316 genes when adjusted for age. Ontology analysis showed that IFN pathway genes were appreciably upregulated, whereas complement and coagulation cascades and extracellular matrix pathway genes were downregulated. RT-qPCR confirmed the changes in genes associated with these pathways in prenatal T21 lungs. Our data demonstrate that specific histological, cellular, and molecular abnormalities occur prenatally in different compartments of human T21 lung, which could be representative of premature stage progression. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

The critical role of collagen VI in lung development and chronic lung disease

Type VI collagen (collagen VI) is an obligate extracellular matrix component found mainly in the basement membrane region of many mammalian tissues and organs, including skeletal muscle and throughout the respiratory system. Collagen VI is probably most recognized in medicine as the genetic cause of a spectrum of muscular dystrophies, including Ullrich Congenital Myopathy and Bethlem Myopathy. Collagen VI is thought to contribute to myopathy, at least in part, by mediating muscle fiber integrity by anchoring myoblasts to the muscle basement membrane. Interestingly, collagen VI myopathies present with restrictive respiratory insufficiency, thought to be due primarily to thoracic muscular weakening. Although it was recently recognized as one of the (if not the) most abundant collagens in the mammalian lung, there is a substantive knowledge gap concerning its role in respiratory system development and function. A few studies have suggested that collagen VI insufficiency is associated with airway epithelial cell survival and altered lung function. Our recent work suggested collagen VI may be a genomic risk factor for chronic lung disease in premature infants. Using this as motivation, we thoroughly assessed the role of collagen VI in lung development and in lung epithelial cell biology. Here, we describe the state-of-the-art for collagen VI cell and developmental biology within the respiratory system, and reveal its essential roles in normal developmental processes and airway epithelial cell phenotype and intracellular signaling.

Repetitive diacetyl vapor exposure promotes ubiquitin proteasome stress and precedes bronchiolitis obliterans pathology

Bronchiolitis obliterans (BO) is a devastating lung disease seen commonly after lung transplant, following severe respiratory tract infection or chemical inhalation exposure. Diacetyl (DA; 2,3-butanedione) is a highly reactive alpha-diketone known to cause BO when inhaled, however, the mechanisms of how inhalation exposure leads to BO development remains poorly understood. In the current work, we combined two clinically relevant models for studying the pathogenesis of DA-induced BO: (1) an in vivo rat model of repetitive DA vapor exposures with recovery and (2) an in vitro model of primary human airway epithelial cells exposed to pure DA vapors. Rats exposed to 5 consecutive days 200 parts-per-million DA 6 h per day had worsening survival, persistent hypoxemia, poor weight gain, and histologic evidence of BO 14 days after DA exposure cessation. At the end of exposure, increased expression of the ubiquitin stress protein ubiquitin-C accumulated within DA-exposed rat lung homogenates and localized primarily to the airway epithelium, the primary site of BO development. Lung proteasome activity increased concurrently with ubiquitin-C expression after DA exposure, supportive of significant proteasome stress. In primary human airway cultures, global proteomics identified 519 significantly modified proteins in DA-exposed samples relative to controls with common pathways of the ubiquitin proteasome system, endosomal reticulum transport, and response to unfolded protein pathways being upregulated and cell-cell adhesion and oxidation-reduction pathways being downregulated. Collectively, these two models suggest that diacetyl inhalation exposure causes abundant protein damage and subsequent ubiquitin proteasome stress prior to the development of chemical-induced BO pathology.

Airway Gene Expression Correlates of Respiratory Syncytial Virus Disease Severity and Microbiome Composition in Infants

BACKGROUND

Respiratory syncytial virus (RSV) is the leading cause of severe respiratory disease in infants. The causes and correlates of severe illness in the majority of infants are poorly defined.

METHODS

We recruited a cohort of RSV-infected infants and simultaneously assayed the molecular status of their airways and the presence of airway microbiota. We used rigorous statistical approaches to identify gene expression patterns associated with disease severity and microbiota composition, separately and in combination.

RESULTS

We measured comprehensive airway gene expression patterns in 106 infants with primary RSV infection. We identified an airway gene expression signature of severe illness dominated by excessive chemokine expression. We also found an association between Haemophilus influenzae, disease severity, and airway lymphocyte accumulation. Exploring the time of onset of clinical symptoms revealed acute activation of interferon signaling following RSV infection in infants with mild or moderate illness, which was absent in subjects with severe illness.

CONCLUSIONS

Our data reveal that airway gene expression patterns distinguish mild/moderate from severe illness. Furthermore, our data identify biomarkers that may be therapeutic targets or useful for measuring efficacy of intervention responses.

Temporal Dysbiosis of Infant Nasal Microbiota Relative to Respiratory Syncytial Virus Infection

BACKGROUND

Respiratory syncytial virus (RSV) is a leading cause of infant respiratory disease. Infant airway microbiota has been associated with respiratory disease risk and severity. The extent to which interactions between RSV and microbiota occur in the airway, and their impact on respiratory disease susceptibility and severity, are unknown.

METHODS

We carried out 16S rRNA microbiota profiling of infants in the first year of life from (1) a cross-sectional cohort of 89 RSV-infected infants sampled during illness and 102 matched healthy controls, and (2) a matched longitudinal cohort of 12 infants who developed RSV infection and 12 who did not, sampled before, during, and after infection.

RESULTS

We identified 12 taxa significantly associated with RSV infection. All 12 taxa were differentially abundant during infection, with 8 associated with disease severity. Nasal microbiota composition was more discriminative of healthy vs infected than of disease severity.

CONCLUSIONS

Our findings elucidate the chronology of nasal microbiota dysbiosis and suggest an altered developmental trajectory associated with RSV infection. Microbial temporal dynamics reveal indicators of disease risk, correlates of illness and severity, and impact of RSV infection on microbiota composition.

Airway gene-expression classifiers for respiratory syncytial virus (RSV) disease severity in infants

Background

A substantial number of infants infected with RSV develop severe symptoms requiring hospitalization. We currently lack accurate biomarkers that are associated with severe illness.

Method

We defined airway gene expression profiles based on RNA sequencing from nasal brush samples from 106 full-tem previously healthy RSV infected subjects during acute infection (day 1–10 of illness) and convalescence stage (day 28 of illness). All subjects were assigned a clinical illness severity score (GRSS). Using AIC-based model selection, we built a sparse linear correlate of GRSS based on 41 genes (NGSS1). We also built an alternate model based upon 13 genes associated with severe infection acutely but displaying stable expression over time (NGSS2).

Results

NGSS1 is strongly correlated with the disease severity, demonstrating a naïve correlation (ρ) of ρ = 0.935 and cross-validated correlation of 0.813. As a binary classifier (mild versus severe), NGSS1 correctly classifies disease severity in 89.6% of the subjects following cross-validation. NGSS2 has slightly less, but comparable, accuracy with a cross-validated correlation of 0.741 and classification accuracy of 84.0%.

Conclusion

Airway gene expression patterns, obtained following a minimally-invasive procedure, have potential utility for development of clinically useful biomarkers that correlate with disease severity in primary RSV infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-00913-2.

Lymphocyte-Specific Biomarkers Associated With Preterm Birth and Bronchopulmonary Dysplasia

Many premature babies who are born with neonatal respiratory distress syndrome (RDS) go on to develop Bronchopulmonary Dysplasia (BPD) and later Post-Prematurity Respiratory Disease (PRD) at one year corrected age, characterized by persistent or recurrent lower respiratory tract symptoms frequently related to inflammation and viral infection. Transcriptomic profiles were generated from sorted peripheral blood CD8+ T cells of preterm and full-term infants enrolled with consent in the NHLBI Prematurity and Respiratory Outcomes Program (PROP) at the University of Rochester and the University at Buffalo. We identified outcome-related gene expression patterns following standard methods to identify markers for oxygen utilization and BPD as outcomes in extremely premature infants. We further identified predictor gene sets for BPD based on transcriptomic data adjusted for gestational age at birth (GAB). RNA-Seq analysis was completed for CD8+ T cells from 145 subjects. Among the subjects with highest risk for BPD (born at <29 weeks gestational age (GA); n=72), 501 genes were associated with oxygen utilization. In the same set of subjects, 571 genes were differentially expressed in subjects with a diagnosis of BPD and 105 genes were different in BPD subjects as defined by physiologic challenge. A set of 92 genes could predict BPD with a moderately high degree of accuracy. We consistently observed dysregulation of TGFB, NRF2, HIPPO, and CD40-associated pathways in BPD. Using gene expression data from both premature and full-term subjects (n=116), we identified a 28 gene set that predicted the PRD status with a moderately high level of accuracy, which also were involved in TGFB signaling. Transcriptomic data from sort-purified peripheral blood CD8+ T cells from 145 preterm and full-term infants identified sets of molecular markers of inflammation associated with independent development of BPD in extremely premature infants at high risk for the disease and of PRD among the preterm and full-term subjects.

Rare Pulmonary Connective Tissue Type Mast Cells Regulate Lung Endothelial Cell Angiogenesis

Within the human lung, mast cells typically reside adjacent to the conducting airway and assume a mucosal phenotype (MC_T). In rare pathologic conditions, connective tissue phenotype mast cells (MC_TCs) can be found in the lung parenchyma. MC_TCs accumulate in the lungs of infants with severe bronchopulmonary dysplasia, a chronic lung disease associated with preterm birth, which is characterized by pulmonary vascular dysmorphia. The human mast cell line (LUVA) was used to model MC_TCs or MC_Ts. The ability of MC_TCs to affect vascular organization during fetal lung development was tested in mouse lung explant cultures. The effect of MC_TCs on in vitro tube formation and barrier function was studied using primary fetal human pulmonary microvascular endothelial cells. The mechanistic role of MC_TC proteases was tested using inhibitors. MC_TCLUVA but not MC_TLUVA was associated with vascular dysmorphia in lung explants. In vitro, the addition of MC_TCLUVA potentiated fetal human pulmonary microvascular endothelial cell interactions, inhibited tube stability, and disrupted endothelial cell junctions. Protease inhibitors ameliorated the ability of MC_TCLUVA to alter endothelial cell angiogenic activities in vitro and ex vivo. These data indicate that MC_TCs may directly contribute to disrupted angiogenesis in bronchopulmonary dysplasia. A better understanding of factors that regulate mast cell subtype and their different effector functions is essential.

Implementing sequence-based antigenic distance calculation into immunological shape space model

BACKGROUND

In 2009, a novel influenza vaccine was distributed worldwide to combat the H1N1 influenza "swine flu" pandemic. However, antibodies induced by the vaccine display differences in their specificity and cross-reactivity dependent on pre-existing immunity. Here, we present a computational model that can capture the effect of pre-existing immunity on influenza vaccine responses. The model predicts the region of the virus hemagglutinin (HA) protein targeted by antibodies after vaccination as well as the level of cross-reactivity induced by the vaccine. We tested our model by simulating a scenario similar to the 2009 pandemic vaccine and compared the results to antibody binding data obtained from human subjects vaccinated with the monovalent 2009 H1N1 influenza vaccine.

RESULTS

We found that both specificity and cross-reactivity of the antibodies induced by the 2009 H1N1 influenza HA protein were affected by the viral strain the individual was originally exposed. Specifically, the level of antigenic relatedness between the original exposure HA antigen and the 2009 HA protein affected antigenic-site immunodominance. Moreover, antibody cross-reactivity was increased when the individual's pre-existing immunity was specific to an HA protein antigenically distinct from the 2009 pandemic strain. Comparison of simulation data with antibody binding data from human serum samples demonstrated qualitative and quantitative similarities between the model and real-life immune responses to the 2009 vaccine.

CONCLUSION

We provide a novel method to evaluate expected outcomes in antibody specificity and cross-reactivity after influenza vaccination in individuals with different influenza HA antigen exposure histories. The model produced similar outcomes as what has been previously reported in humans after receiving the 2009 influenza pandemic vaccine. Our results suggest that differences in cross-reactivity after influenza vaccination should be expected in individuals with different exposure histories.

Airway basal cell injury after acute diacetyl (2,3-butanedione) vapor exposure

RATIONALE

Diacetyl (DA; 2,3-butanedione) is a chemical found commonly in foods and e-cigarettes. When inhaled, DA causes epithelial injury, though the mechanism of repair remain poorly understood. The objective of this study was to evaluate airway basal cell repair after DA vapor exposure.

METHODS

Primary human bronchial epithelial cells were exposed to DA or PBS for 1 h. Lactate dehydrogenase, cleaved caspase 3/7 and trans-epithelial electrical resistance were measured prior to and following exposure. Exposed cultures were analyzed for the airway basal cell markers keratin 5 and p63 as well as ubiquitin and proteasome activity. Cultures were also treated with a proteasome inhibitor (MG132).

RESULTS

DA vapor exposure caused a transient decrease in trans-epithelial electrical resistance in all DA-exposed cultures. Supernatant lactate dehydrogenase and cleaved caspase 3/7 increased significantly at the highest DA concentration but not at lower DA concentrations. Increased keratin 5 ubiquitination occurred after DA exposure but resolved by day 3. Damage to airway basal cells persisted at day 3 in the presence of MG132.

CONCLUSIONS

Diacetyl exposure results in airway basal cell injury with keratin 5 ubiquitination and decreased p63 expression. The ubiquitin-proteasome-pathway partially mediates airway basal cell repair after acute DA exposure.

Transcriptional characterisation of human lung cells identifies novel mesenchymal lineage markers

RATIONALE

The lung mesenchyme gives rise to multiple distinct lineages of cells in the mature respiratory system, including smooth muscle cells of the airway and vasculature. However, a thorough understanding of the specification and mesenchymal cell diversity in the human lung is lacking.

METHODS

We completed single-cell RNA sequencing analysis of fetal human lung tissues. Canonical correlation analysis, clustering, cluster marker gene identification and t-distributed stochastic neighbour embedding representation was performed in Seurat. Cell populations were annotated using ToppFun. Immunohistochemistry and in situ hybridisation were used to validate spatiotemporal gene expression patterns for key marker genes.

RESULTS

We identified molecularly distinct populations representing "committed" fetal human lung endothelial cells, pericytes and smooth muscle cells. Early endothelial lineages expressed "classic" endothelial cell markers (platelet endothelial cell adhesion molecule/CD31 and claudin 5), while pericytes expressed platelet-derived growth factor receptor-β, Thy-1 membrane glycoprotein and basement membrane molecules (collagen IV, laminin and proteoglycans). We observed a large population of "nonspecific" human lung mesenchymal progenitor cells characterised by expression of collagen I and multiple elastin fibre genes (ELN, MFAP2 and FBN1). We closely characterised the diversity of mesenchymal lineages defined by α₂-smooth muscle actin (ACTA2) expression. Two cell populations, with the highest levels of ACTA2 transcriptional activity, expressed unique sets of markers associated with airway or vascular smooth muscle cells. Spatiotemporal analysis of these marker genes confirmed early and persistent spatial specification of airway (HHIP, MYLK and IGF1) and vascular (NTRK3 and MEF2C) smooth muscle cells in the developing human lung.

CONCLUSION

Our data suggest that specification of distinct airway and vascular smooth muscle cell phenotypes is established early in development and can be identified using the markers we provide.

Host Gene Expression in Nose and Blood for the Diagnosis of Viral Respiratory Infection

BACKGROUND

Recently there has been a growing interest in the potential for host transcriptomic analysis to augment the diagnosis of infectious diseases.

METHODS

We compared nasal and blood samples for evaluation of the host transcriptomic response in children with acute respiratory syncytial virus (RSV) infection, symptomatic non-RSV respiratory virus infection, asymptomatic rhinovirus infection, and virus-negative asymptomatic controls. We used nested leave-one-pair-out cross-validation and supervised principal components analysis to define small sets of genes whose expression patterns accurately classified subjects. We validated gene classification scores using an external data set.

RESULTS

Despite lower quality of nasal RNA, the number of genes detected by microarray in each sample type was equivalent. Nasal gene expression signal derived mainly from epithelial cells but also included a variable leukocyte contribution. The number of genes with increased expression in virus-infected children was comparable in nasal and blood samples, while nasal samples also had decreased expression of many genes associated with ciliary function and assembly. Nasal gene expression signatures were as good or better for discriminating between symptomatic, asymptomatic, and uninfected children.

CONCLSUSIONS

Our results support the use of nasal samples to augment pathogen-based tests to diagnose viral respiratory infection.

Collagen VI Deficiency Results in Structural Abnormalities in the Mouse Lung

Collagen VI (COL6) is known for its role in a spectrum of congenital muscular dystrophies, which are often accompanied by respiratory dysfunction. However, little is known regarding the function of COL6 in the lung. We confirmed the presence of COL6 throughout the basement membrane region of mouse lung tissue. Lung structure and organization were studied in a previously described Col6a1^-/- mouse, which does not produce detectable COL6 in the lung. The Col6a1^-/- mouse displayed histopathologic alveolar and airway abnormalities. The airspaces of Col6a1^-/- lungs appeared simplified, with larger (29%; P < 0.01) and fewer (31%; P < 0.001) alveoli. These airspace abnormalities included reduced isolectin B4⁺ alveolar capillaries and surfactant protein C-positive alveolar epithelial type-II cells. Alterations in lung function consistent with these histopathologic changes were evident. Col6a1^-/- mice also displayed multiple airway changes, including increased branching (59%; P < 0.001), increased mucosal thickness (34%; P < 0.001), and increased epithelial cell density (13%; P < 0.001). Comprehensive transcriptome analysis revealed that the loss of COL6 is associated with reductions in integrin-paxillin-phosphatidylinositol 3-kinase signaling in vivo. In vitro, COL6 promoted steady-state phosphorylated paxillin levels and reduced cell density (16% to 28%; P < 0.05) at confluence. Inhibition of phosphatidylinositol 3-kinase, or its downstream effectors, resulted in increased cell density to a level similar to that seen on matrices lacking COL6.

Virus-Specific Antibody, Viral Load, and Disease Severity in Respiratory Syncytial Virus Infection

Background

Maternally derived serum antibody and viral load are thought to influence disease severity in primary respiratory syncytial virus (RSV) infection. As part of the AsPIRES study of RSV pathogenesis, we correlated various serum antibody concentrations and viral load with disease severity.

Methods

Serum neutralizing antibody titers and levels of immunoglobulin G (IgG) to RSV fusion protein (F), attachment proteins of RSV group A and B, the CX3C region of G, and nasal viral load were measured in 139 full-term previously healthy infants with primary RSV infection and correlated with illness severity.

Results

Univariate analysis showed no relationship between measures of serum antibody and severity. However, a multivariate model adjusting for age at the time of infection found a significant 0.56 decrease in severity score for each 2-fold increase in antibody concentration to RSV F. The benefit of antibody was greatest in infants ≤ 2 months of age. Additionally, estimated antibody titer at birth was correlated with age at infection, suggesting that higher antibody titers delay infection. Viral load did not differ by illness severity.

Conclusion

Our data support the concept of maternal immunization with an RSV vaccine during pregnancy as a strategy for reducing the burden of RSV infection in full-term healthy infants exposed to RSV during their first winter.

Mutation of Respiratory Syncytial Virus G Protein's CX3C Motif Attenuates Infection in Cotton Rats and Primary Human Airway Epithelial Cells

Despite being a high priority for vaccine development, no vaccine is yet available for respiratory syncytial virus (RSV). A live virus vaccine is the primary type of vaccine being developed for young children. In this report, we describe our studies of infected cotton rats and primary human airway epithelial cells (pHAECs) using an RSV r19F with a mutation in the CX3C chemokine motif in the RSV G protein (CX4C). Through this CX3C motif, RSV binds to the corresponding chemokine receptor, CX3CR1, and this binding contributes to RSV infection of pHAECs and virus induced host responses that contribute to disease. In both the cotton rat and pHAECs, the CX4C mutation decreased virus replication and disease and/or host responses to infection. Thus, this mutation, or other mutations that block binding to CX3CR1, has the potential to improve a live attenuated RSV vaccine by attenuating both infection and disease pathogenesis.

Integration of transcriptomic and proteomic data identifies biological functions in cell populations from human infant lung

Systems biology uses computational approaches to integrate diverse data types to understand cell and organ behavior. Data derived from complementary technologies, for example transcriptomic and proteomic analyses, are providing new insights into development and disease. We compared mRNA and protein profiles from purified endothelial, epithelial, immune, and mesenchymal cells from normal human infant lung tissue. Signatures for each cell type were identified and compared at both mRNA and protein levels. Cell-specific biological processes and pathways were predicted by analysis of concordant and discordant RNA-protein pairs. Cell clustering and gene set enrichment comparisons identified shared versus unique processes associated with transcriptomic and/or proteomic data. Clear cell-cell correlations between mRNA and protein data were obtained from each cell type. Approximately 40% of RNA-protein pairs were coherently expressed. While the correlation between RNA and their protein products was relatively low (Spearman rank coefficient r_s ~0.4), cell-specific signature genes involved in functional processes characteristic of each cell type were more highly correlated with their protein products. Consistency of cell-specific RNA-protein signatures indicated an essential framework for the function of each cell type. Visualization and reutilization of the protein and RNA profiles are supported by a new web application, "LungProteomics," which is freely accessible to the public.

Aims, Study Design, and Enrollment Results From the Assessing Predictors of Infant Respiratory Syncytial Virus Effects and Severity Study

Background

The majority of infants hospitalized with primary respiratory syncytial virus (RSV) infection have no obvious risk factors for severe disease.

Objective

The aim of this study (Assessing Predictors of Infant RSV Effects and Severity, AsPIRES) was to identify factors associated with severe disease in full-term healthy infants younger than 10 months with primary RSV infection.

Methods

RSV infected infants were enrolled from 3 cohorts during consecutive winters from August 2012 to April 2016 in Rochester, New York. A birth cohort was prospectively enrolled and followed through their first winter for development of RSV infection. An outpatient supplemental cohort was enrolled in the emergency department or pediatric offices, and a hospital cohort was enrolled on admission with RSV infection. RSV was diagnosed by reverse transcriptase-polymerase chain reaction. Demographic and clinical data were recorded and samples collected for assays: buccal swab (cytomegalovirus polymerase chain reaction, PCR), nasal swab (RSV qualitative PCR, complete viral gene sequence, 16S ribosomal ribonucleic acid [RNA] amplicon microbiota analysis), nasal wash (chemokine and cytokine assays), nasal brush (nasal respiratory epithelial cell gene expression using RNA sequencing [RNAseq]), and 2 to 3 ml of heparinized blood (flow cytometry, RNAseq analysis of purified cluster of differentiation [CD]4+, CD8+, B cells and natural killer cells, and RSV-specific antibody). Cord blood (RSV-specific antibody) was also collected for the birth cohort. Univariate and multivariate logistic regression will be used for analysis of data using a continuous Global Respiratory Severity Score (GRSS) as the outcome variable. Novel statistical methods will be developed for integration of the large complex datasets.

Results

A total of 453 infants were enrolled into the 3 cohorts; 226 in the birth cohort, 60 in the supplemental cohort, and 78 in the hospital cohort. A total of 126 birth cohort infants remained in the study and were evaluated for 150 respiratory illnesses. Of the 60 RSV positive infants in the supplemental cohort, 42 completed the study, whereas all 78 of the RSV positive hospital cohort infants completed the study. A GRSS was calculated for each RSV-infected infant and is being used to analyze each of the complex datasets by correlation with disease severity in univariate and multivariate methods.

Conclusions

The AsPIRES study will provide insights into the complex pathogenesis of RSV infection in healthy full-term infants with primary RSV infection. The analysis will allow assessment of multiple factors potentially influencing the severity of RSV infection including the level of RSV specific antibodies, the innate immune response of nasal epithelial cells, the adaptive response by various lymphocyte subsets, the resident airway microbiota, and viral factors. Results of this study will inform disease interventions such as vaccines and antiviral therapies.

Computer-Optimized Immunization Strategy to Boost Highly Cross-Reactive Antibodies to RSV

Respiratory Syncytial Virus (RSV) is a tremendous burden on public health. RSV infects humans within two years of life and causes continual infections throughout their lifetime. RSV is estimated to cause approximately 30 million episodes of acute lower respiratory infection resulting in over 100,000 deaths annually (in children ages 0 – 5 alone). Despite this burden, no vaccine is currently available for RSV. RSV has been classified into two subtypes (A and B) based on genomic and antigenic properties. Evolutionary rates are similar to other RNA viruses with estimates of 1.83–1.95 x 10−3 nucleotide substitutions/site/year. The high mutation rates of RSV produce an evolutionary landscape that has led to genetic and antigenic variation across and within subtypes. The attachment glycoprotein (G-protein) is must abundant on the surface of the virus and contains an almost completely conserved region, the central conserved domain (CCD). We developed a computational model of the humoral immune system that simulates B cell responses to G-protein antigen. We use a sequence-based approach to approximate antigenic differences across (&gt;1000) RSV strains. We then use computer simulations to identify an immunization strategy that can boost antibody levels to the CCD region. Last, we demonstrate that the immunization strategy results in a highly cross-reactive antibody response.

Black Race Is Associated with a Lower Risk of Bronchopulmonary Dysplasia

OBJECTIVE

To use a large current prospective cohort of infants <29 weeks to compare bronchopulmonary dysplasia (BPD) rates in black and white infants.

STUDY DESIGN

The Prematurity and Respiratory Outcome Program (PROP) enrolled 835 infants born in 2011-2013 at <29 weeks of gestation; 728 black or white infants survived to 36 weeks postmenstrual age (PMA). Logistic regression was used to compare BPD outcomes (defined as supplemental oxygen requirement at 36 weeks PMA) between the races, adjusted for gestational age (GA), antenatal steroid use, intubation at birth, and surfactant use at birth.

RESULTS

Of 707 black or white infants with available BPD outcomes, BPD was lower in black infants (38% vs 45%), even though they were of significantly lower GA. At every GA, BPD was more common in white infants. The aOR for BPD was 0.60 (95% CI, 0.42-0.85; P = .004) for black infants compared with white infants after adjusting for GA. Despite the lower rate of BPD, black infants had a higher rate of first-year post-prematurity respiratory disease (black, 79%; white, 63%).

CONCLUSIONS

In this large cohort of recently born preterm infants at <29 weeks GA, compared with white infants, black infants had a lower risk of BPD but an increased risk of persistent respiratory morbidity.

Exome sequencing identifies gene variants and networks associated with extreme respiratory outcomes following preterm birth

Background

Previous studies have identified genetic variants associated with bronchopulmonary dysplasia (BPD) in extremely preterm infants. However, findings with genome-wide significance have been rare, and not replicated. We hypothesized that whole exome sequencing (WES) of premature subjects with extremely divergent phenotypic outcomes could facilitate the identification of genetic variants or gene networks contributing disease risk.

Results

The Prematurity and Respiratory Outcomes Program (PROP) recruited a cohort of > 765 extremely preterm infants for the identification of markers of respiratory morbidity. We completed WES on 146 PROP subjects (85 affected, 61 unaffected) representing extreme phenotypes of early respiratory morbidity. We tested for association between disease status and individual common variants, screened for rare variants exclusive to either affected or unaffected subjects, and tested the combined association of variants across gene loci. Pathway analysis was performed and disease-related expression patterns were assessed. Marginal association with BPD was observed for numerous common and rare variants. We identified 345 genes with variants unique to BPD-affected preterm subjects, and 292 genes with variants unique to our unaffected preterm subjects. Of these unique variants, 28 (19 in the affected cohort and 9 in unaffected cohort) replicate a prior WES study of BPD-associated variants. Pathway analysis of sets of variants, informed by disease-related gene expression, implicated protein kinase A, MAPK and Neuregulin/epidermal growth factor receptor signaling.

Conclusions

We identified novel genes and associated pathways that may play an important role in susceptibility/resilience for the development of lung disease in preterm infants.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0679-7) contains supplementary material, which is available to authorized users.

PPARα-mediated peroxisome induction compensates PPARγ-deficiency in bronchiolar club cells

Despite the important functions of PPARγ in various cell types of the lung, PPARγ-deficiency in club cells induces only mild emphysema. Peroxisomes are distributed in a similar way as PPARγ in the lung and are mainly enriched in club and AECII cells. To date, the effects of PPARγ-deficiency on the overall peroxisomal compartment and its metabolic alterations in pulmonary club cells are unknown. Therefore, we characterized wild-type and club cell-specific PPARγ knockout-mice lungs and used C22 cells to investigate the peroxisomal compartment and its metabolic roles in the distal airway epithelium by means of 1) double-immunofluorescence labelling for peroxisomal proteins, 2) laser-assisted microdissection of the bronchiolar epithelium and subsequent qRT-PCR, 3) siRNA-transfection of PPARγand PPRE dual-luciferase reporter activity in C22 cells, 4) PPARg inhibition by GW9662, 5) GC-MS based lipid analysis. Our results reveal elevated levels of fatty acids, increased expression of PPARα and PPRE activity, a strong overall upregulation of the peroxisomal compartment and its associated gene expression (biogenesis, α-oxidation, β-oxidation, and plasmalogens) in PPARγ-deficient club cells. Interestingly, catalase was significantly increased and mistargeted into the cytoplasm, suggestive for oxidative stress by the PPARγ-deficiency in club cells. Taken together, PPARα-mediated metabolic induction and proliferation of peroxisomes via a PPRE-dependent mechanism could compensate PPARγ-deficiency in club cells.

Cell type-resolved human lung lipidome reveals cellular cooperation in lung function

Cell type-resolved proteome analyses of the brain, heart and liver have been reported, however a similar effort on the lipidome is currently lacking. Here we applied liquid chromatography-tandem mass spectrometry to characterize the lipidome of major lung cell types isolated from human donors, representing the first lipidome map of any organ. We coupled this with cell type-resolved proteomics of the same samples (available at Lungmap.net). Complementary proteomics analyses substantiated the functional identity of the isolated cells. Lipidomics analyses showed significant variations in the lipidome across major human lung cell types, with differences most evident at the subclass and intra-subclass (i.e. total carbon length of the fatty acid chains) level. Further, lipidomic signatures revealed an overarching posture of high cellular cooperation within the human lung to support critical functions. Our complementary cell type-resolved lipid and protein datasets serve as a rich resource for analyses of human lung function.

Dissociation, cellular isolation, and initial molecular characterization of neonatal and pediatric human lung tissues

Human lung morphogenesis begins by embryonic life and continues after birth into early childhood to form a complex organ with numerous morphologically and functionally distinct cell types. Pulmonary organogenesis involves dynamic changes in cell proliferation, differentiation, and migration of specialized cells derived from diverse embryonic lineages. Studying the molecular and cellular processes underlying formation of the fully functional lung requires isolating distinct pulmonary cell populations during development. We now report novel methods to isolate four major pulmonary cell populations from pediatric human lung simultaneously. Cells were dissociated by protease digestion of neonatal and pediatric lung and isolated on the basis of unique cell membrane protein expression patterns. Epithelial, endothelial, nonendothelial mesenchymal, and immune cells were enriched by fluorescence-activated cell sorting. Dead cells and erythrocytes were excluded by 7-aminoactinomycin D uptake and glycophorin-A (CD235a) expression, respectively. Leukocytes were identified by membrane CD45 (protein tyrosine phosphatase, receptor type C), endothelial cells by platelet endothelial cell adhesion molecule-1 (CD31) and vascular endothelial cadherin (CD144), and both were isolated. Thereafter, epithelial cell adhesion molecule (CD326)-expressing cells were isolated from the endothelial- and immune cell-depleted population to enrich epithelial cells. Cells lacking these membrane markers were collected as "nonendothelial mesenchymal" cells. Quantitative RT-PCR and RNA sequencing analyses of population specific transcriptomes demonstrate the purity of the subpopulations of isolated cells. The method efficiently isolates major human lung cell populations that we announce are now available through the National Heart, Lung, and Blood Institute Lung Molecular Atlas Program (LungMAP) for their further study.

T cell developmental arrest in former premature infants increases risk of respiratory morbidity later in infancy

The inverse relationship between gestational age at birth and postviral respiratory morbidity suggests that infants born preterm (PT) may miss a critical developmental window of T cell maturation. Despite a continued increase in younger PT survivors with respiratory complications, we have limited understanding of normal human fetal T cell maturation, how ex utero development in premature infants may interrupt normal T cell development, and whether T cell development has an effect on infant outcomes. In our longitudinal cohort of 157 infants born between 23 and 42 weeks of gestation, we identified differences in T cells present at birth that were dependent on gestational age and differences in postnatal T cell development that predicted respiratory outcome at 1 year of age. We show that naive CD4+ T cells shift from a CD31-TNF-α+ bias in mid gestation to a CD31+IL-8+ predominance by term gestation. Former PT infants discharged with CD31+IL8+CD4+ T cells below a range similar to that of full-term born infants were at an over 3.5-fold higher risk for respiratory complications after NICU discharge. This study is the first to our knowledge to identify a pattern of normal functional T cell development in later gestation and to associate abnormal T cell development with health outcomes in infants.

Association of Dynamic Changes in the CD4 T-Cell Transcriptome With Disease Severity During Primary Respiratory Syncytial Virus Infection in Young Infants

Background

Nearly all children are infected with respiratory syncytial virus (RSV) within the first 2 years of life, with a minority developing severe disease (1%-3% hospitalized). We hypothesized that an assessment of the adaptive immune system, using CD4+ T-lymphocyte transcriptomics, would identify gene expression correlates of disease severity.

Methods

Infants infected with RSV representing extremes of clinical severity were studied. Mild illness (n = 23) was defined as a respiratory rate (RR) < 55 and room air oxygen saturation (SaO2) ≥ 97%, and severe illness (n = 23) was defined as RR ≥ 65 and SaO2 ≤ 92%. RNA from fresh, sort-purified CD4+ T cells was assessed by RNA sequencing.

Results

Gestational age, age at illness onset, exposure to environmental tobacco smoke, bacterial colonization, and breastfeeding were associated (adjusted P < .05) with disease severity. RNA sequencing analysis reliably measured approximately 60% of the genome. Severity of RSV illness had the greatest effect size upon CD4 T-cell gene expression. Pathway analysis identified correlates of severity, including JAK/STAT, prolactin, and interleukin 9 signaling. We also identified genes and pathways associated with timing of symptoms and RSV group (A/B).

Conclusions

These data suggest fundamental changes in adaptive immune cell phenotypes may be associated with RSV clinical severity.

Impact of prematurity and nutrition on the developing gut microbiome and preterm infant growth

BACKGROUND

Identification of factors that influence the neonatal gut microbiome is urgently needed to guide clinical practices that support growth of healthy preterm infants. Here, we examined the influence of nutrition and common practices on the gut microbiota and growth in a cohort of preterm infants.

RESULTS

With weekly gut microbiota samples spanning postmenstrual age (PMA) 24 to 46 weeks, we developed two models to test associations between the microbiota, nutrition and growth: a categorical model with three successive microbiota phases (P1, P2, and P3) and a model with two periods (early and late PMA) defined by microbiota composition and PMA, respectively. The more significant associations with phase led us to use a phase-based framework for the majority of our analyses. Phase transitions were characterized by rapid shifts in the microbiota, with transition out of P1 occurring nearly simultaneously with the change from meconium to normal stool. The rate of phase progression was positively associated with gestational age at birth, and delayed transition to a P3 microbiota was associated with growth failure. We found distinct bacterial metabolic functions in P1-3 and significant associations between nutrition, microbiota phase, and infant growth.

CONCLUSION

The phase-dependent impact of nutrition on infant growth along with phase-specific metabolic functions suggests a pioneering potential for improving growth outcomes by tailoring nutrient intake to microbiota phase.

LungMAP: The Molecular Atlas of Lung Development Program

The National Heart, Lung, and Blood Institute is funding an effort to create a molecular atlas of the developing lung (LungMAP) to serve as a research resource and public education tool. The lung is a complex organ with lengthy development time driven by interactive gene networks and dynamic cross talk among multiple cell types to control and coordinate lineage specification, cell proliferation, differentiation, migration, morphogenesis, and injury repair. A better understanding of the processes that regulate lung development, particularly alveologenesis, will have a significant impact on survival rates for premature infants born with incomplete lung development and will facilitate lung injury repair and regeneration in adults. A consortium of four research centers, a data coordinating center, and a human tissue repository provides high-quality molecular data of developing human and mouse lungs. LungMAP includes mouse and human data for cross correlation of developmental processes across species. LungMAP is generating foundational data and analysis, creating a web portal for presentation of results and public sharing of data sets, establishing a repository of young human lung tissues obtained through organ donor organizations, and developing a comprehensive lung ontology that incorporates the latest findings of the consortium. The LungMAP website (www.lungmap.net) currently contains more than 6,000 high-resolution lung images and transcriptomic, proteomic, and lipidomic human and mouse data and provides scientific information to stimulate interest in research careers for young audiences. This paper presents a brief description of research conducted by the consortium, database, and portal development and upcoming features that will enhance the LungMAP experience for a community of users.

Transcriptomic Biomarkers to Discriminate Bacterial from Nonbacterial Infection in Adults Hospitalized with Respiratory Illness

Lower respiratory tract infection (LRTI) commonly causes hospitalization in adults. Because bacterial diagnostic tests are not accurate, antibiotics are frequently prescribed. Peripheral blood gene expression to identify subjects with bacterial infection is a promising strategy. We evaluated whole blood profiling using RNASeq to discriminate infectious agents in adults with microbiologically defined LRTI. Hospitalized adults with LRTI symptoms were recruited. Clinical data and blood was collected, and comprehensive microbiologic testing performed. Gene expression was measured using RNASeq and qPCR. Genes discriminatory for bacterial infection were identified using the Bonferroni-corrected Wilcoxon test. Constrained logistic models to predict bacterial infection were fit using screened LASSO. We enrolled 94 subjects who were microbiologically classified; 53 as “non-bacterial” and 41 as “bacterial”. RNAseq and qPCR confirmed significant differences in mean expression for 10 genes previously identified as discriminatory for bacterial LRTI. A novel dimension reduction strategy selected three pathways (lymphocyte, α-linoleic acid metabolism, IGF regulation) including eleven genes as optimal markers for discriminating bacterial infection (naïve AUC = 0.94; nested CV-AUC = 0.86). Using these genes, we constructed a classifier for bacterial LRTI with 90% (79% CV) sensitivity and 83% (76% CV) specificity. This novel, pathway-based gene set displays promise as a method to distinguish bacterial from nonbacterial LRTI.

Lung Gene Expression Analysis (LGEA): an integrative web portal for comprehensive gene expression data analysis in lung development

'LungGENS', our previously developed web tool for mapping single-cell gene expression in the developing lung, has been well received by the pulmonary research community. With continued support from the 'LungMAP' consortium, we extended the scope of the LungGENS database to accommodate transcriptomics data from pulmonary tissues and cells from human and mouse at different stages of lung development. Lung Gene Expression Analysis (LGEA) web portal is an extended version of LungGENS useful for the analysis, display and interpretation of gene expression patterns obtained from single cells, sorted cell populations and whole lung tissues. The LGEA web portal is freely available at http://research.cchmc.org/pbge/lunggens/mainportal.html.

The Healthy Infant Nasal Transcriptome: A Benchmark Study

Responses by resident cells are likely to play a key role in determining the severity of respiratory disease. However, sampling of the airways poses a significant challenge, particularly in infants and children. Here, we report a reliable method for obtaining nasal epithelial cell RNA from infants for genome-wide transcriptomic analysis, and describe baseline expression characteristics in an asymptomatic cohort. Nasal epithelial cells were collected by brushing of the inferior turbinates, and gene expression was interrogated by RNA-seq analysis. Reliable recovery of RNA occurred in the absence of adverse events. We observed high expression of epithelial cell markers and similarity to the transcriptome for intrapulmonary airway epithelial cells. We identified genes displaying low and high expression variability, both inherently, and in response to environmental exposures. The greatest gene expression differences in this asymptomatic cohort were associated with the presence of known pathogenic viruses and/or bacteria. Robust bacteria-associated gene expression patterns were significantly associated with the presence of Moraxella. In summary, we have developed a reliable method for interrogating the infant airway transcriptome by sampling the nasal epithelium. Our data demonstrates both the fidelity and feasibility of our methodology, and describes normal gene expression and variation within a healthy infant cohort.