Search

Background: Allergic rhinitis (AR) is strongly associated with a type 2 response, characterized by the cytokines IL-5, IL-4 and IL-13. Several studies have implicated ILC2 and TH2A (CD161+ TH2) but it is not yet entirely clear which subsets are driving the common allergic reactions underlying AR. The objective of this study aims to identify critical pathogenic cell populations associated with AR and to determine their phenotype and functional contribution to disease progression. Methods: We identified integral allergic-specific cell types by transcriptomic sequencing. Whole blood, PBMCs and plasma from a cross-sectional cohort of 216 individuals were analysed by 9-colour flow cytometry and ultra-sensitive cytokine bead arrays using unsupervised clustering algorithms. Clinically active AR cases were further analysed by functional mass cytometry to define phenotype and cytokine secretion (IL-2, IL-3, IL-4, IL-5, IL-9, IL13, IL-17A and IL-22) as well as the expression of the hematopoietic prostaglandin D synthetase (HPGDS). Results: The unbiased analysis revealed that atopy and AR manifestation corelated only with eosinophils, plasma IL-5 and CD161+ TH2 cells. In-depth characterization further revealed that the CD45RB CD161+ TH2 subset were most closely associated with severity. This subset is able to concomitantly secrete multiple cytokines including IL-5, IL-13 and IL-4 and has been previously reported to be associated with other eosinophilic allergies. Conclusion: CD45RB CD161+ TH2 have key roles in driving the allergic response in AR. Neutralizing the CD45RB CD161+ TH2 subset should disrupt the allergic response pathway, thus providing a target for lasting therapeutic interventions.

Human rhinoviruses (HRVs) are frequently associated with asthma exacerbations, and have been found in the airways of asthmatic patients. While HRV-induced acute infection is well-documented, it is less clear whether the nasal epithelium sustains prolonged HRV infections along with the associated activation of host immune responses.To investigate sustainably regulated host responses of human nasal epithelial cells (hNECs) during HRV persistence.Using a time-course study, HRV16 persistence and viral replication dynamics were established using an in vitro infection model of hNECs. RNA sequencing was performed on hNECs in the early and late stages of infection at 3 and 14 days post-infection (dpi), respectively. The functional enrichment of differentially expressed genes (DEGs) was evaluated using gene ontology (GO) and Ingenuity pathway analysis.HRV RNA and protein expression persisted throughout prolonged infections, even after decreased production of infectious virus progeny. GO analysis of unique DEGs indicated altered regulation of pathways related to ciliary function and airway remodeling at 3 dpi and serine-type endopeptidase activity at 14 dpi. The functional enrichment of shared DEGs between the two time-points was related to interferon (IFN) and cytoplasmic pattern recognition receptor (PRR) signaling pathways. Validation of the sustained regulation of candidate genes confirmed the persistent expression of RIG-I and revealed its close co-regulation with interferon-stimulated genes (ISGs) during HRV persistence.The persistence of HRV RNA does not necessarily indicate an active infection during prolonged infection. The sustained expression of RIG-I and ISGs in response to viral RNA persistence highlights the importance of assessing how immune-activating host factors can change during active HRV infection and the immune regulation that persists thereafter.

Background: Some children with asthma experience exacerbations despite long-acting beta2-agonist (LABA) treatment. While this variability is partly caused by genetic variation, no genome-wide study until now has investigated which genetic factors associated with risk of exacerbations despite LABA use in children with asthma. We aimed to assess whether genetic variation was associated with exacerbations in children treated with LABA from a global consortium. Methods: A meta-analysis of genome-wide association studies (meta-GWAS) was performed in 1,425 children and young adults with asthma (age 6-21 years) with reported regular use of LABA from six studies within the PiCA consortium using a random effects model. The primary outcome of each study was defined as any exacerbation within the past 6 or 12 months, including at least one of the following: 1) hospital admissions for asthma, 2) a course of oral corticosteroids or 3) emergency room visits because of asthma. Results: Genome-wide association results for a total of 82 996 common single nucleotide polymorphisms (SNPs, MAF ≥1%) with high imputation quality were meta-analysed. Eight independent variants were suggestively (P-value threshold ≤5 × 10−6) associated with exacerbations despite LABA use. Conclusion: No strong effects of single nucleotide polymorphisms (SNPs) on exacerbations during LABA use were identified. We identified two loci (TBX3 and EPHA7) that were previously implicated in the response to short-acting beta2-agonists (SABA). These loci merit further investigation in response to LABA and SABA use.

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.Analyses of expression profiles from whole blood of 31,684 individuals identify cis-expression quantitative trait loci (eQTL) effects for 88% of genes and trans-eQTL effects for 37% of trait-associated variants.

