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.

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.

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.

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.

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.

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.

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.

MicroRNA-7a overexpression in VMH restores the sympathoadrenal response to hypoglycemia.

It is proposed that the impaired sympathoadrenal response to hypoglycemia induced by recurrent insulin-induced hypoglycemia (RH) is an adaptive phenomenon induced by specific changes in microRNA expression in the ventromedial hypothalamus (VMH). To test this hypothesis, genome-wide microRNAomic profiling of the VMH by RNA-sequencing was performed in control rats and rats treated for RH. Differential expression analysis identified microRNA-7a-5p and microRNA-665 as potential mediators of this phenomenon. To further test this hypothesis, experiments were conducted consisting of targeted lentiviral-mediated overexpression of microRNA-7a-5p and downregulation of microRNA-665 in the VMH. Hyperinsulinemic hypoglycemic clamp experiments demonstrated that targeted overexpression of microRNA-7a-5p (but not downregulation of microRNA-665) in the VMH of RH rats restored the epinephrine response to hypoglycemia. This restored response to hypoglycemia was associated with a restoration of GABAA receptor gene expression. Finally, a direct interaction of microRNA-7a-5p with the 3'-UTR of GABAA receptor α1-subunit (Gabra1) gene was demonstrated in a luciferase assay. These findings indicate that (a) the impaired sympathoadrenal response RH induces is associated with changes in VMH microRNA expression and (b) microRNA-7a-5p, possibly via direct downregulation of GABA receptor gene expression, may serve as a mediator of the altered sympathoadrenal response to hypoglycemia.

The nonlesional skin surface distinguishes atopic dermatitis with food allergy as a unique endotype.

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.

Practical Measures of Clinical Benefit With Ruxolitinib Therapy: An Exploratory Analysis of COMFORT-I.

BACKGROUND: The phase III COMFORT (Controlled Myelofibrosis Study With Oral JAK inhibitor Treatment)-I and COMFORT-II trials in patients with intermediate-2 or high-risk myelofibrosis (MF) showed that ruxolitinib was superior to placebo and best available therapy, respectively, for improvements in spleen volume, MF-related symptoms, and overall survival (OS). However, patients managed in community settings might not have access to the methods used in the COMFORT trials. In this exploratory analysis we summarize efficacy findings of COMFORT-I using practical, community-oriented measures of patient outcomes. PATIENTS AND METHODS: In this post hoc analysis of data from COMFORT-I we evaluated changes from baseline to week 12 in spleen size (palpable length and volume), patient-reported outcomes (Patient Global Impression of Change; Myelofibrosis Symptom Assessment Form; Patient-Reported Outcomes Measurement System Fatigue Scale), body weight, and serum albumin levels in 5 subgroups of ruxolitinib-treated patients on the basis of week 12 spleen length changes from baseline: (1-4) ≥ 50%, 25% to < 50%, 10% to < 25%, or < 10% reduction; and (5) worsening. OS was evaluated in ruxolitinib-treated patients with week 12 spleen length reductions from baseline ≥ 50%, 25% to < 50%, or < 25% (including worsening). RESULTS: In all spleen length subgroups, including patients with worsening spleen length at week 12, ruxolitinib (n = 150) was associated with improvements in spleen volume, patient-reported symptom burden, body weight, and serum albumin levels. Greater reductions in spleen length were associated with prolonged OS. CONCLUSION: A variety of assessment methods beyond palpable spleen length that are easily accessible in the community setting might be useful in evaluating the clinical benefit of ruxolitinib over time in patients with MF.

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.

Measuring Carbon-based Contaminant Mineralization Using Combined CO2 Flux and Radiocarbon Analyses.

A method is described which uses the absence of radiocarbon in industrial chemicals and fuels made from petroleum feedstocks which frequently contaminate the environment. This radiocarbon signal - or rather the absence of signal - is evenly distributed throughout a contaminant source pool (unlike an added tracer) and is not impacted by biological, chemical or physical processes (e.g., the 14C radioactive decay rate is immutable). If the fossil-derived contaminant is fully degraded to CO2, a harmless end-product, that CO2 will contain no radiocarbon. CO2 derived from natural organic matter (NOM) degradation will reflect the NOM radiocarbon content (usually <30,000 years old). Given a known radiocarbon content for NOM (a site background), a two end-member mixing model can be used to determine the CO2 derived from a fossil source in a given soil gas or groundwater sample. Coupling the percent CO2 derived from the contaminant with the CO2 respiration rate provides an estimate for the total amount of contaminant degraded per unit time. Finally, determining a zone of influence (ZOI) representing the volume from which site CO2 is collected allows determining the contaminant degradation per unit time and volume. Along with estimates for total contaminant mass, this can ultimately be used to calculate time-to-remediate or otherwise used by site managers for decision-making.

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.

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.