Search

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

The CD4+ T cell response is critical to host protection against helminth infection. How this response varies across different hosts and tissues remains an important gap in our understanding. Using IL-4-reporter mice to identify responding CD4+ T cells to Nippostrongylus brasiliensis infection, T cell receptor sequencing paired with novel clustering algorithms revealed a broadly reactive and clonally diverse CD4+ T cell response. While the most prevalent clones and clonotypes exhibited some tissue selectivity, most were observed to reside in both the lung and lung-draining lymph nodes. Antigen-reactivity of the broader repertoires was predicted to be shared across both tissues and individual mice. Transcriptome, trajectory, and chromatin accessibility analysis of lung and lymph-node repertoires revealed three unique but related populations of responding IL-4+ CD4+ T cells consistent with T follicular helper, T helper 2, and a transitional population sharing similarity with both populations. The shared antigen reactivity of lymph node and lung repertoires combined with the adoption of tissue-specific gene programs allows for the pairing of cellular and humoral responses critical to the orchestration of anti-helminth immunity., Author summary Using various “omic” approaches, the CD4+ T cell receptor (TCR) repertoire was explored after primary helminth infection. Infection generated a broadly reactive and clonally diverse CD4+ T cell response with the most prevalent clonotypes and predicted antigen specificities residing in both the lung and lung-draining lymph nodes. Tissue-specific programming of responding CD4+ T cells directed the establishment of committed Tfh and Th2 cells, both critical for driving distinct hallmarks of type-2 inflammation. These datasets help to explore the diverse yet tissue-specific nature of anti-helminth immunity.

Air pollution particulate matter 100 genes), which included transcriptional drivers of mucus metaplasia (SPDEF and FOXA3). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the MUC5AC/MUC5B gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., FOXJ1 and MCIDAS). Chronic OE stimulation induced mucus metaplasia–like remodeling characterized by increases in MUC5AC(+) secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM(2.5) exposure.

Rationale: Coronavirus disease (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). ACE2 (angiotensin-converting enzyme 2), and TMPRSS2 (transmembrane protease serine 2) mediate viral infection of host cells. We reasoned that differences in ACE2 or TMPRSS2 gene expression in sputum cells among patients with asthma may identify subgroups at risk for COVID-19 morbidity.Objectives: To determine the relationship between demographic features and sputum ACE2 and TMPRSS2 gene expression in asthma.Methods: We analyzed gene expression for ACE2 and TMPRSS2, and for ICAM-1 (intercellular adhesion molecule 1) (rhinovirus receptor as a comparator) in sputum cells from 330 participants in SARP-3 (Severe Asthma Research Program-3) and 79 healthy control subjects.Measurements and Main Results: Gene expression of ACE2 was lower than TMPRSS2, and expression levels of both genes were similar in asthma and health. Among patients with asthma, male sex, African American race, and history of diabetes mellitus were associated with higher expression of ACE2 and TMPRSS2. Use of inhaled corticosteroids (ICS) was associated with lower expression of ACE2 and TMPRSS2, but treatment with triamcinolone acetonide did not decrease expression of either gene. These findings differed from those for ICAM-1, where gene expression was increased in asthma and less consistent differences were observed related to sex, race, and use of ICS.Conclusions: Higher expression of ACE2 and TMPRSS2 in males, African Americans, and patients with diabetes mellitus provides rationale for monitoring these asthma subgroups for poor COVID-19 outcomes. The lower expression of ACE2 and TMPRSS2 with ICS use warrants prospective study of ICS use as a predictor of decreased susceptibility to SARS-CoV-2 infection and decreased COVID-19 morbidity.

