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

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.