Genome-wide analysis highlights contribution of immune system pathways to the genetic architecture of asthma

Causal role of immune cells in asthma: a Mendelian randomization study

BACKGROUND

Immune cells may have a significant role in the pathophysiology of asthma, according to increasing evidence, although it is yet unclear how immune cells cause asthma. Therefore, we aimed to use Mendelian randomization (MR) methods to investigate this causal relationship.

METHODS

This study explored the causal effects between immune cells and asthma using a two-sample MR technique. Using publicly available genetic data, the causal connection between asthma risk and 731 immune cell phenotypes was investigated. Sensitivity analysis guaranteed the results' stability. To further evaluate the existence of reverse causality, we employed reverse MR analysis.

RESULTS

According to the inverse-variance weighted (IVW) method, five immune cell phenotypes were found to be statistically significantly associated with asthma risk (p < 0.001). Among them, TCRgd %T cell (OR = 0.968, 95%CI = 0.951 - 0.986), TCRgd %lymphocyte (OR = 0.978, 95%CI = 0.965 - 0.991), HLA DR + NK AC (OR = 0.966, 95% CI = 0.947 - 0.986) and CD3 on CD4 Treg (OR = 0.956, 95%CI= 0.931 - 0.981), four phenotypes that resulted in a decreased risk of asthma. CD25 on transitional (OR = 1.033, 95%CI = 1.014 - 1.052) resulted in an increased risk of asthma. Reverse MR analysis revealed that asthma increases HLA DR + NK AC levels (p < 0.05).

CONCLUSION

The results of MR analysis showed a causal relationship between immune cell phenotype and asthma risk, which provides a direction for future asthma treatment.

Genetics of chronic respiratory disease

Chronic respiratory diseases, such as chronic obstructive pulmonary disease (COPD), asthma and interstitial lung diseases are frequently occurring disorders with a polygenic basis that account for a large global burden of morbidity and mortality. Recent large-scale genetic epidemiology studies have identified associations between genetic variation and individual respiratory diseases and linked specific genetic variants to quantitative traits related to lung function. These associations have improved our understanding of the genetic basis and mechanisms underlying common lung diseases. Moreover, examining the overlap between genetic associations of different respiratory conditions, along with evidence for gene-environment interactions, has yielded additional biological insights into affected molecular pathways. This genetic information could inform the assessment of respiratory disease risk and contribute to stratified treatment approaches.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Principled distillation of UK Biobank phenotype data reveals underlying structure in human variation

Data within biobanks capture broad yet detailed indices of human variation, but biobank-wide insights can be difficult to extract due to complexity and scale. Here, using large-scale factor analysis, we distill hundreds of variables (diagnoses, assessments and survey items) into 35 latent constructs, using data from unrelated individuals with predominantly estimated European genetic ancestry in UK Biobank. These factors recapitulate known disease classifications, disentangle elements of socioeconomic status, highlight the relevance of psychiatric constructs to health and improve measurement of pro-health behaviours. We go on to demonstrate the power of this approach to clarify genetic signal, enhance discovery and identify associations between underlying phenotypic structure and health outcomes. In building a deeper understanding of ways in which constructs such as socioeconomic status, trauma, or physical activity are structured in the dataset, we emphasize the importance of considering the interwoven nature of the human phenome when evaluating public health patterns.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

A mutation in Themis contributes to anaphylaxis severity following oral peanut challenge in CC027 mice

BACKGROUND

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized Collaborative Cross strain CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, in contrast to C3H/HeJ (C3H) mice.

OBJECTIVE

This study aimed to determine the genetic basis of orally induced anaphylaxis to peanut in CC027 mice.

METHODS

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 mice and 5 additional Collaborative Cross strains.

RESULTS

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis and 4% having severe anaphylaxis. There were 8 genetic loci associated with variation in response to peanut challenge-6 associated with anaphylaxis (temperature decrease) and 2 associated with peanut-specific IgE levels. There were 2 major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis gene. Consistent with described functions of Themis, we found that CC027 mice have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

CONCLUSIONS

Our results demonstrate a key role for Themis in the orally reactive CC027 mouse model of peanut allergy.

Advances in non-type 2 severe asthma: from molecular insights to novel treatment strategies

Asthma is a prevalent pulmonary disease that affects more than 300 million people worldwide and imposes a substantial economic burden. While medication can effectively control symptoms in some patients, severe asthma attacks, driven by airway inflammation induced by environmental and infectious exposures, continue to be a major cause of asthma-related mortality. Heterogeneous phenotypes of asthma include type 2 (T2) and non-T2 asthma. Non-T2 asthma is often observed in patients with severe and/or steroid-resistant asthma. This review covers the molecular mechanisms, clinical phenotypes, causes and promising treatments of non-T2 severe asthma. Specifically, we discuss the signalling pathways for non-T2 asthma including the activation of inflammasomes, interferon responses and interleukin-17 pathways, and their contributions to the subtypes, progression and severity of non-T2 asthma. Understanding the molecular mechanisms and genetic determinants underlying non-T2 asthma could form the basis for precision medicine in severe asthma treatment.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Shared molecular mechanisms and transdiagnostic potential of neurodevelopmental disorders and immune disorders

The co-occurrence and familial clustering of neurodevelopmental disorders and immune disorders suggest shared genetic risk factors. Based on genome-wide association summary statistics from five neurodevelopmental disorders and four immune disorders, we conducted genome-wide, local genetic correlation and polygenic overlap analysis. We further performed a cross-trait GWAS meta-analysis. Pleotropic loci shared between the two categories of diseases were mapped to candidate genes using multiple algorithms and approaches. Significant genetic correlations were observed between neurodevelopmental disorders and immune disorders, including both positive and negative correlations. Neurodevelopmental disorders exhibited higher polygenicity compared to immune disorders. Around 50%-90% of genetic variants of the immune disorders were shared with neurodevelopmental disorders. The cross-trait meta-analysis revealed 154 genome-wide significant loci, including 8 novel pleiotropic loci. Significant associations were observed for 30 loci with both types of diseases. Pathway analysis on the candidate genes at these loci revealed common pathways shared by the two types of diseases, including neural signaling, inflammatory response, and PI3K-Akt signaling pathway. In addition, 26 of the 30 lead SNPs were associated with blood cell traits. Neurodevelopmental disorders exhibit complex polygenic architecture, with a subset of individuals being at a heightened genetic risk for both neurodevelopmental and immune disorders. The identification of pleiotropic loci has important implications for exploring opportunities for drug repurposing, enabling more accurate patient stratification, and advancing genomics-informed precision in the medical field of neurodevelopmental disorders.

