COVID-19-related Genes in Sputum Cells in Asthma. Relationship to Demographic Features and Corticosteroids.

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.

Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM2.5.

Air pollution particulate matter <2.5 µm (PM2.5) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM2.5-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM2.5 samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM2.5 extracts. We found that PM2.5 organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through CYP1A1, IL1A, and IL1B. This OE exposure also robustly activated a mucus secretory expression program (>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 PM2.5 exposure.

RNA-binding specificity landscape of the pentatricopeptide repeat protein PPR10.

Pentatricopeptide repeat (PPR) proteins comprise a large family of helical repeat proteins that influence gene expression in mitochondria and chloroplasts. PPR tracts can bind RNA via a modular one repeat-one nucleotide mechanism in which the nucleotide is specified by the identities of several amino acids in each repeat. This mode of recognition, the so-called PPR code, offers opportunities for the prediction of native PPR binding sites and the design of proteins to bind specified RNAs. However, a deep understanding of the parameters that dictate the affinity and specificity of PPR-RNA interactions is necessary to realize these goals. We report a comprehensive analysis of the sequence specificity of PPR10, a protein that binds similar RNA sequences of ∼18 nucleotides (nt) near the chloroplast atpH and psaJ genes in maize. We assessed the contribution of each nucleotide in the atpH binding site to PPR10 affinity in vitro by analyzing the effects of single-nucleotide changes at each position. In a complementary approach, the RNAs bound by PPR10 from partially randomized RNA pools were analyzed by deep sequencing. The results revealed three patches in which nucleotide identity has a major impact on binding affinity. These include 5 nt for which protein contacts were not observed in a PPR10-RNA crystal structure and 4 nt that are not explained by current views of the PPR code. These findings highlight aspects of PPR-RNA interactions that pose challenges for binding site prediction and design.

Science Educational Outreach Programs That Benefit Students and Scientists.

Both scientists and the public would benefit from improved communication of basic scientific research and from integrating scientists into education outreach, but opportunities to support these efforts are limited. We have developed two low-cost programs--"Present Your PhD Thesis to a 12-Year-Old" and "Shadow a Scientist"--that combine training in science communication with outreach to area middle schools. We assessed the outcomes of these programs and found a 2-fold benefit: scientists improve their communication skills by explaining basic science research to a general audience, and students' enthusiasm for science and their scientific knowledge are increased. Here we present details about both programs, along with our assessment of them, and discuss the feasibility of exporting these programs to other universities.

Ancient and modern colonization of North America by hemlock woolly adelgid, Adelges tsugae (Hemiptera: Adelgidae), an invasive insect from East Asia.

Hemlock woolly adelgid, Adelges tsugae, is an invasive pest of hemlock trees (Tsuga) in eastern North America. We used 14 microsatellites and mitochondrial COI sequences to assess its worldwide genetic structure and reconstruct its colonization history. The resulting information about its life cycle, biogeography and host specialization could help predict invasion by insect herbivores. We identified eight endemic lineages of hemlock adelgids in central China, western China, Ulleung Island (South Korea), western North America, and two each in Taiwan and Japan, with the Japanese lineages specializing on different Tsuga species. Adelgid life cycles varied at local and continental scales with different sexual, obligately asexual and facultatively asexual lineages. Adelgids in western North America exhibited very high microsatellite heterozygosity, which suggests ancient asexuality. The earliest lineages diverged in Asia during Pleistocene glacial periods, as estimated using approximate Bayesian computation. Colonization of western North America was estimated to have occurred prior to the last glacial period by adelgids directly ancestral to those in southern Japan, perhaps carried by birds. The modern invasion from southern Japan to eastern North America caused an extreme genetic bottleneck with just two closely related clones detected throughout the introduced range. Both colonization events to North America involved host shifts to unrelated hemlock species. These results suggest that genetic diversity, host specialization and host phylogeny are not predictive of adelgid invasion. Monitoring non-native sentinel host trees and focusing on invasion pathways might be more effective methods of preventing invasion than making predictions using species traits or evolutionary history.

Assessing the Ability of a Large Language Model to Score Free-Text Medical Student Clinical Notes: Quantitative Study.

Background: Teaching medical students the skills required to acquire, interpret, apply, and communicate clinical information is an integral part of medical education. A crucial aspect of this process involves providing students with feedback regarding the quality of their free-text clinical notes. Objective: The goal of this study was to assess the ability of ChatGPT 3.5, a large language model, to score medical students' free-text history and physical notes. Methods: This is a single-institution, retrospective study. Standardized patients learned a prespecified clinical case and, acting as the patient, interacted with medical students. Each student wrote a free-text history and physical note of their interaction. The students' notes were scored independently by the standardized patients and ChatGPT using a prespecified scoring rubric that consisted of 85 case elements. The measure of accuracy was percent correct. Results: The study population consisted of 168 first-year medical students. There was a total of 14,280 scores. The ChatGPT incorrect scoring rate was 1.0%, and the standardized patient incorrect scoring rate was 7.2%. The ChatGPT error rate was 86%, lower than the standardized patient error rate. The ChatGPT mean incorrect scoring rate of 12 (SD 11) was significantly lower than the standardized patient mean incorrect scoring rate of 85 (SD 74; P=.002). Conclusions: ChatGPT demonstrated a significantly lower error rate compared to standardized patients. This is the first study to assess the ability of a generative pretrained transformer (GPT) program to score medical students' standardized patient-based free-text clinical notes. It is expected that, in the near future, large language models will provide real-time feedback to practicing physicians regarding their free-text notes. GPT artificial intelligence programs represent an important advance in medical education and medical practice.

Atopic Dermatitis Complicated by Recurrent Eczema Herpeticum Is Characterized by Multiple, Concurrent Epidermal Inflammatory Endotypes.

A subgroup of patients with atopic dermatitis (AD) suffers from recurrent, disseminated herpes simplex virus skin infection, termed eczema herpeticum. To determine the transcriptional mechanisms of the skin and immune system pathobiology that underlie development of AD with eczema herpeticum (ADEH), we performed RNA-sequencing analysis of nonlesional skin (epidermis, dermis) from AD patients with and without a history of ADEH (ADEH+, n = 15; ADEH-, n = 13) along with healthy controls (n = 15). We also performed RNA sequencing on participants' plasmacytoid dendritic cells infected in vitro with herpes simplex virus 1. ADEH+ patients exhibited dysregulated gene expression, limited in the dermis (14 differentially expressed genes) and more widespread in the epidermis (129 differentially expressed genes). ADEH+-upregulated epidermal differentially expressed genes were enriched in type 2 cytokine (IL4R , CCL22, CRLF2, IL7R), interferon (CXCL10, ICAM1, IFI44, IRF7), and IL-36γ (IL36G) inflammatory gene pathways. All ADEH+ participants exhibited type 2 cytokine and inteferon endotypes, and 87% were IL36G-high. In contrast, these endotypes were more variably expressed among ADEH- participants. ADEH+ skin also had dysregulated epidermal differentiation complex gene expression of the late-cornified envelope, S100A, and small proline-rich gene families, which are involved in skin barrier function and antimicrobial activities. Plasmacytoid dendritic cell transcriptional responses to herpes simplex virus 1 infection were unaltered by ADEH status. The study concluded that the pathobiology underlying ADEH+ risk is associated with a unique, multifaceted epidermal inflammation that accompanies dysregulation of epidermal differentiation complex genes. These findings will help direct future studies that define how these inflammatory patterns may drive risk of eczema herpeticum in AD.

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma.

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

Analysis of patient claims data to determine the prevalence of hidradenitis suppurativa in the United States.

BACKGROUND: Recent prevalence estimates for hidradenitis suppurativa (HS), a chronic inflammatory skin condition, are limited by timeliness, population size, and generalizability. OBJECTIVE: We sought to develop prevalence estimates for HS in the United States using large health care claims databases. METHODS: A retrospective analysis used PharMetrics Integrated Database to gather health care claims information for HS among patients with 12 or more months of continuous enrollment in a commercial health care plan throughout 2007. Included patients had: 1 or more diagnoses with International Classification of Diseases, Ninth Revision, Clinical Modification code 705.83 for HS during 2007 without a Current Procedural Terminology code for HS; 1 or more Current Procedural Terminology codes of 11450, 11451, 11462, 11463, 11470, or 11471 during 2007 without International Classification of Diseases, Ninth Revision, Clinical Modification code 705.83; or both. Age- and gender-specific prevalence projections were calculated. RESULTS: Among included patients (n = 7927), mean age (SD) was 38.2 (14.73) years, and 5834 (74%) were women. Most patients (n = 5205; 66%) were aged 30 to 64 years. The overall prevalence estimate was 0.053% (95% confidence interval 0.051-0.054). When adjusted for gender and age, prevalence rates were 0.052% and 0.051%, respectively. The most common procedures for HS were excision of skin and subcutaneous tissue axillary/inguinal simple or intermediate repair. LIMITATIONS: Limitations were a health insured-only population; 12-month enrollment period for 2007; HS-specific procedural codes; and possible HS misclassifications. CONCLUSION: We found a low rate of clinically detected HS (0.053%; approximately 146,000-162,000 patients in the United States in 2007), with affected persons almost 3 times as likely to be female and the highest prevalence in those aged 18 to 44 years.

Atopic dermatitis complicated by recurrent eczema herpeticum is characterized by multiple, concurrent epidermal inflammatory endotypes.

BACKGROUND: A subgroup of atopic dermatitis (AD) patients suffer from recurrent, disseminated herpes simplex virus (HSV) skin infections, termed eczema herpeticum (EH), which can be life-threatening and contribute to AD morbidity. The pathobiology underlying ADEH is unknown. OBJECTIVE: To determine transcriptional mechanisms of skin and immune system pathobiology that underlie ADEH disease. METHODS: We performed whole transcriptome RNA-sequencing of non-lesional skin samples (epidermis, dermis) of AD patients with (ADEH + , n=15) and without (ADEH - , n=13) recurrent EH history, and healthy controls (HC, n=15). We also performed RNA-sequencing on plasmacytoid dendritic cells (pDCs) collected from these participants and infected in vitro with HSV-1. Differential expression, gene set enrichment, and endotyping analyses were performed. RESULTS: ADEH + disease was characterized by dysregulation in skin gene expression, which was limited in dermis (differentially expressed genes [DEGs]=14) and widespread in epidermis (DEGs=129). ADEH + -upregulated epidermal DEGs were enriched in type 2 cytokine (T2) ( IL4R, CCL22, CRLF2, IL7R ), interferon ( CXCL10, ICAM1, IFI44 , and IRF7) , and IL-36γ ( IL36G ) inflammatory pathway genes. At a person-level, all ADEH + participants exhibited T2 and interferon endotypes and 87% were IL36G-high. In contrast, these endotypes were more variably expressed among ADEH - participants. ADEH + patient skin also exhibited dysregulation in epidermal differentiation complex (EDC) genes within the LCE, S100 , and SPRR families, which are involved in skin barrier function, inflammation, and antimicrobial activities. pDC transcriptional responses to HSV-1 infection were not altered by ADEH status. CONCLUSIONS: ADEH + pathobiology is characterized by a unique, multi-faceted epidermal inflammation that accompanies dysregulation in the expression of EDC genes. Key Messages: AD patients with a history of recurrent EH exhibit molecular skin pathobiology that is similar in form, but more severe in degree, than in AD patients without this complication. Non-lesional skin of ADEH + patients concurrently exhibits excessive type 2 cytokine, interferon, and IL-36γ-driven epidermal inflammation. Expression of these inflammatory skin endotypes among ADEH + patients is associated with dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity. Capsule Summary: AD patients with a history of recurrent disseminated HSV-1 skin infections form a unique molecular skin endotype group that concurrently exhibits type 2 cytokine, interferon, and IL-36γ-driven skin inflammation, accompanied by dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity.

Rewriting regulatory DNA to dissect and reprogram gene expression.

Regulatory DNA sequences within enhancers and promoters bind transcription factors to encode cell type-specific patterns of gene expression. However, the regulatory effects and programmability of such DNA sequences remain difficult to map or predict because we have lacked scalable methods to precisely edit regulatory DNA and quantify the effects in an endogenous genomic context. Here we present an approach to measure the quantitative effects of hundreds of designed DNA sequence variants on gene expression, by combining pooled CRISPR prime editing with RNA fluorescence in situ hybridization and cell sorting (Variant-FlowFISH). We apply this method to mutagenize and rewrite regulatory DNA sequences in an enhancer and the promoter of PPIF in two immune cell lines. Of 672 variant-cell type pairs, we identify 497 that affect PPIF expression. These variants appear to act through a variety of mechanisms including disruption or optimization of existing transcription factor binding sites, as well as creation of de novo sites. Disrupting a single endogenous transcription factor binding site often led to large changes in expression (up to -40% in the enhancer, and -50% in the promoter). The same variant often had different effects across cell types and states, demonstrating a highly tunable regulatory landscape. We use these data to benchmark performance of sequence-based predictive models of gene regulation, and find that certain types of variants are not accurately predicted by existing models. Finally, we computationally design 185 small sequence variants (≤10 bp) and optimize them for specific effects on expression in silico. 84% of these rationally designed edits showed the intended direction of effect, and some had dramatic effects on expression (-100% to +202%). Variant-FlowFISH thus provides a powerful tool to map the effects of variants and transcription factor binding sites on gene expression, test and improve computational models of gene regulation, and reprogram regulatory DNA.

Fluorescence tracking of dissolved and particulate organic matter quality in a river-dominated estuary.

Excitation-emission matrix (EEM) fluorescence was combined with parallel factor analysis (PARAFAC) to model base-extracted particulate (POM) and dissolved (DOM) organic matter quality in the Neuse River Estuary (NRE), North Carolina, before and after passage of Hurricane Irene in August 2011. Principle components analysis was used to determine that four of the PARAFAC components (C1-C3 and C6) were terrestrial sources to the NRE. One component (C4), prevalent in DOM of nutrient-impacted streams and estuaries and produced in phytoplankton cultures, was enriched in the POM and in surface sediment pore water DOM. One component (C5) was related to recent autochthonous production. Photoexposure of unfiltered Neuse River water caused an increase in slope ratio values (S(R)) which corresponded to an increase in the ratio C2:C3 for DOM, and the production of C4 fluorescence in both POM and DOM. Changes to the relative abundance of C4 in POM and DOM indicated that advection of pore water DOM from surface sediments into overlying waters could increase the autochthonous quality of DOM in shallow microtidal estuaries. Modeling POM and DOM simultaneously with PARAFAC is an informative technique that is applicable to assessments of estuarine water quality.

IL-13-programmed airway tuft cells produce PGE2, which promotes CFTR-dependent mucociliary function.

Chronic type 2 (T2) inflammatory diseases of the respiratory tract are characterized by mucus overproduction and disordered mucociliary function, which are largely attributed to the effects of IL-13 on common epithelial cell types (mucus secretory and ciliated cells). The role of rare cells in airway T2 inflammation is less clear, though tuft cells have been shown to be critical in the initiation of T2 immunity in the intestine. Using bulk and single-cell RNA sequencing of airway epithelium and mouse modeling, we found that IL-13 expanded and programmed airway tuft cells toward eicosanoid metabolism and that tuft cell deficiency led to a reduction in airway prostaglandin E2 (PGE2) concentration. Allergic airway epithelia bore a signature of PGE2 activation, and PGE2 activation led to cystic fibrosis transmembrane receptor-dependent ion and fluid secretion and accelerated mucociliary transport. These data reveal a role for tuft cells in regulating epithelial mucociliary function in the allergic airway.

Infliximab therapy for patients with moderate to severe hidradenitis suppurativa: a randomized, double-blind, placebo-controlled crossover trial.

BACKGROUND: Biologic therapies with anti-tumor necrosis factor agents are promising treatments for hidradenitis suppurativa (HS). OBJECTIVE: We assessed the efficacy and safety of infliximab (IFX) for the treatment of moderate to severe HS. METHODS: A prospective double-blind treatment phase of 8 weeks where patients received IFX or placebo was followed by an open-label phase where patients taking placebo were given the opportunity to cross over to IFX, and an observational phase. Primary treatment efficacy was based on HS Severity Index. Secondary end points included Dermatology Life Quality Index, visual analog scale, and Physician Global Assessment scores. Inflammatory markers erythrocyte sedimentation rate and C-reactive protein were also assessed. RESULTS: More patients in the IFX than in the placebo group showed a 50% or greater decrease from baseline HS Severity Index score. In addition, statistically and clinically significant improvement from baseline was observed at week 8 in Dermatology Life Quality Index score, visual analog scale score, erythrocyte sedimentation rate, and C-reactive protein compared with placebo. Patients in the placebo group treated with IFX after week 8 (crossover) responded similarly to the original IFX group. Many patients withdrew during the observational phase to continue anti-tumor necrosis factor-alfa therapy. No unexpected serious adverse events were observed. LIMITATIONS: Results are representative of a single center, patients were treated by a single physician, some patients did not return after their last infusion, and the HS Severity Index requires validation. CONCLUSIONS: This clinical study represents the first formal assessment of IFX for treatment of moderate to severe HS. IFX was well tolerated, no unexpected safety issues were identified, and improvements in pain intensity, disease severity, and quality of life were demonstrated with concomitant reduction in clinical markers of inflammation.

PAH mineralization and bacterial organotolerance in surface sediments of the Charleston Harbor estuary.

Semi-volatile organic compounds (SVOCs) in estuarine waters can adversely affect biota but watershed sources can be difficult to identify because these compounds are transient. Natural bacterial assemblages may respond to chronic, episodic exposure to SVOCs through selection of more organotolerant bacterial communities. We measured bacterial production, organotolerance and polycyclic aromatic hydrocarbon (PAH) mineralization in Charleston Harbor and compared surface sediment from stations near a known, permitted SVOC outfall (pulp mill effluent) to that from more pristine stations. Naphthalene additions inhibited an average of 77% of bacterial metabolism in sediments from the more pristine site (Wando River). Production in sediments nearest the outfall was only inhibited an average of 9% and in some cases, was actually stimulated. In general, the stations with the highest rates of bacterial production also were among those with the highest rates of PAH mineralization. This suggests that the capacity to mineralize PAH carbon is a common feature amongst the bacterial assemblage in these estuarine sediments and could account for an average of 5.6% of bacterial carbon demand (in terms of production) in the summer, 3.3% in the spring (April) and only 1.2% in winter (December).

Ecosystem Capacity for Microbial Biodegradation of Munitions Compounds and Phenanthrene in Three Coastal Waterways in North Carolina, United States.

Munitions compounds (i.e., 2,4,6-trinitrotoluene (TNT), octahy-dro-1,3,5,7-tetranitro-1,3,5,7-tetrazocin (HMX), and hexadydro-1,3,5-trinitro-1,3,5-triazin (RDX), also called energetics) were originally believed to be recalcitrant to microbial biodegradation based on historical groundwater chemical attenuation data and laboratory culture work. More recently, it has been established that natural bacterial assemblages in coastal waters and sediment can rapidly metabolize these organic nitrogen sources and even incorporate their carbon and nitrogen into bacterial biomass. Here, we report on the capacity of natural microbial assemblages in three coastal North Carolina (United States) estuaries to metabolize energetics and phenanthrene (PHE), a proxy for terrestrial aromatic compounds. Microbial assemblages generally had the highest ecosystem capacity (mass of the compound mineralized per average estuarine residence time) for HMX (21-5463 kg) > RDX (1.4-5821 kg) ≫ PHE (0.29-660 kg) > TNT (0.25-451 kg). Increasing antecedent precipitation tended to decrease the ecosystem capacity to mineralize TNT in the Newport River Estuary, and PHE and TNT mineralization were often highest with increasing salinity. There was some evidence from the New River Estuary that increased N-demand (due to a phytoplankton bloom) is associated with increased energetic mineralization rates. Using this type of analysis to determine the ecosystem capacity to metabolize energetics can explain why these compounds are rarely detected in seawater and marine sediment, despite the known presence of unexploded ordnance or recent use in military training exercises. Overall, measuring the ecosystem capacity may help predict the effects of climate change (warming and altered precipitation patterns) and other perturbations on exotic compound fate and transport within ecosystems and provide critical information for managers and decision-makers to develop management strategies based on these changes.

Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium.

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.

Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium.

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.

Type 2 and interferon inflammation strongly regulate SARS-CoV-2 related gene expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) outcomes vary from asymptomatic infection to death. This disparity may reflect different airway levels of the SARS-CoV-2 receptor, ACE2, and the spike protein activator, TMPRSS2. 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 (eQTL) for both ACE2 and TMPRSS2, that vary in frequency across world populations. Importantly, we find TMPRSS2 is part of a mucus secretory network, highly upregulated by T2 inflammation through the action of interleukin-13, and that interferon response to respiratory viruses highly upregulates ACE2 expression. Finally, we define airway responses to coronavirus infections in children, finding that these infections upregulate IL6 while also stimulating a more pronounced cytotoxic immune response relative to other respiratory viruses. Our results reveal mechanisms likely influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.