Identification of an alternative short ARID5B isoform associated with B-ALL survival.

Utilizing RNA sequence (RNA-Seq) splice junction data from a cohort of 1841 B-cell acute lymphoblastic leukemia (B-ALL) patients we define transcriptionally distinct isoforms of ARID5B, a risk-associated gene identified in genome wide association studies (GWAS), which associate with disease survival. Short (S) and long (L) ARID5B transcripts, which differ in an encoded BAH-like chromatin interaction domain, show remarkable correlation to the isoform splicing pattern. Testing of the ARID5B proximal promoter of the S & L isoforms indicated that both are functionally independent in luciferase reporter assays. Increased short isoform expression is associated with decreased event-free and overall survival. The abundance of short and long transcripts strongly correlates to B-ALL prognostic stratification, where B-ALL subtypes with poor outcomes express a higher proportion of the S-isoform. These data demonstrate that the analysis of independent promoters and alternative splicing events are essential for improved risk stratification and a more complete understanding of disease pathology.

Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia.

Transcriptome sequencing has identified multiple subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL) of prognostic significance, but a minority of cases lack a known genetic driver. Here, we used integrated whole-genome (WGS) and -transcriptome sequencing (RNA-seq), enhancer mapping, and chromatin topology analysis to identify previously unrecognized genomic drivers in B-ALL. Newly diagnosed (n = 3221) and relapsed (n = 177) B-ALL cases with tumor RNA-seq were studied. WGS was performed to detect mutations, structural variants, and copy number alterations. Integrated analysis of histone 3 lysine 27 acetylation and chromatin looping was performed using HiChIP. We identified a subset of 17 newly diagnosed and 5 relapsed B-ALL cases with a distinct gene expression profile and 2 universal and unique genomic alterations resulting from aberrant recombination-activating gene activation: a focal deletion downstream of PAN3 at 13q12.2 resulting in CDX2 deregulation by the PAN3 enhancer and a focal deletion of exons 18-21 of UBTF at 17q21.31 resulting in a chimeric fusion, UBTF::ATXN7L3. A subset of cases also had rearrangement and increased expression of the PAX5 gene, which is otherwise uncommon in B-ALL. Patients were more commonly female and young adult with median age 35 (range,12-70 years). The immunophenotype was characterized by CD10 negativity and immunoglobulin M positivity. Among 16 patients with known clinical response, 9 (56.3%) had high-risk features including relapse (n = 4) or minimal residual disease >1% at the end of remission induction (n = 5). CDX2-deregulated, UBTF::ATXN7L3 rearranged (CDX2/UBTF) B-ALL is a high-risk subtype of leukemia in young adults for which novel therapeutic approaches are required.

The genetic basis and cell of origin of mixed phenotype acute leukaemia.

Mixed phenotype acute leukaemia (MPAL) is a high-risk subtype of leukaemia with myeloid and lymphoid features, limited genetic characterization, and a lack of consensus regarding appropriate therapy. Here we show that the two principal subtypes of MPAL, T/myeloid (T/M) and B/myeloid (B/M), are genetically distinct. Rearrangement of ZNF384 is common in B/M MPAL, and biallelic WT1 alterations are common in T/M MPAL, which shares genomic features with early T-cell precursor acute lymphoblastic leukaemia. We show that the intratumoral immunophenotypic heterogeneity characteristic of MPAL is independent of somatic genetic variation, that founding lesions arise in primitive haematopoietic progenitors, and that individual phenotypic subpopulations can reconstitute the immunophenotypic diversity in vivo. These findings indicate that the cell of origin and founding lesions, rather than an accumulation of distinct genomic alterations, prime tumour cells for lineage promiscuity. Moreover, these findings position MPAL in the spectrum of immature leukaemias and provide a genetically informed framework for future clinical trials of potential treatments for MPAL.

Tumor-intrinsic and -extrinsic determinants of response to blinatumomab in adults with B-ALL.

