Development and validation of a pharmacogenomics reporting workflow based on the illumina global screening array chip

Background: Microarrays are a well-established and widely adopted technology capable of interrogating hundreds of thousands of loci across the human genome. Combined with imputation to cover common variants not included in the chip design, they offer a cost-effective solution for large-scale genetic studies. Beyond research applications, this technology can be applied for testing pharmacogenomics, nutrigenetics, and complex disease risk prediction. However, establishing clinical reporting workflows requires a thorough evaluation of the assay's performance, which is achieved through validation studies. In this study, we performed pre-clinical validation of a genetic testing workflow based on the Illumina Global Screening Array for 25 pharmacogenomic-related genes. Methods: To evaluate the accuracy of our workflow, we conducted multiple pre-clinical validation studies. Here, we present the results of accuracy and precision assessments, involving a total of 73 cell lines. These assessments encompass reference materials from the Genome-In-A-Bottle (GIAB), the Genetic Testing Reference Material Coordination Program (GeT-RM) projects, as well as additional samples from the 1000 Genomes project (1KGP). We conducted an accuracy assessment of genotype calls for target loci in each indication against established truth sets. Results: In our per-sample analysis, we observed a mean analytical sensitivity of 99.39% and specificity 99.98%. We further assessed the accuracy of star-allele calls by relying on established diplotypes in the GeT-RM catalogue or calls made based on 1KGP genotyping. On average, we detected a diplotype concordance rate of 96.47% across 14 pharmacogenomic-related genes with star allele-calls. Lastly, we evaluated the reproducibility of our findings across replicates and observed 99.48% diplotype and 100% phenotype inter-run concordance. Conclusion: Our comprehensive validation study demonstrates the robustness and reliability of the developed workflow, supporting its readiness for further development for applied testing.

Allopurinol-Induced Stevens–Johnson Syndrome in Javanese Men With Positive HLA‐B*58:01

Background: Allopurinol is the most commonly used drug for the treatment of gout arthritis. However, the use of allopurinol is associated with severe cutaneous adverse reactions (SCARs) and life-threatening immune-mediated reactions that include Stevens–Johnson syndrome (SJS). SJS induced by allopurinol is strongly linked with the presence of HLA-B*58:01 in the Asian population. Such a study has not been conducted in Indonesia. We present two cases with clinical diagnosis of SJS. These patients had Javanese ethnicity, for which evidence on the genetic predisposition of allopurinol-induced SJS/TEN had not been established. Testing for the presence of the HLA-B∗58:01 allele was positive in both cases. Our case report confirms findings from studies in Asian countries that link HLA-B*58:01 and allopurinol-induced SJS/TEN. A larger study is needed to elicit evidence that the HLA-B*58:01 allele can potentially be used as a genetic marker for allopurinol-induced SCARs among different ethnicities in Indonesia.

Implementation of genetic screening test to reduce the incidence of dapsone hypersensitivity syndrome among patients with leprosy in Papua, Indonesia: a study protocol

INTRODUCTION

The mainstay of leprosy treatment is multidrug treatment (MDT), which contains rifampicin, dapsone and clofazimine. The occurrence of dapsone hypersensitivity syndrome (DHS), a sudden, potentially fatal and traumatic adverse reaction due to dapsone, may affect treatment adherence and may result in fatality if untreated. Before MDT administration, screening for HLA-B*13:01 in patients with leprosy can potentially reduce DHS risk. The study aims to assess the effectiveness of using a screening test for HLA-B*13:01 in reducing the incidence of DHS and to evaluate the feasibility of using the quantitative PCR-based screening tool as DHS predictors before dapsone administration using individual patient testing in a referral centralised-lab model.

METHODS AND ANALYSIS

A total of 310 newly diagnosed patients with leprosy will be recruited from health centres in two highly endemic districts in Indonesia. Dried blood will be taken on filter paper as the specimen receptacle to collect DNA from the patients and transported at room temperature to the leprosy referral laboratory before MDT administration. Checking for HLA-B*13:01 from human DNA is performed using the Nala PGx 1301 V.1 kit. The results will be shared with the leprosy health workers on the site via phone call and courier. Patients with a positive test result will be treated with MDT without dapsone, and patients with a negative result will be treated with complete MDT. Physical examination (weight, height, skin, muscle and nerve function examination), complete blood tests (including renal function test) will be carried out at baseline. Follow-up will be performed at the fourth and eighth weeks to observe any development of adverse drug reactions.

ETHICS AND DISSEMINATION

The ethical approval for the study was issued by the Ethical Committee of the National Institute of Health Research and Development, Ministry of Health, Indonesia. Written informed consent will be sought from all participants.

