Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop.

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

Human Epidemiology and RespOnse to SARS-CoV-2 (HEROS): Objectives, Design and Enrollment Results of a 12-City Remote Observational Surveillance Study of Households with Children using Direct-to-Participant Methods.

The Human Epidemiology and Response to SARS-CoV-2 (HEROS) is a prospective multi-city 6-month incidence study which was conducted from May 2020-February 2021. The objectives were to identify risk factors for SARS-CoV-2 infection and household transmission among children and people with asthma and allergic diseases, and to use the host nasal transcriptome sampled longitudinally to understand infection risk and sequelae at the molecular level. To overcome challenges of clinical study implementation due to the coronavirus pandemic, this surveillance study used direct-to-participant methods to remotely enroll and prospectively follow eligible children who are participants in other NIH-funded pediatric research studies and their household members. Households participated in weekly surveys and biweekly nasal sampling regardless of symptoms. The aim of this report is to widely share the methods and study instruments and to describe the rationale, design, execution, logistics and characteristics of a large, observational, household-based, remote cohort study of SARS-CoV-2 infection and transmission in households with children. The study enrolled a total of 5,598 individuals, including 1,913 principal participants (children), 1,913 primary caregivers, 729 secondary caregivers and 1,043 other household children. This study was successfully implemented without necessitating any in-person research visits and provides an approach for rapid execution of clinical research.

IL-13-associated epithelial remodeling correlates with clinical severity in nasal polyposis.

BACKGROUND: Epithelial remodeling is a histopathologic feature of chronic inflammatory airway diseases including chronic rhinosinusitis (CRS). Cell-type shifts and their relationship to CRS endotypes and severity are incompletely described. OBJECTIVE: We sought to understand the relationship of epithelial cell remodeling to inflammatory endotypes and disease outcomes in CRS. METHODS: Using cell-type transcriptional signatures derived from epithelial single-cell sequencing, we analyzed bulk RNA-sequencing data from sinus epithelial brushings obtained from patients with CRS with and without nasal polyps in comparison to healthy controls. RESULTS: The airway epithelium in nasal polyposis displayed increased tuft cell transcripts and decreased ciliated cell transcripts along with an IL-13 activation signature. In contrast, CRS without polyps showed an IL-17 activation signature. IL-13 activation scores were associated with increased tuft cell, goblet cell, and mast cell scores and decreased ciliated cell scores. Furthermore, the IL-13 score was strongly associated with a previously reported activated ("polyp") tuft cell score and a prostaglandin E2 activation signature. The Lund-Mackay score, a computed tomographic metric of sinus opacification, correlated positively with activated tuft cell, mast cell, prostaglandin E2, and IL-13 signatures and negatively with ciliated cell transcriptional signatures. CONCLUSIONS: These results demonstrate that cell-type alterations and prostaglandin E2 stimulation are key components of IL-13-induced epithelial remodeling in nasal polyposis, whereas IL-17 signaling is more prominent in CRS without polyps, and that clinical severity correlates with the degree of IL-13-driven epithelial remodeling.

CCL5 is a potential bridge between type 1 and type 2 inflammation in asthma.

BACKGROUND: Type 1 (T1) inflammation (marked by IFN-γ expression) is now consistently identified in subsets of asthma cohorts, but how it contributes to disease remains unclear. OBJECTIVE: We sought to understand the role of CCL5 in asthmatic T1 inflammation and how it interacts with both T1 and type 2 (T2) inflammation. METHODS: CCL5, CXCL9, and CXCL10 messenger RNA expression from sputum bulk RNA sequencing, as well as clinical and inflammatory data were obtained from the Severe Asthma Research Program III (SARP III). CCL5 and IFNG expression from bronchoalveolar lavage cell bulk RNA sequencing was obtained from the Immune Mechanisms in Severe Asthma (IMSA) cohort and expression related to previously identified immune cell profiles. The role of CCL5 in tissue-resident memory T-cell (TRM) reactivation was evaluated in a T1high murine severe asthma model. RESULTS: Sputum CCL5 expression strongly correlated with T1 chemokines (P < .001 for CXCL9 and CXCL10), consistent with a role in T1 inflammation. CCL5high participants had greater fractional exhaled nitric oxide (P = .009), blood eosinophils (P < .001), and sputum eosinophils (P = .001) in addition to sputum neutrophils (P = .001). Increased CCL5 bronchoalveolar lavage expression was unique to a previously described T1high/T2variable/lymphocytic patient group in the IMSA cohort, with IFNG trending with worsening lung obstruction only in this group (P = .083). In a murine model, high expression of the CCL5 receptor CCR5 was observed in TRMs and was consistent with a T1 signature. A role for CCL5 in TRM activation was supported by the ability of the CCR5 inhibitor maraviroc to blunt reactivation. CONCLUSION: CCL5 appears to contribute to TRM-related T1 neutrophilic inflammation in asthma while paradoxically also correlating with T2 inflammation and with sputum eosinophilia.