SUMMARY The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function., In Brief Using airway epithelial cell cultures, Jackson et al. show that IL-13, a driver of type 2-high asthma, induces emergent mucus secretory expression states for each cell type. This program universally diminishes innate airway defense, produces a pathologic mucus secretome that arrests mucociliary movement, and is recapitulated in type 2 inflamed children., Graphical Abstract

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and non-smokers, we generate a comprehensive atlas of epithelial cell types and states, connect these into lineages, and define cell-specific responses to smoking. Our analysis infers multi-state lineages that develop into surface mucus secretory and ciliated cells and then contrasts these to the unique specification of submucosal gland (SMG) cells. Accompanying knockout studies reveal that tuft-like cells are the likely progenitor of both pulmonary neuroendocrine cells and CFTR-rich ionocytes. Our smoking analysis finds that all cell types, including protected stem and SMG populations, are affected by smoking through both pan-epithelial smoking response networks and hundreds of cell-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium., Chronic lung diseases are characterized by molecular and cellular composition changes. Here the authors use single-cell RNA sequencing to map cell type-specific changes in human tracheal epithelium related to smoking, and to provide evidence for a tuft-like progenitor for pulmonary neuroendocrine cells and ionocytes.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes., ACE2 and TMPRSS2 have received recent attention as entry factors for SARS-CoV-2. Here the authors analyze nasal airway transcriptome data from 695 children determining ACE2 and TMPRSS2 expression is induced by viral and type2 inflammation, respectively, and both exhibit eQTLs that vary across world populations.

High-throughput RNA-sequencing (RNA-seq) technologies provide an unprecedented opportunity to explore the individual transcriptome. Unmapped reads are a large and often overlooked output of standard RNA-seq analyses. Here, we present Read Origin Protocol (ROP), a tool for discovering the source of all reads originating from complex RNA molecules. We apply ROP to samples across 2630 individuals from 54 diverse human tissues. Our approach can account for 99.9% of 1 trillion reads of various read length. Additionally, we use ROP to investigate the functional mechanisms underlying connections between the immune system, microbiome, and disease. ROP is freely available at https://github.com/smangul1/rop/wiki. Electronic supplementary material The online version of this article (10.1186/s13059-018-1403-7) contains supplementary material, which is available to authorized users.

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and nonsmokers, we generated a comprehensive atlas of epithelial cell types and states, connected these into lineages, and defined cell-specific responses to smoking. Our analysis inferred multi-state lineages that develop into surface mucus secretory and ciliated cells and contrasted these to the unique lineage and specialization of submucosal gland (SMG) cells. Our analysis also suggests a lineage relationship between tuft, pulmonary neuroendocrine, and the newly discovered CFTR-rich ionocyte cells. Our smoking analysis found that all cell types, including protected stem and SMG populations, are affected by smoking, through both pan-epithelial smoking response networks and hundreds of cell type-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium.

Skin barrier dysfunction has been reported in both atopic dermatitis (AD) and food allergy (FA). However, only one-third of patients with AD have FA. The purpose of this study was to use a minimally invasive skin tape strip sampling method and a multiomics approach to determine whether children with AD and FA (AD FA+) have stratum corneum (SC) abnormalities that distinguish them from AD without FA (AD FA−) and nonatopic (NA) controls. Transepidermal water loss was found to be increased in AD FA+. Filaggrin and the proportion of ω-hydroxy fatty acid sphingosine ceramide content in nonlesional skin of children with AD FA+ were substantially lower than in AD FA− and NA skin. These abnormalities correlated with morphologic changes in epidermal lamellar bilayer architecture responsible for barrier homeostasis. Shotgun metagenomic studies revealed that the nonlesional skin of AD FA+ had increased abundance of Staphylococcus aureus compared to NA. Increased expression of keratins 5, 14, and 16 indicative of hyperproliferative keratinocytes was observed in the SC of AD FA+. The skin transcriptome of AD FA+ had increased gene expression for dendritic cells and type 2 immune pathways. A network analysis revealed keratins 5, 14, and 16 were positively correlated with AD FA+, whereas filaggrin breakdown products were negatively correlated with AD FA+. These data suggest that the most superficial compartment of nonlesional skin in AD FA+ has unique properties associated with an immature skin barrier and type 2 immune activation.