Variant landscape of the RYR1 gene based on whole genome sequencing of the Singaporean population

The RYR1 gene codes for a ryanodine receptor which is a calcium release channel in the skeletal muscle sarcoplasmic reticulum. It is associated with Malignant Hyperthermia (MH) and congenital myopathies including Central Core Disease (CCD), Multiminicore Disease (MMD) and Congenital Fibre-Type Disproportion (CFTD). There is currently little information on the epidemiology of RYR1 variants in Asians. Our study aims to describe the RYR1 variant landscape in a Singapore cohort unselected for RYR1-associated conditions. Data was retrieved from the SG10K pilot project, where whole genome sequencing was performed on volunteers unselected and undetermined for RYR1-associated conditions. Variants were classified based on pathogenicity using databases ClinVar and InterVar. Allele frequencies of pathogenic variants were compared between Chinese, Indians and Malays. Using databases ExAC, GnomAD and GenomeAsia 100k study, we further compared local allele frequencies to those in Europe, America and Asia. Data was analysed using R Commander. Significant P value was set at p < 0.05. Majority of the RYR1 variants were missense mutations. We identified four pathogenic and four likely pathogenic RYR1 variants. All were related to the aforementioned RYR1-associated conditions. There were 6 carriers of RYR1 pathogenic variants amongst 4810 individuals, corresponding to an allele frequency of 0.06%. The prevalence of pathogenic variants was the highest amongst Indians (4 in 1127 individuals) (p = 0.030). Majority of pathogenic and likely pathogenic mutations were missense and located in mutational hotspots. These variants also occurred at higher frequencies in Asians than globally. This study describes the variant landscape of the RYR1 gene in Singapore. This knowledge will facilitate genetic screening for RYR1-related conditions.

Validation of a multi-gene qPCR-based pharmacogenomics panel across major ethnic groups in Singapore and Indonesia

Aim: The clinical utility of pharmacogenomics (PGx) has been gaining traction alongside growing evidence that adverse drug reactions (ADRs) have significant genetic associations. Nala PGx Core^® is a multi-gene qPCR-based panel of 20 allele variants, comprising 18 SNPs and two CYP2D6 copy number markers across four pharmacogenes - CYP2C9, CYP2C19, CYP2D6 and SLCO1B1. Methods: In this study, we validated the performance of Nala PGx Core^® against benchmark methods, on the Singaporean and Indonesian populations. Results & conclusion: Nala PGx Core^® demonstrated robust and accurate genotyping when compared with other established benchmarks. Furthermore, the panel successfully characterized alleles of clinical relevance, such as CYP2D6*10 and CYP2D6*36, across major ethnic groups present of Singapore and Indonesia, suggesting its potential for adoption in clinical workflows regionally.

Wide Application of Minimally Processed Saliva on Multiple RT-qPCR Kits for SARS-CoV-2 Detection in Indonesia

Saliva as a sample matrix has been an attractive alternative for the detection of SARS-CoV-2. However, due to potential variability in collection and processing steps, evaluating a proposed workflow amongst the local population is recommended. Here, we aim to validate the collection and treatment of human saliva as a direct specimen for RT-qPCR-based detection of SARS-CoV-2 in Indonesia. We demonstrated that SARS-CoV-2 target genes were detected in saliva specimens and remained stable for five days either refrigerated or stored at room temperature. The method of processing saliva specimens described in this report bypasses the need for an RNA-extraction process, thereby reducing the cost, time, and manpower required for processing samples. The developed method was tested across nine commercial kits, including the benchmark, to demonstrate its wide applicability on multiple existing workflows. Our developed method achieved an 86% overall agreement rate compared to paired nasopharyngeal and oropharyngeal swab specimens (NPOP). With the assistance of a saliva sampling device, the collection was found to be more convenient for individuals and improved the overall agreement rate to 97%.

Pharmacogenomics Implementation Training Improves Self-Efficacy and Competency to Drive Adoption in Clinical Practice

Background: Administration of pharmacogenomics (PGx) testing in clinical practice has been suboptimal, presumably due to lack of PGx education. Here, we aim to evaluate the standpoint of PGx testing among a diverse group of healthcare professionals (HCPs) through conducting surveys before and after training. Materials and Methods: Training modules were designed to cover three key learning objectives and deployed in five sections. A pre- and post-training survey questionnaire was used to evaluate participants' self-assessments on employing PGx in clinical practice. Results and Conclusion: Out of all enrollments, 102 survey responses were collected. Overall, respondents agree on the benefits of PGx testing, but have inadequate self-efficacy and competency in utilizing PGx data. Our results show that a 90 min long training significantly improves these, and could lead to greater anticipation of PGx adoption.