Multi-omics in nasal epithelium reveals three axes of dysregulation for asthma risk in the African Diaspora populations

Asthma has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omic signatures of asthma focusing on populations of African ancestry. RNASeq and DNA methylation data are generated from nasal epithelium including cases (current asthma, N = 253) and controls (never-asthma, N = 283) from 7 different geographic sites to identify differentially expressed genes (DEGs) and gene networks. We identify 389 DEGs; the top DEG, FN1, was downregulated in cases (q = 3.26 × 10-9) and encodes fibronectin which plays a role in wound healing. The top three gene expression modules implicate networks related to immune response (CEACAM5; p = 9.62 × 10-16 and CPA3; p = 2.39 × 10-14) and wound healing (FN1; p = 7.63 × 10-9). Multi-omic analysis identifies FKBP5, a co-chaperone of glucocorticoid receptor signaling known to be involved in drug response in asthma, where the association between nasal epithelium gene expression is likely regulated by methylation and is associated with increased use of inhaled corticosteroids. This work reveals molecular dysregulation on three axes - increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response - that may play a critical role in asthma within the African Diaspora.

Breaking down causes, consequences, and mediating effects of telomere length variation on human health

BACKGROUND

Telomeres form repeated DNA sequences at the ends of chromosomes, which shorten with each cell division. Yet, factors modulating telomere attrition and the health consequences thereof are not fully understood. To address this, we leveraged data from 326,363 unrelated UK Biobank participants of European ancestry.

RESULTS

Using linear regression and bidirectional univariable and multivariable Mendelian randomization (MR), we elucidate the relationships between leukocyte telomere length (LTL) and 142 complex traits, including diseases, biomarkers, and lifestyle factors. We confirm that telomeres shorten with age and show a stronger decline in males than in females, with these factors contributing to the majority of the 5.4% of LTL variance explained by the phenome. MR reveals 23 traits modulating LTL. Smoking cessation and high educational attainment associate with longer LTL, while weekly alcohol intake, body mass index, urate levels, and female reproductive events, such as childbirth, associate with shorter LTL. We also identify 24 traits affected by LTL, with risk for cardiovascular, pulmonary, and some autoimmune diseases being increased by short LTL, while longer LTL increased risk for other autoimmune conditions and cancers. Through multivariable MR, we show that LTL may partially mediate the impact of educational attainment, body mass index, and female age at childbirth on proxied lifespan.

CONCLUSIONS

Our study sheds light on the modulators, consequences, and the mediatory role of telomeres, portraying an intricate relationship between LTL, diseases, lifestyle, and socio-economic factors.

Genetic mapping across autoimmune diseases reveals shared associations and mechanisms

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

Blocking CD226 regulates type 2 innate lymphoid cell effector function and alleviates airway hyperreactivity

BACKGROUND

Type 2 innate lymphoid cells (ILC2s) play a pivotal role in type 2 asthma. CD226 is a costimulatory molecule involved in various inflammatory diseases.

OBJECTIVE

We aimed to investigate CD226 expression and function within human and mouse ILC2s, and to assess the impact of targeting CD226 on ILC2-mediated airway hyperreactivity (AHR).

METHODS

We administered IL-33 intranasally to wild-type mice, followed by treatment with anti-CD226 antibody or isotype control. Pulmonary ILC2s were sorted for ex vivo analyses through RNA sequencing and flow cytometry. Next, we evaluated the effects of CD226 on AHR and lung inflammation in wild-type and Rag2-/- mice. Additionally, we compared peripheral ILC2s from healthy donors and asthmatic patients to ascertain the role of CD226 in human ILC2s.

RESULTS

Our findings demonstrated an inducible expression of CD226 in activated ILC2s, enhancing their cytokine secretion and effector functions. Mechanistically, CD226 alters intracellular metabolism and enhances PI3K/AKT and MAPK signal pathways. Blocking CD226 ameliorates ILC2-dependent AHR in IL-33 and Alternaria alternata-induced models. Interestingly, CD226 is expressed and inducible in human ILC2s, and its blocking reduces cytokine production. Finally, we showed that peripheral ILC2s in asthmatic patients exhibited elevated CD226 expression compared to healthy controls.

CONCLUSION

Our findings underscore the potential of CD226 as a novel therapeutic target in ILC2s, presenting a promising avenue for ameliorating AHR and allergic asthma.

Effect of air pollution on asthma

Asthma is a chronic inflammatory airway disease characterized by respiratory symptoms, variable airflow obstruction, bronchial hyperresponsiveness, and airway inflammation. Exposure to air pollution has been linked to an increased risk of asthma development and exacerbation. This review aims to comprehensively summarize recent data on the impact of air pollution on asthma development and exacerbation. Specifically, we reviewed the effects of air pollution on the pathogenic pathways of asthma, including type 2 and non-type 2 inflammatory responses, and airway epithelial barrier dysfunction. Air pollution promotes the release of epithelial cytokines, driving TH2 responses, and induces oxidative stress and the production of proinflammatory cytokines. The enhanced type 2 inflammation, furthered by air pollution-induced dysfunction of the airway epithelial barrier, may be associated with the exacerbation of asthma. Disruption of the TH17/regulatory T cell balance by air pollutants is also related to asthma exacerbation. As the effects of air pollution exposure may accumulate over time, with potentially stronger impacts in the development of asthma during certain sensitive life periods, we also reviewed the effects of air pollution on asthma across the lifespan. Future research is needed to better characterize the sensitive period contributing to the development of air pollution-induced asthma and to map air pollution-associated epigenetic biomarkers contributing to the epigenetic ages onto asthma-related genes.