BackgroundType 2 (T2) inflammation drives airway dysfunction in many patients with asthma; yet, we lack a comprehensive understanding of the airway immune cell types and networks that sustain this inflammation. Moreover, defects in the airway immune system in patients with asthma without T2 inflammation are not established.ObjectivesTo determine the gene networks that sustain T2 airway inflammation in T2-high asthma and to explore the gene networks that characterize T2-low asthma.MethodsNetwork analysis of sputum cell transcriptome expression data from 84 subjects with asthma and 27 healthy control subjects was used to identify immune cell type-enriched networks that underlie asthma subgroups.ResultsSputum T2 gene expression was characterized by an immune cell network derived from multiple innate immune cells, including eosinophils, mast cells/basophils, and inflammatory dendritic cells. Clustering of subjects within this network stratified subjects into T2-high and T2-low groups, but it also revealed a subgroup of T2-high subjects with uniformly higher expression of the T2 network. These "T2-ultrahigh subjects" were characterized clinically by older age and more severe airflow obstruction and pathologically by a second T2 network derived from T2-skewed, CD11b+/CD103-/IRF4+ classical dendritic cells. Subjects with T2-low asthma were differentiated from healthy control subjects by lower expression of a cytotoxic CD8+ T-cell network, which was negatively correlated with body mass index and plasma IL-6 concentrations.ConclusionsPersistent airway T2 inflammation is a complex construct of innate and adaptive immunity gene expression networks that are variable across individuals with asthma and persist despite steroid treatment. Individuals with T2-low asthma exhibit an airway deficiency in cytotoxic T cells associated with obesity-driven inflammation.

The application of conditional reprogramming culture (CRC) methods to nasal airway epithelial cells would allow more wide-spread incorporation of primary airway epithelial culture models into complex lung disease research. In this study, we adapted the CRC method to nasal airway epithelial cells, investigated the growth advantages afforded by this technique over standard culture methods, and determined the cellular and molecular basis of CRC cell culture effects. We found that the CRC method allowed the production of 7.1 × 10(10) cells after 4 passages, approximately 379 times more cells than were generated by the standard bronchial epithelial growth media (BEGM) method. These nasal airway epithelial cells expressed normal basal cell markers and could be induced to form a mucociliary epithelium. Progenitor cell frequency was significantly higher using the CRC method in comparison to the standard culture method, and progenitor cell maintenance was dependent on addition of the Rho-kinase inhibitor Y-27632. Whole-transcriptome sequencing analysis demonstrated widespread gene expression changes in Y-27632-treated basal cells. We found that Y-27632 treatment altered expression of genes fundamental to the formation of the basal cell cytoskeleton, cell-cell junctions, and cell-extracellular matrix (ECM) interactions. Importantly, we found that Y-27632 treatment up-regulated expression of unique basal cell intermediate filament and desmosomal genes. Conversely, Y-27632 down-regulated multiple families of protease/antiprotease genes involved in ECM remodeling. We conclude that Y-27632 fundamentally alters cell-cell and cell-ECM interactions, which preserves basal progenitor cells and allows greater cell amplification.