Validation study of HLA-B*13:01 as a biomarker of dapsone hypersensitivity syndrome in leprosy patients in Indonesia

Leprosy is a stigmatizing, chronic infection which degenerates the nervous system and often leads to incapacitation. Multi-drug therapy which consists of dapsone, rifampicin and clofazimine has been effective to combat this disease. In Indonesia, especially in Papua Island, leprosy is still a problem. Furthermore, there had been higher reports of Dapsone Hypersensitivity Syndrome (DHS) which also challenges leprosy elimination in certain aspects. Globally, DHS has a prevalence rate of 1.4% and a fatality rate up to 13%. The aim of this study is to validate HLA-B*13:01, a previously discovered biomarker for DHS in the Chinese population, as a biomarker for DHS in the Papua population.This is a case-control study of 34 leprosy patients who presented themselves with DHS (case subjects) and 52 leprosy patients without DHS (control subjects). Patients were recruited from 2 provinces: Papua and West Papua. DNA was extracted from 3 ml blood specimens. HLA-B alleles were typed using the gold-standard sequence based typing method. Results were then analysed using logistic regression and risk assessment was carried out. The results of HLA-typing showed that HLA-B*13:01 was the most significant allele associated with DHS, with odds ratio = 233.64 and P-value = 7.11×10⁻⁹, confirming the strong association of HLA-B*13:01 to DHS in the Papua population. The sensitivity of this biomarker is 91.2% and specificity is 96.2%, with an area under the curve of 0.95. HLA-B*13:01 is validated as a biomarker for DHS in leprosy patients in Papua, Indonesia, and can potentially be a good predictor of DHS to help prevent this condition in the future.

Leprosy as a chronic infectious disease that affects the skin and nervous system is treated with a treatment cocktail, including rifampicin, clofazimine and dapsone. Unfortunately, one of these drugs, namely dapsone, may cause the patient to exhibit adverse symptoms that appear as skin hypersensitivity and may potentially lead to death up to 9.9% of the time. In a previous study, it has been established that this adverse drug reaction, which is called dapsone hypersensitivity syndrome (DHS), is associated with an allele in the human leukocyte antigen, HLA-B*13:01. In the light of validating the association of HLA-B*13:01 with DHS in a leprosy endemic area like Indonesia, we conducted a study involving leprosy patients who had DHS during the course of their multi-drug treatment as well as leprosy patients who managed to complete their course of treatment without exhibiting DHS. The results of this study validated the association of HLA-B*13:01 and DHS at a very significant level of evidence and the odds of people carrying at least one allele is 233.3-times the risk of non-carriers. This allele was also able to predict a person at risk of DHS 95% of the time. Thus, screening of HLA-B*13:01 to prevent DHS is warranted.

Analysis of 47 Non-MHC Ankylosing Spondylitis Susceptibility Loci Regarding Associated Variants across Whites and Han Chinese

OBJECTIVE

To present a systematic evaluation of 47 non-MHC ankylosing spondylitis (AS) susceptibility loci that have been initially discovered through white genome-wide association studies in Han Chinese.

METHODS

Originally, 10,743 samples representing north and south Chinese in 4 datasets were obtained. After data quality control and imputation, metaanalysis results of 94,621 variants within 47 loci were extracted. Four ERAP1 single-nucleotide polymorphisms (SNP) and HLA-B27 tag SNP rs13202464 were used for interaction analysis. Population-attributable risk percentages of AS-associated variants were compared. Functional annotations of AS-associated variants were conducted using HaploReg, RegulomeDB, and rVarBase databases.

RESULTS

We revealed 16 AS-associated variants with nominal evidence in Han Chinese, including rs10865331 (p = 6.30 × 10^-10), rs10050860 (p = 4.09 × 10^-5) and rs8070463 (p = 1.03 × 10^-4). Potential susceptible SNP within these 47 loci were also identified, such as rs13024541 (2p15), rs17401719 (5q15), and rs62074054 (17q21). Epistatic interactions between 3 ERAP1 SNP (rs17401719, rs30187, and rs10050860) and HLA-B27 were confirmed. Among the 16 AS-associated variants, rs30187 showed weaker risk effect, while rs10050860 and rs12504282 seemed to attribute more risk in Han Chinese than in whites. Further genomic annotation pinpointed 35 candidate functional SNP, especially in the 2p15, ERAP1, and NPEPPS-TBKBP1 regions.

CONCLUSION

