Diagnostic Yield From a Nationwide Implementation of Precision Medicine for all Children With Cancer.

PURPOSE: Several studies have indicated that broad genomic characterization of childhood cancer provides diagnostically and/or therapeutically relevant information in selected high-risk cases. However, the extent to which such characterization offers clinically actionable data in a prospective broadly inclusive setting remains largely unexplored. METHODS: We implemented prospective whole-genome sequencing (WGS) of tumor and germline, complemented by whole-transcriptome sequencing (RNA-Seq) for all children diagnosed with a primary or relapsed solid malignancy in Sweden. Multidisciplinary molecular tumor boards were set up to integrate genomic data in the clinical decision process along with a medicolegal framework enabling secondary use of sequencing data for research purposes. RESULTS: During the study's first 14 months, 118 solid tumors from 117 patients were subjected to WGS, with complementary RNA-Seq for fusion gene detection in 52 tumors. There was no significant geographic bias in patient enrollment, and the included tumor types reflected the annual national incidence of pediatric solid tumor types. Of the 112 tumors with somatic mutations, 106 (95%) exhibited alterations with a clear clinical correlation. In 46 of 118 tumors (39%), sequencing only corroborated histopathological diagnoses, while in 59 cases (50%), it contributed to additional subclassification or detection of prognostic markers. Potential treatment targets were found in 31 patients (26%), most commonly ALK mutations/fusions (n = 4), RAS/RAF/MEK/ERK pathway mutations (n = 14), FGFR1 mutations/fusions (n = 5), IDH1 mutations (n = 2), and NTRK2 gene fusions (n = 2). In one patient, the tumor diagnosis was revised based on sequencing. Clinically relevant germline variants were detected in 8 of 94 patients (8.5%). CONCLUSION: Up-front, large-scale genomic characterization of pediatric solid malignancies provides diagnostically valuable data in the majority of patients also in a largely unselected cohort.

The complement receptor C3AR constitutes a novel therapeutic target in NPM1-mutated AML.

Mutated nucleophosmin 1 (NPM1) is the most common genetic alteration in acute myeloid leukemia (AML), found in ∼30% of cases. Although mutations in this gene are considered favorable according to current risk stratification guidelines, a large fraction of patients will experience relapse, demonstrating the urgent need for new treatment options. Therefore, we aimed to identify cell surface proteins specifically expressed on NPM1-mutated AML cells, allowing for potential targeting with antibody-based therapies. Herein, we report on an arrayed flow cytometry-based screen directed to 362 cell surface markers. In comparing the cell surface expression on NPM1-mutated AML cells with primitive (CD34+ CD38-) normal bone marrow cells, we identified the complement receptor C3AR as being specifically expressed in NPM1-mutated AML. By flow cytometry and single-cell RNA sequencing, we further show that normal hematopoietic stem and progenitor cells lack detectable C3AR gene and protein expression, making it particularly suitable as a target for antibody therapy. We also demonstrate that C3AR in combination with GPR56 distinguishes the leukemic stem cells (LSCs) in NPM1-mutated AML from the normal hematopoietic stem cells, defining the LSC population, as shown by transplantation into immunodeficient mice. Mechanistically, the stimulation of C3AR-expressing cells with C3a, the ligand of C3AR, leads to the activation of ERK1/2 and increased survival of AML cells, suggesting that this is an important signaling axis in this subtype of AML. Finally, we show that antibodies directed against C3AR efficiently elicit natural killer cell-mediated killing of primary AML cells ex vivo, highlighting C3AR as a candidate therapeutic target in NPM1-mutated AML.

A Study Protocol for Validation and Implementation of Whole-Genome and -Transcriptome Sequencing as a Comprehensive Precision Diagnostic Test in Acute Leukemias.

