Cytokine-induced molecular responses in airway smooth muscle cells inform genome-wide association studies of asthma

Integrative Roles of Pro-Inflammatory Cytokines on Airway Smooth Muscle Structure and Function in Asthma

Asthma has become more appreciated for its heterogeneity with studies identifying type 2 and non-type 2 phenotypes/endotypes that ultimately lead to airflow obstruction, airway hyperresponsiveness, and remodeling. The pro-inflammatory environment in asthma influences airway smooth muscle (ASM) structure and function. ASM has a vast repertoire of inflammatory receptors that, upon activation, contribute to prominent features in asthma, notably immune cell recruitment and activation, hypercontractility, proliferation, migration, and extracellular matrix protein deposition. These pro-inflammatory responses in ASM can be mediated by both type 2 (e.g., IL-4, IL-13, and TSLP) and non-type 2 (e.g., TNFα, IFNγ, IL-17A, and TGFβ) cytokines, highlighting roles for ASM in type 2 and non-type 2 asthma phenotypes/endotypes. In recent years, there has been considerable advances in understanding how pro-inflammatory cytokines promote ASM dysfunction and impair responsiveness to asthma therapy, corticosteroids and long-acting β2-adrenergic receptor agonists (LABAs). Transcriptomic analyses on human ASM cells and tissues have expanded our knowledge in this area but have also raised new questions regarding ASM and its role in asthma. In this review, we discuss how pro-inflammatory cytokines, corticosteroids, and LABAs affect ASM structure and function, with particular focus on changes in gene expression and transcriptional programs in type 2 and non-type 2 asthma.

Targeting cytoskeletal biomechanics to modulate airway smooth muscle contraction in asthma

To contract, to deform, and remodel, the airway smooth muscle cell relies on dynamic changes in the structure of its mechanical force-bearing cytoskeleton. These alternate between a "fluid-like" (relaxed) state characterized by weak contractile protein-protein interactions within the cytoskeletal apparatus and a "solid-like" (contractile) state promoted by strong and highly organized molecular interactions. In this review, we discuss the roles for actin, myosin, factors promoting actin polymerization and depolymerization, adhesome complexes, and cell-cell junctions in these dynamic processes. We describe the relationship between these cytoskeletal factors, extracellular matrix components of bronchial tissue, and mechanical stretch and other changes within the airway wall in the context of the physical mechanisms of cytoskeletal fluidization-resolidification. We also highlight studies that emphasize the distinct processes of cell shortening and force transmission in airway smooth muscle and previously unrecognized roles for actin in cytoskeletal dynamics. Finally, we discuss the implications of these discoveries for understanding and treating airway obstruction in asthma.

Genetic contributions to epigenetic-defined endotypes of allergic phenotypes in children

Background

Asthma is the most common chronic respiratory disease in children, but little is known about genetic contributions to its underlying endotypes. To address this gap, we studied the methylome, transcriptome, and genome from children with extensive phenotyping from birth.

Methods

We performed DNA methylation (DNAm) studies using the Asthma&Allergy array and RNA-sequencing in nasal mucosal cells from 284 children (age 11 years) in the Urban Environment and Childhood Asthma (URECA) birth cohort with genotypes from whole-genome sequencing. Using empirical Bayes matrix factorization on all CpGs on the array, we derived 16 DNAm signatures and tested for associations between phenotypes and gene expression. We then replicated results in two additional cohorts and estimated the heritability of phenotype-associated signatures using single-nucleotide polymorphisms (SNPs) associated with an allergic disease, and with CpGs and genes associated with the signatures.

Findings

Three DNAm signatures were associated with at least one phenotype: allergic asthma, allergic rhinitis, allergic sensitization (atopy), total IgE, exhaled nitric oxide, or blood eosinophils. The genes correlated with each of the three signatures were enriched in networks reflecting inhibited immune response to microbes, impaired epithelial barrier integrity, and activated T2 immune pathways. We replicated the signature-phenotype associations in two additional birth cohorts. The estimated joint SNP heritabilities of the signatures were 0.17 (p=0.0027), 0.30 (p=9.3×10-7), and 0.16 (p=9.0×10-7), respectively.

Interpretation

We identified three significantly heritable DNAm signatures defining asthma and allergy endotypes across diverse populations. Our study demonstrated that epigenetic patterning in airway mucosal cells reflects perturbations in underlying biological processes related to the development of asthma and allergic diseases in childhood.

Funding