Our results provided a detailed spectrum of all the reported non-MHC AS susceptibility loci in Han Chinese, which comprehensively exhibited the ethnic heterogeneity of AS susceptibility and highlighted that 2p15, ERAP1, and NPEPPS-TBKBP1 regions may play a critical role in AS pathogenesis across diverse populations.

Evaluation of Prospective HLA-B*13:01 Screening to Prevent Dapsone Hypersensitivity Syndrome in Patients With Leprosy

Importance

Dapsone hypersensitivity syndrome (DHS) is the most serious adverse reaction associated with dapsone administration and one of the major causes of death in patients with leprosy, whose standard treatment includes multidrug therapy (MDT) with dapsone, rifampicin, and clofazimine. Although the HLA-B*13:01 polymorphism has been identified as the genetic determinant of DHS in the Chinese population, no studies to date have been done to evaluate whether prospective HLA-B*13:01 screening could prevent DHS by identifying patients who should not receive dapsone.

Objective

To evaluate the clinical use of prospective HLA-B*13:01 screening for reduction of the incidence of DHS by excluding dapsone from the treatment for patients with HLA-B*13:01-positive leprosy.

Design, Setting, and Participants

A prospective cohort study was conducted from February 15, 2015, to April 30, 2018, in 21 provinces throughout China. A total of 1539 patients with newly diagnosed leprosy were enrolled who had not received dapsone previously. After excluding patients who had a history of allergy to sulfones or glucose-6-phosphate dehydrogenase deficiency, 1512 individuals underwent HLA-B*13:01 genotyping. All of the patients were followed up weekly for the first 8 weeks after treatment to monitor for adverse events.

Exposures

Patients who were HLA-B*13:01 carriers were instructed to eliminate dapsone from their treatment regimens, and noncarrier patients received standard MDT.

Main Outcomes and Measures

The primary outcome was the incidence of DHS. The historical incidence rate of DHS (1.0%) was used as a control.

Results

Among 1512 patients (1026 [67.9%] men, 486 [32.1%] women; mean [SD] age, 43.1 [16.2] years), 261 (17.3%) were identified as carriers of the HLA-B*13:01 allele. A total of 714 adverse events in 384 patients were observed during the follow-up period. Dapsone hypersensitivity syndrome did not develop in any of the 1251 patients who were HLA-B*13:01-negative who received dapsone, while approximately 13 patients would be expected to experience DHS, based on the historical incidence rate of 1.0% per year (P = 2.05 × 10-5). No significant correlation was found between other adverse events, including dermatologic or other events, and HLA-B*13:01 status.

Conclusions and Relevance

Prospective HLA-B*13:01 screening and subsequent elimination of dapsone from MDT for patients with HLA-B*13:01-positive leprosy may significantly reduce the incidence of DHS in the Chinese population.

Investigation of the predisposing factor of pemphigus and its clinical subtype through a genome-wide association and next generation sequence analysis

BACKGROUND

Pemphigus is an autoimmune blistering disease with pemphigus vulgaris (PV) and foliaceus (PF) as the two major histological subtypes. Associations with HLA molecules have been suggested, but specific HLA risk variants as well as non-HLA risk variants remain to be discovered.

METHODS

We performed a two-stage genome-wide association study in the Chinese Han population through a genome-wide discovery analysis and follow-up validation analysis in a total number of 210 PV, 159 PF and 2493 healthy controls. HLA imputation as well as high coverage next generation sequencing based HLA genotyping was employed to investigate the association of classical HLA alleles and amino acid change.

RESULTS

We have discovered independent novel associations with PF at rs2178077 on 12q24.33, located next to RAN (P_PF = 1.57 × 10^-9 ) and rs3888722 within the MHC region (P = 6.73 × 10^-9 ). For the HLA variants, we confirmed independent genome-wide level risk associations in HLA-DQB1 and HLA-DRB1, with DQB1*05:03 to be the strongest association with PV (P = 8.59 × 10^-68 , OR = 31.16) and PF (P = 4.84 × 10^-17 , OR = 5.64). In addition, DRB1*14 was demonstrated to be a second independent variants (P = 4.2 × 10^-63 , OR = 35.47) for PV, while DRB1*04:06 was demonstrated to be the second independent signal (P = 7.44 × 10^-13 , OR = 5.58) for PF.

CONCLUSIONS

These findings advance our understanding of the genetic basis of pemphigus susceptibility and may offer opportunities for risk prediction and preventive treatment for pemphigus, in particular for PV.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy

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.

A functional SNP associated with atopic dermatitis controls cell type-specific methylation of the VSTM1 gene locus

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.

