Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling

G protein-coupled receptor (GPCR) pharmacogenomics

The field of pharmacogenetics, the investigation of the influence of one or more sequence variants on drug response phenotypes, is a special case of pharmacogenomics, a discipline that takes a genome-wide approach. Massively parallel, next generation sequencing (NGS), has allowed pharmacogenetics to be subsumed by pharmacogenomics with respect to the identification of variants associated with responders and non-responders, optimal drug response, and adverse drug reactions. A plethora of rare and common naturally-occurring GPCR variants must be considered in the context of signals from across the genome. Many fundamentals of pharmacogenetics were established for G protein-coupled receptor (GPCR) genes because they are primary targets for a large number of therapeutic drugs. Functional studies, demonstrating likely-pathogenic and pathogenic GPCR variants, have been integral to establishing models used for in silico analysis. Variants in GPCR genes include both coding and non-coding single nucleotide variants and insertion or deletions (indels) that affect cell surface expression (trafficking, dimerization, and desensitization/downregulation), ligand binding and G protein coupling, and variants that result in alternate splicing encoding isoforms/variable expression. As the breadth of data on the GPCR genome increases, we may expect an increase in the use of drug labels that note variants that significantly impact the clinical use of GPCR-targeting agents. We discuss the implications of GPCR pharmacogenomic data derived from the genomes available from individuals who have been well-phenotyped for receptor structure and function and receptor-ligand interactions, and the potential benefits to patients of optimized drug selection. Examples discussed include the renin-angiotensin system in SARS-CoV-2 (COVID-19) infection, the probable role of chemokine receptors in the cytokine storm, and potential protease activating receptor (PAR) interventions. Resources dedicated to GPCRs, including publicly available computational tools, are also discussed.

Stimulation of platelet P2Y1 receptors by different endogenous nucleotides leads to functional selectivity via biased signalling

BACKGROUND AND PURPOSE

Platelet function during inflammation is dependent on activation by endogenous nucleotides. Non-canonical signalling via the P2Y1 receptor is important for these non-thrombotic functions of platelets. However, apart from ADP, the role of other endogenous nucleotides acting as agonists at P2Y1 receptors is unknown. This study compared the effects of ADP, Ap3A, NAD+ , ADP-ribose, and Up4A on platelet functions contributing to inflammation or haemostasis.

EXPERIMENTAL APPROACH

Platelets obtained from healthy human volunteers were incubated with ADP, Ap3A, NAD+ , ADP-ribose, or Up4A, with aggregation and fibrinogen binding measured (examples of function during haemostasis) or before exposure to fMLP to measure platelet chemotaxis (an inflammatory function). In silico molecular docking of these nucleotides to the binding pocket of P2Y1 receptors was then assessed.

KEY RESULTS

Platelet aggregation and binding to fibrinogen induced by ADP was not mimicked by NAD+ , ADP-ribose, and Up4A. However, these endogenous nucleotides induced P2Y1 -dependent platelet chemotaxis, an effect that required RhoA and Rac-1 activity, but not canonical PLC activity. Analysis of molecular docking of the P2Y1 receptor revealed distinct differences of amino acid interactions and depth of fit within the binding pocket for Ap3A, NAD+ , ADP-ribose, or Up4A compared with ADP.

CONCLUSION AND IMPLICATIONS

Platelet function (aggregation vs motility) can be differentially modulated by biased-agonist activation of P2Y1 receptors. This may be due to the character of the ligand-binding pocket interaction. This has implications for future therapeutic strategies aimed to suppress platelet activation during inflammation without affecting haemostasis as is the requirement of current ant-platelet drugs.

LINKED ARTICLES

This article is part of a themed issue on Platelet purinergic receptor and non-thrombotic disease. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.4/issuetoc.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Identification of a carbohydrate recognition motif of purinergic receptors

As a major class of biomolecules, carbohydrates play indispensable roles in various biological processes. However, it remains largely unknown how carbohydrates directly modulate important drug targets, such as G-protein coupled receptors (GPCRs). Here, we employed P2Y purinoceptor 14 (P2Y14), a drug target for inflammation and immune responses, to uncover the sugar nucleotide activation of GPCRs. Integrating molecular dynamics simulation with functional study, we identified the uridine diphosphate (UDP)-sugar-binding site on P2Y14, and revealed that a UDP-glucose might activate the receptor by bridging the transmembrane (TM) helices 2 and 7. Between TM2 and TM7 of P2Y14, a conserved salt bridging chain (K2.60-D2.64-K7.35-E7.36 [KDKE]) was identified to distinguish different UDP-sugars, including UDP-glucose, UDP-galactose, UDP-glucuronic acid, and UDP-N-acetylglucosamine. We identified the KDKE chain as a conserved functional motif of sugar binding for both P2Y14 and P2Y purinoceptor 12 (P2Y12), and then designed three sugar nucleotides as agonists of P2Y12. These results not only expand our understanding for activation of purinergic receptors but also provide insights for the carbohydrate drug development for GPCRs.

Omics-driven investigation of the biology underlying intrinsic submaximal working capacity and its trainability

Submaximal exercise capacity is an indicator of cardiorespiratory fitness with clinical and public health implications. Submaximal exercise capacity and its response to exercise programs are characterized by heritability levels of about 40%. Using physical working capacity (power output) at a heart rate of 150 beats/min (PWC150) as an indicator of submaximal exercise capacity in subjects of the HERITAGE Family Study, we have undertaken multi-omics and in silico explorations of the underlying biology of PWC150 and its response to 20 wk of endurance training. Our goal was to illuminate the biological processes and identify panels of genes associated with human variability in intrinsic PWC150 (iPWC150) and its trainability (dPWC150). Our bioinformatics approach was based on a combination of genome-wide association, skeletal muscle gene expression, and plasma proteomics and metabolomics experiments. Genes, proteins, and metabolites showing significant associations with iPWC150 or dPWC150 were further queried for the enrichment of biological pathways. We compared genotype-phenotype associations of emerging candidate genes with reported functional consequences of gene knockouts in mouse models. We investigated the associations between DNA variants and multiple muscle and cardiovascular phenotypes measured in HERITAGE subjects. Two panels of prioritized genes of biological relevance to iPWC150 (13 genes) and dPWC150 (6 genes) were identified, supporting the hypothesis that genes and pathways associated with iPWC150 are different from those underlying dPWC150. Finally, the functions of these genes and pathways suggested that human variation in submaximal exercise capacity is mainly driven by skeletal muscle morphology and metabolism and red blood cell oxygen-carrying capacity.NEW & NOTEWORTHY Multi-omics and in silico explorations of the genes and underlying biology of submaximal exercise capacity and its response to 20 wk of endurance training were undertaken. Prioritized genes were identified: 13 genes for variation in submaximal exercise capacity in the sedentary state and 5 genes for the response level to endurance training, with no overlap between them. Genes and pathways associated with submaximal exercise capacity in the sedentary state are different from those underlying trainability.

Sex differences in the genetics of sarcoidosis across European and African ancestry populations

Background

Sex differences in the susceptibility of sarcoidosis are unknown. The study aims to identify sex-dependent genetic variations in two clinical sarcoidosis phenotypes: Löfgren's syndrome (LS) and non-Löfgren's syndrome (non-LS).

Methods