Background Research in transformed immortalized cell lines indicates the cadherin-related family member 3 (CDHR3) protein serves as a receptor for human rhinovirus (HRV)–C. Similar experiments indicate that the CDHR3 coding variant rs6967330 increases CDHR3 protein surface expression. Objective We sought to determine whether CDHR3 is necessary for HRV-C infection of primary airway epithelial cells (AECs) and to identify molecular mechanisms by which CDHR3 variants confer risk for asthma exacerbations. Methods CDHR3 function and influence on HRV-C infection were investigated by using single-cell transcriptomics, CRISPR-Cas9 gene knockout, and genotype-specific donor experiments performed in primary AECs. Nasal airway epithelium cis–expression quantitative trait locus (eQTL) analysis of CDHR3 was performed, followed by association testing for asthma hospitalization in minority children. Results CDHR3 lung expression is exclusive to ciliated AECs and associated with basal bodies during and after motile ciliogenesis. Knockout of CDHR3 in human AECs did not prevent ciliated cell differentiation but was associated with a decrease in transepithelial resistance and an 80% decrease in HRV-C infection of the mucociliary epithelium. AECs from subjects homozygous for the risk-associated rs6967330 single nucleotide polymorphism (SNP) exhibited greater HRV-C infection compared with cells homozygous for the nonrisk allele. AEC cis-eQTL analysis indicated that rs6967330 and other SNPs are eQTLs for CDHR3. Only the eQTL block containing the rs6967330 SNP showed a significant association with childhood asthma hospitalization. Conclusions Genetic deletion and genotype-specific studies in primary AECs indicate CDHR3 is critical to HRV-C infection of ciliated cells. The rs6967330 SNP confers risk of severe childhood asthma exacerbations, likely through increasing HRV-C infection levels and protein surface localization.

The airway epithelium lines the respiratory tract and provides the primary protective barrier against inhalational insults including toxic environmental substances and microorganisms. The airway epithelium also plays a critical role in regulating airway immune responses. The airway epithelial response to the type 2 cytokine, interleukin-13 (IL-13), is critical to airway inflammation, mucus production, and airway hyperresponsiveness present in asthma. Relevant primary cell models of the human airway epithelium are needed to investigate the biology of IL-13-mediated airway epithelial effects. Here, we describe the generation of a differentiated mucociliary human airway epithelium using an in vitro air-liquid interface (ALI) culture model system. We also describe methods to stimulate this culture model with IL-13 and harvest cells and biomolecules to interrogate cellular and molecular aspects of the airway epithelial IL-13 response.

Respiratory illness caused by viral infection is associated with the development and exacerbation of childhood asthma. Little is known about the effects of respiratory viral infections in the absence of illness. Using quantitative PCR (qPCR) for common respiratory viruses and for two genes known to be highly upregulated in viral infections (CCL8/CXCL11), we screened 92 asthmatic and 69 healthy children without illness for respiratory virus infections.We found 21 viral qPCR-positive and 2 suspected virus-infected subjects with high expression of CCL8/CXCL11. We applied a dual RNA-seq workflow to these subjects, together with 25 viral qPCR-negative subjects, to compare qPCR with sequencing-based virus detection and to generate the airway transcriptome for analysis. RNA-seq virus detection achieved 86% sensitivity when compared to qPCR-based screening. We detected additional respiratory viruses in the two CCL8/CXCL11-high subjects and in two of the qPCR-negative subjects. Viral read counts varied widely and were used to stratify subjects into Virus-High and Virus-Low groups. Examination of the host airway transcriptome found that the Virus-High group was characterized by immune cell airway infiltration, downregulation of cilia genes, and dampening of type 2 inflammation. Even the Virus-Low group was differentiated from the No-Virus group by 100 genes, some involved in eIF2 signaling.Respiratory virus infection without illness is not innocuous but may determine the airway function of these subjects by driving immune cell airway infiltration, cellular remodeling, and alteration of asthmogenic gene expression.