Background: Whole-genome sequencing (WGS) and whole-transcriptome sequencing (WTS), with the ability to provide comprehensive genomic information, have become the focal point of research interest as novel techniques that can support precision diagnostics in routine clinical care of patients with various cancer types, including hematological malignancies. This national multi-center study, led by Genomic Medicine Sweden, aims to evaluate whether combined application of WGS and WTS (WGTS) is technically feasible and can be implemented as an efficient diagnostic tool in patients with acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML). In addition to clinical impact assessment, a health-economic evaluation of such strategy will be performed. Methods and Analysis: The study comprises four phases (i.e., retrospective, prospective, real-time validation, and follow-up) including approximately 700 adult and pediatric Swedish AML and ALL patients. Results of WGS for tumor (90×) and normal/germline (30×) samples as well as WTS for tumors only will be compared to current standard of care diagnostics. Primary study endpoints are diagnostic efficiency and improved diagnostic yield. Secondary endpoints are technical and clinical feasibility for routine implementation, clinical utility, and health-economic impact. Discussion: Data from this national multi-center study will be used to evaluate clinical performance of the integrated WGTS diagnostic workflow compared with standard of care. The study will also elucidate clinical and health-economic impacts of a combined WGTS strategy when implemented in routine clinical care. Clinical Trial Registration: [https://doi.org/10.1186/ISRCTN66987142], identifier [ISRCTN66987142].

Transcriptomics paving the way for improved diagnostics and precision medicine of acute leukemia.

Transcriptional profiling of acute leukemia, specifically by RNA-sequencing or whole transcriptome sequencing (WTS), has provided fundamental insights into its underlying disease biology and allows unbiased detection of oncogenic gene fusions, as well as of gene expression signatures that can be used for improved disease classification. While used as a research tool for many years, RNA-sequencing is becoming increasingly used in clinical diagnostics. Here, we highlight key transcriptomic studies of acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML) that have improved our biological understanding of these heterogeneous malignant disorders and have paved the way for translation into clinical diagnostics. Recent single-cell transcriptomic studies of ALL and AML, which provide new insights into the cellular ecosystem of acute leukemia and point to future clinical utility, are also reviewed. Finally, we discuss current challenges that need to be overcome for a more wide-spread adoption of RNA-sequencing in clinical diagnostics and how this technology significantly can aid the identification of genetic alterations in current guidelines and of newly emerging disease entities, some of which are critical to identify because of the availability of targeted therapies, thereby paving the way for improved precision medicine of acute leukemia.

Clinical and genomic characterization of patients diagnosed with the provisional entity acute myeloid leukemia with BCR-ABL1, a Swedish population-based study.

Acute myeloid leukemia (AML) with t(9;22)(q34;q11), also known as AML with BCR-ABL1, is a rare, provisional entity in the WHO 2016 classification and is considered a high-risk disease according to the European LeukemiaNet 2017 risk stratification. We here present a retrospective, population-based study of this disease entity from the Swedish Acute Leukemia Registry. By strict clinical inclusion criteria we aimed to identify genetic markers further distinguishing AML with t(9;22) as a separate entity. Twenty-five patients were identified and next-generation sequencing using a 54-gene panel was performed in 21 cases. Interestingly, no mutations were found in NPM1, FLT3, or DNMT3A, three frequently mutated genes in AML. Instead, RUNX1 was the most commonly mutated gene, with aberrations present in 38% of the cases compared to around 10% in de novo AML. Additional mutations were identified in genes involved in RNA splicing (SRSF2, SF3B1) and chromatin regulation (ASXL1, STAG2, BCOR, BCORL1). Less frequently, mutations were found in IDH2, NRAS, TET2, and TP53. The mutational landscape exhibited a similar pattern as recently described in patients with chronic myeloid leukemia (CML) in myeloid blast crisis (BC). Despite the concomitant presence of BCR-ABL1 and RUNX1 mutations in our cohort, both features of high-risk AML, the RUNX1-mutated cases showed a superior overall survival compared to RUNX1 wildtype cases. Our results suggest that the molecular characteristics of AML with t(9;22)/BCR-ABL1 and CML in myeloid BC are similar and do not support a distinction of the two disease entities based on their underlying molecular alterations.

Clonal competition within complex evolutionary hierarchies shapes AML over time.

Clonal heterogeneity and evolution has major implications for disease progression and relapse in acute myeloid leukemia (AML). To model clonal dynamics in vivo, we serially transplanted 23 AML cases to immunodeficient mice and followed clonal composition for up to 15 months by whole-exome sequencing of 84 xenografts across two generations. We demonstrate vast changes in clonality that both progress and reverse over time, and define five patterns of clonal dynamics: Monoclonal, Stable, Loss, Expansion and Burst. We also show that subclonal expansion in vivo correlates with a more adverse prognosis. Furthermore, clonal expansion enabled detection of very rare clones with AML driver mutations that were undetectable by sequencing at diagnosis, demonstrating that the vast majority of AML cases harbor multiple clones already at diagnosis. Finally, the rise and fall of related clones enabled deconstruction of the complex evolutionary hierarchies of the clones that compete to shape AML over time.