A large-scale genome-wide association and meta-analysis identified four novel susceptibility loci for leprosy

Leprosy, a chronic infectious disease, results from the uncultivable pathogen Mycobacterium leprae (M. leprae), and usually progresses to peripheral neuropathy and permanent progressive deformity if not treated. Previously published genetic studies have identified 18 gene/loci significantly associated with leprosy at the genome-wide significant level. However as a complex disease, only a small proportion of leprosy risk could be explained by those gene/loci. To further identify more susceptibility gene/loci, we hereby performed a three-stage GWAS comprising 8,156 leprosy patients and 15,610 controls of Chinese ancestry. Four novel loci were identified including rs6807915 on 3p25.2 (P=1.94 × 10⁻⁸, OR=0.89), rs4720118 on 7p14.3 (P=3.85 × 10⁻¹⁰, OR=1.16), rs55894533 on 8p23.1 (P=5.07 × 10⁻¹¹, OR=1.15) and rs10100465 on 8q24.11 (P=2.85 × 10⁻¹¹, OR=0.85). Altogether, these findings have provided new insight and significantly expanded our understanding of the genetic basis of leprosy.

Previous studies have shown genetic associations between leprosy and 18 different genes/loci. Here, Wang and colleagues perform genome-wide association study in Han Chinese leprosy patients and describe four novel loci to be associated to the disease.

Cross-Cancer Genome-Wide Analysis of Lung, Ovary, Breast, Prostate, and Colorectal Cancer Reveals Novel Pleiotropic Associations

Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103-14. ©2016 AACR.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

Genome-wide analysis of the genetic regulation of gene expression in human neutrophils

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.

A functional brain-derived neurotrophic factor (BDNF) gene variant increases the risk of moderate-to-severe allergic rhinitis

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.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

MicroRNA related polymorphisms and breast cancer risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Breast cancer risk assessment using genetic variants and risk factors in a Singapore Chinese population

Introduction

Genetic variants for breast cancer risk identified in genome-wide association studies (GWAS) in Western populations require further testing in Asian populations. A risk assessment model incorporating both validated genetic variants and established risk factors may improve its performance in risk prediction of Asian women.

Methods

A nested case-control study of female breast cancer (411 cases and 1,212 controls) within the Singapore Chinese Health Study was conducted to investigate the effects of 51 genetic variants identified in previous GWAS on breast cancer risk. The independent effect of these genetic variants was assessed by creating a summed genetic risk score (GRS) after adjustment for body mass index and the Gail model risk factors for breast cancer.

Results

The GRS was an independent predictor of breast cancer risk in Chinese women. The multivariate-adjusted odds ratios (95% confidence intervals) of breast cancer for the second, third, and fourth quartiles of the GRS were 1.26 (0.90 to 1.76), 1.47 (1.06 to 2.04) and 1.75 (1.27 to 2.41) respectively (P for trend <0.001). In addition to established risk factors, the GRS improved the classification of 6.2% of women for their absolute risk of breast cancer in the next five years.

Conclusions

Genetic variants on top of conventional risk factors can improve the risk prediction of breast cancer in Chinese women.

Genome-wide association study identifies 25 known breast cancer susceptibility loci as risk factors for triple-negative breast cancer

Triple-negative (TN) breast cancer is an aggressive subtype of breast cancer associated with a unique set of epidemiologic and genetic risk factors. We conducted a two-stage genome-wide association study of TN breast cancer (stage 1: 1529 TN cases, 3399 controls; stage 2: 2148 cases, 1309 controls) to identify loci that influence TN breast cancer risk. Variants in the 19p13.1 and PTHLH loci showed genome-wide significant associations (P < 5 × 10(-) (8)) in stage 1 and 2 combined. Results also suggested a substantial enrichment of significantly associated variants among the single nucleotide polymorphisms (SNPs) analyzed in stage 2. Variants from 25 of 74 known breast cancer susceptibility loci were also associated with risk of TN breast cancer (P < 0.05). Associations with TN breast cancer were confirmed for 10 loci (LGR6, MDM4, CASP8, 2q35, 2p24.1, TERT-rs10069690, ESR1, TOX3, 19p13.1, RALY), and we identified associations with TN breast cancer for 15 additional breast cancer loci (P < 0.05: PEX14, 2q24.1, 2q31.1, ADAM29, EBF1, TCF7L2, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, BRCA2, RAD51L1-rs2588809, MKL1). Further, two SNPs independent of previously reported signals in ESR1 [rs12525163 odds ratio (OR) = 1.15, P = 4.9 × 10(-) (4)] and 19p13.1 (rs1864112 OR = 0.84, P = 1.8 × 10(-) (9)) were associated with TN breast cancer. A polygenic risk score (PRS) for TN breast cancer based on known breast cancer risk variants showed a 4-fold difference in risk between the highest and lowest PRS quintiles (OR = 4.03, 95% confidence interval 3.46-4.70, P = 4.8 × 10(-) (69)). This translates to an absolute risk for TN breast cancer ranging from 0.8% to 3.4%, suggesting that genetic variation may be used for TN breast cancer risk prediction.