A meta-analysis of genome-wide association studies was conducted on Europeans and African Americans, totaling 10,103 individuals from three population-based cohorts, Sweden (n = 3,843), Germany (n = 3,342), and the United States (n = 2,918), followed by an SNP lookup in the UK Biobank (UKB, n = 387,945). A genome-wide association study based on Immunochip data consisting of 141,000 single nucleotide polymorphisms (SNPs) was conducted in the sex groups. The association test was based on logistic regression using the additive model in LS and non-LS sex groups independently. Additionally, gene-based analysis, gene expression, expression quantitative trait loci (eQTL) mapping, and pathway analysis were performed to discover functionally relevant mechanisms related to sarcoidosis and biological sex.

Results

We identified sex-dependent genetic variations in LS and non-LS sex groups. Genetic findings in LS sex groups were explicitly located in the extended Major Histocompatibility Complex (xMHC). In non-LS, genetic differences in the sex groups were primarily located in the MHC class II subregion and ANXA11. Gene-based analysis and eQTL enrichment revealed distinct sex-specific gene expression patterns in various tissues and immune cell types. In LS sex groups, a pathway map related to antigen presentation machinery by IFN-gamma. In non-LS, pathway maps related to immune response lectin-induced complement pathway in males and related to maturation and migration of dendritic cells in skin sensitization in females were identified.

Conclusion

Our findings provide new evidence for a sex bias underlying sarcoidosis genetic architecture, particularly in clinical phenotypes LS and non-LS. Biological sex likely plays a role in disease mechanisms in sarcoidosis.

Reduced urine volume and changed renal sphingolipid metabolism in P2ry14-deficient mice

The UDP-glucose receptor P2RY14, a rhodopsin-like G protein-coupled receptor (GPCR), was previously described as receptor expressed in A-intercalated cells of the mouse kidney. Additionally, we found P2RY14 is abundantly expressed in mouse renal collecting duct principal cells of the papilla and epithelial cells lining the renal papilla. To better understand its physiological function in kidney, we took advantage of a P2ry14 reporter and gene-deficient (KO) mouse strain. Morphometric studies showed that the receptor function contributes to kidney morphology. KO mice had a broader cortex relative to the total kidney area than wild-type (WT) mice. In contrast, the area of the outer stripe of the outer medulla was larger in WT compared to KO mice. Transcriptome comparison of the papilla region of WT and KO mice revealed differences in the gene expression of extracellular matrix proteins (e.g., decorin, fibulin-1, fibulin-7) and proteins involved in sphingolipid metabolism (e.g., small subunit b of the serine palmitoyltransferase) and other related GPCRs (e.g., GPR171). Using mass spectrometry, changes in the sphingolipid composition (e.g., chain length) were detected in the renal papilla of KO mice. At the functional level, we found that KO mice had a reduced urine volume but an unchanged glomerular filtration rate under normal chow and salt diets. Our study revealed P2ry14 as a functionally important GPCR in collecting duct principal cells and cells lining the renal papilla and the possible involvement of P2ry14 in nephroprotection by regulation of decorin.

Decoding the genetic and epigenetic basis of asthma

Asthma is a complex and heterogeneous chronic inflammatory disease of the airways. Alongside environmental factors, asthma susceptibility is strongly influenced by genetics. Given its high prevalence and our incomplete understanding of the mechanisms underlying disease susceptibility, asthma is frequently studied in genome-wide association studies (GWAS), which have identified thousands of genetic variants associated with asthma development. Virtually all these genetic variants reside in non-coding genomic regions, which has obscured the functional impact of asthma-associated variants and their translation into disease-relevant mechanisms. Recent advances in genomics technology and epigenetics now offer methods to link genetic variants to gene regulatory elements embedded within non-coding regions, which have started to unravel the molecular mechanisms underlying the complex (epi)genetics of asthma. Here, we provide an integrated overview of (epi)genetic variants associated with asthma, focusing on efforts to link these disease associations to biological insight into asthma pathophysiology using state-of-the-art genomics methodology. Finally, we provide a perspective as to how decoding the genetic and epigenetic basis of asthma has the potential to transform clinical management of asthma and to predict the risk of asthma development.

eQTL Set-Based Association Analysis Identifies Novel Susceptibility Loci for Barrett Esophagus and Esophageal Adenocarcinoma

BACKGROUND

Over 20 susceptibility single-nucleotide polymorphisms (SNP) have been identified for esophageal adenocarcinoma (EAC) and its precursor, Barrett esophagus (BE), explaining a small portion of heritability.

METHODS

Using genetic data from 4,323 BE and 4,116 EAC patients aggregated by international consortia including the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON), we conducted a comprehensive transcriptome-wide association study (TWAS) for BE/EAC, leveraging Genotype Tissue Expression (GTEx) gene-expression data from six tissue types of plausible relevance to EAC etiology: mucosa and muscularis from the esophagus, gastroesophageal (GE) junction, stomach, whole blood, and visceral adipose. Two analytical approaches were taken: standard TWAS using the predicted gene expression from local expression quantitative trait loci (eQTL), and set-based SKAT association using selected eQTLs that predict the gene expression.

RESULTS

Although the standard approach did not identify significant signals, the eQTL set-based approach identified eight novel associations, three of which were validated in independent external data (eQTL SNP sets for EXOC3, ZNF641, and HSP90AA1).

CONCLUSIONS

This study identified novel genetic susceptibility loci for EAC and BE using an eQTL set-based genetic association approach.

IMPACT

This study expanded the pool of genetic susceptibility loci for EAC and BE, suggesting the potential of the eQTL set-based genetic association approach as an alternative method for TWAS analysis.

GPR105-Targeted Therapy Promotes Gout Resolution as a Switch Between NETosis and Apoptosis of Neutrophils

The fate of infiltrating neutrophils in inflamed joints determines the development of acute gouty arthritis (AGA). GPR105 highly expressed in human neutrophils is sensitive to monosodium urate crystals (MSU); nevertheless, the roles of GPR105 in AGA remain unclear. Here, we show that GPR105 is significantly upregulated in peripheral polymorphonuclear neutrophils of AGA patients. GPR105 knockout (GPR105^−/−) prevented NETosis and induced apoptosis of neutrophils under MSU exposure, as well as attenuating inflammatory cascades in AGA. Mechanistically, GPR105 deletion activated cAMP-PKA signals, thereby disrupting Raf-Mek1/2-Erk1/2 pathway-mediated NADPH oxidase activation, contributing to inhibition of NETosis. Whereas, cAMP-PKA activation resulting in GPR105 deficiency modulated PI3K-Akt pathway to regulate apoptosis. More importantly, suppression of cAMP-PKA pathway by SQ22536 and H-89 restored NETosis instead of apoptosis in GPR105^−/− neutrophils, promoting MSU-induced gout flares. Interestingly, lobetyolin was screened out as a potent GPR105 antagonist using molecular docking-based virtual screening and in vitro activity test, which efficiently attenuated MSU-induced inflammatory response interacting with GPR105. Taken together, our study implicated that modulating cell death patterns between NETosis and apoptosis through targeting GPR105 could be a potential therapeutic strategy for the treatment of AGA.

Multiomics Analysis Identifies BIRC3 as a Novel Glucocorticoid Response-Associated Gene

BACKGROUND

Inhaled corticosteroid (ICS) response among patients with asthma is influenced by genetics, but biologically actionable insights based on associations have not been found. Various glucocorticoid response omics datasets are available to interrogate their biological effects.

OBJECTIVE

We sought to identify functionally relevant ICS response genetic associations by integrating complementary multiomics datasets.

METHODS