Single cell analysis of host response to helminth infection reveals the clonal breadth, heterogeneity, and tissue-specific programming of the responding CD4+ T cell repertoire.

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.

Possibilities and promise: Leveraging advances in transcriptomics for clinical decision making in allergic diseases.

Transcriptomics has revolutionized our understanding of the pathobiologic heterogeneity underlying complex allergic diseases, leading to both the discovery of multiple inflammatory allergic disease endotypes and the development of targeted biologic therapies. In addition, transcriptomic endotypes have been associated with disease severity, exacerbation propensity, and responsiveness to nontargeted therapies, suggesting an unrealized potential for transcriptomic assays and endotyping to be used as diagnostic, predictive, and prognostic biomarkers. In this review, we discuss current and emerging transcriptomic technologies and how they have been used to uncover allergic disease endotypes and generate other clinically relevant findings. We then discuss how transcriptomics could be leveraged in the clinic for the delivery of personalized, rational patient care, and we consider strategies for and impediments to their future clinical deployment.

Risk factors for SARS-CoV-2 infection and transmission in households with children with asthma and allergy: A prospective surveillance study.

BACKGROUND: Whether children and people with asthma and allergic diseases are at increased risk for severe acute respiratory syndrome virus 2 (SARS-CoV-2) infection is unknown. OBJECTIVE: Our aims were to determine the incidence of SARS-CoV-2 infection in households with children and to also determine whether self-reported asthma and/or other allergic diseases are associated with infection and household transmission. METHODS: For 6 months, biweekly nasal swabs and weekly surveys were conducted within 1394 households (N = 4142 participants) to identify incident SARS-CoV-2 infections from May 2020 to February 2021, which was the pandemic period largely before a vaccine and before the emergence of SARS-CoV-2 variants. Participant and household infection and household transmission probabilities were calculated by using time-to-event analyses, and factors associated with infection and transmission risk were determined by using regression analyses. RESULTS: In all, 147 households (261 participants) tested positive for SARS-CoV-2. The household SARS-CoV-2 infection probability was 25.8%; the participant infection probability was similar for children (14.0% [95% CI = 8.0%-19.6%]), teenagers (12.1% [95% CI = 8.2%-15.9%]), and adults (14.0% [95% CI = 9.5%-18.4%]). Infections were symptomatic in 24.5% of children, 41.2% of teenagers, and 62.5% of adults. Self-reported doctor-diagnosed asthma was not a risk factor for infection (adjusted hazard ratio [aHR] = 1.04 [95% CI = 0.73-1.46]), nor was upper respiratory allergy or eczema. Self-reported doctor-diagnosed food allergy was associated with lower infection risk (aHR = 0.50 [95% CI = 0.32-0.81]); higher body mass index was associated with increased infection risk (aHR per 10-point increase = 1.09 [95% CI = 1.03-1.15]). The household secondary attack rate was 57.7%. Asthma was not associated with household transmission, but transmission was lower in households with food allergy (adjusted odds ratio = 0.43 [95% CI = 0.19-0.96]; P = .04). CONCLUSION: Asthma does not increase the risk of SARS-CoV-2 infection. Food allergy is associated with lower infection risk, whereas body mass index is associated with increased infection risk. Understanding how these factors modify infection risk may offer new avenues for preventing infection.

Influenza virus infection increases ACE2 expression and shedding in human small airway epithelial cells.

BACKGROUND: Patients with coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demonstrate high rates of co-infection with respiratory viruses, including influenza A (IAV), suggesting pathogenic interactions. METHODS: We investigated how IAV may increase the risk of COVID-19 lung disease, focusing on the receptor angiotensin-converting enzyme (ACE)2 and the protease TMPRSS2, which cooperate in the intracellular uptake of SARS-CoV-2. RESULTS: We found, using single-cell RNA sequencing of distal human nondiseased lung homogenates, that at baseline, ACE2 is minimally expressed in basal, goblet, ciliated and secretory epithelial cells populating small airways. We focused on human small airway epithelial cells (SAECs), central to the pathogenesis of lung injury following viral infections. Primary SAECs from nondiseased donor lungs apically infected (at the air-liquid interface) with IAV (up to 3×105 pfu; ∼1 multiplicity of infection) markedly (eight-fold) boosted the expression of ACE2, paralleling that of STAT1, a transcription factor activated by viruses. IAV increased the apparent electrophoretic mobility of intracellular ACE2 and generated an ACE2 fragment (90 kDa) in apical secretions, suggesting cleavage of this receptor. In addition, IAV increased the expression of two proteases known to cleave ACE2, sheddase ADAM17 (TACE) and TMPRSS2 and increased the TMPRSS2 zymogen and its mature fragments, implicating proteolytic autoactivation. CONCLUSION: These results indicate that IAV amplifies the expression of molecules necessary for SARS-CoV-2 infection of the distal lung. Furthermore, post-translational changes in ACE2 by IAV may increase vulnerability to lung injury such as acute respiratory distress syndrome during viral co-infections. These findings support efforts in the prevention and treatment of influenza infections during the COVID-19 pandemic.

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.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

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.

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy.