Blinatumomab, a bispecific antibody that directs CD3+ T cells to CD19+ tumor cells, shows variable efficacy in B-progenitor acute lymphoblastic leukemia (B-ALL). To determine tumor-intrinsic and -extrinsic determinants of response, we studied 44 adults with relapsed or refractory B-ALL (including 2 minimal residual disease positive) treated with blinatumomab using bulk tumor and single-cell sequencing. The overall response rate in patients with hematological disease was 55%, with a high response rate in those with CRLF2-rearranged Philadelphia chromosome-like ALL (12 [75%] of 16). Pretreatment samples of responders exhibited a tumor-intrinsic transcriptomic signature of heightened immune response. Multiple mechanisms resulted in loss of CD19 expression, including CD19 mutations, CD19-mutant allele-specific expression, low CD19 RNA expression, and mutations in CD19 signaling complex member CD81. Patients with low hypodiploid ALL were prone to CD19- relapse resulting from aneuploidy-mediated loss of the nonmutated CD19 allele. Increased expression of a CD19 isoform with intraexonic splicing of exon 2, CD19 ex2part, at baseline or during therapy was associated with treatment failure. These analyses demonstrate both tumor-intrinsic and -extrinsic factors influence blinatumomab response. We show that CD19 mutations are commonly detected in CD19- relapse during blinatumomab treatment. Identification of the CD19 ex2part splice variant represents a new biomarker predictive of blinatumomab therapy failure.

The age of the bone marrow microenvironment influences B-cell acute lymphoblastic leukemia progression via CXCR5-CXCL13.

B-cell acute lymphoblastic leukemia (B-ALL) occurs most commonly in children, whereas chronic myeloid leukemia is more frequent in adults. The myeloid bias of hematopoiesis in elderly individuals has been considered causative, but the age of the bone marrow microenvironment (BMM) may be contributory. Using various murine models of B-ALL in young vs old mice, we recapitulated B-ALL preponderance in children vs adults. We showed differential effects of young vs old BM macrophages on B-ALL cell function. Molecular profiling using RNA- and ATAC-sequencing revealed pronounced differences in young vs old BMM-derived macrophages and enrichment for gene sets associated with inflammation. In concordance with the role of C-X-C motif chemokine (CXCL) 13 for disease-associated B-cell chemoattraction, we found CXCL13 to be highly expressed in young macrophages on a translational compared with a transcriptional level. Inhibition of CXCL13 in BM macrophages impaired leukemia cell migration and decreased the proliferation of cocultured B-ALL cells, whereas recombinant CXCL13 increased pAKT and B-ALL cell expansion. Pretreatment of B-ALL-initiating cells with CXCL13 accelerated B-ALL progression. Deficiency of Cxcr5, the receptor for CXCL13, on B-ALL-initiating cells prolonged murine survival, whereas high expression of CXCR5 in pediatric B-ALL may predict central nervous system relapse. CXCL13 staining was increased in bone sections from pediatric compared with adult patients with B-ALL. Taken together, our study shows that the age of the BMM and, in particular, BM macrophages influence the leukemia phenotype. The CXCR5-CXCL13 axis may act as prognostic marker and an attractive novel target for the treatment of B-ALL.

Molecular classification improves risk assessment in adult BCR-ABL1-negative B-ALL.

