The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth.

BACKGROUND: The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences. METHODS: DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants. RESULTS: 128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study. CONCLUSIONS: We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture.

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Admixture mapping of severe asthma exacerbations in Hispanic/Latino children and youth.

BACKGROUND: In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations. OBJECTIVE: We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles. METHODS: We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases. RESULTS: Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels. CONCLUSIONS: Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2.

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 (ORT allele ) = 0.82, p = 9.05 × 10-6 and replication: ORT allele  = 0.89, p = 5.35 × 10-3 ) and rs943126 from pantothenate kinase 1 (PANK1) (discovery: ORC allele  = 0.85, p = 3.10 × 10-5 and replication: ORC 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.

Race- and Ethnicity-Based Spirometry Reference Equations: Are They Accurate for Genetically Admixed Children?

BACKGROUND: Variation in genetic ancestry among genetically admixed racial and ethnic groups may influence the fit of guideline-recommended spirometry reference equations, which rely on self-identified race and ethnicity. RESEARCH QUESTION: What is the influence of genetic ancestry on the fit of race- and ethnicity-based spirometry reference equations in populations of genetically admixed children? STUDY DESIGN AND METHODS: Cross-sectional fit of guideline-recommended race- and ethnicity-based spirometry reference equations was evaluated in healthy control participants from case-control studies of asthma. Anthropometry, blood samples, and spirometric measurements were obtained for 599 genetically admixed children 8 to 21 years of age. Genetic ancestry was estimated using genome-wide genotype data. Equation fit, measured as a mean z score, was assessed in self-identified African American (n = 275) and Puerto Rican (n = 324) children as well as genetic ancestry-defined strata of each population. RESULTS: For African American children, African American-derived equations fit for predicting FEV1 and FVC in those with an African ancestry more than the median (81.4%-100.0%), whereas composite equations for "other/mixed" populations fit for predicting FEV1 and FVC in those with African ancestry at or less than the median (30.7%-81.3%). For Puerto Rican children with African ancestry at or less than the median (6.4%-21.3%), White-derived equations fit both FEV1 and FVC, whereas for those with African ancestry more than the median (21.4%-87.5%), White-derived equations fit the FEV1 and the composite equations fit the FVC. INTERPRETATION: Guideline-recommended spirometry reference equations yielded biased estimates of lung function in genetically admixed children with high variation of African ancestry. Spirometry could benefit from reference equations that incorporate genetic ancestry, either for more precise application of the current equations or the derivation and use of new equations.

Role of Sex on the Genetic Susceptibility to Childhood Asthma in Latinos and African Americans.

Asthma is a respiratory disease whose prevalence changes throughout the lifespan and differs by sex, being more prevalent in males during childhood and in females after puberty. In this study, we assessed the influence of sex on the genetic susceptibility to childhood asthma in admixed populations. Sex-interaction and sex-stratified genome-wide association studies (GWAS) were performed in 4291 Latinos and 1730 African Americans separately, and results were meta-analyzed. Genome-wide (p ≤ 9.35 × 10-8) and suggestive (p ≤ 1.87 × 10-6) population-specific significance thresholds were calculated based on 1000 Genomes Project data. Additionally, protein quantitative trait locus (pQTL) information was gathered for the suggestively associated variants, and enrichment analyses of the proteins identified were carried out. Four independent loci showed interaction with sex at a suggestive level. The stratified GWAS highlighted the 17q12-21 asthma locus as a contributor to asthma susceptibility in both sexes but reached genome-wide significance only in females (p-females < 9.2 × 10-8; p-males < 1.25 × 10-2). Conversely, genetic variants upstream of ligand-dependent nuclear receptor corepressor-like gene (LCORL), previously involved in height determination and spermatogenesis, were associated with asthma only in males (minimum p = 5.31 × 10-8 for rs4593128). Enrichment analyses revealed an overrepresentation of processes related to the immune system and highlighted differences between sexes. In conclusion, we identified sex-specific polymorphisms that could contribute to the differences in the prevalence of childhood asthma between males and females.

Lymph node-resident dendritic cells drive TH2 cell development involving MARCH1.

Type 2 T helper (TH2) cells are protective against parasitic worm infections but also aggravate allergic inflammation. Although the role of dendritic cells (DCs) in TH2 cell differentiation is well established, the underlying mechanisms are largely unknown. Here, we show that DC induction of TH2 cells depends on membrane-associated RING-CH-1 (MARCH1) ubiquitin ligase. The pro-TH2 effect of MARCH1 relied on lymph node (LN)­resident DCs, which triggered T cell receptor (TCR) signaling and induced GATA-3 expression from naïve CD4+ T cells independent of tissue-driven migratory DCs. Mice with mutations in the ubiquitin acceptor sites of MHCII and CD86, the two substrates of MARCH1, failed to develop TH2 cells. These findings suggest that TH2 cell development depends on ubiquitin-mediated clearance of antigen-presenting and costimulatory molecules by LN-resident DCs and consequent control of TCR signaling.

Racial/ethnic differences in eligibility for asthma biologics among pediatric populations.

BACKGROUND: Asthma is a heterogeneous disease. Clinical blood parameters differ by race/ethnicity and are used to distinguish asthma subtypes and inform therapies. Differences in subtypes may explain population-specific trends in asthma outcomes. However, these differences in racial/ethnic minority pediatric populations are unclear. OBJECTIVE: We investigated the association of blood parameters and asthma subtypes with asthma outcomes and examined population-specific eligibility for biologic therapies in minority pediatric populations. METHODS: Using data from 2 asthma case-control studies of pediatric minority populations, we performed case-control (N = 3738) and case-only (N = 2743) logistic regressions to quantify the association of blood parameters and asthma subtypes with asthma outcomes. Heterogeneity of these associations was tested using an interaction term between race/ethnicity and each exposure. Differences in therapeutic eligibility were investigated using chi-square tests. RESULTS: Race/ethnicity modified the association between total IgE and asthma exacerbations. Elevated IgE level was associated with worse asthma outcomes in Puerto Ricans. Allergic asthma was associated with worse outcomes in Mexican Americans, whereas eosinophilic asthma was associated with worse outcomes in Puerto Ricans. A lower proportion of Puerto Ricans met dosing criteria for allergic asthma-directed biologic therapy than other groups. A higher proportion of Puerto Ricans qualified for eosinophilic asthma-directed biologic therapy than African Americans. CONCLUSIONS: We found population-specific associations between blood parameters and asthma subtypes with asthma outcomes. Our findings suggest that eligibility for asthma biologic therapies differs across pediatric racial/ethnic populations. These findings call for more studies in diverse populations for equitable treatment of minority patients with asthma.

A genome-wide study of DNA methylation in white blood cells and asthma in Latino children and youth.

Latinos are heavily affected with childhood asthma. Little is known about epigenetic mechanisms of asthma in Latino youth. We conducted a meta-analysis of two epigenome-wide association studies (EWAS) of asthma, using DNA from white blood cells (WBCs) from 1,136 Latino children and youth aged 6 to 20 years. Genes near the top CpG sites in this EWAS were examined in a pathway enrichment analysis, and we then assessed whether our results replicated those from publicly available data from three independent EWAS conducted in non-Latino populations. We found that DNA methylation profiles differed between subjects with and without asthma. After adjustment for covariates and multiple testing, two CpGs were differentially methylated at a false discovery rate (FDR)-adjusted P < 0.1, and 193 CpG sites were differentially methylated at FDR-adjusted P < 0.2. The two top CpGs are near genes relevant to inflammatory signalling, including CAMK1D (Calcium/Calmodulin Dependent Protein Kinase ID) and TIGIT (T Cell Immunoreceptor With Ig And ITIM Domains). Moreover, 25 genomic regions were differentially methylated between subjects with and without asthma, at Sidák-corrected P < 0.10. An enrichment analysis then identified the TGF-beta pathway as most relevant to asthma in our analysis, and we replicated some of the top signals from publicly available EWAS datasets in non-Hispanic populations. In conclusion, we have identified novel epigenetic markers of asthma in WBCs from Latino children and youth, while also replicating previous results from studies conducted in non-Latinos.

A genome-wide association study of asthma hospitalizations in adults.

BACKGROUND: Little is known about the genetic determinants of severe asthma exacerbations. OBJECTIVES: We aimed to identify genetic variants associated with asthma hospitalizations. METHODS: We conducted a genome-wide association study of asthma hospitalizations in 34,167 white British adults with asthma, 1,658 of whom had at least 1 asthma-related hospitalization. This analysis was conducted by using logistic regression under an additive genetic model with adjustment for age, sex, body mass index, smoking status, and the first 5 principal components derived from genotypic data. We then analyzed data from 2 cohorts of Latino children and adolescents for replication and conducted quantitative trait locus and functional annotation analyses. RESULTS: At the chromosome 6p21.3 locus, the single-nucleotide polymorphism (SNP) rs56151658 (8 kb from the promoter of HLA-DQB1) was most significantly associated with asthma hospitalizations (for test allele A, odds ratio = 1.36 [95% CI = 1.22-1.52]; P = 3.11 × 10-8); 21 additional SNPs in this locus were associated with asthma hospitalizations at a P value less than 1 × 10-6. In the replication cohorts, multiple SNPs in strong linkage disequilibrium with rs56151658 were associated with severe asthma exacerbations at a P value of .01 or less in the same direction of association as in the discovery cohort. Three HLA genes (HLA-DQA2, HLA-DRB6, and HLA-DOB) were also shown to mediate the estimated effects of the SNPs associated with asthma hospitalizations through effects on gene expression in lung tissue. CONCLUSIONS: We identified strong candidate genes for asthma hospitalizations in adults in the region for class II HLA genes through genomic, quantitative trait locus, and summary data-based mendelian randomization analyses.

Genome-wide association study reveals a novel locus for asthma with severe exacerbations in diverse populations.

BACKGROUND: Severe asthma exacerbations are a major cause of asthma morbidity and increased healthcare costs. Several studies have shown racial and ethnic differences in asthma exacerbation rates. We aimed to identify genetic variants associated with severe exacerbations in two high-risk populations for asthma. METHODS: A genome-wide association study of asthma in children and youth with severe exacerbations was performed in 1283 exacerbators and 2027 controls without asthma of Latino ancestry. Independent suggestive variants (P ≤ 5 × 10-6 ) were selected for replication in 448 African Americans exacerbators and 595 controls. Case-only analyses were performed comparing the exacerbators with additional 898 Latinos and 524 African Americans asthma patients without exacerbations, while adjusting by treatment category as a proxy of asthma severity. We analyzed the functionality of associated variants with in silico methods and by correlating genotypes with methylation levels in whole blood in a subset of 473 Latinos. RESULTS: We identified two genome-wide significant associations for susceptibility to asthma with severe exacerbations, including a novel locus located at chromosome 2p21 (rs4952375, odds ratio = 1.39, P = 3.8 × 10-8 ), which was also associated with asthma exacerbations in a case-only analysis (odds ratio = 1.25, P = 1.95 × 10-3 ). This polymorphism is an expression quantitative trait locus of the long intergenic non-protein coding RNA 1913 (LINC01913) in lung tissues (P = 1.3 × 10-7 ) and influences methylation levels of the protein kinase domain-containing cytoplasmic (PKDCC) gene in whole-blood cells (P = 9.8 × 10-5 ). CONCLUSION: We identified a novel susceptibility locus for severe asthma exacerbations in Hispanic/Latino and African American youths with functional effects in gene expression and methylation status of neighboring genes.

Differential asthma odds following respiratory infection in children from three minority populations.

RATIONALE: Severe early-life respiratory illnesses, particularly those caused by respiratory syncytial virus (RSV) and human rhinovirus (HRV), are strongly associated with the development of asthma in children. Puerto Rican children in particular have a strikingly high asthma burden. However, prior studies of the potential associations between early-life respiratory illnesses and asthma in Puerto Rican and other minority populations have been limited. OBJECTIVES: We sought to determine whether early-life respiratory illness was associated with asthma in Puerto Rican, Mexican American, and African American children. METHODS: Using a logistic regression analysis, we examined the association between early-life respiratory illnesses (report of upper respiratory infection (URI), pneumonia, bronchitis, and bronchiolitis/RSV) within the first two years of life and physician-diagnosed asthma after the age of two in a large cohort of Puerto Rican, Mexican American, and African American children. MEASUREMENTS AND MAIN RESULTS: While early-life respiratory illnesses were associated with greater asthma odds in Puerto Ricans, Mexican Americans, and African Americans, these associations were stronger among Puerto Rican children. Specifically, in Puerto Ricans, the odds was 6.15 (95% CI: 4.21-9.05) if the child reported at least one of the following respiratory illness: URI, pneumonia, bronchitis or bronchiolitis. The odds were also higher in Puerto Ricans when considering these conditions separately. CONCLUSIONS: We observed population-specific associations between early-life respiratory illnesses and asthma, which were especially significant and stronger in Puerto Ricans. Taken together with the known high burden of RSV in Puerto Rico, our results may help explain the high burden of asthma in Puerto Ricans.

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.

In utero tobacco smoke exposure, DNA methylation, and asthma in Latino children.

BACKGROUND: Maternal smoking during pregnancy is a risk factor for chronic disease later in life and has been associated with variability of DNA methylation at specific cytosine-phosphate-guanine (CpG) loci. We assessed the role of DNA methylation as a potential mediator of adverse effects of in utero tobacco smoke exposures on asthma outcomes in Latino children from the US mainland and Puerto Rico. METHODS: Relationships between self-reported exposure and DNA methylation at CpG loci previously reported to be associated with maternal smoking were assessed in a subsample consisting of 572 children aged 8-21 years (310 cases with asthma, 262 healthy controls), sampled from a larger asthma case-control study. Subsequently, we assessed associations between top loci and asthma-related outcomes, followed by mediation analysis for loci for which associations with outcomes were observed. RESULTS: Self-reported maternal smoking was associated with a -1.5% (95% confidence interval (CI) = -2.4%, -0.6%) lower methylation at CpG locus cg05575921 on the AHRR gene; a 1% increase in DNA methylation at the same locus resulted in an odds ratio (OR) of 0.90 (95% CI = 0.83, 0.96) for the odds of asthma. The OR for the indirect effect of maternal smoking on asthma mediated through methylation at the cg05575921 locus was 1.18 (95% CI = 1.07, 1.68), compared to the OR for the total effect of exposure in the parent study of 1.48 (95% CI = 1.03, 2.11). CONCLUSIONS: Our findings suggest potential mediation by DNA methylation in the association between maternal smoking during pregnancy and asthma status.

Racial/Ethnic-Specific Differences in the Effects of Inhaled Corticosteroid Use on Bronchodilator Response in Patients With Asthma.

American Thoracic Society guidelines recommend inhaled corticosteroid (ICS) therapy, plus a short-acting bronchodilator, in patients with persistent asthma. However, few prior studies have examined the efficacy of this combination in children of all racial/ethnic groups. We evaluated the association between ICS use and bronchodilator response (BDR) in three pediatric populations with persistent asthma (656 African American, 916 Puerto Rican, and 398 Mexican American children). The association was assessed using multivariable quantile regression. After adjusting for baseline forced expiratory volume in one second and use of controller medications, ICS use was significantly associated with increased BDR only among Mexican Americans (1.56%, P = 0.028) but not African Americans (0.49%, P = 0.426) or Puerto Ricans (0.16%, P = 0.813). Our results demonstrate that ICS augmentation is disproportionate across racial/ethnic groups, where improved BDR is observed in Mexican Americans only. This study highlights the complexities of treating asthma in children, and reinforces the importance of investigating the influence of race/ethnicity on pharmacological response.

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.