Variants with p-values<10^-4 from a previous ICS response genome-wide association study were re-ranked based on integrative scores determined from: i) glucocorticoid receptor (GR)- and ii) RNA polymerase II (RNAP II)-binding regions inferred from ChIP-Seq data for three airway cell types, iii) glucocorticoid response element (GRE) motifs, iv) differentially expressed genes in response to glucocorticoid exposure according to 20 transcriptomic datasets, and v) expression quantitative trait loci (eQTLs) from GTEx. Candidate variants were tested for association with ICS response and asthma in six independent studies.

RESULTS

Four variants had significant (q-value<0.05) multiomics integrative scores. These variants were in a locus consisting of 52 variants in high LD (r²≥0.8) near GR-binding sites by the gene BIRC3. Variants were also BIRC3 eQTLs in lung, and two were within/near putative GRE motifs. BIRC3 had increased RNAP II occupancy and gene expression with glucocorticoid exposure in two ChIP-Seq and 13 transcriptomic datasets. Some BIRC3 variants in the 52-variant locus were associated (p-value<0.05) with ICS response in three independent studies and others with asthma in one study.

CONCLUSION

BIRC3 should be prioritized for further functional studies of ICS response.

CLINICAL IMPLICATION

Genetic variation near BIRC3 may influence ICS response in people with asthma.

Targeting the P2Y13 Receptor Suppresses IL-33 and HMGB1 Release and Ameliorates Experimental Asthma

RATIONALE

The alarmins IL-33 and HMGB1 (high mobility group box 1) contribute to type-2 inflammation and asthma pathogenesis.

OBJECTIVES

To determine whether P2Y13 receptor (P2Y13-R), a purinergic G protein-coupled receptor (GPCR) and risk allele for asthma, regulates the release of IL-33 and HMGB1.

METHODS

Bronchial biopsies were obtained from healthy and asthmatic subjects. Primary human airway epithelial cells (AECs), primary mouse (m)AECs, or C57Bl/6 mice were inoculated with various aeroallergens or respiratory viruses, and the nuclear-to-cytoplasmic translocation and release of alarmins measured by immunohistochemistry and ELISA. The role of P2Y13-R in AEC function and in the onset, progression, and an exacerbation of experimental asthma, was assessed using pharmacological antagonists and P2Y13-R gene-deleted mice.

MEASUREMENTS AND MAIN RESULTS

Aeroallergen-exposure induced the extracellular release of ADP and ATP, nucleotides that activate P2Y13-R. ATP, ADP, aeroallergen (house dust mite, cockroach or Alternaria) or virus exposure induced the nuclear-to-cytoplasmic translocation and subsequent release of IL-33 and HMGB1, and this response was ablated by genetic deletion or pharmacological antagonism of P2Y13. In mice, prophylactic or therapeutic P2Y13-R blockade attenuated asthma onset, and critically, ablated the severity of a rhinovirus-associated exacerbation in a high-fidelity experimental model of chronic asthma. Moreover, P2Y13-R antagonism derepressed antiviral immunity, increasing IFN-λ production and decreasing viral copies in the lung.

CONCLUSIONS

We identify P2Y13-R as a novel gatekeeper of the nuclear alarmins IL-33 and HMGB1, and demonstrate that the targeting of this GPCR via genetic deletion or treatment with a small-molecule antagonist protects against the onset and exacerbations of experimental asthma.

Leveraging auxiliary data from arbitrary distributions to boost GWAS discovery with Flexible cFDR

Genome-wide association studies (GWAS) have identified thousands of genetic variants that are associated with complex traits. However, a stringent significance threshold is required to identify robust genetic associations. Leveraging relevant auxiliary covariates has the potential to boost statistical power to exceed the significance threshold. Particularly, abundant pleiotropy and the non-random distribution of SNPs across various functional categories suggests that leveraging GWAS test statistics from related traits and/or functional genomic data may boost GWAS discovery. While type 1 error rate control has become standard in GWAS, control of the false discovery rate can be a more powerful approach. The conditional false discovery rate (cFDR) extends the standard FDR framework by conditioning on auxiliary data to call significant associations, but current implementations are restricted to auxiliary data satisfying specific parametric distributions, typically GWAS p-values for related traits. We relax these distributional assumptions, enabling an extension of the cFDR framework that supports auxiliary covariates from arbitrary continuous distributions ("Flexible cFDR"). Our method can be applied iteratively, thereby supporting multi-dimensional covariate data. Through simulations we show that Flexible cFDR increases sensitivity whilst controlling FDR after one or several iterations. We further demonstrate its practical potential through application to an asthma GWAS, leveraging various functional genomic data to find additional genetic associations for asthma, which we validate in the larger, independent, UK Biobank data resource.

Novel pharmacological therapies for the treatment of bronchial asthma

Asthma has long been recognised as a chronic inflammatory disease of the airways, often in response to inhaled allergens prompting inappropriate activation of the immune response. involving a range of cells including mast cells, Th2 lymphocytes and eosinophils and a wide range of inflammatory mediators. First-line therapy for treatment of persistent asthma involves the use of inhaled corticosteroids (ICS) in combination with inhaled β2-agonists enabling both the control of the underlying airways inflammation and a reduction of airway hyperresponsiveness. However, many patients remain symptomatic despite high-dose therapy. There is therefore a continued unmet clinical need to develop specifically new anti-inflammatory therapies for patients with asthma, either as an add-on therapy to ICS or as replacement monotherapies. The success of fixed dose combination inhalers containing both a bronchodilator and an anti-inflammatory drug has also led to the development of "bifunctional" drugs which are molecules specifically designed to have two distinct pharmacological actions based on distinct pharmacophores. In this review we will discuss these different pharmacological approaches under development for the treatment of bronchial asthma and the available pre-clinical and clinical data.

Purinergic Receptors in the Airways: Potential Therapeutic Targets for Asthma?

Extracellular ATP functions as a signaling messenger through its actions on purinergic receptors, and is known to be involved in numerous physiological and pathophysiological processes throughout the body, including in the lungs and airways. Consequently, purinergic receptors are considered to be promising therapeutic targets for many respiratory diseases, including asthma. This review explores how online bioinformatics resources combined with recently generated datasets can be utilized to investigate purinergic receptor gene expression in tissues and cell types of interest in respiratory disease to identify potential therapeutic targets, which can then be investigated further. These approaches show that different purinergic receptors are expressed at different levels in lung tissue, and that purinergic receptors tend to be expressed at higher levels in immune cells and at more moderate levels in airway structural cells. Notably, P2RX1, P2RX4, P2RX7, P2RY1, P2RY11, and P2RY14 were revealed as the most highly expressed purinergic receptors in lung tissue, therefore suggesting that these receptors have good potential as therapeutic targets for asthma and other respiratory diseases.

The central role of IL-33/IL-1RL1 pathway in asthma: From pathogenesis to intervention

Interleukin-33 (IL-33), a member of the IL-1 family, and its cognate receptor, Interleukin-1 receptor like-1 (IL-1RL1 or ST2), are susceptibility genes for childhood asthma. In response to cellular damage, IL-33 is released from barrier tissues as an 'alarmin' to activate the innate immune response. IL-33 drives type 2 responses by inducing signalling through its receptor IL-1RL1 in several immune and structural cells, thereby leading to type 2 cytokine and chemokine production. IL-1RL1 gene transcript encodes different isoforms generated through alternative splicing. Its soluble isoform, IL-1RL1-a or sST2, acts as a decoy receptor by sequestering IL-33, thereby inhibiting IL1RL1-b/IL-33 signalling. IL-33 and its receptor IL-1RL1 are therefore considered as putative biomarkers or targets for pharmacological intervention in asthma. This review will provide an overview of the genetics and biology of the IL-33/IL-1RL1 pathway in the context of asthma pathogenesis. It will discuss the potential and complexities of targeting the cytokine or its receptor, how genetics or biomarkers may inform precision medicine for asthma targeting this pathway, and the possible positioning of therapeutics targeting IL-33 or its receptor in the expanding landscape of novel biologicals applied in asthma management.