Lipids in the stratum corneum of atopic dermatitis (AD) patients differ substantially in composition from healthy subjects. We hypothesized that hyperactivated type 2 immune response alters AD skin lipid metabolism. We have analyzed stratum corneum lipids from nonlesional and lesional skin of AD subjects and IL-13 skin-specific Tg mice. We also directly examined the effects of IL-4/IL-13 on human keratinocytes in vitro. Mass spectrometric analysis of lesional stratum corneum from AD subjects and IL-13 Tg mice revealed an increased proportion of short-chain (N-14:0 to N-24:0) NS ceramides, sphingomyelins, and 14:0–22:0 lysophosphatidylcholines (14:0–22:0 LPC) with a simultaneous decline in the proportion of corresponding long-chain species (N-26:0 to N-32:0 sphingolipids and 24:0–30:0 LPC) when compared with healthy controls. An increase in short-chain LPC species was also observed in nonlesional AD skin. Similar changes were observed in IL-4/IL-13–driven responses in Ca2+-differentiated human keratinocytes in vitro, all being blocked by STAT6 silencing with siRNA. RNA sequencing analysis performed on stratum corneum of AD as compared with healthy subjects identified decreased expression of fatty acid elongases ELOVL3 and ELOVL6 that contributed to observed changes in atopic skin lipids. IL-4/IL-13 also inhibited ELOVL3 and ELOVL6 expression in keratinocyte cultures in a STAT6-dependent manner. Downregulation of ELOVL3/ELOVL6 expression in keratinocytes by siRNA decreased the proportion of long-chain fatty acids globally and in sphingolipids. Thus, our data strongly support the pathogenic role of type 2 immune activation in AD skin lipid metabolism.

The adaptation of the clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR associated endonuclease 9 (CRISPR-Cas9) machinery from prokaryotic organisms has resulted in a gene editing system that is highly versatile, easily constructed, and can be leveraged to generate human cells knocked out (KO) for a specific gene. While standard transfection techniques can be used for the introduction of CRISPR-Cas9 expression cassettes to many cell types, delivery by this method is not efficient in many primary cell types, including primary human airway epithelial cells (AECs). More efficient delivery in AECs can be achieved through lentiviral-mediated transduction, allowing the CRISPR-Cas9 system to be integrated into the genome of the cell, resulting in stable expression of the nuclease machinery and increasing editing rates. In parallel, advancements have been made in the culture, expansion, selection, and differentiation of AECs, which allow the robust generation of a bulk edited AEC population from transduced cells. Applying these methods, we detail here our latest protocol to generate mucociliary epithelial cultures knocked out for a specific gene from donor-isolated primary human basal airway epithelial cells. This protocol includes methods to: (1) design and generate lentivirus which targets a specific gene for KO with CRISPR-Cas9 machinery, (2) efficiently transduce AECs, (3) culture and select for a bulk edited AEC population, (4) molecularly screen AECs for Cas9 cutting and specific sequence edits, and (5) further expand and differentiate edited cells to a mucociliary airway epithelial culture. The AEC knockouts generated using this protocol provide an excellent primary cell model system with which to characterize the function of genes involved in airway dysfunction and disease.