Risk of atopic dermatitis and the atopic march paradigm in children of mothers with atopic illnesses: A birth cohort study from the United Kingdom

BACKGROUND

Atopic dermatitis (AD) is thought to precede the onset of other allergic illness (OAI) in a temporal progression (ie, atopic march), yet the timing and progression has been questioned. It is also unclear how parental allergic illness impacts the development of these illnesses in offspring.

OBJECTIVE

(1) Explore risk of incident AD and (2) timing of allergic disease onset in children of mothers with AD compared with mothers without AD from the United Kingdom.

METHODS

We created a birth cohort of mother-child pairs using IQVIA Medical Research Data database and developed Cox proportional models to examine the above associations (hazard ratio, HR [95% confidence interval, CI]).

RESULTS

Among 1,224,243 child-mother pairs, mean child (standard deviation) follow-up time was 10.8 (8.3) years and 50.1% were males (N = 600,905). Children were 59% (HR = 1.59 [1.57, 1.60]) more likely to have AD if their mothers had AD compared with no AD with mean age of first AD diagnosis at 3.3 (4.8) years. Most children with any diagnosis of AD present with AD first (91.0%); however, in those with asthma, only 67.8% developed AD first.

CONCLUSION

Children born to mothers with AD are more prone to develop AD and some develop OAI first, suggesting that not all follow the same sequential pathway.

Sexual dimorphism in the gut microbiota and sexual dimorphism in chronic diseases: Association or causation?

Understanding the sexual dimorphism in diseases is essential to investigate the pathogenesis of some chronic diseases (e.g., autoimmune diseases, etc). The gut microbiota has been found to show a notable impact on the pathology of several chronic diseases in recent years. Intriguingly, the composition of the gut microbiota varies between sexes. Here, we reviewed 'facts and fiction' regarding sexual dimorphism in chronic diseases and sexual dimorphism in the gut microbiota respectively. The association and causative relationship between them aiming to elucidate the pathological mechanisms of sexual dimorphism in chronic diseases were further explored. The development of gender-special food products based on the sexual dimorphism in the gut microbiota were recommended, which would be beneficial to facilitating the personalized treatment.

Omics approaches in asthma research: Challenges and opportunities

Asthma, a chronic respiratory disease with a global prevalence of approximately 300 million individuals, presents a significant societal and economic burden. This multifaceted syndrome exhibits diverse clinical phenotypes and pathogenic endotypes influenced by various factors. The advent of omics technologies has revolutionized asthma research by delving into the molecular foundation of the disease to unravel its underlying mechanisms. Omics technologies are employed to systematically screen for potential biomarkers, encompassing genes, transcripts, methylation sites, proteins, and even the microbiome components. This review provides an insightful overview of omics applications in asthma research, with a special emphasis on genetics, transcriptomics, epigenomics, and the microbiome. We explore the cutting-edge methods, discoveries, challenges, and potential future directions in the realm of asthma omics research. By integrating multi-omics and non-omics data through advanced statistical techniques, we aspire to advance precision medicine in asthma, guiding diagnosis, risk assessment, and personalized treatment strategies for this heterogeneous condition.

Plasma protein signatures of adult asthma

BACKGROUND

Adult asthma is complex and incompletely understood. Plasma proteomics is an evolving technique that can both generate biomarkers and provide insights into disease mechanisms. We aimed to identify plasma proteomic signatures of adult asthma.

METHODS

Protein abundance in plasma was measured in individuals from the Agricultural Lung Health Study (ALHS) (761 asthma, 1095 non-case) and the Atherosclerosis Risk in Communities study (470 asthma, 10,669 non-case) using the SOMAScan 5K array. Associations with asthma were estimated using covariate adjusted logistic regression and meta-analyzed using inverse-variance weighting. Additionally, in ALHS, we examined phenotypes based on both asthma and seroatopy (asthma with atopy (n = 207), asthma without atopy (n = 554), atopy without asthma (n = 147), compared to neither (n = 948)).

RESULTS

Meta-analysis of 4860 proteins identified 115 significantly (FDR<0.05) associated with asthma. Multiple signaling pathways related to airway inflammation and pulmonary injury were enriched (FDR<0.05) among these proteins. A proteomic score generated using machine learning provided predictive value for asthma (AUC = 0.77, 95% CI = 0.75-0.79 in training set; AUC = 0.72, 95% CI = 0.69-0.75 in validation set). Twenty proteins are targeted by approved or investigational drugs for asthma or other conditions, suggesting potential drug repurposing. The combined asthma-atopy phenotype showed significant associations with 20 proteins, including five not identified in the overall asthma analysis.

CONCLUSION

This first large-scale proteomics study identified over 100 plasma proteins associated with current asthma in adults. In addition to validating previous associations, we identified many novel proteins that could inform development of diagnostic biomarkers and therapeutic targets in asthma management.

Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease

The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.Communicated by Ramaswamy H. Sarma.

TGFβ prevents IgE-mediated allergic disease by restraining T follicular helper 2 differentiation

