Exploring proteomic plasma biomarkers in eosinophilic and neutrophilic asthma

BACKGROUND

Few biomarkers identify eosinophilic and neutrophilic asthma beyond cell concentrations in blood or sputum. Finding novel biomarkers for asthma endotypes could give insight about disease mechanisms and guide tailored treatment. Our aim was to investigate clinical characteristics and inflammation-related plasma proteins in relation to blood eosinophil and neutrophil concentrations in subjects with and without asthma.

METHODS

We included 24-26-year-old subjects (n = 2063) from the Swedish population-based cohort BAMSE. Subjects with asthma (n = 239) and without asthma (n = 1824) were subdivided based on blood eosinophil and neutrophil concentrations (cut-offs 0.3 × 10⁹ /L and 5.0 × 10⁹ /L, respectively). We measured the levels of 92 plasma proteins using Olink Proseek Multiplex Inflammation Panel Assay. Group statistics tests were used to analyse the data, as well as adjusted multiple logistic regression models.

RESULTS

Among subjects with asthma, 21.8% had eosinophilic asthma and 20.5% neutrophilic asthma. Eosinophilic asthma, but not neutrophilic asthma, was associated with a distinct clinical phenotype with, for example, higher proportions of eczema and sensitization. Most plasma proteins that associated with high eosinophil and/or neutrophil blood concentrations in subjects with asthma showed similar associations in subjects without asthma. However, out of these proteins, MMP10 levels were associated with eosinophilic asthma and were significantly higher as compared to controls with high eosinophilic concentration, while CCL4 levels associated with high neutrophil concentration only in subjects with asthma.

CONCLUSIONS

Eosinophilic asthma was associated with a clear clinical phenotype. With our definitions, we identified MMP10 as a possible plasma biomarker for eosinophilic asthma and CCL4 was linked to neutrophilic asthma. These proteins should be evaluated further in clinical settings and using sputum granulocytes to define the asthma endotypes.

Multi-ancestry genome-wide association study of asthma exacerbations

BACKGROUND

Asthma exacerbations are a serious public health concern due to high healthcare resource utilization, work/school productivity loss, impact on quality of life, and risk of mortality. The genetic basis of asthma exacerbations has been studied in several populations, but no prior study has performed a multi-ancestry meta-analysis of genome-wide association studies (meta-GWAS) for this trait. We aimed to identify common genetic loci associated with asthma exacerbations across diverse populations and to assess their functional role in regulating DNA methylation and gene expression.

METHODS

A meta-GWAS of asthma exacerbations in 4989 Europeans, 2181 Hispanics/Latinos, 1250 Singaporean Chinese, and 972 African Americans analyzed 9.6 million genetic variants. Suggestively associated variants (p ≤ 5 × 10^-5 ) were assessed for replication in 36,477 European and 1078 non-European asthma patients. Functional effects on DNA methylation were assessed in 595 Hispanic/Latino and African American asthma patients and in publicly available databases. The effect on gene expression was evaluated in silico.

RESULTS

One hundred and twenty-six independent variants were suggestively associated with asthma exacerbations in the discovery phase. Two variants independently replicated: rs12091010 located at vascular cell adhesion molecule-1/exostosin like glycosyltransferase-2 (VCAM1/EXTL2) (discovery: odds ratio (OR_{T allele} ) = 0.82, p = 9.05 × 10^-6 and replication: OR_{T allele}  = 0.89, p = 5.35 × 10^-3 ) and rs943126 from pantothenate kinase 1 (PANK1) (discovery: OR_{C allele}  = 0.85, p = 3.10 × 10^-5 and replication: OR_{C allele}  = 0.89, p = 1.30 × 10^-2 ). Both variants regulate gene expression of genes where they locate and DNA methylation levels of nearby genes in whole blood.

CONCLUSIONS

This multi-ancestry study revealed novel suggestive regulatory loci for asthma exacerbations located in genomic regions participating in inflammation and host defense.

Gene-environment interactions in childhood asthma revisited; expanding the interaction concept

