Genome-wide association study of asthma, total IgE, and lung function in a cohort of Peruvian children

Adjusting for principal components can induce collider bias in genome-wide association studies

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Asthma-Genomic Advances Toward Risk Prediction

Asthma is a common complex airway disease whose prediction of disease risk and most severe outcomes is crucial in clinical practice for adequate clinical management. This review discusses the latest findings in asthma genomics and current obstacles faced in moving forward to translational medicine. While genome-wide association studies have provided valuable insights into the genetic basis of asthma, there are challenges that must be addressed to improve disease prediction, such as the need for diverse representation, the functional characterization of genetic variants identified, variant selection for genetic testing, and refining prediction models using polygenic risk scores.

Genome-wide association study based on clustering by obesity-related variables uncovers a genetic architecture of obesity in the Japanese and the UK populations

Whether all obesity-related variants contribute to the onset of obesity or one or a few variants cause obesity in genetically heterogeneous populations remains obscure. Here, we investigated the genetic architecture of obesity by clustering the Japanese and British populations with obesity using obesity-related factors. In Step-1, we conducted a genome-wide association study (GWAS) with body mass index (BMI) as the outcome for eligible participants. In Step-2, we assigned participants with obesity (BMI ≥25 kg/m2) to five clusters based on obesity-related factors. Subsequently, participants from each cluster and those with a BMI <25 kg/m2 were combined. A GWAS was conducted for each cluster. Several previously identified obesity-related genes were verified in Step-1. Of the genes detected in Step-1, unique obesity-related genes were detected separately for each cluster in Step-2. Our novel findings suggest that a smaller sample size with increased homogeneity may provide insights into the genetic architecture of obesity.

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Nasal epithelial gene expression and total IgE in children and adolescents with asthma

BACKGROUND

Little is known about nasal epithelial gene expression and total IgE in youth.

OBJECTIVE

We aimed to identify genes whose nasal epithelial expression differs by total IgE in youth, and group them into modules that could be mapped to airway epithelial cell types.

METHODS

We conducted a transcriptome-wide association study of total IgE in 469 Puerto Ricans aged 9 to 20 years who participated in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study, separately in all subjects and in those with asthma. We then attempted to replicate top findings for each analysis using data from 3 cohorts. Genes with a Benjamini-Hochberg-adjusted P value of less than .05 in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study and a P value of less than .05 in the same direction of association in 1 or more replication cohort were considered differentially expressed genes (DEGs). DEGs for total IgE in subjects with asthma were further dissected into gene modules using coexpression analysis, and such modules were mapped to specific cell types in airway epithelia using public single-cell RNA-sequencing data.

RESULTS

A higher number of DEGs for total IgE were identified in subjects with asthma (n = 1179 DEGs) than in all subjects (n = 631 DEGs). In subjects with asthma, DEGs were mapped to 11 gene modules. The top module for positive correlation with total IgE was mapped to myoepithelial and mucus secretory cells in lower airway epithelia and was regulated by IL-4, IL5, IL-13, and IL-33. Within this module, hub genes included CDH26, FETUB, NTRK2, CCBL1, CST1, and CST2. Furthermore, an enrichment analysis showed overrepresentation of genes in signaling pathways for synaptogenesis, IL-13, and ferroptosis, supporting interactions between interleukin- and acetylcholine-induced responses.

CONCLUSIONS

Our findings for nasal epithelial gene expression support neuroimmune coregulation of total IgE in youth with asthma.

From gene identifications to therapeutic targets for asthma: Focus on great potentials of TSLP, ORMDL3, and GSDMB

Asthma is a chronic respiratory disease, and clinically, asthma exacerbations remain difficult to treat. The disease is caused by combinations of and interactions between genetic and environmental factors. Genomic and genetic approaches identified many novel genes to treat asthma and brought new insights into the disease. The products of the genes have functional roles in regulating physiological or pathophysiological processes in airway structural cells and immune system cells. Genetic factors also interact with environmental factors such as air pollutants, and bacterial and viral infections to trigger the disease. Thymic stromal lymphopoietin (TSLP), orosomucoid-like 3 (ORMDL3), and gasdermin B (GSDMB) are three genes identified by genetic studies to have a great potential as therapeutic targets of asthma. TSLP is an important driver of type 2 inflammation. ORMDL3 mediates cell stress, sphingolipid synthesis, and viral and bacterial infections. GSDMB regulates cell pyroptosis through its N and C terminals and can bind sulfatides to influence inflammatory response. Investigating inhibitors or modulators for these pathways would bring a new landscape for therapeutics of asthma in future.