Arrayed molecular barcoding identifies TNFSF13 as a positive regulator of acute myeloid leukemia-initiating cells.

Dysregulation of cytokines in the bone marrow (BM) microenvironment promotes acute myeloid leukemia (AML) cell growth. Due to the complexity and low throughput of in vivo stem-cell based assays, studying the role of cytokines in the BM niche in a screening setting is challenging. Here, we developed an ex vivo cytokine screen using 11 arrayed molecular barcodes, allowing for a competitive in vivo readout of leukemia-initiating capacity. With this approach, we assessed the effect of 114 murine cytokines on MLL-AF9 AML mouse cells and identified the tumor necrosis factor ligand superfamily member 13 (TNFSF13) as a positive regulator of leukemia-initiating cells. By using Tnfsf13-/- recipient mice, we confirmed that TNFSF13 supports leukemia initiation also under physiological conditions. TNFSF13 was secreted by normal myeloid cells but not by leukemia mouse cells, suggesting that mature myeloid BM cells support leukemia cells by secreting TNFSF13. TNFSF13 supported leukemia cell proliferation in an NF-κB-dependent manner by binding TNFRSF17 and suppressed apoptosis. Moreover, TNFSF13 supported the growth and survival of several human myeloid leukemia cell lines, demonstrating that our findings translate to human disease. Taken together, using arrayed molecular barcoding, we identified a previously unrecognized role of TNFSF13 as a positive regulator of AML-initiating cells. The arrayed barcoded screening methodology is not limited to cytokines and leukemia, but can be extended to other types of ex vivo screens, where a multiplexed in vivo read-out of stem cell functionality is needed.

Microenvironmental control of breast cancer subtype elicited through paracrine platelet-derived growth factor-CC signaling.

Breast tumors of the basal-like, hormone receptor-negative subtype remain an unmet clinical challenge, as there is high rate of recurrence and poor survival in patients following treatment. Coevolution of the malignant mammary epithelium and its underlying stroma instigates cancer-associated fibroblasts (CAFs) to support most, if not all, hallmarks of cancer progression. Here we delineate a previously unappreciated role for CAFs as determinants of the molecular subtype of breast cancer. We identified paracrine crosstalk between cancer cells expressing platelet-derived growth factor (PDGF)-CC and CAFs expressing the cognate receptors in human basal-like mammary carcinomas. Genetic or pharmacological intervention of PDGF-CC activity in mouse models of cancer resulted in conversion of basal-like breast cancers into a hormone receptor-positive state that enhanced sensitivity to endocrine therapy in previously resistant tumors. We conclude that specification of breast cancer to the basal-like subtype is under microenvironmental control and is therapeutically actionable.

Agonistic targeting of TLR1/TLR2 induces p38 MAPK-dependent apoptosis and NFκB-dependent differentiation of AML cells.

Acute myeloid leukemia (AML) is associated with poor survival, and there is a strong need to identify disease vulnerabilities that might reveal new treatment opportunities. Here, we found that Toll-like receptor 1 (TLR1) and TLR2 are upregulated on primary AML CD34+CD38- cells relative to corresponding normal bone marrow cells. Activating the TLR1/TLR2 complex by the agonist Pam3CSK4 in MLL-AF9-driven human AML resulted in induction of apoptosis by p38 MAPK-dependent activation of Caspase 3 and myeloid differentiation in a NFκB-dependent manner. By using murine Trp53-/-MLL-AF9 AML cells, we demonstrate that p53 is dispensable for Pam3CSK4-induced apoptosis and differentiation. Moreover, murine AML1-ETO9a-driven AML cells also were forced into apoptosis and differentiation on TLR1/TLR2 activation, demonstrating that the antileukemic effects observed were not confined to MLL-rearranged AML. We further evaluated whether Pam3CSK4 would exhibit selective antileukemic effects. Ex vivo Pam3CSK4 treatment inhibited murine and human leukemia-initiating cells, whereas murine normal hematopoietic stem and progenitor cells (HSPCs) were relatively less affected. Consistent with these findings, primary human AML cells across several genetic subtypes of AML were more vulnerable for TLR1/TLR2 activation relative to normal human HSPCs. In the MLL-AF9 AML mouse model, treatment with Pam3CSK4 provided proof of concept for in vivo therapeutic efficacy. Our results demonstrate that TLR1 and TLR2 are upregulated on primitive AML cells and that agonistic targeting of TLR1/TLR2 forces AML cells into apoptosis by p38 MAPK-dependent activation of Caspase 3, and differentiation by activating NFκB, thus revealing a new putative strategy for therapeutically targeting AML cells.