UDP-glucose and P2Y14 receptor amplify allergen-induced airway eosinophilia

Airway eosinophilia is a hallmark of allergic asthma and is associated with mucus production, airway hyperresponsiveness, and shortness of breath. Although glucocorticoids are widely used to treat asthma, their prolonged use is associated with several side effects. Furthermore, many individuals with eosinophilic asthma are resistant to glucocorticoid treatment, and they have an unmet need for novel therapies. Here, we show that UDP-glucose (UDP-G), a nucleotide sugar, is selectively released into the airways of allergen-sensitized mice upon their subsequent challenge with that same allergen. Mice lacking P2Y14R, the receptor for UDP-G, had decreased airway eosinophilia and airway hyperresponsiveness compared with wild-type mice in a protease-mediated model of asthma. P2Y14R was dispensable for allergic sensitization and for the production of type 2 cytokines in the lung after challenge. However, UDP-G increased chemokinesis in eosinophils and enhanced their response to the eosinophil chemoattractant, CCL24. In turn, eosinophils triggered the release of UDP-G into the airway, thereby amplifying eosinophilic recruitment. This positive feedback loop was sensitive to therapeutic intervention, as a small molecule antagonist of P2Y14R inhibited airway eosinophilia. These findings thus reveal a pathway that can be therapeutically targeted to treat asthma exacerbations and glucocorticoid-resistant forms of this disease.

Uridine diphosphate-glucose/P2Y14R axis is a nonchemokine pathway that selectively promotes eosinophil accumulation

Allergic asthma is a chronic inflammatory airway disease characterized by dysregulated type 2 immune responses, including degranulating airway eosinophils that induce tissue damage and airway hyperresponsiveness (AHR). The type 2 cytokines interleukin 5 (IL-5) and IL-13 and the eosinophil-specific chemokine CCL11/CCL24/CCL26 axis recruit, activate, and regulate eosinophils in the airways. In this issue of the JCI, Karcz et al. identified a mechanism involving the nucleotide sugar UDP-glucose (UDP-G) and the purinergic receptor P2Y14R in amplifying eosinophil accumulation in the lung. During type 2 inflammation, UDP-G activates P2Y14R on eosinophils, inducing the cells to move and migrate into the lung. Pharmacologically or genetically inhibiting P2Y14R on eosinophils attenuated eosinophil infiltration and AHR. Future experiments, including identifying additional type 2 factors regulating P2Y14R expression on lung eosinophils, are necessary to ascertain the impact of targeting P2Y14R as an alternative or adjunctive therapy to current type 2 biologics for the treatment of asthma.

The influence of human genetic variation on Epstein-Barr virus sequence diversity

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

E-MAGMA: an eQTL-informed method to identify risk genes using genome-wide association study summary statistics

MOTIVATION

Genome-wide association studies have successfully identified multiple independent genetic loci that harbour variants associated with human traits and diseases, but the exact causal genes are largely unknown. Common genetic risk variants are enriched in non-protein-coding regions of the genome and often affect gene expression (expression quantitative trait loci, eQTL) in a tissue-specific manner. To address this challenge, we developed a methodological framework, E-MAGMA, which converts genome-wide association summary statistics into gene-level statistics by assigning risk variants to their putative genes based on tissue-specific eQTL information.

RESULTS

We compared E-MAGMA to three eQTL informed gene-based approaches using simulated phenotype data. Phenotypes were simulated based on eQTL reference data using GCTA for all genes with at least one eQTL at chromosome 1. We performed 10 simulations per gene. The eQTL-h2 (i.e., the proportion of variation explained by the eQTLs) was set at 1%, 2%, and 5%. We found E-MAGMA outperforms other gene-based approaches across a range of simulated parameters (e.g. the number of identified causal genes). When applied to genome-wide association summary statistics for five neuropsychiatric disorders, E-MAGMA identified more putative candidate causal genes compared to other eQTL-based approaches. By integrating tissue-specific eQTL information, these results show E-MAGMA will help to identify novel candidate causal genes from genome-wide association summary statistics and thereby improve the understanding of the biological basis of complex disorders.

AVAILABILITY

A tutorial and input files are made available in a github repository: https://github.com/eskederks/eMAGMA-tutorial.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

P2Y14 Receptor as a Target for Neutrophilia Attenuation in Severe COVID-19 Cases: From Hematopoietic Stem Cell Recruitment and Chemotaxis to Thrombo-inflammation

The global SARS-CoV-2 pandemic starting in 2019 has already reached more than 2.3 million deaths. Despite the scientific community’s efforts to investigate the COVID-19 disease, a drug for effectively treating or curing patients yet needs to be discovered. Hematopoietic stem cells (HSC) differentiating into immune cells for defense express COVID-19 entry receptors, and COVID-19 infection hinders their differentiation. The importance of purinergic signaling in HSC differentiation and innate immunity has been recognized. The metabotropic P2Y14 receptor subtype, activated by UDP-glucose, controls HSC differentiation and mobilization. Thereon, the exacerbated activation of blood immune cells amplifies the inflammatory state observed in COVID-19 patients, specially through the continuous release of reactive oxygen species and extracellular neutrophil traps (NETs). Further, the P2Y14 subtype, robustly inhibits the infiltration of neutrophils into various epithelial tissues, including lungs and kidneys. Here we discuss findings suggesting that antagonism of the P2Y14 receptor could prevent the progression of COVID-19-induced systemic inflammation, which often leads to severe illness and death cases. Considering the modulation of neutrophil recruitment of extreme relevance for respiratory distress and lung failure prevention, we propose that P2Y14 receptor inhibition by its selective antagonist PPTN could limit neutrophil recruitment and NETosis, hence limiting excessive formation of oxygen reactive species and proteolytic activation of the kallikrein-kinin system and subsequent bradykinin storm in the alveolar septa of COVID-19 patients.

Biochemical transformation of bacterial lipopolysaccharides by acyloxyacyl hydrolase reduces host injury and promotes recovery

Animals can sense the presence of microbes in their tissues and mobilize their own defenses by recognizing and responding to conserved microbial structures (often called microbe-associated molecular patterns (MAMPs)). Successful host defenses may kill the invaders, yet the host animal may fail to restore homeostasis if the stimulatory microbial structures are not silenced. Although mice have many mechanisms for limiting their responses to lipopolysaccharide (LPS), a major Gram-negative bacterial MAMP, a highly conserved host lipase is required to extinguish LPS sensing in tissues and restore homeostasis. We review recent progress in understanding how this enzyme, acyloxyacyl hydrolase (AOAH), transforms LPS from stimulus to inhibitor, reduces tissue injury and death from infection, prevents prolonged post-infection immunosuppression, and keeps stimulatory LPS from entering the bloodstream. We also discuss how AOAH may increase sensitivity to pulmonary allergens. Better appreciation of how host enzymes modify LPS and other MAMPs may help prevent tissue injury and hasten recovery from infection.