Genomic classification has improved risk assignment of pediatric, but not adult B-lineage acute lymphoblastic leukemia (B-ALL). The international UKALLXII/ECOG-ACRIN E2993 (#NCT00002514) trial accrued 1229 adolescent/adult patients with BCR-ABL1- B-ALL (aged 14 to 65 years). Although 93% of patients achieved remission, 41% relapsed at a median of 13 months (range, 28 days to 12 years). Five-year overall survival (OS) was 42% (95% confidence interval, 39, 44). Transcriptome sequencing, gene expression profiling, cytogenetics, and fusion polymerase chain reaction enabled genomic subtyping of 282 patient samples, of which 264 were eligible for trial, accounting for 64.5% of E2993 patients. Among patients with outcome data, 29.5% with favorable outcomes (5-year OS 65% to 80%) were deemed standard risk (DUX4-rearranged [9.2%], ETV6-RUNX1/-like [2.3%], TCF3-PBX1 [6.9%], PAX5 P80R [4.1%], high-hyperdiploid [6.9%]); 50.2% had high-risk genotypes with 5-year OS of 0% to 27% (Ph-like [21.2%], KMT2A-AFF1 [12%], low-hypodiploid/near-haploid [14.3%], BCL2/MYC-rearranged [2.8%]); 20.3% had intermediate-risk genotypes with 5-year OS of 33% to 45% (PAX5alt [12.4%], ZNF384/-like [5.1%], MEF2D-rearranged [2.8%]). IKZF1 alterations occurred in 86% of Ph-like, and TP53 mutations in patients who were low-hypodiploid (54%) and BCL2/MYC-rearranged (33%) but were not independently associated with outcome. Of patients considered high risk based on presenting age and white blood cell count, 40% harbored subtype-defining genetic alterations associated with standard- or intermediate-risk outcomes. We identified distinct immunophenotypic features for DUX4-rearranged, PAX5 P80R, ZNF384-R/-like, and Ph-like genotypes. These data in a large adult B-ALL cohort treated with a non-risk-adapted approach on a single trial show the prognostic importance of genomic analyses, which may translate into future therapeutic benefits.

MD-ALL: an integrative platform for molecular diagnosis of B-acute lymphoblastic leukemia.

B-acute lymphoblastic leukemia (B-ALL) consists of dozens of subtypes defined by distinct gene expression profiles (GEPs) and various genetic lesions. With the application of transcriptome sequencing (RNA-seq), multiple novel subtypes have been identified, which lead to an advanced B-ALL classification and risk-stratification system. However, the complexity of analyzing RNA-seq data for B-ALL classification hinders the implementation of the new B-ALL taxonomy. Here, we introduce MD-ALL (Molecular Diagnosis of ALL), an integrative platform featuring sensitive and accurate B-ALL classification based on GEPs and sentinel genetic alterations from RNA-seq data. In this study, we systematically analyzed 2,955 B-ALL RNA-seq samples and generated a reference dataset representing all the reported B-ALL subtypes. Using multiple machine learning algorithms, we identified the feature genes and then established highly sensitive and accurate models for B-ALL classification using either bulk or single-cell RNA-seq data. Importantly, this platform integrates multiple aspects of key genetic lesions acquired from RNAseq data, which include sequence mutations, large-scale copy number variations, and gene rearrangements, to perform comprehensive and definitive B-ALL classification. Through validation in a hold-out cohort of 974 samples, our models demonstrated superior performance for B-ALL classification compared with alternative tools. Moreover, to ensure accessibility and user-friendly navigation even for users with limited or no programming background, we developed an interactive graphical user interface for this MD-ALL platform, using the R Shiny package. In summary, MD-ALL is a user-friendly B-ALL classification platform designed to enable integrative, accurate, and comprehensive B-ALL subtype classification. MD-ALL is available from https://github.com/gu-lab20/MD-ALL.

High response rates and transition to transplant after novel targeted and cellular therapies in adults with relapsed/refractory acute lymphoblastic leukemia with Philadelphia-like fusions.

Philadelphia (Ph)-like acute lymphoblastic leukemia (ALL) is associated with a poor response to standard chemotherapy. However, outcomes with novel antibody and cellular therapies in relapsed/refractory (r/r) Ph-like ALL are largely unknown. We conducted a single-center retrospective analysis of adult patients (n = 96) with r/r B-ALL and fusions associated with Ph-like who received novel salvage therapies. Patients were treated with 149 individual novel regimens (blinatumomab = 83, inotuzumab ozogamicin [InO] = 36, and CD19CAR T cells = 30). The median age at first novel salvage therapy was 36 years (range; 18-71). Ph-like fusions were IGH::CRLF2 (n = 48), P2RY8::CRLF2 (n = 26), JAK2 (n = 9), ABL-class (n = 8), EPOR::IGH (n = 4) and ETV6::NTRK2 (n = 1). CD19CAR T cells were administered later in the course of therapy compared to blinatumomab and InO (p < .001) and more frequently in recipients who relapsed after allogeneic hematopoietic cell transplantation (alloHCT) (p = .002). Blinatumomab was administered at an older age compared to InO and CAR T-cells (p = .004). The complete remission (CR)/CR with incomplete hematologic recovery (CRi) rates were 63%, 72%, and 90% following blinatumomab, InO and CD19CAR, respectively, among which 50%, 50%, and 44% of responders underwent consolidation with alloHCT, respectively. In multivariable analysis, the type of novel therapy (p = .044) and pretreatment marrow blasts (p = .006) predicted the CR/CRi rate, while the Ph-like fusion subtype (p = .016), pretreatment marrow blasts (p = .022) and post-response consolidation with alloHCT (p < .001) influenced event-free survival. In conclusion, novel therapies are effective in inducing high remission rates in patients with r/r Ph-like ALL and successfully transitioning the responders to alloHCT.

Recurrent retinal detachment after pars plana vitrectomy with silicone oil tamponade for rhegmatogenous retinal detachment.

BACKGROUND: The recurrence of retinal detachment following rhegmatogenous retinal detachment (RRD) is a relatively common complication that can lead to reduced visual acuity and requires further surgery. The purpose of this study was to investigate the risk factors and visual outcomes of recurrent RRD following pars plana vitrectomy (PPV) with silicone oil tamponade for primary RRD. METHODS: This was a retrospective follow-up study of 343 eyes that underwent initial PPV surgery with silicone oil tamponade for primary RRD. Patients were divided into a recurrence group and a reattachment group. The main outcome measures included causative factors, visual outcomes related to the recurrence of RRD, and the perioperative factors most affecting the recurrence of RRD. RESULTS: After retinal reattachment, we observed RRD recurrence after PPV for primary RRD in 42 out of 343 eyes (12.2%) during the follow-up period. Most causes of recurrence (69%) occurred within 6 months of surgery. Multivariate logistic regression analysis showed that a PVR ≥ Grade C (odds ratio [OR]: 4.015; 95% confidence interval [CI] 1.721-9.367; P = 0.001) was a significant predictor for the development of recurrent RRD. Compared with the reattachment group, the recurrence group exhibited a significant decline in best-corrected visual acuity (BCVA) at the last follow-up visit (P = 0.000). Eyes with PVR prior to primary surgery, or at the diagnosis of re-detachment, showed a worse final BCVA. CONCLUSIONS: Our analysis shows that the predominant risk factor for the recurrence of RRD is a PVR ≥ Grade C. PVR prior to primary surgery, or at the diagnosis of re-detachment, was also shown to limit the recovery of final visual acuity.

ShinyCNV: a Shiny/R application to view and annotate DNA copy number variations.

Motivation: Single nucleotide polymorphism (SNP) array is the most widely used platform to assess somatic copy number variations (CNVs) in cancer studies. Many SNP data-based CNV callers are available, however, the false positive rates from automated calling are commonly high, and reported breakpoints can be inaccurate. Manual review for each reported CNV by visualizing the SNP data is important, but is challenging for users lacking computational experience. To address this, we present a Shiny/R application ShinyCNV, an interactive graphical user interface to view and annotate CNVs. Results: With this application, normalized SNP data, which includes log R ratio (LRR) and B allele frequency, can be plotted against the reported CNVs, and users can visually check the reliability of CNVs per se or adjust the incorrectly assigned breakpoints. Further, the interactive LRR spectrum panel within ShinyCNV can facilitate the process to identify commonly affected CNV regions from a group of samples, and to visually check if important focal gains/losses are missing from reported CNVs. ShinyCNV is designed to be intuitive for cancer researchers and can be easily installed for either personal use or deployed on servers to provide online service. Availability and implementation: ShinyCNV and the tutorial are freely available from https://github.com/gzhmat/ShinyCNV. Supplementary information: Supplementary data are available at Bioinformatics online.

Genomic and outcome analyses of Ph-like ALL in NCI standard-risk patients: a report from the Children's Oncology Group.

Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL; BCR-ABL1-like ALL) in children with National Cancer Institute (NCI) intermediate- or high-risk (HR) ALL is associated with poor outcome. Ph-like ALL is characterized by genetic alterations that activate cytokine receptor and kinase signaling and may be amenable to treatment with tyrosine kinase inhibitors. The prevalence, outcome, and potential for targeted therapy of Ph-like ALL in standard-risk (SR) ALL is less clear. We retrospectively analyzed a cohort of 1023 SR childhood B-ALL consecutively enrolled in the Children's Oncology Group AALL0331 clinical trial. The Ph-like ALL gene expression profile was identified in 206 patients, and 67 patients with either BCR-ABL1 (n = 6) or ETV6-RUNX1 (n = 61) were excluded from downstream analysis, leaving 139 of 1023 (13.6%) as Ph-like. Targeted reverse transcription polymerase chain reaction assays and RNA-sequencing identified kinase-activating alterations in 38.8% of SR Ph-like cases, including CRLF2 rearrangements (29.5% of Ph-like), ABL-class fusions (1.4%), JAK2 fusions (1.4%), an NTRK3 fusion (0.7%), and other sequence mutations (IL7R, KRAS, NRAS; 5.6%). Patients with Ph-like ALL had inferior 7-year event-free survival compared with non-Ph-like ALL (82.4 ± 3.6% vs 90.7 ± 1.0%, P = .0022), with no difference in overall survival (93.2 ± 2.4% vs 95.8 ± 0.7%, P = .14). These findings illustrate the significant differences in the spectrum of kinase alterations and clinical outcome of Ph-like ALL based on presenting clinical features and establish that genomic alterations potentially targetable with approved kinase inhibitors are less frequent in SR than in HR ALL.