Allergic diseases are common, affecting more than 20% of the population. Genetic variants in the TGFβ pathway are strongly associated with atopy. To interrogate the mechanisms underlying this association, we examined patients and mice with Loeys-Dietz syndrome (LDS) who harbor missense mutations in the kinase domain of TGFΒR1/2. We demonstrate that LDS mutations lead to reduced TGFβ signaling and elevated total and allergen-specific IgE, despite the presence of wild-type T regulatory cells in a chimera model. Germinal center activity was enhanced in LDS and characterized by a selective increase in type 2 follicular helper T cells (TFH2). Expression of Pik3cg was increased in LDS TFH cells and associated with reduced levels of the transcriptional repressor SnoN. PI3Kγ/mTOR signaling in LDS naïve CD4+ T cells was elevated after T cell receptor cross-linking, and pharmacologic inhibition of PI3Kγ or mTOR prevented exaggerated TFH2 and antigen-specific IgE responses after oral antigen exposure in an adoptive transfer model. Naïve CD4+ T cells from nonsyndromic allergic individuals also displayed decreased TGFβ signaling, suggesting that our mechanistic discoveries may be broadly relevant to allergic patients in general. Thus, TGFβ plays a conserved, T cell-intrinsic, and nonredundant role in restraining TFH2 development via the PI3Kγ/mTOR pathway and thereby protects against allergic disease.

Identification of differences in CD4+ T-cell gene expression between people with asthma and healthy controls

Functional enrichment analysis of genome-wide association study (GWAS)-summary statistics has suggested that CD4+ T-cells play an important role in asthma pathogenesis. Despite this, CD4+ T-cells are under-represented in asthma transcriptome studies. To fill the gap, 3'-RNA-Seq was used to generate gene expression data on CD4+ T-cells (isolated within 2 h from collection) from peripheral blood from participants with well-controlled asthma (n = 32) and healthy controls (n = 11). Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify sets of co-expressed genes (modules) associated with the asthma phenotype. We identified three modules associated with asthma, which are strongly enriched for GWAS-identified asthma genes, antigen processing/presentation and immune response to viral infections. Through integration of publicly available eQTL and GWAS summary statistics (colocalisation), and protein-protein interaction (PPI) data, we identified PTPRC, a potential druggable target, as a putative master regulator of the asthma gene-expression profiles. Using a co-expression network approach, with integration of external genetic and PPI data, we showed that CD4+ T-cells from peripheral blood from asthmatics have different expression profiles, albeit small in magnitude, compared to healthy controls, for sets of genes involved in immune response to viral infections (upregulated) and antigen processing/presentation (downregulated).

Genetic relationship between the immune system and autism

Autism spectrum disorder (ASD) is a common and complex neurodevelopmental condition. The pathophysiology of ASD is poorly defined; however, it includes a strong genetic component and there is increasing evidence to support a role of immune dysregulation. Nonetheless, it is unclear which immune phenotypes link to ASD through genetics. Hence, we investigated the genetic correlation between ASD and diverse classes of immune conditions and markers; and if these immune-related genetic factors link to specific autistic-like traits in the population. We estimated global and local genetic correlations between ASD (n = 55,420) and 11 immune phenotypes (n = 14,256-755,406) using genome-wide association study summary statistics. Subsequently, polygenic scores (PGS) for these immune phenotypes were calculated in a population-based sample (n = 2487) and associated to five autistic-like traits (i.e., attention to detail, childhood behaviour, imagination, rigidity, social skills), and a total autistic-like traits score. Sex-stratified PGS analyses were also performed. At the genome-wide level, ASD was positively correlated with allergic diseases (ALG), and negatively correlated with lymphocyte count, rheumatoid arthritis (RA), and systemic lupus erythematosus (SLE) (FDR-p = 0.01-0.02). At the local genetic level, ASD was correlated with RA, C-reactive protein, and granulocytes and lymphocyte counts (p = 5.8 × 10-6-0.002). In the general population sample, increased genetic liability for SLE, RA, ALG, and lymphocyte levels, captured by PGS, was associated with the total autistic score and with rigidity and childhood behaviour (FDR-p = 0.03). In conclusion, we demonstrated a genetic relationship between ASD and immunity that depends on the type of immune phenotype considered; some increase likelihood whereas others may potentially help build resilience. Also, this relationship may be restricted to specific genetic loci and link to specific autistic dimensions (e.g., rigidity).

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease

BACKGROUND

Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features.

OBJECTIVE

We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO).

METHODS

We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry.

RESULTS

In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants.

CONCLUSIONS

PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Perinatal, obstetric and parental risk factors for asthma in the offspring throughout childhood: a longitudinal cohort study

OBJECTIVES

Asthma is a common chronic and burdensome disease which typically begins in childhood. The aim of this study was to assess perinatal and obstetric factors which may increase the risk of developing asthma in the offspring.

METHODS

Data from five consecutive waves (n=7,073 children, from birth to 15 years old) from a nationally-representative birth cohort of people born in the United Kingdom between 2000 and 2002, the Millennium Cohort Study (MCS), were used. The Kaplan-Meier survival curve was used to graphically display the risk of developing asthma from early childhood to adolescence. The Z-based Wald test was used to prove significant covariate loading.

RESULTS

Cox regression analyzing the influence of covariates on asthma development risk showed a significant likelihood ratio test, χ2(18)=899.30, p<0.01. A parent with asthma (OR=2.02, p<0.01), a younger maternal age at delivery (OR=0.98, p<0.05), and the use of assisted reproductive technology (OR=1.43, p<0.05) were associated with an increased risk of developing asthma in the offspring.

CONCLUSIONS

Perinatal factors (a younger maternal age, assisted reproductive technology) and a parental factor (a parent with asthma) increased the risk for developing asthma in the offspring.

Can polygenic risk scores help explain disease prevalence differences around the world? A worldwide investigation

Complex disorders are caused by a combination of genetic, environmental and lifestyle factors, and their prevalence can vary greatly across different populations. The extent to which genetic risk, as identified by Genome Wide Association Study (GWAS), correlates to disease prevalence in different populations has not been investigated systematically. Here, we studied 14 different complex disorders and explored whether polygenic risk scores (PRS) based on current GWAS correlate to disease prevalence within Europe and around the world. A clear variation in GWAS-based genetic risk was observed based on ancestry and we identified populations that have a higher genetic liability for developing certain disorders. We found that for four out of the 14 studied disorders, PRS significantly correlates to disease prevalence within Europe. We also found significant correlations between worldwide disease prevalence and PRS for eight of the studied disorders with Multiple Sclerosis genetic risk having the highest correlation to disease prevalence. Based on current GWAS results, the across population differences in genetic risk for certain disorders can potentially be used to understand differences in disease prevalence and identify populations with the highest genetic liability. The study highlights both the limitations of PRS based on current GWAS but also the fact that in some cases, PRS may already have high predictive power. This could be due to the genetic architecture of specific disorders or increased GWAS power in some cases.