Investigation of gene-environment interactions (GxE) may provide important insights into the gene regulatory framework in response to environmental factors of relevance for childhood asthma. Over the years, different methodological strategies have been applied, more recently using genome-wide approaches. The best example to date is the major asthma locus on the 17q12-21 chromosome region, viral infections, and airway epithelium processes where recent studies have shed much light on mechanisms in childhood asthma. However, there are challenges with the traditional single variant-single exposure interaction models, as they do not encompass the complexity and cumulative effects of multiple exposures or multiple genetic variants. As such, we need to redefine our traditional GxE thinking, and we propose in this review to expand the GxE concept by also evaluating other omics layers, such as epigenetics, transcriptomics, metabolomics, and proteomics. In addition, host factors such as age, gender, and other exposures are very likely to influence GxE effects and need firmly to be considered in future studies.

ADRB2 haplotypes and asthma exacerbations in children and young adults: An individual participant data meta-analysis

BACKGROUND

The polymorphism Arg16 in β2 -adrenergic receptor (ADRB2) gene has been associated with an increased risk of exacerbations in asthmatic children treated with long-acting β2 -agonists (LABA). However, it remains unclear whether this increased risk is mainly attributed to this single variant or the combined effect of the haplotypes of polymorphisms at codons 16 and 27.

OBJECTIVE

We assessed whether the haplotype analysis could explain the association between the polymorphisms at codons 16 (Arg16Gly) and 27 (Gln27Glu) in ADRB2 and risk of asthma exacerbations in patients treated with inhaled corticosteroids (ICS) plus LABA.

METHODS

The study was undertaken using data from 10 independent studies (n = 5903) participating in the multi-ethnic Pharmacogenomics in Childhood Asthma (PiCA) consortium. Asthma exacerbations were defined as asthma-related use of oral corticosteroids or hospitalizations/emergency department visits in the past 6 or 12 months prior to the study visit/enrolment. The association between the haplotypes and the risk of asthma exacerbations was performed per study using haplo.stats package adjusted for age and sex. Results were meta-analysed using the inverse variance weighting method assuming random-effects.

RESULTS

In subjects treated with ICS and LABA (n = 832, age: 3-21 years), Arg16/Gln27 versus Gly16/Glu27 (OR: 1.40, 95% CI: 1.05-1.87, I2  = 0.0%) and Arg16/Gln27 versus Gly16/Gln27 (OR: 1.43, 95% CI: 1.05-1.94, I2  = 0.0%), but not Gly16/Gln27 versus Gly16/Glu27 (OR: 0.99, 95% CI: 0.71-1.39, I2  = 0.0%), were significantly associated with an increased risk of asthma exacerbations. The sensitivity analyses indicated no significant association between the ADRB2 haplotypes and asthma exacerbations in the other treatment categories, namely as-required short-acting β2 -agonists (n = 973), ICS monotherapy (n = 2623), ICS plus leukotriene receptor antagonists (LTRA; n = 338), or ICS plus LABA plus LTRA (n = 686).

CONCLUSION AND CLINICAL RELEVANCE

The ADRB2 Arg16 haplotype, presumably mainly driven by the Arg16, increased the risk of asthma exacerbations in patients treated with ICS plus LABA. This finding could be beneficial in ADRB2 genotype-guided treatment which might improve clinical outcomes in asthmatic patients.

Genome-wide association studies of exacerbations in children using long-acting beta2-agonists

BACKGROUND

Some children with asthma experience exacerbations despite long-acting beta2-agonist (LABA) treatment. While this variability is partly caused by genetic variation, no genome-wide study until now has investigated which genetic factors associated with risk of exacerbations despite LABA use in children with asthma. We aimed to assess whether genetic variation was associated with exacerbations in children treated with LABA from a global consortium.

METHODS

A meta-analysis of genome-wide association studies (meta-GWAS) was performed in 1,425 children and young adults with asthma (age 6-21 years) with reported regular use of LABA from six studies within the PiCA consortium using a random effects model. The primary outcome of each study was defined as any exacerbation within the past 6 or 12 months, including at least one of the following: 1) hospital admissions for asthma, 2) a course of oral corticosteroids or 3) emergency room visits because of asthma.

RESULTS

Genome-wide association results for a total of 82 996 common single nucleotide polymorphisms (SNPs, MAF ≥1%) with high imputation quality were meta-analysed. Eight independent variants were suggestively (P-value threshold ≤5 × 10^-6 ) associated with exacerbations despite LABA use.