Targetable kinase gene fusions in high-risk B-ALL: a study from the Children's Oncology Group.

Philadelphia chromosome-like (Ph-like) acute lymphoblastic leukemia (ALL) is a high-risk subtype characterized by genomic alterations that activate cytokine receptor and kinase signaling. We examined the frequency and spectrum of targetable genetic lesions in a retrospective cohort of 1389 consecutively diagnosed patients with childhood B-lineage ALL with high-risk clinical features and/or elevated minimal residual disease at the end of remission induction therapy. The Ph-like gene expression profile was identified in 341 of 1389 patients, 57 of whom were excluded from additional analyses because of the presence of BCR-ABL1 (n = 46) or ETV6-RUNX1 (n = 11). Among the remaining 284 patients (20.4%), overexpression and rearrangement of CRLF2 (IGH-CRLF2 or P2RY8-CRLF2) were identified in 124 (43.7%), with concomitant genomic alterations activating the JAK-STAT pathway (JAK1, JAK2, IL7R) identified in 63 patients (50.8% of those with CRLF2 rearrangement). Among the remaining patients, using reverse transcriptase polymerase chain reaction or transcriptome sequencing, we identified targetable ABL-class fusions (ABL1, ABL2, CSF1R, and PDGFRB) in 14.1%, EPOR rearrangements or JAK2 fusions in 8.8%, alterations activating other JAK-STAT signaling genes (IL7R, SH2B3, JAK1) in 6.3% or other kinases (FLT3, NTRK3, LYN) in 4.6%, and mutations involving the Ras pathway (KRAS, NRAS, NF1, PTPN11) in 6% of those with Ph-like ALL. We identified 8 new rearrangement partners for 4 kinase genes previously reported to be rearranged in Ph-like ALL. The current findings provide support for the precision-medicine testing and treatment approach for Ph-like ALL implemented in Children's Oncology Group ALL trials.