Identifying genetic variants that influence the abundance of cell states in single-cell data

To understand genetic mechanisms driving disease, it is essential but difficult to map how risk alleles affect the composition of cells present in the body. Single-cell profiling quantifies granular information about tissues, but variant-associated cell states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce GeNA (Genotype-Neighborhood Associations), a statistical tool to identify cell state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of scRNA-seq peripheral blood profiling from 969 individuals,1GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (p=1.96×10-11) associates with increased abundance of NK cells expressing TNF-α response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-TNF treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

From preschool wheezing to asthma: Environmental determinants

Wheezing is common among preschool children, representing a group of highly heterogeneous conditions with varying natural history. Several phenotypes of wheezing have been proposed to facilitate the identification of young children who are at risk of subsequent development of asthma. Epidemiological and immunological studies across different populations have revealed the key role of environmental factors in influencing the progression from preschool wheezing to childhood asthma. Significant risk factors include severe respiratory infections, allergic sensitization, and exposure to tobacco smoke. In contrast, a farming/rural environment has been linked to asthma protection in both human and animal studies. Early and intense exposures to microorganisms and microbial metabolites have been demonstrated to alter host immune responses to allergens and viruses, thereby driving the trajectory away from wheezing illness and asthma. Ongoing clinical trials of candidate microbes and microbial products have shown promise in shaping the immune function to reduce episodes of viral-induced wheezing. Moreover, restoring immune training may be especially important for young children who had reduced microbial exposure due to pandemic restrictions. A comprehensive understanding of the role of modifiable environmental factors will pave the way for developing targeted prevention strategies for preschool wheezing and asthma.

Whole Genome Sequencing Identifies Novel Common and Low-Frequency Variants Associated With Age-Related Macular Degeneration

Purpose

To identify associations of common, low-frequency, and rare variants with advanced age-related macular degeneration (AMD) using whole genome sequencing (WGS).

Methods