CONCLUSION

No strong effects of single nucleotide polymorphisms (SNPs) on exacerbations during LABA use were identified. We identified two loci (TBX3 and EPHA7) that were previously implicated in the response to short-acting beta2-agonists (SABA). These loci merit further investigation in response to LABA and SABA use.

Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma

Inhaled corticosteroids (ICS) are the most common asthma controller medication. An important contribution of genetic factors in ICS response has been evidenced. Here, we aimed to identify novel genetic markers involved in ICS response in asthma. A genome-wide association study (GWAS) of the change in lung function after 6 weeks of ICS treatment was performed in 166 asthma patients from the SLOVENIA study. Patients with an improvement in lung function ≥8% were considered as ICS responders. Suggestively associated variants (p-value ≤ 5 × 10⁻⁶) were evaluated in an independent study (n = 175). Validation of the association with asthma exacerbations despite ICS use was attempted in European (n = 2681) and admixed (n = 1347) populations. Variants previously associated with ICS response were also assessed for replication. As a result, the SNP rs1166980 from the ROBO2 gene was suggestively associated with the change in lung function (OR for G allele: 7.01, 95% CI: 3.29–14.93, p = 4.61 × 10⁻⁷), although this was not validated in CAMP. ROBO2 showed gene-level evidence of replication with asthma exacerbations despite ICS use in Europeans (minimum p-value = 1.44 × 10⁻⁵), but not in admixed individuals. The association of PDE10A-T with ICS response described by a previous study was validated. This study suggests that ROBO2 could be a potential novel locus for ICS response in Europeans.

Genome-wide association study of asthma exacerbations despite inhaled corticosteroid use

RATIONALE

Substantial variability in response to asthma treatment with inhaled corticosteroids (ICS) has been described among individuals and populations, suggesting the contribution of genetic factors. Nonetheless, only a few genes have been identified to date. We aimed to identify genetic variants associated with asthma exacerbations despite ICS use in European children and young adults and to validate the findings in non-Europeans. Moreover, we explored whether a gene-set enrichment analysis could suggest potential novel asthma therapies.

METHODS

A genome-wide association study (GWAS) of asthma exacerbations was tested in 2681 children of European descent treated with ICS from eight studies. Suggestive association signals were followed up for replication in 538 European asthma patients. Further evaluation was performed in 1773 non-Europeans. Variants revealed by published GWAS were assessed for replication. Additionally, gene-set enrichment analysis focused on drugs was performed.

RESULTS

10 independent variants were associated with asthma exacerbations despite ICS treatment in the discovery phase (p≤5×10-6). Of those, one variant at the CACNA2D3-WNT5A locus was nominally replicated in Europeans (rs67026078; p=0.010), but this was not validated in non-European populations. Five other genes associated with ICS response in previous studies were replicated. Additionally, an enrichment of associations in genes regulated by trichostatin A treatment was found.

CONCLUSIONS

The intergenic region of CACNA2D3 and WNT5A was revealed as a novel locus for asthma exacerbations despite ICS treatment in European populations. Genes associated were related to trichostatin A, suggesting that this drug could regulate the molecular mechanisms involved in treatment response.

Pharmacogenomic associations of adverse drug reactions in asthma: systematic review and research prioritisation

A systematic review of pharmacogenomic studies capturing adverse drug reactions (ADRs) related to asthma medications was undertaken, and a survey of Pharmacogenomics in Childhood Asthma (PiCA) consortia members was conducted. Studies were eligible if genetic polymorphisms were compared with suspected ADR(s) in a patient with asthma, as either a primary or secondary outcome. Five studies met the inclusion criteria. The ADRs and polymorphisms identified were change in lung function tests (rs1042713), adrenal suppression (rs591118), and decreased bone mineral density (rs6461639) and accretion (rs9896933, rs2074439). Two of these polymorphisms were replicated within the paper, but none had external replication. Priorities from PiCA consortia members (representing 15 institution in eight countries) for future studies were tachycardia (SABA/LABA), adrenal suppression/crisis and growth suppression (corticosteroids), sleep/behaviour disturbances (leukotriene receptor antagonists), and nausea and vomiting (theophylline). Future pharmacogenomic studies in asthma should collect relevant ADR data as well as markers of efficacy.