Identification of ETV6-RUNX1-like and DUX4-rearranged subtypes in paediatric B-cell precursor acute lymphoblastic leukaemia.

Fusion genes are potent driver mutations in cancer. In this study, we delineate the fusion gene landscape in a consecutive series of 195 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP ALL). Using RNA sequencing, we find in-frame fusion genes in 127 (65%) cases, including 27 novel fusions. We describe a subtype characterized by recurrent IGH-DUX4 or ERG-DUX4 fusions, representing 4% of cases, leading to overexpression of DUX4 and frequently co-occurring with intragenic ERG deletions. Furthermore, we identify a subtype characterized by an ETV6-RUNX1-like gene-expression profile and coexisting ETV6 and IKZF1 alterations. Thus, this study provides a detailed overview of fusion genes in paediatric BCP ALL and adds new pathogenetic insights, which may improve risk stratification and provide therapeutic options for this disease.

Transcriptome analysis of controlled and therapy-resistant childhood asthma reveals distinct gene expression profiles.

BACKGROUND: Children with problematic severe asthma have poor disease control despite high doses of inhaled corticosteroids and additional therapy, leading to personal suffering, early deterioration of lung function, and significant consumption of health care resources. If no exacerbating factors, such as smoking or allergies, are found after extensive investigation, these children are given a diagnosis of therapy-resistant (or therapy-refractory) asthma (SA). OBJECTIVE: We sought to deepen our understanding of childhood SA by analyzing gene expression and modeling the underlying regulatory transcription factor networks in peripheral blood leukocytes. METHODS: Gene expression was analyzed by using Cap Analysis of Gene Expression in children with SA (n = 13), children with controlled persistent asthma (n = 15), and age-matched healthy control subjects (n = 9). Cap Analysis of Gene Expression sequencing detects the transcription start sites of known and novel mRNAs and noncoding RNAs. RESULTS: Sample groups could be separated by hierarchical clustering on 1305 differentially expressed transcription start sites, including 816 known genes and several novel transcripts. Ten of 13 tested novel transcripts were validated by means of RT-PCR and Sanger sequencing. Expression of RAR-related orphan receptor A (RORA), which has been linked to asthma in genome-wide association studies, was significantly upregulated in patients with SA. Gene network modeling revealed decreased glucocorticoid receptor signaling and increased activity of the mitogen-activated protein kinase and Jun kinase cascades in patients with SA. CONCLUSION: Circulating leukocytes from children with controlled asthma and those with SA have distinct gene expression profiles, demonstrating the possible development of specific molecular biomarkers and supporting the need for novel therapeutic approaches.

Risk of childhood asthma is associated with CpG-site polymorphisms, regional DNA methylation and mRNA levels at the GSDMB/ORMDL3 locus.

Single-nucleotide polymorphisms (SNPs) in GSDMB (Gasdermin B) and ORMDL3 (ORMDL sphingolipid biosynthesis regulator 3) are strongly associated with childhood asthma, but the molecular alterations contributing to disease remain unknown. We investigated the effects of asthma-associated SNPs on DNA methylation and mRNA levels of GSDMB and ORMDL3. Genetic association between GSDMB/ORMDL3 and physician-diagnosed childhood asthma was confirmed in the Swedish birth-cohort BAMSE. CpG-site SNPs (rs7216389 and rs4065275) showed differences in DNA methylation depending on carrier status of the risk alleles, and were significantly associated with methylation levels in two CpG sites in the 5' UTR (untranslated region) of ORMDL3. In the Swedish Search study, we found significant differences in DNA methylation between asthmatics and controls in five CpG sites; after adjusting for lymphocyte and neutrophil cell counts, three remained significant: one in IKZF3 [IKAROS family zinc finger 3 (Aiolos); cg16293631] and two in the CpG island (CGI) of ORMDL3 (cg02305874 and cg16638648). Also, cg16293631 and cg02305874 correlated with mRNA levels of ORMDL3. The association between methylation and asthma was independent of the genotype in rs7216389, rs4065275 and rs12603332. Both SNPs and CpG sites showed significant associations with ORMDL3 mRNA levels. SNPs influenced expression independently of methylation, and the residual association between methylation and expression was not mediated by these SNPs. We found a differentially methylated region in the CGI shore of ORMDL3 with six CpG sites less methylated in CD8(+) T-cells. In summary, this study supports that there are differences in DNA methylation at this locus between asthmatics and controls; and both SNPs and CpG sites are independently associated with ORMDL3 expression.