WGS data were obtained for 2123 advanced AMD patients (participants of clinical trials for advanced AMD) and 2704 controls (participants of clinical trials for asthma [N = 2518] and Alzheimer's disease [N = 186]), and joint genotype calling was performed, followed by quality control of the dataset. Single variant association analyses were performed for all identified common, low-frequency, and rare variants. Gene-based tests were executed for rare and low-frequency variants using SKAT-O and three groups of variants based on putative impact information: (1) all variants, (2) modifier impact variants, and (3) high- and moderate-impact variants. To ascertain independence of the identified associations from previously reported AMD and asthma loci, conditional analyses were performed.

Results

Previously identified AMD variants at the CFH, ARMS2/HTRA1, APOE, and C3 loci were associated with AMD at a genome-wide significance level. We identified new single variant associations for common variants near the PARK7 gene and in the long non-coding RNA AC103876.1, and for a rare variant near the TENM3 gene. In addition, gene-based association analyses identified a burden of modifier variants in eight intergenic and gene-spanning regions and of high- and moderate-impact variants in the C3, CFHR5, SLC16A8, and CFI genes.

Conclusions

We describe the largest WGS study in AMD to date. We confirmed previously identified associations and identified several novel associations that are worth exploring in further follow-up studies.

Early-immune development in asthma: A review of the literature

This review presents a comprehensive examination of the various factors contributing to the immunopathogenesis of asthma from the prenatal to preschool period. We focus on the contributions of genetic and environmental components as well as the role of the nasal and gut microbiome on immune development. Predisposing genetic factors, including inherited genes associated with increased susceptibility to asthma, are discussed alongside environmental factors such as respiratory viruses and pollutant exposure, which can trigger or exacerbate asthma symptoms. Furthermore, the intricate interplay between the nasal and gut microbiome and the immune system is explored, emphasizing their influence on allergic immune development and response to environmental stimuli. This body of literature underscores the necessity of a comprehensive approach to comprehend and manage asthma, as it emphasizes the interactions of multiple factors in immune development and disease progression.

Identification of asthma-related genes using asthmatic blood eQTLs of Korean patients

BACKGROUND

More than 200 asthma-associated genetic variants have been identified in genome-wide association studies (GWASs). Expression quantitative trait loci (eQTL) data resources can help identify causal genes of the GWAS signals, but it can be difficult to find an eQTL that reflects the disease state because most eQTL data are obtained from normal healthy subjects.

METHODS

We performed a blood eQTL analysis using transcriptomic and genotypic data from 433 Korean asthma patients. To identify asthma-related genes, we carried out colocalization, Summary-based Mendelian Randomization (SMR) analysis, and Transcriptome-Wide Association Study (TWAS) using the results of asthma GWASs and eQTL data. In addition, we compared the results of disease eQTL data and asthma-related genes with two normal blood eQTL data from Genotype-Tissue Expression (GTEx) project and a Japanese study.

RESULTS

We identified 340,274 cis-eQTL and 2,875 eGenes from asthmatic eQTL analysis. We compared the disease eQTL results with GTEx and a Japanese study and found that 64.1% of the 2,875 eGenes overlapped with the GTEx eGenes and 39.0% with the Japanese eGenes. Following the integrated analysis of the asthmatic eQTL data with asthma GWASs, using colocalization and SMR methods, we identified 15 asthma-related genes specific to the Korean asthmatic eQTL data.

CONCLUSIONS

We provided Korean asthmatic cis-eQTL data and identified asthma-related genes by integrating them with GWAS data. In addition, we suggested these asthma-related genes as therapeutic targets for asthma. We envisage that our findings will contribute to understanding the etiological mechanisms of asthma and provide novel therapeutic targets.

Serum Allergen-Specific Immunoglobulin E in Cats with Inflammatory Bronchial Disease

The etiology of feline inflammatory bronchial disease is poorly understood. This study compares the degree of allergen-specific serum IgE responses between cats with feline asthma, chronic bronchitis, mixed inflammation, and clinically healthy cats (HCs). The retrospective case-control study used serum from eighteen cats with eosinophilic inflammation (EI), ten with neutrophilic inflammation (NI), six with mixed inflammation (MI), and fourteen HCs. Affected cats were categorized into groups based on bronchoalveolar lavage cytology. The measurement of IgE for 34 different allergens including fungal organisms, weeds, grasses, trees, mites, and insects was performed using an indirect ELISA. Positive reactions to allergens were detected in the serum of 17/18 cats with EI, 8/10 with NI, 6/6 with MI, and 11/14 HCs (p = 0.364). When overall positive reactions were compared between groups, cats with MI (p = <0.01) had significantly more positive reactions against mite allergens than HCs. Blood eosinophils inversely correlated with the absolute amount of allergen-specific serum IgE expressed in ELISA absorbance units (EAs) (p = 0.014). Sensitization against dust mites seems to be more prevalent in cats with MI. However, positive IgE reactions can be observed in healthy and diseased cats, and, therefore, need to be interpreted in the light of clinical findings and environmental conditions of individual patients.

Revealing polygenic pleiotropy using genetic risk scores for asthma

In this study we examined how genetic risk for asthma associates with different features of the disease and with other medical conditions and traits. Using summary statistics from two multi-ancestry genome-wide association studies of asthma, we modeled polygenic risk scores (PRSs) and validated their predictive performance in the UK Biobank. We then performed phenome-wide association studies of the asthma PRSs with 371 heritable traits in the UK Biobank. We identified 228 total significant associations across a variety of organ systems, including associations that varied by PRS model, sex, age of asthma onset, ancestry, and human leukocyte antigen region alleles. Our results highlight pervasive pleiotropy between asthma and numerous other traits and conditions and elucidate pathways that contribute to asthma and its comorbidities.

Circulating CC16 and Asthma: A Population-based, Multicohort Study from Early Childhood through Adult Life

Rationale: Club cell secretory protein (CC16) is an antiinflammatory protein highly expressed in the airways. CC16 deficiency has been associated with lung function deficits, but its role in asthma has not been established conclusively. Objectives: To determine 1) the longitudinal association of circulating CC16 with the presence of active asthma from early childhood through adult life and 2) whether CC16 in early childhood predicts the clinical course of childhood asthma into adult life. Methods: We assessed the association of circulating CC16 and asthma in three population-based birth cohorts: the Tucson Children's Respiratory Study (years 6-36; total participants, 814; total observations, 3,042), the Swedish Barn/Children, Allergy, Milieu, Stockholm, Epidemiological survey (years 8-24; total participants, 2,547; total observations, 3,438), and the UK Manchester Asthma and Allergy Study (years 5-18; total participants, 745; total observations, 1,626). Among 233 children who had asthma at the first survey in any of the cohorts, baseline CC16 was also tested for association with persistence of symptoms. Measurements and Main Results: After adjusting for covariates, CC16 deficits were associated with increased risk for the presence of asthma in all cohorts (meta-analyzed adjusted odds ratio per 1-SD CC16 decrease, 1.20; 95% confidence interval [CI], 1.12-1.28; P < 0.0001). The association was particularly strong for asthma with frequent symptoms (meta-analyzed adjusted relative risk ratio, 1.40; 95% CI, 1.24-1.57; P < 0.0001), was confirmed for both atopic and nonatopic asthma, and was independent of lung function impairment. After adjustment for known predictors of persistent asthma, children with asthma in the lowest CC16 tertile had a nearly fourfold increased risk for having frequent symptoms persisting into adult life compared with children with asthma in the other two CC16 tertiles (meta-analyzed adjusted odds ratio, 3.72; 95% CI, 1.78-7.76; P < 0.0001). Conclusions: Circulating CC16 deficits are associated with the presence of asthma with frequent symptoms from childhood through midadult life and predict the persistence of asthma symptoms into adulthood. These findings support a possible protective role of CC16 in asthma and its potential use for risk stratification.

Genetic analyses of chr11p15.5 region identify MUC5AC-MUC5B associated with asthma-related phenotypes

OBJECTIVE

Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) in chr11p15.5 region associated with asthma and idiopathic interstitial pneumonias (IIPs). We sought to identify functional genes for asthma by combining SNPs and mRNA expression in bronchial epithelial cells (BEC) in the Severe Asthma Research Program (SARP).

METHODS

Correlation analyses of mRNA expression of six candidate genes (AP2A2, MUC6, MUC2, MUC5AC, MUC5B, and TOLLIP) and asthma phenotypes were performed in the longitudinal cohort (n = 156) with RNAseq in BEC, and replicated in the cross-sectional cohort (n = 155). eQTL (n = 114) and genetic association analysis of asthma severity (426 severe vs. 531 non-severe asthma) were performed, and compared with previously published GWASs of IIPs and asthma.

RESULTS

Higher expression of AP2A2 and MUC5AC and lower expression of MUC5B in BEC were correlated with asthma, asthma exacerbations, and T2 biomarkers (P < 0.01). SNPs associated with asthma and IIPs in previous GWASs were eQTL SNPs for MUC5AC, MUC5B, or TOLLIP, however, they were not in strong linkage disequilibrium. The risk alleles for asthma or protective alleles for IIPs were associated with higher expression of MUC5AC and lower expression of MUC5B. rs11603634, rs12788104, and rs28415845 associated with moderate-to-severe asthma or adult onset asthma in previous GWASs were not associated with asthma severity (P > 0.8).

CONCLUSIONS

SNPs associated with asthma in chr11p15.5 region are not associated with asthma severity neither with IIPs. Higher expression of MUC5AC and lower expression of MUC5B are risk for asthma but protective for IIPs.