An Automated Functional Annotation Pipeline That Rapidly Prioritizes Clinically Relevant Genes for Autism Spectrum Disorder

Human genetic studies have implicated more than a hundred genes in Autism Spectrum Disorder (ASD). Understanding how variation in implicated genes influence expression of co-occurring conditions and drug response can inform more effective, personalized approaches for treatment of individuals with ASD. Rapidly translating this information into the clinic requires efficient algorithms to sort through the myriad of genes implicated by rare gene-damaging single nucleotide and copy number variants, and common variation detected in genome-wide association studies (GWAS). To pinpoint genes that are more likely to have clinically relevant variants, we developed a functional annotation pipeline. We defined clinical relevance in this project as any ASD associated gene with evidence indicating a patient may have a complex, co-occurring condition that requires direct intervention (e.g., sleep and gastrointestinal disturbances, attention deficit hyperactivity, anxiety, seizures, depression), or is relevant to drug development and/or approaches to maximizing efficacy and minimizing adverse events (i.e., pharmacogenomics). Starting with a list of all candidate genes implicated in all manifestations of ASD (i.e., idiopathic and syndromic), this pipeline uses databases that represent multiple lines of evidence to identify genes: (1) expressed in the human brain, (2) involved in ASD-relevant biological processes and resulting in analogous phenotypes in mice, (3) whose products are targeted by approved pharmaceutical compounds or possessing pharmacogenetic variation and (4) whose products directly interact with those of genes with variants recommended to be tested for by the American College of Medical Genetics (ACMG). Compared with 1000 gene sets, each with a random selection of human protein coding genes, more genes in the ASD set were annotated for each category evaluated (p ≤ 1.99 × 10⁻²). Of the 956 ASD-implicated genes in the full set, 18 were flagged based on evidence in all categories. Fewer genes from randomly drawn sets were annotated in all categories (x = 8.02, sd = 2.56, p = 7.75 × 10⁻⁴). Notably, none of the prioritized genes are represented among the 59 genes compiled by the ACMG, and 78% had a pathogenic or likely pathogenic variant in ClinVar. Results from this work should rapidly prioritize potentially actionable results from genetic studies and, in turn, inform future work toward clinical decision support for personalized care based on genetic testing.

Purinergic Signaling: Impact of GPCR Structures on Rational Drug Design

The purinergic signaling system includes membrane-bound receptors for extracellular purines and pyrimidines, and enzymes/transporters that regulate receptor activation by endogenous agonists. Receptors include: adenosine (A1 , A2A , A2B, and A3 ) and P2Y (P2Y1 , P2Y2 , P2Y4 , P2Y6 , P2Y11 , P2Y12 , P2Y13 , and P2Y14 ) receptors (all GPCRs), as well as P2X receptors (ion channels). Receptor activation, especially accompanying physiological stress or damage, creates a temporal sequence of signaling to counteract this stress and either mobilize (P2Rs) or suppress (ARs) immune responses. Thus, modulation of this large signaling family has broad potential for treating chronic diseases. Experimentally determined structures represent each of the three receptor families. We focus on selective purinergic agonists (A1 , A3 ), antagonists (A3 , P2Y14 ), and allosteric modulators (P2Y1 , A3 ). Examples of applying structure-based design, including the rational modification of known ligands, are presented for antithrombotic P2Y1 R antagonists and anti-inflammatory P2Y14 R antagonists and A3 AR agonists. A3 AR agonists are a potential, nonaddictive treatment for chronic neuropathic pain.

Exploration of Alternative Scaffolds for P2Y14 Receptor Antagonists Containing a Biaryl Core

Various heteroaryl and bicyclo-aliphatic analogues of zwitterionic biaryl P2Y₁₄ receptor (P2Y₁₄R) antagonists were synthesized, and affinity was measured in P2Y₁₄R-expressing Chinese hamster ovary cells by flow cytometry. Given this series' low water solubility, various polyethylene glycol derivatives of the distally binding piperidin-4-yl moiety of moderate affinity were synthesized. Rotation of previously identified 1,2,3-triazole attached to the central m-benzoic acid core (25) provided moderate affinity but not indole and benzimidazole substitution of the aryl-triazole. The corresponding P2Y₁₄R region is predicted by homology modeling as a deep, sterically limited hydrophobic pocket, with the outward pointing piperidine moiety being the most flexible. Bicyclic-substituted piperidine ring derivatives of naphthalene antagonist 1, e.g., quinuclidine 17 (MRS4608, IC₅₀ ≈ 20 nM at hP2Y₁₄R/mP2Y₁₄R), or of triazole 2, preserved affinity. Potent antagonists 1, 7a, 17, and 23 (10 mg/kg) protected in an ovalbumin/Aspergillus mouse asthma model, and PEG conjugate 12 reduced chronic pain. Thus, we expanded P2Y₁₄R antagonist structure-activity relationship, introducing diverse physical-chemical properties.

On the cross-population generalizability of gene expression prediction models

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

Differential DNA methylation profiles of peripheral blood mononuclear cells in allergic asthmatic children following dust mite immunotherapy

BACKGROUND/PURPOSE

Allergen-specific immunotherapy (SIT) is now considered curative to allergic diseases such as asthma. Mechanistically, our previous work showed DNA hypermethylation of cytokine genes, in T-helper cells, in allergic asthmatic children treated with allergen-SIT. In this study, we extended to work to assess possible changes in the DNA methylomes of peripheral blood mononuclear cells (PBMCs), isolated from mite allergen-SIT asthmatic children, to explore further the underlying methylation changes.

METHODS

Thirteen allergic asthmatic children who received Der p-SIT, 12 non-SIT allergic asthmatic controls, and 12 healthy controls were enrolled. Bisulfite-converted DNA from Der p-stimulated PBMCs was analyzed using Human Methylation 450 k BeadChip. Pyrosequencing and quantitative real-time PCR were used to validate the DNA methylation levels and the gene expression of individual samples.

RESULTS

We identified 108 significantly differentially methylated regions (DMRs) unique to Der p-treated PBMCs, with 53 probes linked to demethylated DMRs, and 55 probes linked to methylated DMRs. Three associated genes (BCL6, HSPG2, and HSP90AA1), of selected DMRs, were subjected to bisulfite pyrosequencing. Of these, BCL6 showed significant hypomethylation, while HSPG2 and HSP90AA1 were hypermethylated in SIT group, compared to the AA group. Furthermore, SIT group had significantly higher gene expression of BCL6 and lower gene expression of HSPG2. KEGG pathway analysis further revealed DMR genes involved in ECM-receptor interactions, asthma, and antigen processing and presentation pathways.

CONCLUSIONS

Several DNA regions showed DNA methylation altered by Der p specific immunotherapy, indicating desensitization-associated methylomes. Genes belonging to these SIT-altered pathways may represent therapeutic targets for better clinical management of asthma.

Proteolytic, lipidergic and polysaccharide molecular recognition shape innate responses to house dust mite allergens

House dust mites (HDMs) are sources of an extensive repertoire of allergens responsible for a range of allergic conditions. Technological advances have accelerated the identification of these allergens and characterized their putative roles within HDMs. Understanding their functional bioactivities is illuminating how they interact with the immune system to cause disease and how interrelations between them are essential to maximize allergic responses. Two types of allergen bioactivity, namely proteolysis and peptidolipid/lipid binding, elicit IgE and stimulate bystander responses to unrelated allergens. Much of this influence arises from Toll-like receptor (TLR) 4 or TLR2 signalling and, in the case of protease allergens, the activation of additional pleiotropic effectors with strong disease linkage. Of related interest is the interaction of HDM allergens with common components of the house dust matrix, through either their binding to allergens or their autonomous modulation of immune receptors. Herein, we provide a contemporary view of how proteolysis, lipid-binding activity and interactions with polysaccharides and polysaccharide molecular recognition systems coordinate the principal responses which underlie allergy. The power of the catalytically competent group 1 HDM protease allergen component is demonstrated by a review of disclosures surrounding the efficacy of novel inhibitors produced by structure-based design.