Bacterial salivary microbiome associates with asthma among african american children and young adults

Several studies have shown that the airways of asthma patients contain higher diversity of bacteria and are enriched in pathogenic species. However, sampling the airways in children is challenging. Here we aimed to identify differences in the salivary bacterial composition between African Americans children with and without asthma. Saliva samples from 57 asthma cases and 57 healthy controls were analyzed by means of 16S ribosomal RNA amplicon profiling. Measurements of bacterial diversity and genus relative abundance were compared between cases and controls using the nonparametric Wilcoxon test and multivariate regression models. A total of five phyla and a mean of 56 genera were identified. Among them, 15 genera had a relative abundance greater than 1%, being Prevotella, Haemophilus, Streptococcus, and Veillonella the most abundant genera. Differences between cases and controls were found in terms of diversity, as well as in relative abundance for Streptococcus genus (13.0% in cases vs 18.3% in controls; P = .003) and Veillonella genus (11.1% in cases vs 8.0% in controls; P = .002). These differences remained significant after correction for multiple comparisons and when potential confounders were taken into account in logistic regression models. In conclusion, we identified changes in the salivary microbiota associated with asthma among African Americans.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations

BACKGROUND

Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations.

RESULTS

We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response.

CONCLUSIONS

Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

Genome-wide association study of inhaled corticosteroid response in admixed children with asthma

BACKGROUND

Inhaled corticosteroids (ICS) are the most widely prescribed and effective medication to control asthma symptoms and exacerbations. However, many children still have asthma exacerbations despite treatment, particularly in admixed populations, such as Puerto Ricans and African Americans. A few genome-wide association studies (GWAS) have been performed in European and Asian populations, and they have demonstrated the importance of the genetic component in ICS response.

OBJECTIVE

We aimed to identify genetic variants associated with asthma exacerbations in admixed children treated with ICS and to validate previous GWAS findings.

METHODS

A meta-analysis of two GWAS of asthma exacerbations was performed in 1347 admixed children treated with ICS (Hispanics/Latinos and African Americans), analysing 8.7 million genetic variants. Those with P ≤ 5 × 10^-6 were followed up for replication in 1697 asthmatic patients from six European studies. Associations of ICS response described in published GWAS were followed up for replication in the admixed populations.

RESULTS

A total of 15 independent variants were suggestively associated with asthma exacerbations in admixed populations (P ≤ 5 × 10^-6 ). One of them, located in the intergenic region of APOBEC3B and APOBEC3C, showed evidence of replication in Europeans (rs5995653, P = 7.52 × 10^-3 ) and was also associated with change in lung function after treatment with ICS (P = 4.91 × 10^-3 ). Additionally, the reported association of the L3MBTL4-ARHGAP28 genomic region was confirmed in admixed populations, although a different variant was identified.

CONCLUSIONS AND CLINICAL RELEVANCE

This study revealed the novel association of APOBEC3B and APOBEC3C with asthma exacerbations in children treated with ICS and replicated previously identified genomic regions. This contributes to the current knowledge about the multiple genetic markers determining responsiveness to ICS which could lead in the future the clinical identification of those asthma patients who are not able to respond to such treatment.

Genomic Predictors of Asthma Phenotypes and Treatment Response

Asthma is a complex respiratory disease considered as the most common chronic condition in children. A large genetic contribution to asthma susceptibility is predicted by the clustering of asthma and allergy symptoms among relatives and the large disease heritability estimated from twin studies, ranging from 55 to 90%. Genetic basis of asthma has been extensively investigated in the past 40 years using linkage analysis and candidate-gene association studies. However, the development of dense arrays for polymorphism genotyping has enabled the transition toward genome-wide association studies (GWAS), which have led the discovery of several unanticipated asthma genes in the last 11 years. Despite this, currently known risk variants identified using many thousand samples from distinct ethnicities only explain a small proportion of asthma heritability. This review examines the main findings of the last 2 years in genomic studies of asthma using GWAS and admixture mapping studies, as well as the direction of studies fostering integrative perspectives involving omics data. Additionally, we discuss the need for assessing the whole spectrum of genetic variation in association studies of asthma susceptibility, severity, and treatment response in order to further improve our knowledge of asthma genes and predictive biomarkers. Leveraging the individual's genetic information will allow a better understanding of asthma pathogenesis and will facilitate the transition toward a more precise diagnosis and treatment.