Large-scale plasma proteomics comparisons through genetics and disease associations

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Machine learning reveals genetic modifiers of the immune microenvironment of cancer

Heritability in the immune tumor microenvironment (iTME) has been widely observed yet remains largely uncharacterized. Here, we developed a machine learning approach to map iTME modifiers within loci from genome-wide association studies (GWASs) for breast cancer (BrCa) incidence. A random forest model was trained on a positive set of immune-oncology (I-O) targets, and then used to assign I-O target probability scores to 1,362 candidate genes in linkage disequilibrium with 155 BrCa GWAS loci. Cluster analysis of the most probable candidates revealed two subfamilies of genes related to effector functions and adaptive immune responses, suggesting that iTME modifiers impact multiple aspects of anticancer immunity. Two of the top ranking BrCa candidates, LSP1 and TLR1, were orthogonally validated as iTME modifiers using BrCa patient biopsies and comparative mapping studies, respectively. Collectively, these data demonstrate a robust and flexible framework for functionally fine-mapping GWAS risk loci to identify translatable therapeutic targets.

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.

Genetic Variants Associated With Hidradenitis Suppurativa

Importance

Hidradenitis suppurativa (HS) is a common and severely morbid chronic inflammatory skin disease that is reported to be highly heritable. However, the genetic understanding of HS is insufficient, and limited genome-wide association studies (GWASs) have been performed for HS, which have not identified significant risk loci.

Objective

To identify genetic variants associated with HS and to shed light on the underlying genes and genetic mechanisms.

Design, Setting, and Participants

This genetic association study recruited 753 patients with HS in the HS Program for Research and Care Excellence (HS ProCARE) at the University of North Carolina Department of Dermatology from August 2018 to July 2021. A GWAS was performed for 720 patients (after quality control) with controls from the Add Health study and then meta-analyzed with 2 large biobanks, UK Biobank (247 cases) and FinnGen (673 cases). Variants at 3 loci were tested for replication in the BioVU biobank (290 cases). Data analysis was performed from September 2021 to December 2022.

Main Outcomes and Measures

Main outcome measures are loci identified, with association of P < 1 × 10-8 considered significant.

Results

A total of 753 patients were recruited, with 720 included in the analysis. Mean (SD) age at symptom onset was 20.3 (10.57) years and at enrollment was 35.3 (13.52) years; 360 (50.0%) patients were Black, and 575 (79.7%) were female. In a meta-analysis of the 4 studies, 2 HS-associated loci were identified and replicated, with lead variants rs10512572 (P = 2.3 × 10-11) and rs17090189 (P = 2.1 × 10-8) near the SOX9 and KLF5 genes, respectively. Variants at these loci are located in enhancer regulatory elements detected in skin tissue.

Conclusions and Relevance

In this genetic association study, common variants associated with HS located near the SOX9 and KLF5 genes were associated with risk of HS. These or other nearby genes may be associated with genetic risk of disease and the development of clinical features, such as cysts, comedones, and inflammatory tunnels, that are unique to HS. New insights into disease pathogenesis related to these genes may help predict disease progression and novel treatment approaches in the future.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Priority index for asthma (PIA): In silico discovery of shared and distinct drug targets for adult- and childhood-onset disease

Asthma is a chronic disease that is caused by a combination of genetic risks and environmental triggers and can affect both adults and children. Genome-wide association studies have revealed partly distinct genetic architectures for its two age-of-onset subtypes (namely, adult-onset and childhood-onset). We reason that identifying shared and distinct drug targets between these subtypes may inform the development of subtype-specific therapeutic strategies. In attempting this, we here introduce Priority Index for Asthma or PIA, a genetics-led and network-driven drug target prioritisation tool for asthma. We demonstrate the validity of the tool in improving drug target prioritisation for asthma compared to the status quo methods, as well as in capturing the underlying etiology and existing therapeutics for the disease. We also illustrate how PIA can be used to prioritise drug targets for adult- and childhood-onset asthma, as well as to identify shared and distinct pathway crosstalk genes. Shared crosstalk genes are mostly involved in JAK-STAT signaling, with clinical evidence supporting that targeting this pathway may be a promising drug repurposing opportunity for both subtypes. Crosstalk genes specific to childhood-onset asthma are enriched for PI3K-AKT-mTOR signaling, and we identify genes that are already targeted by licensed medications as repurposed drug candidates for this subtype. We make all our results accessible and reproducible at http://www.genetictargets.com/PIA. Collectively, our study has significant implications for asthma computational medicine research and can guide the future development of subtype-specific therapeutic strategies for the disease.

The Impact of CC16 on Pulmonary Epithelial-Driven Host Responses during Mycoplasma pneumoniae Infection in Mouse Tracheal Epithelial Cells

Club Cell Secretory Protein (CC16) plays many protective roles within the lung; however, the complete biological functions, especially regarding the pulmonary epithelium during infection, remain undefined. We have previously shown that CC16-deficient (CC16-/-) mouse tracheal epithelial cells (MTECs) have enhanced Mp burden compared to CC16-sufficient (WT) MTECs; therefore, in this study, we wanted to further define how the pulmonary epithelium responds to infection in the context of CC16 deficiency. Using mass spectrometry and quantitative proteomics to analyze proteins secreted apically from MTECs grown at an air-liquid interface, we investigated the protective effects that CC16 elicits within the pulmonary epithelium during Mycoplasma pneumoniae (Mp) infection. When challenged with Mp, WT MTECs have an overall reduction in apical protein secretion, whereas CC16-/- MTECs have increased apical protein secretion compared to their unchallenged controls. Following Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) assessment, many of the proteins upregulated from CC16-/- MTECS (unchallenged and during Mp infection) were related to airway remodeling, which were not observed by WT MTECs. These findings suggest that CC16 may be important in providing protection within the pulmonary epithelium during respiratory infection with Mp, which is the major causative agent of community-acquired pneumoniae.