BACKGROUND: Expression profiling of skin biopsies has established molecular features of the skin in atopic dermatitis (AD). Invasiveness of biopsies has prevented their use in defining individual level AD pathobiological mechanisms (endotypes) in large research studies. OBJECTIVE: To determine if minimally invasive skin tape strip transcriptome analysis identifies gene expression dysregulation in AD and molecular disease endotypes. METHODS: We sampled non-lesional and lesional skin tape strips and biopsies from adult Caucasian subjects AD patients (18 males, 12 females; age (Mean±SE) 36.3±2.2 yrs) and healthy controls (9 males, 16 females; age (Mean±SE) 34.8±2.2 yrs). Ampliseq whole transcriptome sequencing was performed on extracted RNA. Differential expression, clustering/pathway analyses, immunostaining of skin biopsies, and clinical trait correlations were performed. RESULTS: Skin tape expression profiles were distinct from skin biopsy profiles and better sampled epidermal differentiation complex genes. Skin tape expression of 29 immune and epidermis-related genes (FDR

Background Ectopic olfactory receptors (ORs) are found in the skin, but their expression and biological function in normal skin and skin form patients with atopic dermatitis (AD) are unknown. Objectives We sought to characterize the expression of ORs in the skin and assess OR-mediated biological responses of primary human keratinocytes in the presence of odorant ligands. Methods OR expression was examined by using whole-transcriptome sequencing of skin tape strips collected from patients with AD and healthy control (HC) subjects. OR10G7 and filaggrin 1 (FLG-1) expression was analyzed by using RT-PCR and immunostaining in skin biopsy specimens and primary human keratinocytes from patients with AD and HC subjects. ATP and cyclic AMP production by control and OR10G7 small interfering RNA–transfected keratinocytes in response to odorant stimulation with acetophenone and eugenol was assessed. Results A total of 381 OR gene transcripts were detected in the skin samples, with the greatest OR expression detected in the skin tape strips corresponding to the upper granular layer of the skin. OR10G7 expression was significantly increased in skin biopsy specimens from patients with AD compared with those from HC subjects (P = .01) and inversely correlated with FLG-1 expression (P = .009). OR10G7 expression was greatest in undifferentiated keratinocytes from patients with AD and was downregulated with progressive differentiation. Primary human keratinocytes produced ATP, an essential neurotransmitter in sensory pathways, in response to acetophenone and eugenol, odorants previously identified as potential ligands for this receptor. This response was abolished in OR10G7 small interfering RNA–transfected keratinocytes. Conclusions OR10G7 is expressed at significantly greater levels in undifferentiated keratinocytes from patients with AD compared with HC subjects. OR10G7 is likely involved in transmission of skin-induced chemosensory responses to odorant stimulation, which might modulate differential nociceptive responses in AD skin.

Background Respiratory illness caused by viral infection is associated with the development and exacerbation of childhood asthma. Little is known about the effects of respiratory viral infections in the absence of illness. Using quantitative PCR (qPCR) for common respiratory viruses and for two genes known to be highly upregulated in viral infections (CCL8/CXCL11), we screened 92 asthmatic and 69 healthy children without illness for respiratory virus infections. Results We found 21 viral qPCR-positive and 2 suspected virus-infected subjects with high expression of CCL8/CXCL11. We applied a dual RNA-seq workflow to these subjects, together with 25 viral qPCR-negative subjects, to compare qPCR with sequencing-based virus detection and to generate the airway transcriptome for analysis. RNA-seq virus detection achieved 86% sensitivity when compared to qPCR-based screening. We detected additional respiratory viruses in the two CCL8/CXCL11-high subjects and in two of the qPCR-negative subjects. Viral read counts varied widely and were used to stratify subjects into Virus-High and Virus-Low groups. Examination of the host airway transcriptome found that the Virus-High group was characterized by immune cell airway infiltration, downregulation of cilia genes, and dampening of type 2 inflammation. Even the Virus-Low group was differentiated from the No-Virus group by 100 genes, some involved in eIF2 signaling. Conclusions Respiratory virus infection without illness is not innocuous but may determine the airway function of these subjects by driving immune cell airway infiltration, cellular remodeling, and alteration of asthmogenic gene expression. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-1140-8) contains supplementary material, which is available to authorized users.

Genome-wide association studies of asthma have identified genetic variants in the IL1RL1 gene, but the molecular mechanisms conferring risk are unknown. IL1RL1 encodes the ST2 receptor (ST2L) for IL-33 and an inhibitory decoy receptor (sST2). IL-33 promotes type 2 inflammation, which is present in some but not all asthmatics. We find that two single nucleotide polymorphisms (SNPs) in IL1RL1 - rs1420101 and rs11685480 - are strongly associated with plasma sST2 levels, though neither is an expression quantitative trait locus (eQTL) in whole blood. Rather, rs1420101 and rs11685480 mark eQTLs in airway epithelial cells and distal lung parenchyma, respectively. We find that the genetically determined plasma sST2 reservoir, derived from the lung, neutralizes IL-33 activity, and these eQTL SNPs additively increase the risk of airway type 2 inflammation among asthmatics. These risk variants define a population of asthmatics at risk of IL-33-driven type 2 inflammation.