Rule-based models of the interplay between genetic and environmental factors in childhood allergy.

Both genetic and environmental factors are important for the development of allergic diseases. However, a detailed understanding of how such factors act together is lacking. To elucidate the interplay between genetic and environmental factors in allergic diseases, we used a novel bioinformatics approach that combines feature selection and machine learning. In two materials, PARSIFAL (a European cross-sectional study of 3113 children) and BAMSE (a Swedish birth-cohort including 2033 children), genetic variants as well as environmental and lifestyle factors were evaluated for their contribution to allergic phenotypes. Monte Carlo feature selection and rule based models were used to identify and rank rules describing how combinations of genetic and environmental factors affect the risk of allergic diseases. Novel interactions between genes were suggested and replicated, such as between ORMDL3 and RORA, where certain genotype combinations gave odds ratios for current asthma of 2.1 (95% CI 1.2-3.6) and 3.2 (95% CI 2.0-5.0) in the BAMSE and PARSIFAL children, respectively. Several combinations of environmental factors appeared to be important for the development of allergic disease in children. For example, use of baby formula and antibiotics early in life was associated with an odds ratio of 7.4 (95% CI 4.5-12.0) of developing asthma. Furthermore, genetic variants together with environmental factors seemed to play a role for allergic diseases, such as the use of antibiotics early in life and COL29A1 variants for asthma, and farm living and NPSR1 variants for allergic eczema. Overall, combinations of environmental and life style factors appeared more frequently in the models than combinations solely involving genes. In conclusion, a new bioinformatics approach is described for analyzing complex data, including extensive genetic and environmental information. Interactions identified with this approach could provide useful hints for further in-depth studies of etiological mechanisms and may also strengthen the basis for risk assessment and prevention.

Transcriptome analysis reveals upregulation of bitter taste receptors in severe asthmatics.

The causes of severe childhood asthma are poorly understood. Our aim was to define global patterns of gene expression in children with severe therapy-resistant and controlled asthma. White blood cells were isolated and the global transcriptome profile was characterised using the Affymetrix Human Gene ST 1.0 chip in children with severe, therapy-resistant asthma (n=17), controlled asthma (n=19) and healthy controls (n=18). Receptor expression was studied in separated leukocyte fractions from asthmatic adults (n=12). Overall, 1378 genes were differentially expressed between children with severe/controlled asthma and controls. Three significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways were represented: natural killer cell-mediated cytotoxicity (upregulated in controlled asthma); N-glycan biosynthesis (downregulated in severe asthma); and bitter taste transduction receptors (TAS2Rs) (upregulated in severe asthma). Quantitative PCR experiments confirmed upregulation of TAS2Rs in severe asthmatics. TAS2R expression was replicated in leukocytes from adult asthmatics, in which TAS2R agonists also inhibited LPS-induced cytokine release. Significant correlations between expression of TAS2Rs and clinical markers of asthma severity were found in both adults and children. In conclusion, specific gene expression patterns were observed in children with severe, therapy-resistant asthma. The increased expression of bronchodilatory TAS2Rs suggests a new target for the treatment of asthma.

Assessment of the neuropeptide S system in anxiety disorders.