ZRANB3 is an African-specific type 2 diabetes locus associated with beta-cell mass and insulin response

Genome analysis of diverse human populations has contributed to the identification of novel genomic loci for diseases of major clinical and public health impact. Here, we report a genome-wide analysis of type 2 diabetes (T2D) in sub-Saharan Africans, an understudied ancestral group. We analyze ~18 million autosomal SNPs in 5,231 individuals from Nigeria, Ghana and Kenya. We identify a previously-unreported genome-wide significant locus: ZRANB3 (Zinc Finger RANBP2-Type Containing 3, lead SNP p = 2.831 × 10-9). Knockdown or genomic knockout of the zebrafish ortholog results in reduction in pancreatic β-cell number which we demonstrate to be due to increased apoptosis in islets. siRNA transfection of murine Zranb3 in MIN6 β-cells results in impaired insulin secretion in response to high glucose, implicating Zranb3 in β-cell functional response to high glucose conditions. We also show transferability in our study of 32 established T2D loci. Our findings advance understanding of the genetics of T2D in non-European ancestry populations.

The Purinergic System as a Pharmacological Target for the Treatment of Immune-Mediated Inflammatory Diseases

Immune-mediated inflammatory diseases (IMIDs) encompass a wide range of seemingly unrelated conditions, such as multiple sclerosis, rheumatoid arthritis, psoriasis, inflammatory bowel diseases, asthma, chronic obstructive pulmonary disease, and systemic lupus erythematosus. Despite differing etiologies, these diseases share common inflammatory pathways, which lead to damage in primary target organs and frequently to a plethora of systemic effects as well. The purinergic signaling complex comprising extracellular nucleotides and nucleosides and their receptors, the P2 and P1 purinergic receptors, respectively, as well as catabolic enzymes and nucleoside transporters is a major regulatory system in the body. The purinergic signaling complex can regulate the development and course of IMIDs. Here we provide a comprehensive review on the role of purinergic signaling in controlling immunity, inflammation, and organ function in IMIDs. In addition, we discuss the possible therapeutic applications of drugs acting on purinergic pathways, which have been entering clinical development, to manage patients suffering from IMIDs.

Identification of Novel Susceptibility Loci and Genes for Prostate Cancer Risk: A Transcriptome-Wide Association Study in Over 140,000 European Descendants

Genome-wide association study-identified prostate cancer risk variants explain only a relatively small fraction of its familial relative risk, and the genes responsible for many of these identified associations remain unknown. To discover novel prostate cancer genetic loci and possible causal genes at previously identified risk loci, we performed a transcriptome-wide association study in 79,194 cases and 61,112 controls of European ancestry. Using data from the Genotype-Tissue Expression Project, we established genetic models to predict gene expression across the transcriptome for both prostate models and cross-tissue models and evaluated model performance using two independent datasets. We identified significant associations for 137 genes at P < 2.61 × 10-6, a Bonferroni-corrected threshold, including nine genes that remained significant at P < 2.61 × 10-6 after adjusting for all known prostate cancer risk variants in nearby regions. Of the 128 remaining associated genes, 94 have not yet been reported as potential target genes at known loci. We silenced 14 genes and many showed a consistent effect on viability and colony-forming efficiency in three cell lines. Our study provides substantial new information to advance our understanding of prostate cancer genetics and biology. SIGNIFICANCE: This study identifies novel prostate cancer genetic loci and possible causal genes, advancing our understanding of the molecular mechanisms that drive prostate cancer.

Recent progress on the discovery of P2Y14 receptor antagonists

The P2Y₁₄ receptor (P2Y₁₄R), a G protein-coupled receptor (GPCR), is activated by extracellular nucleotides. P2Y₁₄R is involved in inflammatory, diabetes, immune processes and other related complications, and is therefore an attractive therapeutic target. As the three-dimensional structure of the P2Y₁₄R has not yet been elucidated, homology modeling based on the crystallography of the closely related P2Y₁₂R have been used in the structure-based design of P2Y₁₄R ligands. Several P2Y₁₄R antagonists with excellent potency and high subtype-selectivity have been discovered in recent years. In this review, development of novel small molecules as antagonists of P2Y₁₄R was described.

Pharmacological strategies for targeting platelet activation in asthma

The activation of platelets during host defence and inflammatory disorders has become increasingly documented. Clinical studies of patients with asthma reveal heightened platelet activation and accumulation into lung tissue. Accompanying studies in animal models of allergic lung inflammation, using protocols of experimentally induced thrombocytopenia proclaim an important role for platelets during the leukocyte recruitment cascade, tissue integrity, and lung function. The functions of platelets during these inflammatory events are clearly distinct to platelet functions during haemostasis and clot formation, and have led to the concept that a dichotomy (or polytomy, depending on what else platelets do) in platelet activation exists. The platelet, therefore, presents us with novel opportunities for modulating these inflammatory responses. This review discusses the rationale and effectiveness of current anti-platelet drugs in their use to supress inflammation with regard to asthma, and the need to consider novel possibilities for pharmacological modulation of platelet function associated with inflammation that are pharmacologically distinct to current anti-platelet therapies.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Association study in African-admixed populations across the Americas recapitulates asthma risk loci in non-African populations

Asthma is a complex disease with striking disparities across racial and ethnic groups. Despite its relatively high burden, representation of individuals of African ancestry in asthma genome-wide association studies (GWAS) has been inadequate, and true associations in these underrepresented minority groups have been inconclusive. We report the results of a genome-wide meta-analysis from the Consortium on Asthma among African Ancestry Populations (CAAPA; 7009 asthma cases, 7645 controls). We find strong evidence for association at four previously reported asthma loci whose discovery was driven largely by non-African populations, including the chromosome 17q12-q21 locus and the chr12q13 region, a novel (and not previously replicated) asthma locus recently identified by the Trans-National Asthma Genetic Consortium (TAGC). An additional seven loci reported by TAGC show marginal evidence for association in CAAPA. We also identify two novel loci (8p23 and 8q24) that may be specific to asthma risk in African ancestry populations.

Eleven loci with new reproducible genetic associations with allergic disease risk

BACKGROUND

A recent genome-wide association study (GWAS) identified 99 loci that contain genetic risk variants shared between asthma, hay fever, and eczema. Many more risk loci shared between these common allergic diseases remain to be discovered, which could point to new therapeutic opportunities.

OBJECTIVE

We sought to identify novel risk loci shared between asthma, hay fever, and eczema by applying a gene-based test of association to results from a published GWAS that included data from 360,838 subjects.

METHODS

We used approximate conditional analysis to adjust the results from the published GWAS for the effects of the top risk variants identified in that study. We then analyzed the adjusted GWAS results with the EUGENE gene-based approach, which combines evidence for association with disease risk across regulatory variants identified in different tissues. Novel gene-based associations were followed up in an independent sample of 233,898 subjects from the UK Biobank study.

RESULTS