Background Gene expression patterns (particularly, cilia-associated genes) of nasal mucosa, the first-line defense system, in allergic rhinitis (AR) are not well understood. Objective We sought to screen for AR-associated genes in inferior turbinate (IT) from patients with AR, and to validate the expression of common cilia-related genes and ciliary shedding. Methods Prime View™ Human Gene Expression Array, which consisted of more than 530 000 probes covering more than 36 000 transcripts and variants, was employed to compare individual gene expression of ITs from control subjects (n = 11) and patients with AR (n = 19). Gene ontology (GO) analysis was performed with Cytoscape software. Eight of the common cilia-related genes were validated with quantitative polymerase chain reaction. We applied a semiquantitative scoring system for immunofluorescence assay to demonstrate ciliary shedding in 5 areas per paraffin section, with individual sections being scored between 0 (normal ciliary distribution) and 1 (ciliary shedding). Results Compared with control subjects, 160 (38 upregulated and 122 downregulated) genes were differentially expressed for at least 2 folds (all P Conclusion The significant aberrant cilia-related gene expression, revealed by microarray assays, might be the critical driver of AR where ciliary shedding is prominent.

Key immune signatures of SARS-CoV-2 infection may associate with either adverse immune reactions (severity) or simply an ongoing anti-viral response (temporality); how immune signatures contribute to severe manifestations and/or temporal progression of disease and whether longer disease duration correlates with severity remain unknown. Patient blood was comprehensively immunophenotyped via mass cytometry and multiplex cytokine arrays, leading to the identification of 327 basic subsets that were further stratified into more than 5000 immunotypes and correlated with 28 plasma cytokines. Low-density neutrophil abundance was closely correlated with hepatocyte growth factor levels, which in turn correlated with disease severity. Deep analysis also revealed additional players, namely conventional type 2 dendritic cells, natural killer T cells, plasmablasts and CD16+ monocytes, that can influence COVID-19 severity independent of temporal progression. Herein, we provide interactive network analysis and data visualization tools to facilitate data mining and hypothesis generation for elucidating COVID-19 pathogenesis.

SARS-CoV-2 is the novel coronavirus responsible for the current COVID-19 pandemic. Severe complications are observed only in a small proportion of infected patients but the cellular mechanisms underlying this progression are still unknown. Comprehensive flow cytometry of whole blood samples from 54 COVID-19 patients reveals a dramatic increase in the number of immature neutrophils. This increase strongly correlates with disease severity and is associated with elevated IL-6 and IP-10 levels, two key players in the cytokine storm. The most pronounced decrease in cell counts is observed for CD8 T-cells and VD2 γδ T-cells, which both exhibit increased differentiation and activation. ROC analysis reveals that the count ratio of immature neutrophils to VD2 (or CD8) T-cells predicts pneumonia onset (0.9071) as well as hypoxia onset (0.8908) with high sensitivity and specificity. It would thus be a useful prognostic marker for preventive patient management and improved healthcare resource management., COVID-19 severity is associated with cytokine levels and lymphopenia, but the role of immune cell subsets is not well understood. Here the authors immunophenotype whole blood samples from 54 COVID-19 patients and find that the immature neutrophil-to-VD2 T-cell ratio is associated with severe COVID-19.

SARS-CoV-2 is the novel coronavirus responsible for the current COVID-19 pandemic. Severe complications are observed only in a small proportion of infected patients but the cellular mechanisms underlying this progression are still unknown. Comprehensive flow cytometry of whole blood samples from 54 COVID-19 patients revealed a dramatic increase in the number of immature neutrophils. This increase strongly correlated with disease severity and was associated with elevated IL-6 and IP-10 levels, two key players in the cytokine storm. The most pronounced decrease in cell counts was observed for CD8 T-cells and VD2 γδ T-cells, which both exhibited increased differentiation and activation. ROC analysis revealed that the count ratio of immature neutrophils to CD8 or VD2 T-cells predicts pneumonia onset (0.9071) as well as hypoxia onset (0.8908) with high sensitivity and specificity. It would thus be a useful prognostic marker for preventive patient management and improved healthcare resource management.

Background Atopic dermatitis (AD) is a common skin disease affecting up to 20% of the global population, with significant clinical heterogeneity and limited information about molecular subtypes and actionable biomarkers. Although alterations in the skin microbiome have been described in subjects with AD during progression to flare state, the prognostic value of baseline microbiome configurations has not been explored. Objective Our aim was to identify microbial signatures on AD skin that are predictive of disease fate. Methods Nonlesional skin of patients with AD and healthy control subjects were sampled at 2 time points separated by at least 4 weeks. Using whole metagenome analysis of skin microbiomes of patients with AD and control subjects (n = 49 and 189 samples), we identified distinct microbiome configurations (dermotypes A and B). Blood was collected for immunophenotyping, and skin surface samples were analyzed for correlations with natural moisturizing factors and antimicrobial peptides. Results Dermotypes were robust and validated across 2 additional cohorts (63 individuals), with strong enrichment of subjects with AD in dermotype B. Dermotype B was characterized by reduced microbial richness, depletion of Cutibacterium acnes, Dermacoccus and Methylobacterium species, individual-specific outlier abundance of Staphylococcus species (eg, S epidermidis, S capitis, S aureus), and enrichment in metabolic pathways (eg, branched chain amino acids and arginine biosynthesis) and virulence genes (eg, β-toxin, δ-toxin) that defined a pathogenic ecology. Skin surface and circulating host biomarkers exhibited a distinct microbial-associated signature that was further reflected in more severe itching, frequent flares, and increased disease severity in patients harboring the dermotype B microbiome. Conclusion We report distinct clusters of microbial profiles that delineate the role of microbiome configurations in AD heterogeneity, highlight a mechanism for ongoing inflammation, and provide prognostic utility toward microbiome-based disease stratification.

Novel wide array blood biomarkers of multisystem dysregulation, compared to conventional clinical and blood biomarkers, are potentially able to provide more accurate prognostic information of long-term mortality risks. We identified biomarker signatures of all-cause and disease-specific mortality from a comprehensive range of analytes related to six major physiological functions: cytokine, chemokine, and growth factors; glucose metabolism regulators and adipokines; adhesion molecules; acute phase response; pathogen-specific antibodies; and bone remodeling. A total of 144 elderly were prospectively followed up on mortality for median 136 months. Plasma levels of 93 biomarkers measured by enzyme-linked immunosorbent assay, Luminex, FlowCytomix, DNA quantification, and recombinant DNA technology at baseline were compared among deceased and surviving elderly and in a referent group of 72 healthy young adults. The elderly, and especially deceased elderly, exhibited differential profiles of the composite index of each physiological function from young adults. In Cox regression, we identified and validated in an independent cohort of 357 elderly the specific and common biomarkers predicting all-cause, cardiovascular disease-related, neoplasm-related, and respiratory disease-related mortalities after controlling age, sex, and major comorbidities. These biomarkers had 74.3% correct classification for deceased elderly from surviving elderly. We reported for the first time, stem cell growth factor-β and gastric inhibitory polypeptide as specific biomarkers of mortality risk.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD ( P = 1.71 × 10 –9 , odds ratio [OR] = 4.35) and rs149308743 in CARD9 ( P = 2.09 × 10 –8 , OR = 4.75); three low-frequency variants: rs76418789 in IL23R ( P = 1.03 × 10 –10 , OR = 1.36), rs146466242 in FLG ( P = 3.39 × 10 –12 , OR = 1.45), and rs55882956 in TYK2 ( P = 1.04 × 10 –6 , OR = 1.30); and two common variants: rs780668 in SLC29A3 ( P = 2.17 × 10 –9 , OR = 1.14) and rs181206 in IL27 ( P = 1.08 × 10 –7 , OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

Aim: International collaboration is needed to enable large-scale pharmacogenomics studies in childhood asthma. Here, we describe the design of the Pharmacogenomics in Childhood Asthma (PiCA) consortium. Materials & methods: Investigators of each study participating in PiCA provided data on the study characteristics by answering an online questionnaire. Results: A total of 21 studies, including 14,227 children/young persons (58% male), from 12 different countries are currently enrolled in the PiCA consortium. Fifty six percent of the patients are Caucasians. In total, 7619 were inhaled corticosteroid users. Among patients from 13 studies with available data on asthma exacerbations, a third reported exacerbations despite inhaled corticosteroid use. In the future pharmacogenomics studies within the consortium, the pharmacogenomics analyses will be performed separately in each center and the results will be meta-analyzed. Conclusion: PiCA is a valuable platform to perform pharmacogenetics studies within a multiethnic pediatric asthma population.

Neuropeptide Y, a widely circulating neurotransmitter, plays a pivotal role in energy balance, immunomodulation and asthma, and several NPY polymorphisms are promising genetic risk factors for asthma and obesity. We explored the associations of candidate NPY gene polymorphisms with prevalent asthma and its relationship with obesity in young adult asthma patients free of other chronic medical morbidity.Five common gene variants of NPY (rs16147 (-399T/C), rs17149106 (-602G/T), rs16140 (+1000C/G), rs5573 (+1201A/G), rs5574 (+5327C/T)) previously validated to account for most of the NPY expression in vitro and in vivo were investigated in 126 physician-diagnosed asthma patients without other chronic medical morbidity and 182 healthy controls (21-35years). Plasma levels of NPY, adiponectin, and CRP were determined using ELISA, and IL-6 was measured by Luminex in a subgroup of 70 patients and 69 age- and sex-matched healthy controls.In logistic regression models controlling for gender and obesity, the CT genotype of rs5574 (OR=0.54, 95%CI: 0.30-0.89) and the GT genotype of rs17149106 (OR=5.58, 95%CI: 1.09-28.54) were significantly associated with asthma. No significant interaction between NPY SNP polymorphisms and obesity were detected. Plasma NPY level was correlated with adiponectin levels (p0.05). Compared with the healthy controls, patients with asthma had higher BMI (p0.001), adiponectin (p0.05), IL-6 (p=0.001) and CRP (p0.001), and lower NPY levels (p0.01).The CT genotype of rs5574 and the GT genotype of rs17149106 are significantly associated with prevalent asthma.

The human nasal epithelium is the primary site of exposure to influenza virus, the initiator of host responses to influenza and the resultant pathologies. Influenza virus may cause serious respiratory infection resulting in major complications, as well as severe impairment of the airways. Here, we elucidated the global transcriptomic changes during H3N2 infection of human nasal epithelial cells from multiple individuals. Using RNA sequencing, we characterized the differentially-expressed genes and pathways associated with changes occurring at the nasal epithelium following infection. We used in vitro differentiated human nasal epithelial cell culture model derived from seven different donors who had no concurrent history of viral infections. Statistical analysis highlighted strong transcriptomic signatures significantly associated with 24 and 48 h after infection, but not at the earlier 8-h time point. In particular, we found that the influenza infection induced in the nasal epithelium early and altered responses in interferon gamma signaling, B-cell signaling, apoptosis, necrosis, smooth muscle proliferation, and metabolic alterations. These molecular events initiated at the infected nasal epithelium may potentially adversely impact the airway, and thus the genes we identified could serve as potential diagnostic biomarkers or therapeutic targets for influenza infection and associated disease management.

IntroductionThe pathways underlying chronic rhinosinusitis with nasal polyps (CRSwNP) are unclear. We conducted genome-wide gene expression analysis to determine pathways and candidate gene sets associated with CRSwNP.MethodsWe performed whole-transcriptome RNA sequencing on 42 polyp (CRSwNP-NP) and 33 paired nonpolyp inferior turbinate (CRSwNP-IT) tissues from patients with CRSwNP and 28 inferior turbinate samples from non-CRS controls (CS-IT). We analysed the differentially expressed genes (DEGs) and the gene sets that were enriched in functional pathways.ResultsPrincipal component-informed analysis revealed cilium function and immune regulation as the two main Gene Ontology (GO) categories differentiating CRSwNP patients from controls. We detected 6182 and 1592 DEGs between CRSwNP-NP versus CS-IT and between CRSwNP-NP versus CRSwNP-IT tissues, respectively. Atopy status did not have a major impact on gene expression in various tissues. GO analysis on these DEGs implicated extracellular matrix (ECM) disassembly, O-glycan processing, angiogenesis and host viral response in CRSwNP pathogenesis. Ingenuity Pathway Analysis identified significant enrichment of type 1 interferon signalling and axonal guidance canonical pathways, angiogenesis, and collagen and fibrotic changes in CRSwNP (CRSwNP-NP and CRSwNP-IT) tissues compared with CS-IT. Finally, gene set enrichment analysis implicated sets of genes co-regulated in processes associated with inflammatory response and aberrant cell differentiation in polyp formation.ConclusionsGene signatures involved in defective host defences (including cilia dysfunction and immune dysregulation), inflammation and abnormal metabolism of ECM are implicated in CRSwNP. Functional validation of these gene expression patterns will open opportunities for CRSwNP therapeutic interventions such as biologics and immunomodulators.

Background: Atopic dermatitis (AD) is a chronic allergic disease typically accompanied by atopy and thus, a tendency to develop allergic diseases such as allergic rhinitis, asthma or food allergies. Currently, individuals with AD are classified into those presenting with AD alone and those presenting with AD along with other allergic diseases (AD+). It is important to identify the various endophenotypes of AD using anthropometric, environmental, socio-economic, and disease history data in order to improve disease management. To characterize the phenotypic differences among Singaporean Chinese individuals with AD alone and AD+, and identify the socioeconomic, lifestyle, and environmental factors associated with these different presentations. Methods: Based on data collected via a standardized/validated questionnaire, 4604 participants (mean age: 22.1 years) were classified into three groups: 1) AD alone group; 2) AD with other allergic diseases group (AD+); and 3) Control group. Results: Participants were less sensitized to common inhalant allergens in the AD alone group versus the Control group (67% vs. 72%, respectively; p 23) was associated with the disease and the difference was more pronounced in the AD alone group compared to the AD+ group (Odds Ratio: 1.38; 95% Confidence Interval: 1.4–1.67; p

Resistin is a key cytokine associated with metabolic and inflammatory diseases. Especially in East Asian populations, the expression levels are strongly influenced by genetic polymorphisms. Mechanisms and functional implications of this genetic control are still unknown. By employing reporter assays, EMSA, inhibition studies, bisulphite sequencing, ChIP-Seq and gene-editing we show that the p50/p50 homodimer known to act as repressor for a number of pro-inflammatory genes plays a central role in the genetic regulation of resistin in monocytes along with promoter methylation. In the common RETN haplotype p50/p50 constitutively dampens the expression by binding to the promoter. In an Asian haplotype variant however this interaction is disrupted by the A allele of rs3219175. The SNP is in very close linkage to rs34861192, a CpG SNP, located 280 bp upstream which provides an allele-specific C-methylation site. rs34861192 is located in a 100 bp region found to be methylated in the common but not in the Asian haplotype, resulting in the latter having a higher basal expression, which also associates with elevated histone acetylation (H3K27ac). Genotype associations within cohort data of 200 East Asian individuals revealed significant associations between this haplotype and the plasma levels of factors such as TGF-b, S100B, sRAGE and IL-8 as well as with myeloid DC counts. Thus, the common RETN haplotype is tightly regulated by the epigenetic mechanism linked to p50/p50-binding. This control is lost in the Asian haplotype, which may have evolved to balance the antagonistic RETN effects on pathogen protection vs. metabolic and inflammatory disease induction.

SummaryWhile many disease-associated variants have been identified through genome-wide association studies, their downstream molecular consequences remain unclear.To identify these effects, we performedcis-andtrans-expressionquantitative trait locus (eQTL) analysis in blood from 31,684 individuals through the eQTLGen Consortium.We observed thatcis-eQTLs can be detected for 88% of the studied genes, but that they have a different genetic architecture compared to disease-associated variants, limiting our ability to usecis-eQTLs to pinpoint causal genes within susceptibility loci.In contrast, trans-eQTLs (detected for 37% of 10,317 studied trait-associated variants) were more informative. Multiple unlinked variants, associated to the same complex trait, often converged on trans-genes that are known to play central roles in disease etiology.We observed the same when ascertaining the effect of polygenic scores calculated for 1,263 genome-wide association study (GWAS) traits. Expression levels of 13% of the studied genes correlated with polygenic scores, and many resulting genes are known to drive these traits.

Background The IL-13 receptor α2 (IL-13Rα2) is a receptor for IL-13 which has conflicting roles in mediating IL-13 responses in the lower airway, with little known about its impact on upper airway diseases. We sought to investigate the expression of IL-13 receptors, IL-13Rα1 and IL-13Rα2, in chronically inflamed nasal epithelium, and explore IL-13-induced signaling pathways in an in vitro model of human nasal epithelial cells (hNECs). Methods The protein and mRNA expression levels of IL-13 and its receptors in nasal biopsies of patients with nasal polyps (NP) and healthy controls were evaluated. We investigated goblet cell stimulation with mucus hypersecretion induced by IL-13 (10 ng/mL, 72 hours) treatment in hNECs using a pseudostratified epithelium in air-liquid interface (ALI) culture. Results There were significant increases in IL-13, IL-13Rα1, and IL-13Rα2 mRNA and protein levels in NP epithelium with healthy controls as baseline. MUC5AC mRNA positively correlated with IL-13Rα2 (r = .5886, P = .002) but not with IL-13Rα1 in primary hNECs. IL-13 treatment resulted in a significant increase in mRNA and protein levels of IL-13Rα2 only in hNECs. IL-13 treatment induced an activation of extracellular signal-regulated kinases (ERK)1/2 and an upregulation of C-JUN, where the IL-13-induced effects on hNECs could be attenuated by ERK1/2 inhibitor (50 μmol/L) or dexamethasone (10-4 -10-7 mol/L) treatment. Conclusions IL-13Rα2 has a potential role in IL-13-induced MUC5AC and ciliary changes through ERK1/2 signal pathway in the nasal epithelium. IL-13Rα2 may contribute to airway inflammation and aberrant remodeling which are the main pathological features of CRSwNP.

Meiotic “point recombination” refers to homologue recombination events affecting only individual SNPs. Driven mostly by gene conversion, it is common process that allows for a gradual adaptation and maturation of haplotypes during genetic evolution. In contrast to crossover recombination it is not tied to predetermined recombination sites and therefore assumed to occur largely randomly. Our analysis of archaic human haplotypes however revealed striking differences in the local point recombination rate. A linkage-study of 1.9 million SNPs defined by the sequence of denisovan hominids revealed low rates in introns and quiescent intergenic regions but high rates in splice sites, exons, 5’- and 3’-UTRs, and CpG islands. Correlations with ChIP-Seq tracks from ENCODE and other public sources identified a number of epigenetic modifications, that associated directly with these recombination events. A particularly tight association was observed for 5-hydroxymethylcytosine marks (5hmC). The mark was enriched in virtually all of the functional regions associated with elevated point recombination rates, including CpG islands and ‘poised’ bivalent regions. As intermediate of oxidative demethylation, 5hmC is also a marker of recently opened gene loci. The data, thus, supports a model of ‘guided’ evolution, in which point recombination is directed by 5hmC marks towards the functionally relevant regions.

Genome-wide association studies (GWAS) have led to the discovery of several susceptibility loci for leprosy with robust evidence, providing biological insight into the role of host genetic factors in mycobacterial infection. However, the identified loci only partially explain disease heritability, and additional genetic risk factors remain to be discovered. We performed a 3-stage GWAS of leprosy in the Chinese population using 8,313 cases and 16,017 controls. Besides confirming all previously published loci, we discovered six new susceptibility loci, and further gene prioritization analysis of these loci implicated BATF3, CCDC88B and CIITA-SOCS1 as new susceptibility genes for leprosy. A systematic evaluation of pleiotropic effects demonstrated a high tendency for leprosy susceptibility loci to show association with autoimmunity and inflammatory diseases. Further analysis suggests that molecular sensing of infection might have a similar pathogenic role across these diseases, whereas immune responses have discordant roles in infectious and inflammatory diseases.

Background: Chronic obstructive pulmonary disease (COPD) is currently the third leading cause of death and there is a huge unmet clinical need to identify disease biomarkers in peripheral blood. Compared to gene level differential expression approaches to identify gene signatures, network analyses provide a biologically intuitive approach which leverages the co-expression patterns in the transcriptome to identify modules of co-expressed genes. Methods: A weighted gene co-expression network analysis (WGCNA) was applied to peripheral blood transcriptome from 238 COPD subjects to discover co-expressed gene modules. We then determined the relationship between these modules and forced expiratory volume in 1 s (FEV1). In a second, independent cohort of 381 subjects, we determined the preservation of these modules and their relationship with FEV1. For those modules that were significantly related to FEV1, we determined the biological processes as well as the blood cell-specific gene expression that were over-represented using additional external datasets. Results: Using WGCNA, we identified 17 modules of co-expressed genes in the discovery cohort. Three of these modules were significantly correlated with FEV1 (FDR

Background Expression quantitative trait loci (eQTL) databases represent a valuable resource to link disease-associated SNPs to specific candidate genes whose gene expression is significantly modulated by the SNP under investigation. We previously identified signal inhibitory receptor on leukocytes-1 (SIRL-1) as a powerful regulator of human innate immune cell function. While it is constitutively high expressed on neutrophils, on monocytes the SIRL-1 surface expression varies strongly between individuals. The underlying mechanism of regulation, its genetic control as well as potential clinical implications had not been explored yet. Methods Whole blood eQTL data of a Chinese cohort was used to identify SNPs regulating the expression of VSTM1, the gene encoding SIRL-1. The genotype effect was validated by flow cytometry (cell surface expression), correlated with electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) and bisulfite sequencing (C-methylation) and its functional impact studied the inhibition of reactive oxygen species (ROS). Results We found a significant association of a single CpG-SNP, rs612529T/C, located in the promoter of VSTM1. Through flow cytometry analysis we confirmed that primarily in the monocytes the protein level of SIRL-1 is strongly associated with genotype of this SNP. In monocytes, the T allele of this SNP facilitates binding of the transcription factors YY1 and PU.1, of which the latter has been recently shown to act as docking site for modifiers of DNA methylation. In line with this notion rs612529T associates with a complete demethylation of the VSTM1 promoter correlating with the allele-specific upregulation of SIRL-1 expression. In monocytes, this upregulation strongly impacts the IgA-induced production of ROS by these cells. Through targeted association analysis we found a significant Meta P value of 1.14 × 10–6 for rs612529 for association to atopic dermatitis (AD). Conclusion Low expression of SIRL-1 on monocytes is associated with an increased risk for the manifestation of an inflammatory skin disease. It thus underlines the role of both the cell subset and this inhibitory immune receptor in maintaining immune homeostasis in the skin. Notably, the genetic regulation is achieved by a single CpG-SNP, which controls the overall methylation state of the promoter gene segment. Electronic supplementary material The online version of this article (doi:10.1186/s13073-017-0404-6) contains supplementary material, which is available to authorized users.

Background: Genetic variation may partly explain heterogeneity in response to asthma treatment . Recently, the P harmacogenomics i n C hildhood A sthma (PiCA) consortium was initiated to perform large-scale pharmacogenomics studies. Aims: To describe the characteristics of the asthmatic children currently included in PiCA and to assess treatment response definitions within this growing consortium. Methods: An online survey including questions about data collection, patient characteristics and outcome phenotypes was send to the 20 studies included in the PiCA consortium. Results: In total, 14,634 children/young adults (up to 25 years) from 11 different countries are enrolled in the PiCA consortium. Preliminary analysis of patients showed that 60% of the patients are male. PiCA is ethnically diverse: 47% Caucasians, 25% Latinos, 14% Asians, 9% African-Americans, 2% have mixed/other-ethnicity and 3% is unknown. PiCA currently includes data of 2,056 severe asthmatics. Furthermore, PiCA has data of 8,907 short acting β 2 -adrenergic agonist (SABA) users, 7,908 ICS users, 2,054 leukotriene modifier (LTM) users and 1,543 users of long acting β 2 -adrenergic agonists (LABA). Ninety-five percent of the studies have information regarding asthma exacerbations and lung function measurements and 80% have data on asthma symptoms. Discussion: The children within the PiCA consortium are a reflection of the global heterogeneous pediatric asthma population. Different outcome measures reflect different dimensions of asthma. By classifying three commonly used outcome measures using data that are available within the PiCA consortium, we will be able to study distinct response phenotypes in relation to genetics and clinical factors.

Genetic polymorphism within the 9q32 locus is linked with increased risk of several diseases, including Crohn’s disease (CD), primary biliary cholangitis (PBC) and leprosy. The most likely disease-causing gene within 9q32 is TNFSF15, which encodes the pro-inflammatory cytokine TNF super-family member 15, but it was unknown whether these disparate diseases were associated with the same genetic variance in 9q32, and how variance within this locus might contribute to pathology. Using genetic data from published studies on CD, PBC and leprosy we revealed that bearing a T allele at rs6478108/rs6478109 (r2 = 1) or rs4979462 was significantly associated with increased risk of CD and decreased risk of leprosy, while the T allele at rs4979462 was associated with significantly increased risk of PBC. In vitro analyses showed that the rs6478109 genotype significantly affected TNFSF15 expression in cells from whole blood of controls, while functional annotation using publicly-available data revealed the broad cell type/tissue-specific regulatory potential of variance at rs6478109 or rs4979462. In summary, we provide evidence that variance within TNFSF15 has the potential to affect cytokine expression across a range of tissues and thereby contribute to protection from infectious diseases such as leprosy, while increasing the risk of immune-mediated diseases including CD and PBC.

Background Cohort studies indicated that in certain individuals the basophils do not respond toward allergens due to a desensitization of their Fc epsilon receptor pathway. Cause and functional role as well as the implications on allergic reactions, however, are not clear yet. Methods A cross-sectional study was carried out in the tropical urban environment of Singapore, where the allergic response is dominated by a single allergen (house dust mite; HDM). Blood samples were collected from 476 individuals and analyzed comprehensively to correlate the functional state of their basophils with the clinical state as well as the composition of the cellular and soluble plasma components. Results Inactivation of basophils ('basophil anergy') was observed in about 10% of the cohort. It was associated with a downregulation of basophil Syk and an apparent reduction in the incidence of allergic rhinitis. Correlations on the cohort level suggest that it represents a transitional state to be passed through during the interconversion of responder and nonresponder state. Conclusions Basophil anergy thus seems to function as activation barrier to prevent unwanted reactions against minor allergens. It may therefore be relevant for diagnostic purposes or therapeutic interventions of allergic diseases.

Neutrophils are an abundant immune cell type involved in both antimicrobial defence and autoimmunity. The regulation of their gene expression, however, is still largely unknown. Here we report an eQTL study on isolated neutrophils from 114 healthy individuals of Chinese ethnicity, identifying 21,210 eQTLs on 832 unique genes. Unsupervised clustering analysis of these eQTLs confirms their role in inflammatory responses and immunological diseases but also indicates strong involvement in dermatological pathologies. One of the strongest eQTL identified (rs2058660) is also the tagSNP of a linkage block reported to affect leprosy and Crohn's disease in opposite directions. In a functional study, we can link the C allele with low expression of the β-chain of IL18-receptor (IL18RAP). In neutrophils, this results in a reduced responsiveness to IL-18, detected both on the RNA and protein level. Thus, the polymorphic regulation of human neutrophils can impact beneficial as well as pathological inflammatory responses., Neutrophils are abundant immune cells important for antimicrobial defence and in autoimmunity. Here, by mapping expression quantitative trait loci (eQTL) in neutrophils of Chinese ethnicity from Singapore, Andiappan et al. provide a resource for understanding immune-related trait associated genetic variants.

IgA nephropathy (IgAN) is one of the most common primary glomerulonephritis. Previously identified genome-wide association study (GWAS) loci explain only a fraction of disease risk. To identify novel susceptibility loci in Han Chinese, we conduct a four-stage GWAS comprising 8,313 cases and 19,680 controls. Here, we show novel associations at ST6GAL1 on 3q27.3 (rs7634389, odds ratio (OR)=1.13, P=7.27 × 10−10), ACCS on 11p11.2 (rs2074038, OR=1.14, P=3.93 × 10−9) and ODF1-KLF10 on 8q22.3 (rs2033562, OR=1.13, P=1.41 × 10−9), validate a recently reported association at ITGAX-ITGAM on 16p11.2 (rs7190997, OR=1.22, P=2.26 × 10−19), and identify three independent signals within the DEFA locus (rs2738058, P=1.15 × 10−19; rs12716641, P=9.53 × 10−9; rs9314614, P=4.25 × 10−9, multivariate association). The risk variants on 3q27.3 and 11p11.2 show strong association with mRNA expression levels in blood cells while allele frequencies of the risk variants within ST6GAL1, ACCS and DEFA correlate with geographical variation in IgAN prevalence. Our findings expand our understanding on IgAN genetic susceptibility and provide novel biological insights into molecular mechanisms underlying IgAN., IgA nephropathy is a major cause of end-stage renal disease in China, occurring at a high frequency in Asian populations. Here Li and colleagues conduct a four-stage genome-wide association study in a Chinese population, identifying novel loci and variants associated with disease risk.

Background Brain-derived neurotrophic factor (BDNF) is a secretory protein that has been implicated in the pathogenesis of allergic rhinitis (AR), atopic asthma, and eczema, but it is currently unknown whether BDNF polymorphisms influence susceptibility to moderate-to-severe AR. Objective We sought to identify disease associations and the functional effect of BDNF genetic variants in patients with moderate-to-severe AR. Methods Tagging single nucleotide polymorphisms (SNPs) spanning the BDNF gene were selected from the human HapMap Han Chinese from Beijing (CHB) data set, and associations with moderate-to-severe AR were assessed in 2 independent cohorts of Chinese patients (2216 from Shandong province and 1239 living in Singapore). The functional effects of the BDNF genetic variants were determined by using both in vitro and ex vivo assays. Results The tagging SNP rs10767664 was significantly associated with the risk of moderate-to-severe AR in both Singapore Chinese ( P = .0017; odds ratio, 1.324) and Shandong Chinese populations ( P = .039; odds ratio, 1.180). The coding nonsynonymous SNP rs6265 was in perfect linkage with rs10767664 and conferred increased BDNF protein secretion by a human cell line in vitro . Subjects bearing the AA genotype of rs10767664 exhibited increased risk of moderate-to-severe AR and displayed increased BDNF protein and total IgE levels in plasma. Using a large-scale expression quantitative trait locus study, we demonstrated that BDNF SNPs are significantly associated with altered BDNF concentrations in peripheral blood. Conclusion A common genetic variant of the BDNF gene is associated with increased risk of moderate-to-severe AR, and the AA genotype is associated with increased BDNF mRNA levels in peripheral blood. Together, these data indicate that functional BDNF gene variants increase the risk of moderate-to-severe AR.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn’s disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus., Author Summary Many variants in the genome, including variants associated with disease, affect the expression of genes. These so-called expression quantitative trait loci (eQTL) can be used to gain insight in the downstream consequences of disease. While it has been shown that many disease-associated variants alter gene expression in a cell-type dependent manner, eQTL datasets for specific cell types may not always be available and their sample size is often limited. We present a method that is able to detect cell type specific effects within eQTL datasets that have been generated from whole tissues (which may be composed of many cell types), in our case whole blood. By combining numerous whole blood datasets through meta-analysis, we show that we are able to detect eQTL effects that are specific for neutrophils and lymphocytes (two blood cell types). Additionally, we show that the variants associated with some diseases may preferentially alter the gene expression in one of these cell types. We conclude that our method is an alternative method to detect cell type specific eQTL effects, that may complement generating cell type specific eQTL datasets and that may be applied on other cell types and tissues as well.

Background Allergic rhinitis (AR) and asthma are common allergic conditions with a shared genetic component to their cause. The 17q12-21 locus includes several genes that have been linked to asthma susceptibility, but the role of this locus in AR is unclear. Asthma and AR in adults of Chinese ethnicity in Singapore are predominately caused by sensitization against house dust mites with a nearly complete penetrance of the allergen, which presents a unique opportunity for accurately identifying genetic associations with allergic diseases. Objective We sought to define the functional role of 17q12-21 in patients with AR and allergic asthma. Methods We asked whether single nucleotide polymorphisms (SNPs) in the 17q12-21 locus were associated with AR or asthma in a cohort of 3460 ethnic Chinese subjects residing in Singapore (1435 in the discovery phase and 2025 in the validation phase). Full-blood mRNA gene expression data, plasma IgE levels, and immune cell frequencies in peripheral blood were tested against the tag SNP genotypes. Luciferase assays were used to measure the effect of putative promoter SNPs on expression of the asthma-associated orosomucoid-like 3 gene (ORMDL3) . Results Within 17q12-21, only the tag SNP rs8076131 was significantly associated with asthma ( P = 8.53 × 10 −10 ; odds ratio, 0.6715), and AR status was independent of SNPs in this region. C-A alleles at rs8076131 resulted in significantly increased ORMDL3 expression in HEK293 cells in vitro relative to T-G alleles. Moreover, subjects with the risk genotype AA exhibited significantly higher total IgE levels and higher blood eosinophil counts than those with the lower-risk genotypes. Conclusion The 17q12-21 locus has a strong genetic association with allergic asthma but not with AR. The polymorphic effect of this locus is attributed to the linkage set tagged by rs8076131, which affects the expression of ORMDL3 , protein phosphatase 1, regulatory inhibitor subunit 1B ( PPP1R1B ), zona pellucida binding protein 2 ( ZPBP2 ), and gasdermin B (GSDMB) and is correlated with high IgE levels and eosinophil counts in subjects bearing the risk genotype.