BACKGROUND: Global gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context. METHODS: We used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age. RESULTS: The SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10- 9), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10- 8). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells. CONCLUSION: Focusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids.

A Transcriptomic Method to Determine Airway Immune Dysfunction in T2-High and T2-Low Asthma.

BACKGROUND: Type 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. OBJECTIVES: To determine the gene networks that sustain T2 airway inflammation in T2-high asthma and to explore the gene networks that characterize T2-low asthma. METHODS: Network 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. RESULTS: Sputum 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. CONCLUSIONS: Persistent 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.

Expression and function of the ectopic olfactory receptor OR10G7 in patients with atopic dermatitis.

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.

Functional genomics of CDHR3 confirms its role in HRV-C infection and childhood asthma exacerbations.

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.

An admixture mapping meta-analysis implicates genetic variation at 18q21 with asthma susceptibility in Latinos.

BACKGROUND: Asthma is a common but complex disease with racial/ethnic differences in prevalence, morbidity, and response to therapies. OBJECTIVE: We sought to perform an analysis of genetic ancestry to identify new loci that contribute to asthma susceptibility. METHODS: We leveraged the mixed ancestry of 3902 Latinos and performed an admixture mapping meta-analysis for asthma susceptibility. We replicated associations in an independent study of 3774 Latinos, performed targeted sequencing for fine mapping, and tested for disease correlations with gene expression in the whole blood of more than 500 subjects from 3 racial/ethnic groups. RESULTS: We identified a genome-wide significant admixture mapping peak at 18q21 in Latinos (P = 6.8 × 10-6), where Native American ancestry was associated with increased risk of asthma (odds ratio [OR], 1.20; 95% CI, 1.07-1.34; P = .002) and European ancestry was associated with protection (OR, 0.86; 95% CI, 0.77-0.96; P = .008). Our findings were replicated in an independent childhood asthma study in Latinos (P = 5.3 × 10-3, combined P = 2.6 × 10-7). Fine mapping of 18q21 in 1978 Latinos identified a significant association with multiple variants 5' of SMAD family member 2 (SMAD2) in Mexicans, whereas a single rare variant in the same window was the top association in Puerto Ricans. Low versus high SMAD2 blood expression was correlated with case status (13.4% lower expression; OR, 3.93; 95% CI, 2.12-7.28; P < .001). In addition, lower expression of SMAD2 was associated with more frequent exacerbations among Puerto Ricans with asthma. CONCLUSION: Ancestry at 18q21 was significantly associated with asthma in Latinos and implicated multiple ancestry-informative noncoding variants upstream of SMAD2 with asthma susceptibility. Furthermore, decreased SMAD2 expression in blood was strongly associated with increased asthma risk and increased exacerbations.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma.

Short-acting ß2-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the United States. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African-American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single-nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1830 Latinos (total = 2779). Finally, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p = 7.69 × 10-9). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958 and rs7081864, p ≤ 5 × 10-8). Our results failed to replicate in three additional populations of 416 Latinos and 1615 African Americans. Our findings indicate that both population-specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

Alternative splicing of interleukin-33 and type 2 inflammation in asthma.

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.

Minimally invasive skin tape strip RNA sequencing identifies novel characteristics of the type 2-high atopic dermatitis disease endotype.

BACKGROUND: Expression profiling of skin biopsy specimens has established molecular features of the skin in patients with atopic dermatitis (AD). The invasiveness of biopsies has prevented their use in defining individual-level AD pathobiological mechanisms (endotypes) in large research studies. OBJECTIVE: We sought to determine whether minimally invasive skin tape strip transcriptome analysis identifies gene expression dysregulation in AD and molecular disease endotypes. METHODS: We sampled nonlesional and lesional skin tape strips and biopsy specimens from white adult patients with AD (18 male and 12 female patients; age [mean ± SE], 36.3 ± 2.2 years) and healthy control subjects (9 male and 16 female subjects; age [mean ± SE], 34.8 ± 2.2 years). AmpliSeq whole-transcriptome sequencing was performed on extracted RNA. Differential expression, clustering/pathway analyses, immunostaining of skin biopsy specimens, 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 (false discovery rate < 5%) separated patients with AD from healthy subjects. Agnostic gene set analyses and clustering revealed 50% of patients with AD exhibited a type 2 inflammatory signature (type 2-high endotype) characterized by differential expression of 656 genes, including overexpression of IL13, IL4R, CCL22, CCR4 (log2 fold change = 5.5, 2.0, 4.0, and 4.1, respectively) and at a pathway level by TH2/dendritic cell activation. Both expression and immunostaining of skin biopsy specimens indicated this type 2-high group was enriched for inflammatory, type 2-skewed dendritic cells expressing FcεRI. The type 2-high endotype group exhibited more severe disease by using both the Eczema Area and Severity Index score and body surface area covered by lesions. CONCLUSION: Minimally invasive expression profiling of nonlesional skin reveals stratification in AD molecular pathology by type 2 inflammation that correlates with disease severity.