A genome-wide association study of early-onset breast cancer identifies PFKM as a novel breast cancer gene and supports a common genetic spectrum for breast cancer at any age

Early-onset breast cancer (EOBC) causes substantial loss of life and productivity, creating a major burden among women worldwide. We analyzed 1,265,548 Hapmap3 single-nucleotide polymorphisms (SNP) among a discovery set of 3,523 EOBC incident cases and 2,702 population control women ages ≤ 51 years. The SNPs with smallest P values were examined in a replication set of 3,470 EOBC cases and 5,475 control women. We also tested EOBC association with 19,684 genes by annotating each gene with putative functional SNPs, and then combining their P values to obtain a gene-based P value. We examined the gene with smallest P value for replication in 1,145 breast cancer cases and 1,142 control women. The combined discovery and replication sets identified 72 new SNPs associated with EOBC (P < 4 × 10(-8)) located in six genomic regions previously reported to contain SNPs associated largely with later-onset breast cancer (LOBC). SNP rs2229882 and 10 other SNPs on chromosome 5q11.2 remained associated (P < 6 × 10(-4)) after adjustment for the strongest published SNPs in the region. Thirty-two of the 82 currently known LOBC SNPs were associated with EOBC (P < 0.05). Low power is likely responsible for the remaining 50 unassociated known LOBC SNPs. The gene-based analysis identified an association between breast cancer and the phosphofructokinase-muscle (PFKM) gene on chromosome 12q13.11 that met the genome-wide gene-based threshold of 2.5 × 10(-6). In conclusion, EOBC and LOBC seem to have similar genetic etiologies; the 5q11.2 region may contain multiple distinct breast cancer loci; and the PFKM gene region is worthy of further investigation. These findings should enhance our understanding of the etiology of breast cancer.

HLA-B*13:01 and the dapsone hypersensitivity syndrome

BACKGROUND

Dapsone is used in the treatment of infections and inflammatory diseases. The dapsone hypersensitivity syndrome, which is associated with a reported mortality of 9.9%, develops in about 0.5 to 3.6% of persons treated with the drug. Currently, no tests are available to predict the risk of the dapsone hypersensitivity syndrome.

METHODS

We performed a genomewide association study involving 872 participants who had received dapsone as part of multidrug therapy for leprosy (39 participants with the dapsone hypersensitivity syndrome and 833 controls), using log-additive tests of single-nucleotide polymorphisms (SNPs) and imputed HLA molecules. For a replication analysis, we genotyped 24 SNPs in an additional 31 participants with the dapsone hypersensitivity syndrome and 1089 controls and performed next-generation sequencing for HLA-B and HLA-C typing at four-digit resolution in an independent series of 37 participants with the dapsone hypersensitivity syndrome and 201 controls.

RESULTS

Genomewide association analysis showed that SNP rs2844573, located between the HLA-B and MICA loci, was significantly associated with the dapsone hypersensitivity syndrome among patients with leprosy (odds ratio, 6.18; P=3.84×10(-13)). HLA-B*13:01 was confirmed to be a risk factor for the dapsone hypersensitivity syndrome (odds ratio, 20.53; P=6.84×10(-25)). The presence of HLA-B*13:01 had a sensitivity of 85.5% and a specificity of 85.7% as a predictor of the dapsone hypersensitivity syndrome, and its absence was associated with a reduction in risk by a factor of 7 (from 1.4% to 0.2%). HLA-B*13:01 is present in about 2 to 20% of Chinese persons, 1.5% of Japanese persons, 1 to 12% of Indians, and 2 to 4% of Southeast Asians but is largely absent in Europeans and Africans.

CONCLUSIONS

HLA-B*13:01 was associated with the development of the dapsone hypersensitivity syndrome among patients with leprosy. (Funded by the National Natural Science Foundation of China and others.).

An association study of TOLL and CARD with leprosy susceptibility in Chinese population