Vitamin D status, sleep patterns, genetic susceptibility, and the risk of incident adult-onset asthma: a large prospective cohort study

Introduction

Vitamin D has been known to be associated with asthma, particularly in children, while the evidence among adults is limited and inconclusive. This study aimed to investigate the association between serum, vitamin D concentrations, and the incidence of adult-onset asthma and also the modified effect caused by sleep patterns and genetic risks.

Methods

A prospective cohort study with 307,872 participants aged between 37 and 73 years was conducted based on the UK Biobank, with a median follow-up of 12 years. The Cox proportional hazard model was applied to evaluate the association between vitamin D status and incident adult-onset asthma, and the modified effect was investigated by conducting stratified analysis according to sleep pattern score and genetic risk score, and subgroup analyses were performed by sex, age, BMI, and smoking status as well.

Results

Individuals with optimal vitamin D concentration were associated with 11.1% reduced risk of incident asthma compared to those participants with deficient vitamin D (HR = 0.889; 95% CI: 0.820-0.964; p = 0.005). Moreover, stratification analysis demonstrated that the protective effect of vitamin D on asthma risk was modified by sleep patterns or genetic susceptibility, with the strongest protective effect being observed in the subpopulation with a moderate sleep pattern (HR = 0.883; 95% CI: 0.797-0.977; p = 0.016) and a moderate genetic risk (HR = 0.817; 95% CI: 0.711-0.938; p = 0.004). In subgroup analyses, the protective effect of optimal vitamin D levels was only significant among men, individuals younger than 60 years of age, overweight individuals, and current or previous smokers.

Conclusion

Increased serum vitamin D levels were associated with a lower risk of incident adult-onset asthma, and this association was modified by sleep patterns and genetic predisposition to some extent.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

An integrated cell atlas of the lung in health and disease

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

Mendelian Randomization Analyses of Chronic Immune-Mediated Diseases, Circulating Inflammatory Biomarkers, and Cytokines in Relation to Liver Cancer

Liver cancer is closely linked to chronic inflammation. While observational studies have reported positive associations between extrahepatic immune-mediated diseases and systemic inflammatory biomarkers and liver cancer, the genetic association between these inflammatory traits and liver cancer remains elusive and merits further investigation. We conducted a two-sample Mendelian randomization (MR) analysis, using inflammatory traits as exposures and liver cancer as the outcome. The genetic summary data of both exposures and outcome were retrieved from previous genome-wide association studies (GWAS). Four MR methods, including inverse-variance-weighted (IVW), MR-Egger regression, weighted-median, and weighted-mode methods, were employed to examine the genetic association between inflammatory traits and liver cancer. Nine extrahepatic immune-mediated diseases, seven circulating inflammatory biomarkers, and 187 inflammatory cytokines were analyzed in this study. The IVW method suggested that none of the nine immune-mediated diseases were associated with the risk of liver cancer, with odds ratios of 1.08 (95% CI 0.87-1.35) for asthma, 0.98 (95% CI 0.91-1.06) for rheumatoid arthritis, 1.01 (95% CI 0.96-1.07) for type 1 diabetes, 1.01 (95% CI 0.98-1.03) for psoriasis, 0.98 (95% CI 0.89-1.08) for Crohn's disease, 1.02 (95% CI 0.91-1.13) for ulcerative colitis, 0.91 (95% CI 0.74-1.11) for celiac disease, 0.93 (95% CI 0.84-1.05) for multiple sclerosis, and 1.05 (95% CI 0.97-1.13) for systemic lupus erythematosus. Similarly, no significant association was found between circulating inflammatory biomarkers and cytokines and liver cancer after correcting for multiple testing. The findings were consistent across all four MR methods used in this study. Our findings do not support a genetic association between extrahepatic inflammatory traits and liver cancer. However, larger-scale GWAS summary data and more genetic instruments are needed to confirm these findings.

Cross-disorder genetic analysis of immune diseases reveals distinct gene associations that converge on common pathways

Genome-wide association studies (GWAS) have mapped thousands of susceptibility loci associated with immune-mediated diseases. To assess the extent of the genetic sharing across nine immune-mediated diseases we apply genomic structural equation modelling to GWAS data from European populations. We identify three disease groups: gastrointestinal tract diseases, rheumatic and systemic diseases, and allergic diseases. Although loci associated with the disease groups are highly specific, they converge on perturbing the same pathways. Finally, we test for colocalization between loci and single-cell eQTLs derived from peripheral blood mononuclear cells. We identify the causal route by which 46 loci predispose to three disease groups and find evidence for eight genes being candidates for drug repurposing. Taken together, here we show that different constellations of diseases have distinct patterns of genetic associations, but that associated loci converge on perturbing different nodes in T cell activation and signalling pathways.

Human germline heterozygous gain-of-function STAT6 variants cause severe allergic disease

STAT6 (signal transducer and activator of transcription 6) is a transcription factor that plays a central role in the pathophysiology of allergic inflammation. We have identified 16 patients from 10 families spanning three continents with a profound phenotype of early-life onset allergic immune dysregulation, widespread treatment-resistant atopic dermatitis, hypereosinophilia with esosinophilic gastrointestinal disease, asthma, elevated serum IgE, IgE-mediated food allergies, and anaphylaxis. The cases were either sporadic (seven kindreds) or followed an autosomal dominant inheritance pattern (three kindreds). All patients carried monoallelic rare variants in STAT6 and functional studies established their gain-of-function (GOF) phenotype with sustained STAT6 phosphorylation, increased STAT6 target gene expression, and TH2 skewing. Precision treatment with the anti-IL-4Rα antibody, dupilumab, was highly effective improving both clinical manifestations and immunological biomarkers. This study identifies heterozygous GOF variants in STAT6 as a novel autosomal dominant allergic disorder. We anticipate that our discovery of multiple kindreds with germline STAT6 GOF variants will facilitate the recognition of more affected individuals and the full definition of this new primary atopic disorder.