National Institute of Allergy and Infectious Diseases and the National Institutes of Health, Office of the Director.

Development of an epigenetic clock resistant to changes in immune cell composition

Epigenetic clocks are age predictors that use machine-learning models trained on DNA CpG methylation values to predict chronological or biological age. Increases in predicted epigenetic age relative to chronological age (epigenetic age acceleration) are connected to aging-associated pathologies, and changes in epigenetic age are linked to canonical aging hallmarks. However, epigenetic clocks rely on training data from bulk tissues whose cellular composition changes with age. Here, we found that human naive CD8+ T cells, which decrease in frequency during aging, exhibit an epigenetic age 15-20 years younger than effector memory CD8+ T cells from the same individual. Importantly, homogenous naive T cells isolated from individuals of different ages show a progressive increase in epigenetic age, indicating that current epigenetic clocks measure two independent variables, aging and immune cell composition. To isolate the age-associated cell intrinsic changes, we created an epigenetic clock, the IntrinClock, that did not change among 10 immune cell types tested. IntrinClock shows a robust predicted epigenetic age increase in a model of replicative senescence in vitro and age reversal during OSKM-mediated reprogramming.

Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents

Several single nucleotide polymorphisms (SNPs) in multiple interleukin receptor genes could be associated with asthma risk and/or phenotype. Interleukin-17 (IL-17) has been implicated in tissue inflammation and autoimmune diseases. As no previous studies have uncovered the potential role of IL17 receptor A (RA) gene variants in asthma risk, we aimed to explore the association of four IL17RA SNPs (i.e., rs4819554A/G, rs879577C/T, rs41323645G/A, and rs4819555C/T) with asthma susceptibility/phenotype in our region. TaqMan allelic discrimination analysis was used to genotype 192 individuals. We found that the rs4819554 G/G genotype significantly reduced disease risk in the codominant (OR = 0.15, 95%CI = 0.05-0.45, p < 0.001), dominant (OR = 0.49, 95%CI = 0.26-0.93, p = 0.028), and recessive (OR = 0.18, 95%CI = 0.07-0.52, p < 0.001) models. Similarly, rs879577 showed reduced disease risk associated with the T allele across all genetic models. However, the A allele of rs41323645 was associated with increased disease risk in all models. The G/A and A/A genotypes have higher ORs of 2.47 (95%CI = 1.19-5.14) and 3.86 (95%CI = 1.62-9.18), respectively. Similar trends are observed in the dominant 2.89 (95%CI = 1.47-5.68, p = 0.002) and recessive 2.34 (95%CI = 1.10-4.98, p = 0.025) models. For the rs4819555 variant, although there was no significant association identified under any models, carriers of the rs4819554*A demonstrated an association with a positive family history of asthma (71.4% in carriers vs. 27% in non-carriers; p = 0.025) and the use of relievers for >2 weeks (52.2% of carriers vs. 28.8% of non-carriers; p = 0.047). Meanwhile, the rs4819555*C carriers displayed a significant divergence in the asthma phenotype, specifically atopic asthma (83.3% vs. 61.1%; p = 0.007), showed a higher prevalence of chest tightness (88.9% vs. 61.5%; p = 0.029), and were more likely to report comorbidities (57.7% vs. 16.7%, p = 0.003). The most frequent haplotype in the asthma group was ACAC, with a frequency of 22.87% vs. 1.36% in the controls (p < 0.001). In conclusion, the studied IL17RA variants could be essential in asthma susceptibility and phenotype in children and adolescents.

A functional genomics pipeline to identify high-value asthma and allergy CpGs in the human methylome

BACKGROUND

DNA methylation of cytosines at cytosine-phosphate-guanine (CpG) dinucleotides (CpGs) is a widespread epigenetic mark, but genome-wide variation has been relatively unexplored due to the limited representation of variable CpGs on commercial high-throughput arrays.

OBJECTIVES

To explore this hidden portion of the epigenome, this study combined whole-genome bisulfite sequencing with in silico evidence of gene regulatory regions to design a custom array of high-value CpGs. This study focused on airway epithelial cells from children with and without allergic asthma because these cells mediate the effects of inhaled microbes, pollution, and allergens on asthma and allergic disease risk.

METHODS