Type 2 inflammation occurs in a large subgroup of asthmatics, and novel cytokine-directed therapies are being developed to treat this population. In mouse models, interleukin-33 (IL-33) activates lung resident innate lymphoid type 2 cells (ILC2s) to initiate airway type 2 inflammation. In human asthma, which is chronic and difficult to model, the role of IL-33 and the target cells responsible for persistent type 2 inflammation remain undefined. Full-length IL-33 is a nuclear protein and may function as an "alarmin" during cell death, a process that is uncommon in chronic stable asthma. We demonstrate a previously unidentified mechanism of IL-33 activity that involves alternative transcript splicing, which may operate in stable asthma. In human airway epithelial cells, alternative splicing of the IL-33 transcript is consistently present, and the deletion of exons 3 and 4 (Δ exon 3,4) confers cytoplasmic localization and facilitates extracellular secretion, while retaining signaling capacity. In nonexacerbating asthmatics, the expression of Δ exon 3,4 is strongly associated with airway type 2 inflammation, whereas full-length IL-33 is not. To further define the extracellular role of IL-33 in stable asthma, we sought to determine the cellular targets of its activity. Comprehensive flow cytometry and RNA sequencing of sputum cells suggest basophils and mast cells, not ILC2s, are the cellular sources of type 2 cytokines in chronic asthma. We conclude that IL-33 isoforms activate basophils and mast cells to drive type 2 inflammation in chronic stable asthma, and novel IL-33 inhibitors will need to block all biologically active isoforms.

High throughput RNA sequencing technologies have provided invaluable research opportunities across distinct scientific domains by producing quantitative readouts of the transcriptional activity of both entire cellular populations and single cells. The majority of RNA-Seq analyses begin by mapping each experimentally produced sequence (i.e., read) to a set of annotated reference sequences for the organism of interest. For both biological and technical reasons, a significant fraction of reads remains unmapped. In this work, we develop Read Origin Protocol (ROP) to discover the source of all reads originating from complex RNA molecules, recombinant T and B cell receptors, and microbial communities. We applied ROP to 8,641 samples across 630 individuals from 54 tissues. A fraction of RNA-Seq data (n=86) was obtained in-house; the remaining data was obtained from the Genotype-Tissue Expression (GTEx v6) project. To generalize the reported number of accounted reads, we also performed ROP analysis on thousands of different, randomly selected, and publicly available RNA-Seq samples in the Sequence Read Archive (SRA). Our approach can account for 99.9% of 1 trillion reads of various read length across the merged dataset (n=10641). Using in-house RNA-Seq data, we show that immune profiles of asthmatic individuals are significantly different from the profiles of control individuals, with decreased average per sample T and B cell receptor diversity. We also show that immune diversity is inversely correlated with microbial load. Our results demonstrate the potential of ROP to exploit unmapped reads in order to better understand the functional mechanisms underlying connections between the immune system, microbiome, human gene expression, and disease etiology. ROP is freely available athttps://github.com/smangul1/ropand currently supports human and mouse RNA-Seq reads.

Targeted knockout of genes in primary human cells using CRISPR-Cas9 mediated genome-editing represents a powerful approach to study gene function and to discern molecular mechanisms underlying complex human diseases. We used lentiviral delivery of CRISPR-Cas9 machinery and conditional reprogramming culture methods to knockout the MUC18 gene in human primary nasal airway epithelial cells (AECs). Massively parallel sequencing technology was used to confirm that the genome of essentially all cells in the edited AEC populations contained coding region insertions and deletions (indels). Correspondingly, we found mRNA expression of MUC18 was greatly reduced and protein expression was absent. Characterization of MUC18 knockout cell populations stimulated with TLR2, 3 and 4 agonists revealed that IL-8 (a pro-inflammatory chemokine) responses of AECs were greatly reduced in the absence of functional MUC18 protein. Our results show the feasibility of CRISPR-Cas9 mediated gene knockouts in AEC culture (both submerged and polarized), and suggest a pro-inflammatory role for MUC18 in airway epithelial response to bacterial and viral stimuli.