No impact of helminth coinfection in patients with smear positive tuberculosis on immunoglobulin levels using a novel method measuring Mycobacterium tuberculosis-specific antibodies

Helminth/tuberculosis (TB)-coinfection can reduce cell-mediated immunity against Mycobacterium tuberculosis (Mtb) and increase disease severity, although the effects are highly helminth species dependent. Mtb have long been ranked as the number one single infectious agent claiming the most lives. The only licensed vaccine for TB (BCG) offers highly variable protection against TB, and almost no protection against transmission of Mtb. In recent few years the identification of naturally occurring antibodies in humans that are protective during Mtb infection has reignited the interest in adaptive humoral immunity against TB and its possible implementation in novel TB vaccine design. The effects of helminth/TB coinfection on the humoral response against Mtb during active pulmonary TB are however still unclear, and specifically the effect by globally prevalent helminth species such as Ascaris lumbricoides, Strongyloides stercoralis, Ancylostoma duodenale, Trichuris trichiura. Plasma samples from smear positive TB patients were used to measure both total and Mtb-specific antibody responses in a Peruvian endemic setting where these helminths are dominating. Mtb-specific antibodies were detected by a novel approach coating ELISA-plates with a Mtb cell-membrane fraction (CDC1551) that contains a broad range of Mtb surface proteins. Compared to controls without helminths or TB, helminth/TB coinfected patients had high levels of Mtb-specific IgG (including an IgG1 and IgG2 subclass response) and IgM, which were similarly increased in TB patients without helminth infection. These data, indicate that helminth/TB coinfected have a sustained humoral response against Mtb at the level of active TB only. More studies on the species-specific impact of helminths on the adaptive humoral response against Mtb using a larger sample size, and in relation to TB disease severity, are needed.

Understanding and controlling asthma in Latin America: A review of recent research informed by the SCAALA programme

Asthma is an important health concern in Latin America (LA) where it is associated with variable prevalence and disease burden between countries. High prevalence and morbidity have been observed in some regions, particularly marginalized urban populations. Research over the past 10 years from LA has shown that childhood disease is primarily non-atopic. The attenuation of atopy may be explained by enhanced immune regulation induced by intense exposures to environmental factors such as childhood infections and poor environmental conditions of the urban poor. Non-atopic symptoms are associated with environmental and lifestyle factors including poor living conditions, respiratory infections, psychosocial stress, obesity, and a diet of highly processed foods. Ancestry (particularly African) and genetic factors increase asthma risk, and some of these factors may be specific to LA settings. Asthma in LA tends to be poorly controlled and depends on access to health care and medications. There is a need to improve management and access to medication through primary health care. Future research should consider the heterogeneity of asthma to identify relevant endotypes and underlying causes. The outcome of such research will need to focus on implementable strategies relevant to populations living in resource-poor settings where the disease burden is greatest.

New considerations of health disparities within allergy and immunology

The pandemic, political upheavals, and social justice efforts in our society have resulted in attention to persistent health disparities and the urgent need to address them. Using a scoping review, we describe published updates to address disparities and targets for interventions to improve gaps in care within allergy and immunology. These disparities-related studies provide a broad view of our current understanding of how social determinants of health threaten patient outcomes and our ability to advance health equity efforts in our field. We outline next steps to improve access to care and advance health equity for patients with allergic/immunologic diseases through actions taken at the individual, community, and policy levels, which could be applied outside of our field. Key among these are efforts to increase the diversity among our trainees, providers, and scientific teams and enhancing efforts to participate in advocacy work and public health interventions. Addressing health disparities requires advancing our understanding of the interplay between social and structural barriers to care and enacting the needed interventions in various key areas to effect change.

Recent progress in the genetic and epigenetic underpinnings of atopy

In the past 2 years, there continue to be advances in our understanding of the genetic and epigenetic underpinnings of atopy pertaining to disease risk and disease severity. The joint role of genetics and the environment has been emphasized in multiple studies. Combining genetics with family history, biomarkers, and comorbidities is further refining our ability to predict the development of individual atopic diseases as well as the advancement of the atopic march. Polygenic risk scores will be an important next step for the field moving toward clinical translation of the genetic findings thus far. A systems biology approach, as illustrated by studies of the microbiome and epigenome, will be necessary to fully understand disease development and to develop increasingly targeted therapeutics.

Zinc finger protein 33B demonstrates sex interaction with atopy-related markers in childhood asthma