Previous genome-wide association studies (GWASs) identified multiple susceptibility loci that have highlighted the important role of TLR (Toll-like receptor) and CARD (caspase recruitment domain) genes in leprosy. A large three-stage candidate gene-based association study of 30 TLR and 47 CARD genes was performed in the leprosy samples of Chinese Han. Of 4363 SNPs investigated, eight SNPs showed suggestive association (P < 0.01) in our previously published GWAS datasets (Stage 1). Of the eight SNPs, rs2735591 and rs4889841 showed significant association (P < 0.001) in an independent series of 1504 cases and 1502 controls (Stage 2), but only rs2735591 (next to BCL10) showed significant association in the second independent series of 938 cases and 5827 controls (Stage 3). Rs2735591 showed consistent association across the three stages (P > 0.05 for heterogeneity test), significant association in the combined validation samples (P_corrected = 5.54 × 10⁻⁴ after correction for 4363 SNPs tested) and genome-wide significance in the whole GWAS and validation samples (P = 1.03 × 10⁻⁹, OR = 1.24). In addition, we demonstrated the lower expression of BCL10 in leprosy lesions than normal skins and a significant gene connection between BCL10 and the eight previously identified leprosy loci that are associated with NFκB, a major regulator of downstream inflammatory responses, which provides further biological evidence for the association. We have discovered a novel susceptibility locus on 1p22, which implicates BCL10 as a new susceptibility gene for leprosy. Our finding highlights the important role of both innate and adaptive immune responses in leprosy.

Germline variation in TP53 regulatory network genes associates with breast cancer survival and treatment outcome

Germline variation in the TP53 network genes PRKAG2, PPP2R2B, CCNG1, PIAS1 and YWHAQ was previously suggested to have an impact on drug response in vitro. Here, we investigated the effect on breast cancer survival of germline variation in these genes in 925 Finnish breast cancer patients and further analyzed five single nucleotide polymorphisms (SNPs) in PRKAG2 (rs1029946, rs4726050, rs6464153, rs7789699) and PPP2R2B (rs10477313) for 10-year survival in breast cancer patients, interaction with TP53 R72P and MDM2-SNP309, outcome after specific adjuvant therapy and correlation to tumor characteristics in 4,701 invasive cases from four data sets. We found evidence for carriers of PRKAG2-rs1029946 and PRKAG2-rs4726050 having improved survival in the pooled data (HR 0.53, 95% CI 0.3-0.9; p = 0.023 for homozygous carriers of the rare G-allele and HR 0.85, 95% CI 0.7-0.9; p = 0.049 for carriers of the rare G allele, respectively). PRKAG2-rs4726050 showed a significant interaction with MDM2-SNP309, with PRKAG2-rs4726050 rare G-allele having a dose-dependent effect for better breast cancer survival confined only to MDM2 SNP309 rare G-allele carriers (HR 0.45, 95% CI 0.2-0.7; p = 0.001). This interaction also emerged as an independent predictor of better survival (p = 0.047). PPP2R2B-rs10477313 rare A-allele was found to predict better survival (HR 0.82, 95% CI 0.6-0.9; p = 0.018), especially after hormonal therapy (HR 0.66, 95% CI 0.5-0.9; p = 0.048). These findings warrant further studies and suggest that genetic markers in TP53 network genes such as PRKAG2 and PPP2R2B might affect prognosis and treatment outcome in breast cancer patients.

Large-scale genotyping identifies 41 new loci associated with breast cancer risk

Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ∼9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10(-8)). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.

Genome-wide association studies identify four ER negative-specific breast cancer risk loci

Estrogen receptor (ER)-negative tumors represent 20-30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry. The etiology and clinical behavior of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10(-12) and LGR6, P = 1.4 × 10(-8)), 2p24.1 (P = 4.6 × 10(-8)) and 16q12.2 (FTO, P = 4.0 × 10(-8)), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.

A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11