Of the 19,432 genes tested, 30 had a significant gene-based association at a Bonferroni-corrected P value of 2.5 × 10-6. Of these, 20 were also significantly associated (P < .05/30 = .0016) with disease risk in the replication sample, including 19 that were located in 11 loci not reported to contain allergy risk variants in previous GWASs. Among these were 9 genes with a known function that is directly relevant to allergic disease: FOSL2, VPRBP, IPCEF1, PRR5L, NCF4, APOBR, IL27, ATXN2L, and LAT. For 4 genes (eg, ATXN2L), a genetically determined decrease in gene expression was associated with decreased allergy risk, and therefore drugs that inhibit gene expression or function are predicted to ameliorate disease symptoms. The opposite directional effect was observed for 14 genes, including IL27, a cytokine known to suppress TH2 responses.

CONCLUSION

Using a gene-based approach, we identified 11 risk loci for allergic disease that were not reported in previous GWASs. Functional studies that investigate the contribution of the 19 associated genes to the pathophysiology of allergic disease and assess their therapeutic potential are warranted.

Modulation of blood inflammatory markers by benralizumab in patients with eosinophilic airway diseases

Background

Benralizumab, a humanized, afucosylated, monoclonal antibody that targets interleukin-5 receptor α, depletes eosinophils and basophils by enhanced antibody-dependent cell-mediated cytotoxicity. It demonstrated efficacy for patients with moderate to severe asthma and, in a Phase IIa trial, for chronic obstructive pulmonary disease (COPD) with eosinophilic inflammation. We investigated effects of benralizumab 100 mg every 8 weeks (first three doses every 4 weeks) subcutaneous on blood inflammatory markers through proteomic and gene-expression analyses collected during two Phase II studies of patients with eosinophilic asthma and eosinophilic COPD.

Methods

Serum samples for proteomic analysis and whole blood for gene expression analysis were collected at baseline and 52 weeks (asthma study) or 32 weeks (COPD study) post-treatment. Proteomic analyses were conducted on a custom set of 90 and 147 Rules-Based Medicine analytes for asthma and COPD, respectively. Gene expression was profiled by Affymetrix Human Genome U133 plus 2 arrays (~ 54 K probes). Gene set variation analysis (GSVA) was used to determine transcriptomic activity of immune signatures. Treatment-related differences between analytes, genes, and gene signatures were analyzed for the overall population and for patient subgroups stratified by baseline blood eosinophil count (eosinophil-high [≥300 cells/μL] and eosinophil-low [< 300 cells/μL]) via t-test and repeated measures analysis of variance.

Results

Eosinophil chemokines eotaxin-1 and eotaxin-2 were significantly upregulated (false discovery rate [FDR] < 0.05) by approximately 2.1- and 1.4-fold in the asthma study and by 2.3- and 1.7-fold in the COPD study following benralizumab treatment. Magnitude of upregulation of these two chemokines was greater for eosinophil-high patients than eosinophil-low patients in both studies. Benralizumab was associated with significant reductions (FDR < 0.05) in expression of genes associated with eosinophils and basophils, such as CLC, IL-5Rα, and PRSS33; immune-signaling complex genes (FCER1A); G-protein–coupled receptor genes (HRH4, ADORA3, P2RY14); and further immune-related genes (ALOX15 and OLIG2). The magnitude of downregulation of gene expression was greater for eosinophil-high than eosinophil-low patients. GSVA on immune signatures indicated significant treatment reductions (FDR < 0.05) in eosinophil-associated signatures.

Conclusions

Benralizumab is highly selective, modulating blood proteins or genes associated with eosinophils or basophils. Modulated protein and gene expression patterns are most prominently altered in eosinophil-high vs. eosinophil-low patients.

Trial registration

NCT01227278 and NCT01238861.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0968-8) contains supplementary material, which is available to authorized users.

Lessons Learned From GWAS of Asthma

Asthma is a common complex disease of the airways. Genome-wide association studies (GWASs) of asthma have identified many risk alleles and loci that have been replicated in worldwide populations. Although the risk alleles identified by GWAS have small effects and explain only a small portion of prevalence, the discovery of asthma loci can provide an understanding of its genetic architecture and the molecular pathways involved in disease pathogenesis. These discoveries can translate into advances in clinical care by identifying therapeutic targets, preventive strategies and ultimately approaches for personalized medicine. In this review, we summarize results from GWAS of asthma from the past 10 years and the insights gleaned from these discoveries.

High-Throughput Sequencing in Respiratory, Critical Care, and Sleep Medicine Research. An Official American Thoracic Society Workshop Report

High-throughput, "next-generation" sequencing methods are now being broadly applied across all fields of biomedical research, including respiratory disease, critical care, and sleep medicine. Although there are numerous review articles and best practice guidelines related to sequencing methods and data analysis, there are fewer resources summarizing issues related to study design and interpretation, especially as applied to common, complex, nonmalignant diseases. To address these gaps, a single-day workshop was held at the American Thoracic Society meeting in May 2017, led by the American Thoracic Society Section on Genetics and Genomics. The aim of this workshop was to review the design, analysis, interpretation, and functional follow-up of high-throughput sequencing studies in respiratory, critical care, and sleep medicine research. This workshop brought together experts in multiple fields, including genetic epidemiology, biobanking, bioinformatics, and research ethics, along with physician-scientists with expertise in a range of relevant diseases. The workshop focused on application of DNA and RNA sequencing research in common chronic diseases and did not cover sequencing studies in lung cancer, monogenic diseases (e.g., cystic fibrosis), or microbiome sequencing. Participants reviewed and discussed study design, data analysis and presentation, interpretation, functional follow-up, and reporting of results. This report summarizes the main conclusions of the workshop, specifically addressing the application of these methods in respiratory, critical care, and sleep medicine research. This workshop report may serve as a resource for our research community as well as for journal editors and reviewers of sequencing-based manuscript submissions in our research field.

Genomic Predictors of Asthma Phenotypes and Treatment Response

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

Interleukin 33 Selectively Augments Rhinovirus-Induced Type 2 Immune Responses in Asthmatic but not Healthy People

Interleukin- 33 (IL-33) is an epithelial-derived cytokine that initiates type 2 immune responses to allergens, though whether IL-33 has the ability to modify responses to respiratory viral infections remains unclear. This study aimed to investigate the effects of IL-33 on rhinovirus (RV)-induced immune responses by circulating leukocytes from people with allergic asthma, and how this response may differ from non-allergic controls. Our experimental approach involved co-exposing peripheral blood mononuclear cells to IL-33 and RV in order to model how the functions of virus-responsive lymphocytes could be modified after recruitment to an airway environment enriched in IL-33. In the current study, IL-33 enhanced RV-induced IL-5 and IL-13 release by cells from people with allergic asthma, but had no effect on IL-5 and IL-13 production by cells from healthy donors. In asthmatic individuals, IL-33 also enhanced mRNA and surface protein expression of ST2 (the IL-33 receptor IL1RL1), while soluble ST2 concentrations were low. In contrast, IL-33 had no effect on mRNA and surface expression of ST2 in healthy individuals. In people with allergic asthma, RV-activated ST2⁺ innate lymphoid cells (ST2⁺ILC) were the predominant source of IL-33 augmented IL-13 release. In contrast, RV-activated natural killer cells (NK cells) were the predominant source of IL-33 augmented IFNγ release in healthy individuals. This suggests that the effects of IL-33 on the cellular immune response to RV differ between asthmatic and healthy individuals. These findings provide a mechanism by which RV infections and IL-33 might interact in asthmatic individuals to exacerbate type 2 immune responses and allergic airway inflammation.