BACKGROUND

Sex differences related to immune responses can influence atopic manifestations in childhood asthma. While genome-wide association studies have investigated a sex-specific genetic architecture of the immune response, gene-by-sex interactions have not been extensively analysed for atopy-related markers including allergy skin tests, IgE and eosinophils in asthmatic children.

METHODS

We performed a genome-wide gene-by-sex interaction analysis for atopy-related markers using whole-genome sequencing data based on 889 trios from the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) and 284 trios from the Childhood Asthma Management Program (CAMP). We also tested the findings in UK Biobank participants with self-reported childhood asthma. Furthermore, downstream analyses in GACRS integrated gene expression to disentangle observed associations.

RESULTS

Single nucleotide polymorphism (SNP) rs1255383 at 10q11.21 demonstrated a genome-wide significant gene-by-sex interaction (pinteraction=9.08×10-10) for atopy (positive skin test) with opposite direction of effects between females and males. In the UK Biobank participants with a history of childhood asthma, the signal was consistently observed with the same sex-specific effect directions for high eosinophil count (pinteraction=0.0058). Gene expression of ZNF33B (zinc finger protein 33B), located at 10q11.21, was moderately associated with atopy in girls, but not in boys.

CONCLUSIONS

We report SNPs in/near a zinc finger gene as novel sex-differential loci for atopy-related markers with opposite effect directions in females and males. A potential role for ZNF33B should be studied further as an important driver of sex-divergent features of atopy in childhood asthma.

Multi-omic association study identifies DNA methylation-mediated genotype and smoking exposure effects on lung function in children living in urban settings

Impaired lung function in early life is associated with the subsequent development of chronic respiratory disease. Most genetic associations with lung function have been identified in adults of European descent and therefore may not represent those most relevant to pediatric populations and populations of different ancestries. In this study, we performed genome-wide association analyses of lung function in a multiethnic cohort of children (n = 1,035) living in low-income urban neighborhoods. We identified one novel locus at the TDRD9 gene in chromosome 14q32.33 associated with percent predicted forced expiratory volume in one second (FEV1) (p = 2.4x10-9; βz = -0.31, 95% CI = -0.41- -0.21). Mendelian randomization and mediation analyses revealed that this genetic effect on FEV1 was partially mediated by DNA methylation levels at this locus in airway epithelial cells, which were also associated with environmental tobacco smoke exposure (p = 0.015). Promoter-enhancer interactions in airway epithelial cells revealed chromatin interaction loops between FEV1-associated variants in TDRD9 and the promoter region of the PPP1R13B gene, a stimulator of p53-mediated apoptosis. Expression of PPP1R13B in airway epithelial cells was significantly associated the FEV1 risk alleles (p = 1.3x10-5; β = 0.12, 95% CI = 0.06-0.17). These combined results highlight a potential novel mechanism for reduced lung function in urban youth resulting from both genetics and smoking exposure.

GWAS in the southern African context

Researchers would generally adjust for the possible confounding effect of population structure by considering global ancestry proportions or top principle components. Alternatively, researchers would conduct admixture mapping to increase the power to detect variants with an ancestry effect. This is sufficient in simple admixture scenarios, however, populations from southern Africa can be complex multi-way admixed populations. Duan et al. (2018) first described local ancestry adjusted allelic (LAAA) analysis as a robust method for discovering association signals, while producing minimal false positive hits. Their simulation study, however, was limited to a two-way admixed population. Realizing that their findings might not translate to other admixture scenarios, we simulated a three- and five-way admixed population to compare the LAAA model to other models commonly used in genome-wide association studies (GWAS). We found that, given our admixture scenarios, the LAAA model identifies the most causal variants in most of the phenotypes we tested across both the three-way and five-way admixed populations. The LAAA model also produced a high number of false positive hits which was potentially caused by the ancestry effect size that we assumed. Considering the extent to which the various models tested differed in their results and considering that the source of a given association is unknown, we recommend that researchers use multiple GWAS models when analysing populations with complex ancestry.

[Genetic aspects involved in asthma]

Asthma is an etiologically heterogeneous disease resulting from a complex interaction between genetic. The genetic aspects involved in asthma, which were analyzed from the perspective of the traditional model of multifactorial inheritance, were susceptibility, host factors, and environmental exposures. In the present paper, studies on their family aggregation, concordance in twins, and heritability were analized; as well as the current knowledge about candidate genes, genome wide association studies, and epigenomics contributions and other omic studies that have increased our knowledge about their pathophysiology and environmental interactions.