Multiregion sequencing reveals the intratumor heterogeneity of driver mutations in TP53-driven non-small cell lung cancer.

Intratumor heterogeneity (ITH) in non-small cell lung cancer (NSCLC) may account for resistance after a period of targeted therapies because drugs destroy only a portion of tumor cells. The recognition of ITH helps identify high-risk patients to make effective treatment decisions. However, ITH studies are confounded by interpatient heterogeneity in NSCLC and a large amount of passenger mutations. To address these issues, we recruited NSCLC patients carrying TP53 mutations and selected driver mutations within recurrently mutated genes in NSCLC. A total of 12-paired normal-tumor tissues were subjected to whole-genome/whole-exome sequencing. From these, 367 non-silent mutations were selected as driver mutations and deeply sequenced in 61 intratumoral microdissections. We identified a universal prevalence of heterogeneity in all 12 tumors, indicating branched evolution. Although TP53 mutations were observed in single biopsy of all 12 tumors, most tumors consist of both TP53 mutated and non-mutated cells in separate regions within the same tumor. This suggests the late molecular timing of the acquisition of TP53 mutations; therefore, the detection of TP53 mutations in a single biopsy may simply not reflect the early malignant potential. In addition, we identified regions of loss of heterozygosity surrounding TP53 and CDKN2A mutations in tumor 711, which also exhibited heterogeneity in different regional samples. Because the ITH of driver mutations likely has clinical consequences, further efforts are needed to limit the impact of ITH and to improve therapeutic efficiency, which will benefit NSCLC patients receiving targeted treatments.