Genome-wide association studies (GWAS) of breast cancer defined by hormone receptor status have revealed loci contributing to susceptibility of estrogen receptor (ER)-negative subtypes. To identify additional genetic variants for ER-negative breast cancer, we conducted the largest meta-analysis of ER-negative disease to date, comprising 4754 ER-negative cases and 31 663 controls from three GWAS: NCI Breast and Prostate Cancer Cohort Consortium (BPC3) (2188 ER-negative cases; 25 519 controls of European ancestry), Triple Negative Breast Cancer Consortium (TNBCC) (1562 triple negative cases; 3399 controls of European ancestry) and African American Breast Cancer Consortium (AABC) (1004 ER-negative cases; 2745 controls). We performed in silico replication of 86 SNPs at P ≤ 1 × 10(-5) in an additional 11 209 breast cancer cases (946 with ER-negative disease) and 16 057 controls of Japanese, Latino and European ancestry. We identified two novel loci for breast cancer at 20q11 and 6q14. SNP rs2284378 at 20q11 was associated with ER-negative breast cancer (combined two-stage OR = 1.16; P = 1.1 × 10(-8)) but showed a weaker association with overall breast cancer (OR = 1.08, P = 1.3 × 10(-6)) based on 17 869 cases and 43 745 controls and no association with ER-positive disease (OR = 1.01, P = 0.67) based on 9965 cases and 22 902 controls. Similarly, rs17530068 at 6q14 was associated with breast cancer (OR = 1.12; P = 1.1 × 10(-9)), and with both ER-positive (OR = 1.09; P = 1.5 × 10(-5)) and ER-negative (OR = 1.16, P = 2.5 × 10(-7)) disease. We also confirmed three known loci associated with ER-negative (19p13) and both ER-negative and ER-positive breast cancer (6q25 and 12p11). Our results highlight the value of large-scale collaborative studies to identify novel breast cancer risk loci.

Genome-wide association analysis identifies three new breast cancer susceptibility loci

Breast cancer is the most common cancer among women. To date, 22 common breast cancer susceptibility loci have been identified accounting for ∼8% of the heritability of the disease. We attempted to replicate 72 promising associations from two independent genome-wide association studies (GWAS) in ∼70,000 cases and ∼68,000 controls from 41 case-control studies and 9 breast cancer GWAS. We identified three new breast cancer risk loci at 12p11 (rs10771399; P = 2.7 × 10(-35)), 12q24 (rs1292011; P = 4.3 × 10(-19)) and 21q21 (rs2823093; P = 1.1 × 10(-12)). rs10771399 was associated with similar relative risks for both estrogen receptor (ER)-negative and ER-positive breast cancer, whereas the other two loci were associated only with ER-positive disease. Two of the loci lie in regions that contain strong plausible candidate genes: PTHLH (12p11) has a crucial role in mammary gland development and the establishment of bone metastasis in breast cancer, and NRIP1 (21q21) encodes an ER cofactor and has a role in the regulation of breast cancer cell growth.

Common breast cancer susceptibility loci are associated with triple-negative breast cancer

Triple-negative breast cancers are an aggressive subtype of breast cancer with poor survival, but there remains little known about the etiologic factors that promote its initiation and development. Commonly inherited breast cancer risk factors identified through genome-wide association studies display heterogeneity of effect among breast cancer subtypes as defined by the status of estrogen and progesterone receptors. In the Triple Negative Breast Cancer Consortium (TNBCC), 22 common breast cancer susceptibility variants were investigated in 2,980 Caucasian women with triple-negative breast cancer and 4,978 healthy controls. We identified six single-nucleotide polymorphisms, including rs2046210 (ESR1), rs12662670 (ESR1), rs3803662 (TOX3), rs999737 (RAD51L1), rs8170 (19p13.1), and rs8100241 (19p13.1), significantly associated with the risk of triple-negative breast cancer. Together, our results provide convincing evidence of genetic susceptibility for triple-negative breast cancer.

A genome-wide association scan on estrogen receptor-negative breast cancer

Introduction

Breast cancer is a heterogeneous disease and may be characterized on the basis of whether estrogen receptors (ER) are expressed in the tumour cells. ER status of breast cancer is important clinically, and is used both as a prognostic indicator and treatment predictor. In this study, we focused on identifying genetic markers associated with ER-negative breast cancer risk.

Methods

We conducted a genome-wide association analysis of 285,984 single nucleotide polymorphisms (SNPs) genotyped in 617 ER-negative breast cancer cases and 4,583 controls. We also conducted a genome-wide pathway analysis on the discovery dataset using permutation-based tests on pre-defined pathways. The extent of shared polygenic variation between ER-negative and ER-positive breast cancers was assessed by relating risk scores, derived using ER-positive breast cancer samples, to disease state in independent, ER-negative breast cancer cases.

Results

Association with ER-negative breast cancer was not validated for any of the five most strongly associated SNPs followed up in independent studies (1,011 ER-negative breast cancer cases, 7,604 controls). However, an excess of small P-values for SNPs with known regulatory functions in cancer-related pathways was found (global P = 0.052). We found no evidence to suggest that ER-negative breast cancer shares a polygenic basis to disease with ER-positive breast cancer.

Conclusions

ER-negative breast cancer is a distinct breast cancer subtype that merits independent analyses. Given the clinical importance of this phenotype and the likelihood that genetic effect sizes are small, greater sample sizes and further studies are required to understand the etiology of ER-negative breast cancers.