LPS inactivation by a host lipase allows lung epithelial cell sensitization for allergic asthma

This study provides strong evidence that intestinal commensal LPS desensitizes lung epithelial cells and therefore diminishes allergic responses to inhaled allergens. A host lipase, acyloxyacyl hydrolase (AOAH), prevents the desensitization by inactivating commensal LPS.

Allergic asthma is a chronic inflammatory disease primarily mediated by Th2 immune mechanisms. Numerous studies have suggested that early life exposure to lipopolysaccharide (LPS) is negatively associated with allergic asthma. One proposed mechanism invokes desensitization of lung epithelial cells by LPS. We report here that acyloxyacyl hydrolase (AOAH), a host lipase that degrades and inactivates LPS, renders mice more susceptible to house dust mite (HDM)–induced allergic asthma. Lung epithelial cells from Aoah^−/− mice are refractory to HDM stimulation, decreasing dendritic cell activation and Th2 responses. Antibiotic treatment that diminished commensal LPS-producing bacteria normalized Aoah^−/− responses to HDM, while giving LPS intrarectally ameliorated asthma. Aoah^−/− mouse feces, plasma, and lungs contained more bioactive LPS than did those of Aoah^+/+ mice. By inactivating commensal LPS, AOAH thus prevents desensitization of lung epithelial cells. An enzyme that prevents severe lung inflammation/injury in Gram-negative bacterial pneumonia has the seemingly paradoxical effect of predisposing to a Th2-mediated airway disease.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Proteomic and Phosphoproteomic Changes Induced by Prolonged Activation of Human Eosinophils with IL-3

Purified human eosinophils treated for 18-24 h with IL-3 adopt a unique activated phenotype marked by increased reactivity to aggregated immunoglobulin-G (IgG). To characterize this phenotype, we quantified protein abundance and phosphorylation by multiplexed isobaric labeling combined with high-resolution mass spectrometry. Purified blood eosinophils of five individuals were treated with IL-3 or no cytokine for 20 h, and comparative data were obtained on abundance of 5385 proteins and phosphorylation at 7330 sites. The 1150 proteins that were significantly up-regulated ( q < 0.05, pairwise t test with Benjamini-Hochberg correction) by IL-3 included the IL3RA and CSF2RB subunits of the IL-3 receptor, the low-affinity receptor for IgG (FCGR2B), 96 proteins involved in protein translation, and 55 proteins involved in cytoskeleton organization. Among the 703 proteins that decreased were 78 mitochondrial proteins. Dynamic regulation of protein phosphorylation was detected at 4218 sites. These included multiple serines in CSF2RB; Y694 of STAT5, a key site of activating phosphorylation downstream of IL3RA/CSF2RB; and multiple sites in RPS6KA1, RPS6, and EIF4B, which are responsible for translational initiation. We conclude that IL-3 up-regulates overall protein synthesis and targets specific proteins for up-regulation, including its own receptor.

Embracing polygenicity: a review of methods and tools for psychiatric genetics research

The availability of genome-wide genetic data on hundreds of thousands of people has led to an equally rapid growth in methodologies available to analyse these data. While the motivation for undertaking genome-wide association studies (GWAS) is identification of genetic markers associated with complex traits, once generated these data can be used for many other analyses. GWAS have demonstrated that complex traits exhibit a highly polygenic genetic architecture, often with shared genetic risk factors across traits. New methods to analyse data from GWAS are increasingly being used to address a diverse set of questions about the aetiology of complex traits and diseases, including psychiatric disorders. Here, we give an overview of some of these methods and present examples of how they have contributed to our understanding of psychiatric disorders. We consider: (i) estimation of the extent of genetic influence on traits, (ii) uncovering of shared genetic control between traits, (iii) predictions of genetic risk for individuals, (iv) uncovering of causal relationships between traits, (v) identifying causal single-nucleotide polymorphisms and genes or (vi) the detection of genetic heterogeneity. This classification helps organise the large number of recently developed methods, although some could be placed in more than one category. While some methods require GWAS data on individual people, others simply use GWAS summary statistics data, allowing novel well-powered analyses to be conducted at a low computational burden.

Tools and drugs for uracil nucleotide-activated P2Y receptors

P2Y receptors (P2YRs) are a family of G protein-coupled receptors activated by extracellular nucleotides. Physiological P2YR agonists include purine and pyrimidine nucleoside di- and triphosphates, such as ATP, ADP, UTP, UDP, nucleotide sugars, and dinucleotides. Eight subtypes exist, P2Y₁, P2Y₂, P2Y₄, P2Y₆, P2Y₁₁, P2Y₁₂, P2Y₁₃, and P2Y₁₄, which represent current or potential future drug targets. Here we provide a comprehensive overview of ligands for the subgroup of the P2YR family that is activated by uracil nucleotides: P2Y₂ (UTP, also ATP and dinucleotides), P2Y₄ (UTP), P2Y₆ (UDP), and P2Y₁₄ (UDP, UDP-glucose, UDP-galactose). The physiological agonists are metabolically unstable due to their fast hydrolysis by ectonucleotidases. A number of agonists with increased potency, subtype-selectivity and/or enzymatic stability have been developed in recent years. Useful P2Y₂R agonists include MRS2698 (6-01, highly selective) and PSB-1114 (6-05, increased metabolic stability). A potent and selective P2Y₂R antagonist is AR-C118925 (10-01). For studies of the P2Y₄R, MRS4062 (3-15) may be used as a selective agonist, while PSB-16133 (10-06) is a selective antagonist. Several potent P2Y₆R agonists have been developed including 5-methoxyuridine 5'-O-((R_p)α-boranodiphosphate) (6-12), PSB-0474 (3-11), and MRS2693 (3-26). The isocyanate MRS2578 (10-08) is used as a selective P2Y₆R antagonist, although its reactivity and low water-solubility are limiting. With MRS2905 (6-08), a potent and metabolically stable P2Y₁₄R agonist is available, while PPTN (10-14) represents a potent and selective P2Y₁₄R antagonist. The radioligand [³H]UDP can be used to label P2Y₁₄Rs. In addition, several fluorescent probes have been developed. Uracil nucleotide-activated P2YRs show great potential as drug targets, especially in inflammation, cancer, cardiovascular and neurodegenerative diseases.

Analysis combining correlated glaucoma traits identifies five new risk loci for open-angle glaucoma

Open-angle glaucoma (OAG) is a major cause of blindness worldwide. To identify new risk loci for OAG, we performed a genome-wide association study in 3,071 OAG cases and 6,750 unscreened controls, and meta-analysed the results with GWAS data for intraocular pressure (IOP) and optic disc parameters (the overall meta-analysis sample size varying between 32,000 to 48,000 participants), which are glaucoma-related traits. We identified and independently validated four novel genome-wide significant associations within or near MYOF and CYP26A1, LINC02052 and CRYGS, LMX1B, and LMO7 using single variant tests, one additional locus (C9) using gene-based tests, and two genetic pathways - "response to fluid shear stress" and "abnormal retina morphology" - in pathway-based tests. Interestingly, some of the new risk loci contribute to risk of other genetically-correlated eye diseases including myopia and age-related macular degeneration. To our knowledge, this study is the first integrative study to combine genetic data from OAG and its correlated traits to identify new risk variants and genetic pathways, highlighting the future potential of combining genetic data from genetically-correlated eye traits for the purpose of gene discovery and mapping.