DNMT3A Arg882 mutation drives chronic myelomonocytic leukemia through disturbing gene expression/DNA methylation in hematopoietic cells.

The gene encoding DNA methyltransferase 3A (DNMT3A) is mutated in ∼20% of acute myeloid leukemia cases, with Arg882 (R882) as the hotspot. Here, we addressed the transformation ability of the DNMT3A-Arg882His (R882H) mutant by using a retroviral transduction and bone marrow transplantation (BMT) approach and found that the mutant gene can induce aberrant proliferation of hematopoietic stem/progenitor cells. At 12 mo post-BMT, all mice developed chronic myelomonocytic leukemia with thrombocytosis. RNA microarray analysis revealed abnormal expressions of some hematopoiesis-related genes, and the DNA methylation assay identified corresponding changes in methylation patterns in gene body regions. Moreover, DNMT3A-R882H increased the CDK1 protein level and enhanced cell-cycle activity, thereby contributing to leukemogenesis.

Genomic landscape of CD34+ hematopoietic cells in myelodysplastic syndrome and gene mutation profiles as prognostic markers.

Myelodysplastic syndrome (MDS) includes a group of diseases characterized by dysplasia of bone marrow myeloid lineages with ineffective hematopoiesis and frequent evolution to acute myeloid leukemia (AML). Whole-genome sequencing was performed in CD34(+) hematopoietic stem/progenitor cells (HSPCs) from eight cases of refractory anemia with excess blasts (RAEB), the high-risk subtype of MDS. The nucleotide substitution patterns were found similar to those reported in AML, and mutations of 96 protein-coding genes were identified. Clonal architecture analysis revealed the presence of subclones in six of eight cases, whereas mutation detection of CD34(+) versus CD34(-) cells revealed heterogeneity of HSPC expansion status. With 39 marker genes belonging to eight functional categories, mutations were analyzed in 196 MDS cases including mostly RAEB (n = 89) and refractory cytopenia with multilineage dysplasia (RCMD) (n = 95). At least one gene mutation was detected in 91.0% of RAEB, contrary to that in RCMD (55.8%), suggesting a higher mutational burden in the former group. Gene abnormality patterns differed between MDS and AML, with mutations of activated signaling molecules and NPM1 being rare, whereas those of spliceosome more common, in MDS. Finally, gene mutation profiles also bore prognostic value in terms of overall survival and progression free survival.

Genome-wide association study identifies a novel susceptibility gene for serum TSH levels in Chinese populations.