This study identified differentially methylated regions from whole-genome bisulfite sequencing in nasal epithelial cell DNA from a total of 39 children with and without allergic asthma of both European and African ancestries. This study selected CpGs from differentially methylated regions, previous allergy or asthma epigenome-wide association studies (EWAS), or genome-wide association study loci, and overlapped them with functional annotations for inclusion on a custom Asthma&Allergy array. This study used both the custom and EPIC arrays to perform EWAS of allergic sensitization (AS) in nasal epithelial cell DNA from children in the URECA (Urban Environment and Childhood Asthma) birth cohort and using the custom array in the INSPIRE [Infant Susceptibility to Pulmonary Infections and Asthma Following RSV Exposure] birth cohort. Each CpG on the arrays was assigned to its nearest gene and its promotor capture Hi-C interacting gene and performed expression quantitative trait methylation (eQTM) studies for both sets of genes.

RESULTS

Custom array CpGs were enriched for intermediate methylation levels compared to EPIC CpGs. Intermediate methylation CpGs were further enriched among those associated with AS and for eQTMs on both arrays.

CONCLUSIONS

This study revealed signature features of high-value CpGs and evidence for epigenetic regulation of genes at AS EWAS loci that are robust to race/ethnicity, ascertainment, age, and geography.

MicroRNA-15b-5p inhibits tumor necrosis factor alpha-induced proliferation, migration, and extracellular matrix production of airway smooth muscle cells via targeting yes-associated protein 1

The excessive proliferation and the deposition of extracellular matrix (ECM) of airway smooth muscle (ASM) cells facilitates airway remodeling in asthma. This study explores how microRNA-15b-5p (miR-15b-5p) functions in modulating the proliferation, migration, inflammatory response, and ECM deposition of ASM cells. MiR-15b-5p and yes-associated protein 1 (YAP1) mRNA expression levels in tumor necrosis factor alpha (TNF-α)-induced ASM cells were, respectively, examined by real-time quantitative polymerase-chain reaction. Besides, the proliferative ability and migrative potential of ASM cells were examined by cell counting kit-8 assay, 5-bromo-2 ‘-deoxyuridine assay, and transwell assays, respectively. Interleukin-6 and interleukin-8 levels in ASM cells were detected by enzyme-linked immunosorbent assay. YAP1, collagen I, and collagen III expressions in ASM cells were detected by Western blot. With dual-luciferase reporter gene assay, the relations between miR-15b-5p and YAP1 3ʹUTR in ASM cells was examined. MiR-15b-5p expression level was reduced in ASM cells treated with TNF-α. MiR-15b-5p repressed TNF-α-initiated growth and migration of ASM cells and also suppressed IL-6 and IL-8 secretion, and inhibited collagen I and collagen III expressions in ASM cells. Furthermore, it was validated that YAP1 was a downstream target of miR-15b-5p in ASM cells. Notably, YAP1 overexpression attenuated the inhibitory effects of miR-15b-5p up-regulation on the proliferation, migration, and inflammatory response, as well as ECM deposition of TNF-α-induced ASM cells. In conclusion, miR-15b-5p/YAP1 axis modulates the growth, migration, inflammatory response, and ECM deposition of ASM cells, thus participating in the pathogenesis of asthma.

Multi-omics colocalization with genome-wide association studies reveals a context-specific genetic mechanism at a childhood onset asthma risk locus

BACKGROUND

Genome-wide association studies (GWASs) have identified thousands of variants associated with asthma and other complex diseases. However, the functional effects of most of these variants are unknown. Moreover, GWASs do not provide context-specific information on cell types or environmental factors that affect specific disease risks and outcomes. To address these limitations, we used an upper airway epithelial cell (AEC) culture model to assess transcriptional and epigenetic responses to rhinovirus (RV), an asthma-promoting pathogen, and provide context-specific functional annotations to variants discovered in GWASs of asthma.

METHODS

Genome-wide genetic, gene expression, and DNA methylation data in vehicle- and RV-treated upper AECs were collected from 104 individuals who had a diagnosis of airway disease (n=66) or were healthy participants (n=38). We mapped cis expression and methylation quantitative trait loci (cis-eQTLs and cis-meQTLs, respectively) in each treatment condition (RV and vehicle) in AECs from these individuals. A Bayesian test for colocalization between AEC molecular QTLs and adult onset asthma and childhood onset asthma GWAS SNPs, and a multi-ethnic GWAS of asthma, was used to assign the function to variants associated with asthma. We used Mendelian randomization to demonstrate DNA methylation effects on gene expression at asthma colocalized loci.

RESULTS