Author Correction: Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

In the version of this article initially published, in Fig. 3, the y-axis numbering did not match the log scale indicated in the axis label. The error has been corrected in the HTML and PDF version of the article.

Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

Allergic rhinitis is the most common clinical presentation of allergy, affecting 400 million people worldwide, with increasing incidence in westernized countries^1,2. To elucidate the genetic architecture and understand the underlying disease mechanisms, we carried out a meta-analysis of allergic rhinitis in 59,762 cases and 152,358 controls of European ancestry and identified a total of 41 risk loci for allergic rhinitis, including 20 loci not previously associated with allergic rhinitis, which were confirmed in a replication phase of 60,720 cases and 618,527 controls. Functional annotation implicated genes involved in various immune pathways, and fine mapping of the HLA region suggested amino acid variants important for antigen binding. We further performed genome-wide association study (GWAS) analyses of allergic sensitization against inhalant allergens and nonallergic rhinitis, which suggested shared genetic mechanisms across rhinitis-related traits. Future studies of the identified loci and genes might identify novel targets for treatment and prevention of allergic rhinitis.

A vascular endothelial growth factor receptor gene variant is associated with susceptibility to acute respiratory distress syndrome

BACKGROUND

The acute respiratory distress syndrome (ARDS) is one of the main causes of mortality in adults admitted to intensive care units. Previous studies have demonstrated the existence of genetic variants involved in the susceptibility and outcomes of this syndrome. We aimed to identify novel genes implicated in sepsis-induced ARDS susceptibility.

METHODS

We first performed a prioritization of candidate genes by integrating our own genomic data from a transcriptomic study in an animal model of ARDS and from the only published genome-wide association study of ARDS study in humans. Then, we selected single nucleotide polymorphisms (SNPs) from prioritized genes to conduct a case-control discovery association study in patients with sepsis-induced ARDS (n = 225) and population-based controls (n = 899). Finally, we validated our findings in an independent sample of 661 sepsis-induced ARDS cases and 234 at-risk controls.

RESULTS

Three candidate genes were prioritized: dynein cytoplasmic-2 heavy chain-1, fms-related tyrosine kinase 1 (FLT1), and integrin alpha-1. Of those, a SNP from FLT1 gene (rs9513106) was associated with ARDS in the discovery study, with an odds ratio (OR) for the C allele of 0.76, 95% confidence interval (CI) 0.58-0.98 (p = 0.037). This result was replicated in an independent study (OR = 0.78, 95% CI = 0.62-0.98, p = 0.039), showing consistent direction of effects in a meta-analysis (OR = 0.77, 95% CI = 0.65-0.92, p = 0.003).

CONCLUSIONS

We identified FLT1 as a novel ARDS susceptibility gene and demonstrated that integration of genomic data can be a valid procedure to identify novel susceptibility genes. These results contribute to previous firm associations and functional evidences implicating FLT1 gene in other complex traits that are mechanistically linked, through the key role of endothelium, to the pathophysiology of ARDS.

Rationale and design of the multiethnic Pharmacogenomics in Childhood Asthma consortium

AIM

International collaboration is needed to enable large-scale pharmacogenomics studies in childhood asthma. Here, we describe the design of the Pharmacogenomics in Childhood Asthma (PiCA) consortium.

MATERIALS & METHODS

Investigators of each study participating in PiCA provided data on the study characteristics by answering an online questionnaire.

RESULTS

A total of 21 studies, including 14,227 children/young persons (58% male), from 12 different countries are currently enrolled in the PiCA consortium. Fifty six percent of the patients are Caucasians. In total, 7619 were inhaled corticosteroid users. Among patients from 13 studies with available data on asthma exacerbations, a third reported exacerbations despite inhaled corticosteroid use. In the future pharmacogenomics studies within the consortium, the pharmacogenomics analyses will be performed separately in each center and the results will be meta-analyzed.

CONCLUSION

PiCA is a valuable platform to perform pharmacogenetics studies within a multiethnic pediatric asthma population.