BACKGROUND: The G protein-coupled receptor neuropeptide S receptor 1 (NPSR1) and its ligand neuropeptide S (NPS) form a signaling system mainly implicated in susceptibility to asthma and inflammatory disorders in humans and regulation of anxiety and arousal in rodents. We addressed here the role of NPS and NPSR1 as susceptibility genes for human anxiety disorders. METHODS: We performed comprehensive association analysis of genetic variants in NPS and NPSR1 in three independent study samples. We first studied a population-based sample (Health 2000, Finland) of 321 anxiety disorder patients and 1317 control subjects and subsequently a Spanish clinical panic disorder sample consisting of 188 cases and 315 control subjects. In addition, we examined a birth cohort of 2020 children (Barn Allergi Miljö Stockholm Epidemiologi [BAMSE], Sweden). We then tested whether alleles of the most significantly associated single nucleotide polymorphisms alter DNA-protein complex formation in electrophoretic mobility shift assays. Finally, we compared acute stress responses on the gene expression level in wild-type and Npsr1(-/-) mice. RESULTS: We confirmed previously observed epidemiological association between anxiety and asthma in two population-based cohorts. Single nucleotide polymorphisms within NPS and NPSR1 associated with panic disorder diagnosis in the Finnish and Spanish samples and with parent-reported anxiety/depression in the BAMSE sample. Moreover, some of the implicated single nucleotide polymorphisms potentially affect transcription factor binding. Expression of neurotrophin-3, a neurotrophic factor connected to stress and panic reaction, was significantly downregulated in brain regions of stressed Npsr1(-/-) mice, whereas interleukin-1 beta, an active stress-related immunotransmitter, was upregulated. CONCLUSIONS: Our results suggest that NPS-NPSR1 signaling is likely involved in anxiety.

Evidence for genetic association and interaction between the TYK2 and IRF5 genes in systemic lupus erythematosus.

OBJECTIVE: Several candidate genes have been implicated in susceptibility for systemic lupus erythematosus (SLE), a complex autoimmune disease. The proposed genes include members of the type I interferon (IFN) pathway and genes involved in immunological defense functions. Our aim was to systematically replicate 6 such genes, TYK2, IRF5, CTLA4, PDCD1, FCGR2A, and NOD2. METHODS: Single-nucleotide polymorphisms in TYK2, IRF5, CTLA4, PDCD1, FCGR2A, and NOD2 were genotyped in 277 SLE patients and 356 healthy controls from Finland, giving a power of 42%-70% for different genes at published allele frequencies. RESULTS: Significant association was seen for rs2304256 (p = 0.0001) and rs12720270 (p = 0.0031) in TYK2 and rs10954213 (p = 0.0043) in IRF5 in our samples, but not for the other genes. We found evidence for genetic interaction (p = 0.014) between rs2304256 in TYK2 and rs10954213 in IRF5, both members of the type I IFN pathway, strengthening the role of the type I IFN pathway in the pathogenesis of SLE. CONCLUSION: The IFN pathway genes IRF5 and TYK2 may act epistatically in increasing risk for SLE, but our lack of replication does not exclude effects of the other genes studied.

Biological and genetic interaction between tenascin C and neuropeptide S receptor 1 in allergic diseases.

Neuropeptide S receptor 1 (NPSR1, GPRA 154, GPRA) has been verified as a susceptibility gene for asthma and related phenotypes. The ligand for NPSR1, Neuropeptide S (NPS), activates signalling through NPSR1 and microarray analysis has identified Tenascin C (TNC) as a target gene of NPS-NPSR1 signalling. TNC has previously been implicated as a risk gene for asthma. We aimed therefore to study the genetic association of TNC in asthma- and allergy-related disorders as well as the biological and genetic interactions between NPSR1 and TNC. Regulation of TNC was investigated using NPS stimulated NPSR1 transfected cells. We genotyped 12 TNC SNPs in the cross-sectional PARSIFAL study (3113 children) and performed single SNP association, haplotype association and TNC and NPSR1 gene-gene interaction analyses. Our experimental results show NPS-dependent upregulation of TNC-mRNA. The genotyping results indicate single SNP and haplotype associations for several SNPs in TNC with the most significant association to rhinoconjunctivitis for a haplotype, with a frequency of 29% in cases (P = 0.0005). In asthma and atopic sensitization significant gene-gene interactions were found between TNC and NPSR1 SNPs, indicating that depending on the NPSR1 genotype, TNC can be associated with either an increased or a decreased risk of disease. We conclude that variations in TNC modifies, not only risk for asthma, but also for rhinoconjunctivitis. Furthermore, we show epistasis based on both a direct suggested regulatory effect and a genetic interaction between NPSR1 and TNC. These results suggest merging of previously independent pathways of importance in the development of asthma- and allergy-related traits.