Asthma and allergic disease-associated GWAS SNPs were specifically enriched among molecular QTLs in AECs, but not in GWASs from non-immune diseases, and in AEC eQTLs, but not among eQTLs from other tissues. Colocalization analyses of AEC QTLs with asthma GWAS variants revealed potential molecular mechanisms of asthma, including QTLs at the TSLP locus that were common to both the RV and vehicle treatments and to both childhood onset and adult onset asthma, as well as QTLs at the 17q12-21 asthma locus that were specific to RV exposure and childhood onset asthma, consistent with clinical and epidemiological studies of these loci.

CONCLUSIONS

This study provides evidence of functional effects for asthma risk variants in AECs and insight into RV-mediated transcriptional and epigenetic response mechanisms that modulate genetic effects in the airway and risk for asthma.

Using "-omics" Data to Inform Genome-wide Association Studies (GWASs) in the Osteoporosis Field

PURPOSE OF REVIEW

Osteoporosis constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 loci influencing bone mineral density (BMD); however, few of the causal genes have been identified. Here, we review approaches that use "-omics" data and genetic- and systems genetics-based analytical strategies to facilitate causal gene discovery.

RECENT FINDINGS

The bone field is beginning to adopt approaches that are commonplace in other disease disciplines. The slower progress has been due in part to the lack of large-scale "omics" data on bone and bone cells. This is however changing, and approaches such as eQTL colocalization, transcriptome-wide association studies (TWASs), network, and integrative approaches are beginning to provide significant insight into the genes responsible for BMD GWAS associations. The use of "-omics" data to inform BMD GWASs has increased in recent years, leading to the identification of novel regulators of BMD in humans. The ultimate goal will be to use this information to develop more effective therapies to treat and ultimately prevent osteoporosis.

Psychosocial experiences modulate asthma-associated genes through gene-environment interactions

Social interactions and the overall psychosocial environment have a demonstrated impact on health, particularly for people living in disadvantaged urban areas. Here, we investigated the effect of psychosocial experiences on gene expression in peripheral blood immune cells of children with asthma in Metro Detroit. Using RNA-sequencing and a new machine learning approach, we identified transcriptional signatures of 19 variables including psychosocial factors, blood cell composition, and asthma symptoms. Importantly, we found 169 genes associated with asthma or allergic disease that are regulated by psychosocial factors and 344 significant gene-environment interactions for gene expression levels. These results demonstrate that immune gene expression mediates the link between negative psychosocial experiences and asthma risk.

Airway smooth muscle pathophysiology in asthma

The airway smooth muscle (ASM) cell plays a central role in the pathogenesis of asthma and constitutes an important target for treatment. These cells control muscle tone and thus regulate the opening of the airway lumen and air passage. Evidence indicates that ASM cells participate in the airway hyperresponsiveness as well as the inflammatory and remodeling processes observed in asthmatic subjects. Therapeutic approaches require a comprehensive understanding of the structure and function of the ASM in both the normal and disease states. This review updates current knowledge about ASM and its effects on airway narrowing, remodeling, and inflammation in asthma.

Cytokine-induced molecular responses in airway smooth muscle cells inform genome-wide association studies of asthma

BACKGROUND

A challenge in the post-GWAS era is to assign function to disease-associated variants. However, available resources do not include all tissues or environmental exposures that are relevant to all diseases. For example, exaggerated bronchoconstriction of airway smooth muscle cells (ASMCs) defines airway hyperresponsiveness (AHR), a cardinal feature of asthma. However, the contribution of ASMC to genetic and genomic studies has largely been overlooked. Our study aimed to address the gap in data availability from a critical tissue in genomic studies of asthma.

METHODS

We developed a cell model of AHR to discover variants associated with transcriptional, epigenetic, and cellular responses to two AHR promoting cytokines, IL-13 and IL-17A, and performed a GWAS of bronchial responsiveness (BRI) in humans.

RESULTS

Our study revealed significant response differences between ASMCs from asthma cases and controls, including genes implicated in asthma susceptibility. We defined molecular quantitative trait loci (QTLs) for expression (eQTLs) and methylation (meQTLs), and cellular QTLs for contractility (coQTLs) and performed a GWAS of BRI in human subjects. Variants in asthma GWAS were significantly enriched for ASM QTLs and BRI-associated SNPs, and near genes enriched for ASM function, many with small P values that did not reach stringent thresholds of significance in GWAS.

CONCLUSIONS

Our study identified significant differences between ASMCs from asthma cases and controls, potentially reflecting trained tolerance in these cells, as well as a set of variants, overlooked in previous GWAS, which reflect the AHR component of asthma.