Thyroid-stimulating hormone (TSH) is a sensitive indicator of thyroid function. High and low TSH levels reflect hypothyroidism and hyperthyroidism, respectively. Even within the normal range, small differences in TSH levels, on the order of 0.5-1.0 mU/l, are associated with significant differences in blood pressure, BMI, dyslipidemia, risk of atrial fibrillation and atherosclerosis. Most of the variance in TSH levels is thought to be genetically influenced. We conducted a genome-wide association study of TSH levels in 1346 Chinese Han individuals. In the replication study, we genotyped four candidate SNPs with the top association signals in an independent isolated Chinese She cohort (n = 3235). We identified a novel serum TSH susceptibility locus within XKR4 at 8q12.1 (rs2622590, Pcombined = 2.21 × 10(-10)), and we confirmed two previously reported TSH susceptibility loci near FOXE1 at 9q22.33 and near CAPZB at 1p36.13, respectively. The rs2622590_T allele at XKR4 and the rs925489_C allele near FOXE1 were correlated with low TSH levels and were found to be nominally associated to patients with papillary thyroid carcinoma (PTC) (OR = 1.41, P= 0.014 for rs2622590_T, and OR = 1.61, P= 0.030 for rs925489_C). The rs2622590 and rs925489 genotypes were also correlated with the expression levels of FOXE1 and XKR4, respectively, in PTC tissues (P = 2.41 × 10(-4) and P= 0.02). Our findings suggest that the SNPs in XKR4 and near FOXE1 are involved in the regulation of TSH levels.

Robust evidence for five new Graves' disease risk loci from a staged genome-wide association analysis.

Graves' disease (GD), characterized by autoantibodies targeting antigens specifically expressed in thyroid tissues causing hyperthyroidism, is triggered by a combination of genetic and environmental factors. However, only a few loci for GD risk were confirmed in the various ethnic groups, and additional genetic determinants have to be detected. In this study, we carried out a three-stage study in 9529 patients with GD and 9984 controls to identify new risk loci for GD and found genome-wide significant associations in the overall populations for five novel susceptibility loci: the GPR174-ITM2A at Xq21.1, C1QTNF6-RAC2 at 22q12.3-13.1, SLAMF6 at 1q23.2, ABO at 9q34.2 and an intergenic region harboring two non-coding RNAs at 14q32.2 and one previous indefinite locus, TG at 8q24.22 (Pcombined < 5 × 10(-8)). The genotypes of corresponding variants at 14q32.2 and 8q24.22 were correlated with the expression levels of C14orf64 and a TG transcript skipping exon 46, respectively. This study increased the number of GD loci with compelling evidence and indicated that non-coding RNAs might be potentially involved in the pathogenesis of GD.

Identification of BACH2 as a susceptibility gene for Graves' disease in the Chinese Han population based on a three-stage genome-wide association study.

The BACH2 gene regulates B cell differentiation and function and has been reported to be a shared susceptibility gene for several autoimmune diseases. Our previous genome-wide association study (GWAS) indicated that several single nucleotide polymorphisms (SNPs) in the BACH2 gene are associated with Graves' disease (GD) in the Chinese Han population; however, the association did not achieve genome-wide significance levels. Recently, this association of BACH2 with GD was confirmed in Caucasians in the UK population, but fine mapping in this region has not yet been reported. Here, we provide a refined analysis of a 331-kb region in the BACH2 gene, which harbors 359 SNPs, using GWAS data from 1,442 GD patients and 1,468 controls. The SNPs rs2474619 and rs9344996 were implied as the independent variants associated with GD by forward and two-locus logistic regression analysis. We genotyped eight out of 10 tagSNPs with P < 1 × 10(-3) in 3,508 GD patients and 3,209 controls, the results also showed that rs2474619 was independently associated with GD in the combined population from GWAS and the second stage (P = 1.81 × 10(-5)). The rs2474619 and rs9344996 were further genotyped in the third stage cohorts, and rs2474619 showed evidence of association with GD at genome-wide significance levels in the combined population (P = 3.28 × 10(-8), odds ratio = 1.13). The association of rs9344996 with GD can be explained by its linkage to rs2474619 in the combined population. Our study clearly demonstrated that BACH2 is a susceptibility gene for GD in the Chinese Han population and further supported rs2474619, in intron 2 of BACH2, is the best association signal with GD. However, the mechanism by which BACH2 confers increased risk of GD requires further study.