Asthma in the Americas: An Update: A Joint Perspective from the Brazilian Thoracic Society, Canadian Thoracic Society, Latin American Thoracic Society, and American Thoracic Society

Asthma affects a large number of people living in the Americas, a vast and diverse geographic region comprising 35 nations in the Caribbean and North, Central, and South America. The marked variability in the prevalence, morbidity, and mortality from asthma across and within nations in the Americas offers a unique opportunity to improve our understanding of the risk factors and management of asthma phenotypes and endotypes in children and adults. Moreover, a better assessment of the causes and treatment of asthma in less economically developed regions in the Americas would help diagnose and treat individuals migrating from those areas to Canada and the United States. In this focused review, we first assess the epidemiology of asthma, review known and potential risk factors, and examine commonalities and differences in asthma management across the Americas. We then discuss future directions in research and health policies to improve the prevention, diagnosis, and management of pediatric and adult asthma in the Americas, including standardized and periodic assessment of asthma burden across the region; large-scale longitudinal studies including omics and comprehensive environmental data on racially and ethnically diverse populations; and dissemination and implementation of guidelines for asthma management across the spectrum of disease severity. New initiatives should recognize differences in socioeconomic development and health care systems across the region while paying particular attention to novel or more impactful risk factors for asthma in the Americas, including indoor pollutants such as biomass fuel, tobacco use, infectious agents and the microbiome, and psychosocial stressor and chronic stress.

Unlocking the secrets of long non-coding RNAs in asthma

Asthma is a very heterozygous disease, divided in subtypes, such as eosinophilic and neutrophilic asthma. Phenotyping and endotyping of patients, especially patients with severe asthma who are refractory to standard treatment, are crucial in asthma management and are based on a combination of clinical and biological features. Nevertheless, the quest remains to find better biomarkers that distinguish asthma subtypes in a more clear and objective manner and to find new therapeutic targets to treat people with therapy-resistant asthma. In the past, research to identify asthma subtypes mainly focused on expression profiles of protein-coding genes. However, advances in RNA-sequencing technologies and the discovery of non-coding RNAs as important post-transcriptional regulators have provided an entire new field of research opportunities in asthma. This review focusses on long non-coding RNAs (lncRNAs) in asthma; these are non-coding RNAs with a length of more than 200 nucleotides. Many lncRNAs are differentially expressed in asthma, and several have been associated with asthma severity or inflammatory phenotype. Moreover, in vivo and in vitro functional studies have identified the mechanisms of action of specific lncRNAs. Although lncRNAs remain not widely studied in asthma, the current studies show the potential of lncRNAs as biomarkers and therapeutic targets as well as the need for further research.

Targeted analysis of genomic regions enriched in African ancestry reveals novel classical HLA alleles associated with asthma in Southwestern Europeans

Despite asthma has a considerable genetic component, an important proportion of genetic risks remain unknown, especially for non-European populations. Canary Islanders have the largest African genetic ancestry observed among Southwestern Europeans and the highest asthma prevalence in Spain. Here we examined broad chromosomal regions previously associated with an excess of African genetic ancestry in Canary Islanders, with the aim of identifying novel risk variants associated with asthma susceptibility. In a two-stage cases-control study, we revealed a variant within HLA-DQB1 significantly associated with asthma risk (rs1049213, meta-analysis p = 1.30 × 10^–7, OR [95% CI] = 1.74 [1.41–2.13]) previously associated with asthma and broad allergic phenotype. Subsequent fine-mapping analyses of classical HLA alleles revealed a novel allele significantly associated with asthma protection (HLA-DQA1*01:02, meta-analysis p = 3.98 × 10^–4, OR [95% CI] = 0.64 [0.50–0.82]) that had been linked to infectious and autoimmune diseases, and peanut allergy. HLA haplotype analyses revealed a novel haplotype DQA1*01:02-DQB1*06:04 conferring asthma protection (meta-analysis p = 4.71 × 10^–4, OR [95% CI] = 0.47 [0.29– 0.73]).

Multiethnic genome-wide and HLA association study of total serum IgE level

BACKGROUND

Total serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE.

OBJECTIVE

We aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901).

METHODS

We performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data.

RESULTS

We identified 6 loci at genome-wide significance (P < 5 × 10^-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (P_discovery = 2 × 10^-4; P_replication = 5 × 10^-4; P_{discovery+replication} = 4 × 10^-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (P_{Hispanic/Latino discovery+replication} = 8 × 10^-8).

CONCLUSION

We performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.