Precise immunofluorescence canceling for highly multiplexed imaging to capture specific cell states.

Cell states are regulated by the response of signaling pathways to receptor ligand-binding and intercellular interactions. High-resolution imaging has been attempted to explore the dynamics of these processes and, recently, multiplexed imaging has profiled cell states by achieving a comprehensive acquisition of spatial protein information from cells. However, the specificity of antibodies is still compromised when visualizing activated signals. Here, we develop Precise Emission Canceling Antibodies (PECAbs) that have cleavable fluorescent labeling. PECAbs enable high-specificity sequential imaging using hundreds of antibodies, allowing for reconstruction of the spatiotemporal dynamics of signaling pathways. Additionally, combining this approach with seq-smFISH can effectively classify cells and identify their signal activation states in human tissue. Overall, the PECAb system can serve as a comprehensive platform for analyzing complex cell processes.

Identification of the hybrid gene LILRB5-3 by long-read sequencing and implication of its novel signaling function.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) on human chromosome 19q13.4 encode 11 immunoglobulin superfamily receptors, exhibiting genetic diversity within and between human populations. Among the LILR genes, the genomic region surrounding LILRB3 and LILRA6 has yet to be fully characterized due to their significant sequence homology, which makes it difficult to differentiate between them. To examine the LILRB3 and LILRA6 genomic region, a tool named JoGo-LILR CN Caller, which can call copy number from short-read whole genome sequencing (srWGS) data, was applied to an extensive international srWGS dataset comprising 2,504 samples. During this process, a previously unreported loss of both LILRB3 and LILRA6 was detected in three samples. Using long-read sequencing of these samples, we have discovered a novel large deletion (33,692 bp) in the LILRB3 and LILRA6 genomic regions in the Japanese population. This deletion spanned three genes, LILRB3, LILRA6, and LILRB5, resulting in LILRB3 exons 12-13 being located immediately downstream of LILRB5 exons 1-12 with the loss of LILRA6, suggesting the potential expression of a hybrid gene between LILRB5 and LILRB3 (LILRB5-3). Transcription and subsequent translation of the LILRB5-3 hybrid gene were also verified. The hybrid junction was located within the intracellular domain, resulting in an LILRB5 extracellular domain fused to a partial LILRB3 intracellular domain with three immunoreceptor tyrosine-based inhibitory motifs (ITIMs), suggesting that LILRB5-3 acquired a novel signaling function. Further application of the JoGo-LILR tool to srWGS samples suggested the presence of the LILRB5-3 hybrid gene in the CEU population. Our findings provide insight into the genetic and functional diversity of the LILR family.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases.

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Analysis of Genetic Polymorphism of Bitter Taste Receptor TAS2R38 and TAS2R46, and Its Relationship with Eating and Drinking Habits in Japanese ToMMo Subjects.

Human type 2 taste receptor (TAS2R) genes encode bitter-taste receptors that are activated by various bitter ligands. It has been said that TAS2R38 may detect bitter substances and then suppress their intake by controlling gustatory or digestive responses. The major haplotypes of TAS2R38 involve three non-synonymous, closely-linked single-nucleotide polymorphisms (SNPs), leading to three amino acid substitutions (A49P, V262A and I296V) and resulting in a PAV or AVI allele. The allele frequency of AVI/PAV was 0.42/0.58 in this study. The genotype frequency distributions of TAS2R38 were 18.32%, 46.95% and 33.95% for AVI/AVI, AVI/PAV and PAV/PAV, respectively, and were in Hardy-Weinberg equilibrium. Five haplotype combinations of minor alleles were identified: AVI/AAV, AVI/AVV, AAI/PAV, AVI/PVV, AVI/AAI, with corresponding frequencies of 0.49%, 0.10%, 0.10%, 0.05%, 0.05%, respectively, in 2,047 Japanese Tohoku Medical Megabank Organization (ToMMo) subjects (2KJPN). The 16 subjects with these minor alleles were excluded from the questionnaire analysis, which found no significant differences among the major TAS2R38 genotypes (AVI/AVI, AVI/PAV and PAV/PAV) in the intake frequency of cruciferous vegetables or in the frequency of drinking alcohol. This result differs from previous data using American and European subjects. This is the first study to analyze the relationship between TAS2R38 genotype and the eating and drinking habits of Japanese subjects. It was also shown that there were no relationships at all between the genetic polymorphism of TAS2R46 and the phenotypes such as clinical BMI, eating and drinking habits among the 3 genotypes of TAS2R46 (∗/∗, ∗/W, W/W) at position W250∗ (∗stop codon).

HLA-DQA1*05 and upstream variants of PPARGC1B are associated with infliximab persistence in Japanese Crohn's disease patients.

Recently, the HLA-DQA1*05 (rs2097432) genetic variation has been reported to be linked to early infliximab (IFX) treatment failure in the Caucasian Crohn's disease (CD) population, but that evidence is scarce in the Asian population. This study aimed to investigate the relationship between rs2097432 and the cumulative discontinuation-free time of IFX (IFX persistence) in 189 Japanese biologics-naive CD patients. We also performed a genome-wide association study (GWAS) to discover novel genetic predictors for IFX persistence. The C allele of rs2097432 significantly increased the risk of early discontinuation of IFX [Hazard ratio (HR) = 2.23 and P-value = 0.026]. In GWAS, one locus tagged by rs73277969, located upstream of PPARGC1B which attenuates macrophage-mediated inflammation, reached genome-wide significance (HR = 6.04 and P-value = 7.93E-9). Pathway analysis suggested association of signaling by PDGF and FCGR activation signaling with IFX persistence (P-value = 8.56E-5 and 5.80E-4, respectively).

Novel HLA allele associations with susceptibility, staging, symptomatic state, autoimmune hepatitis and hepatocellular carcinoma events for primary biliary cholangitis in the Japanese population.

Primary biliary cholangitis (PBC) is a rare autoimmune disease with a clear predisposition for human leukocyte antigen (HLA)-DR/DQ-associated loss of immune tolerance for the E2 component of the pyruvate dehydrogenase complex. Three-field-resolution HLA imputation of 1,670 Japanese PBC patients and 2,328 healthy controls was conducted using Japanese population-specific HLA reference panels. Eighteen previously reported Japanese PBC-associated HLA alleles were confirmed and extended to 3-field-resolution, including HLA-DRB1*08:03 to HLA-DRB1*08:03:02, HLA-DQB1*03:01 to HLA-DQB1*03:01:01, HLA-DQB1*04:01 to HLA-DQB1*04:01:01 and HLA-DQB1*06:04 to HLA-DQB1*06:04:01. In addition, additional significant novel HLA alleles were identified, including 3 novel susceptible HLA-DQA1 alleles: HLA-DQA1*03:03:01, HLA-DQA1*04:01:01, HLA-DQA1*01:04:01 and 1 novel protective HLA-DQA1 allele, HLA-DQA1*05:05:01. In addition, PBC patients carrying HLA-DRB1*15:01:01 and HLA-DQA1*03:03:01 would have a higher predisposition toward developing concomitant autoimmune hepatitis (AIH). Further, late-stage and symptomatic PBC shared the same susceptible HLA alleles of HLA-A*26:01:01, HLA-DRB1*09:01:02 and HLA-DQB1*03:03:02. Lastly, HLA-DPB1*05:01:01 was identified as a potential risk HLA allele for development of hepatocellular carcinoma (HCC) in PBC patients. In conclusion, we have extended the current knowledge of HLA allele associations to 3-field resolution and identified novel HLA allele associations with predisposition risk, staging, symptomatic state, and AIH and HCC events for Japanese PBC patients.

Genetic architecture of the inflammatory bowel diseases across East Asian and European ancestries.

Inflammatory bowel diseases (IBDs) are chronic disorders of the gastrointestinal tract with the following two subtypes: Crohn's disease (CD) and ulcerative colitis (UC). To date, most IBD genetic associations were derived from individuals of European (EUR) ancestries. Here we report the largest IBD study of individuals of East Asian (EAS) ancestries, including 14,393 cases and 15,456 controls. We found 80 IBD loci in EAS alone and 320 when meta-analyzed with ~370,000 EUR individuals (~30,000 cases), among which 81 are new. EAS-enriched coding variants implicate many new IBD genes, including ADAP1 and GIT2. Although IBD genetic effects are generally consistent across ancestries, genetics underlying CD appears more ancestry dependent than UC, driven by allele frequency (NOD2) and effect (TNFSF15). We extended the IBD polygenic risk score (PRS) by incorporating both ancestries, greatly improving its accuracy and highlighting the importance of diversity for the equitable deployment of PRS.

Prediction Model with HLA-A*33:03 Reveals Number of Days to Develop Liver Cancer from Blood Test.

The development of liver cancer in patients with hepatitis B is a major problem, and several models have been reported to predict the development of liver cancer. However, no predictive model involving human genetic factors has been reported to date. For the items incorporated in the prediction model reported so far, we selected items that were significant in predicting liver carcinogenesis in Japanese patients with hepatitis B and constructed a prediction model of liver carcinogenesis by the Cox proportional hazard model with the addition of Human Leukocyte Antigen (HLA) genotypes. The model, which included four items-sex, age at the time of examination, alpha-fetoprotein level (log10AFP) and presence or absence of HLA-A*33:03-revealed an area under the receiver operating characteristic curve (AUROC) of 0.862 for HCC prediction within 1 year and an AUROC of 0.863 within 3 years. A 1000 repeated validation test resulted in a C-index of 0.75 or higher, or sensitivity of 0.70 or higher, indicating that this predictive model can distinguish those at high risk of developing liver cancer within a few years with high accuracy. The prediction model constructed in this study, which can distinguish between chronic hepatitis B patients who develop hepatocellular carcinoma (HCC) early and those who develop HCC late or not, is clinically meaningful.

Genome-Wide Association Study and Transcriptome of Japanese Patients with Developmental Dysplasia of the Hip Demonstrates an Association with the Ferroptosis Signaling Pathway.

This study examined the association between developmental dysplasia of the hip (DDH) and disease-associated loci in a Japanese cohort. A genome-wide association study (GWAS) of 238 Japanese patients with DDH and 2044 healthy individuals was performed. As a replicate, GWAS was also conducted on the UK Biobank data with 3315 cases and matched 74,038 controls. Gene set enrichment analyses (GSEAs) of both the genetics and transcriptome of DDH were performed. Transcriptome analysis of cartilage specimens from DDH-associated osteoarthritis and femoral neck fractures was performed as a control. Most of the lead variants were very low-frequency ones in the UK, and variants in the Japanese GWAS could not be replicated with the UK GWAS. We assigned DDH-related candidate variants to 42 and 81 genes from the Japanese and UK GWASs, respectively, using functional mapping and annotation. GSEA of gene ontology, disease ontology, and canonical pathways identified the most enriched pathway to be the ferroptosis signaling pathway, both in the Japanese gene set as well as the Japanese and UK merged set. Transcriptome GSEA also identified significant downregulation of genes in the ferroptosis signaling pathway. Thus, the ferroptosis signaling pathway may be associated with the pathogenic mechanism of DDH.

Design and implementation of a hybrid cloud system for large-scale human genomic research.

In the field of genomic medical research, the amount of large-scale information continues to increase due to advances in measurement technologies, such as high-performance sequencing and spatial omics, as well as the progress made in genomic cohort studies involving more than one million individuals. Therefore, researchers require more computational resources to analyze this information. Here, we introduce a hybrid cloud system consisting of an on-premise supercomputer, science cloud, and public cloud at the Kyoto University Center for Genomic Medicine in Japan as a solution. This system can flexibly handle various heterogeneous computational resource-demanding bioinformatics tools while scaling the computational capacity. In the hybrid cloud system, we demonstrate the way to properly perform joint genotyping of whole-genome sequencing data for a large population of 11,238, which can be a bottleneck in sequencing data analysis. This system can be one of the reference implementations when dealing with large amounts of genomic medical data in research centers and organizations.

CFH-CFHR1 hybrid genes in two cases of atypical hemolytic uremic syndrome.

Atypical hemolytic uremic syndrome (aHUS) is a rare complement-mediated disease that manifests as the triad of thrombotic microangiopathy. We identified two aHUS patients with neither anti-complement factor H (CFH) antibodies nor causative variants of seven aHUS-related genes (CFH, CFI, CFB, C3, MCP, THBD, and DGKE); however, their plasma showed increased levels of hemolysis by hemolytic assay, which strongly suggests CFH-related abnormalities. Using a copy number variation (CNV) analysis of the CFH/CFHR gene cluster, we identified CFH-CFHR1 hybrid genes in these patients. We verified the absence of aHUS-related abnormal CNVs of the CFH gene in control genomes of 2036 individuals in the general population, which suggests that pathogenicity is related to these hybrid genes. Our study emphasizes that, for patients suspected of having aHUS, it is important to perform an integrated analysis based on a clinical examination, functional analysis, and detailed genetic investigation.

Sex-Specific Differences in the Transcriptome of the Human Dorsolateral Prefrontal Cortex in Schizophrenia.

Schizophrenia presents clinical and biological differences between males and females. This study investigated transcriptional profiles in the dorsolateral prefrontal cortex (DLPFC) using postmortem data from the largest RNA-sequencing (RNA-seq) database on schizophrenic cases and controls. Data for 154 male and 113 female controls and 160 male and 93 female schizophrenic cases were obtained from the CommonMind Consortium. In the RNA-seq database, the principal component analysis showed that sex effects were small in schizophrenia. After we analyzed the impact of sex-specific differences on gene expression, the female group showed more significantly changed genes compared with the male group. Based on the gene ontology analysis, the female sex-specific genes that changed were overrepresented in the mitochondrion, ATP (phosphocreatine and adenosine triphosphate)-, and metal ion-binding relevant biological processes. An ingenuity pathway analysis revealed that the differentially expressed genes related to schizophrenia in the female group were involved in midbrain dopaminergic and γ-aminobutyric acid (GABA)-ergic neurons and microglia. We used methylated DNA-binding domain-sequencing analyses and microarray to investigate the DNA methylation that potentially impacts the sex differences in gene transcription using a maternal immune activation (MIA) murine model. Among the sex-specific positional genes related to schizophrenia in the PFC of female offspring from MIA, the changes in the methylation and transcriptional expression of loci ACSBG1 were validated in the females with schizophrenia in independent postmortem samples by real-time PCR and pyrosequencing. Our results reveal potential genetic risks in the DLPFC for the sex-dependent prevalence and symptomology of schizophrenia.

Secure secondary utilization system of genomic data using quantum secure cloud.

Secure storage and secondary use of individual human genome data is increasingly important for genome research and personalized medicine. Currently, it is necessary to store the whole genome sequencing information (FASTQ data), which enables detections of de novo mutations and structural variations in the analysis of hereditary diseases and cancer. Furthermore, bioinformatics tools to analyze FASTQ data are frequently updated to improve the precision and recall of detected variants. However, existing secure secondary use of data, such as multi-party computation or homomorphic encryption, can handle only a limited algorithms and usually requires huge computational resources. Here, we developed a high-performance one-stop system for large-scale genome data analysis with secure secondary use of the data by the data owner and multiple users with different levels of data access control. Our quantum secure cloud system is a distributed secure genomic data analysis system (DSGD) with a "trusted server" built on a quantum secure cloud, the information-theoretically secure Tokyo QKD Network. The trusted server will be capable of deploying and running a variety of sequencing analysis hardware, such as GPUs and FPGAs, as well as CPU-based software. We demonstrated that DSGD achieved comparable throughput with and without encryption on the trusted server Therefore, our system is ready to be installed at research institutes and hospitals that make diagnoses based on whole genome sequencing on a daily basis.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms.

BACKGROUND: Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. RESULTS: Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms. CONCLUSION: The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Genome-wide Survival Analysis for Macular Neovascularization Development in Central Serous Chorioretinopathy Revealed Shared Genetic Susceptibility with Polypoidal Choroidal Vasculopathy.

PURPOSE: To identify susceptibility genes for macular neovascularization (MNV) development in central serous chorioretinopathy (CSC). DESIGN: Genome-wide survival analysis using a longitudinal cohort study. PARTICIPANTS: We included 402 and 137 patients with CSC but without MNV at their first visit from the Kyoto CSC Cohort and Kobe CSC dataset, respectively. All patients underwent detailed ophthalmic examinations, including multimodal imaging, such as fundus autofluorescence, spectral-domain OCT, and fluorescein angiography/indocyanine green angiography or OCT angiography. METHODS: We conducted a genome-wide survival analysis using the Kyoto CSC Cohort. We applied the Cox proportional hazard model to adjust for age, sex, and the first principal component. Single nucleotide polymorphisms (SNPs) with P values < 1.0 × 10-5 were carried forward to the replication in the Kobe CSC dataset. Moreover, we evaluated the contribution of previously reported age-related macular degeneration (AMD) susceptibility loci. We used FUMA and ToppFun for the functional enrichment analysis. MAIN OUTCOME MEASURES: The association between SNPs and MNV development in patients with CSC. RESULTS: Rs370974631 near ARMS2 displayed a genome-wide significant association in the meta-analysis of discovery and replication result (hazard ratio [HR]meta, 3.63; Pmeta = 5.76 × 10-9). Among previously reported AMD susceptibility loci, we additionally identified CFH rs800292 (HR, 0.39, P = 2.55 × 10-4), COL4A3 rs4276018 (HR, 0.26, P = 1.56 × 10-3), and B3GALTL rs9564692 (HR, 0.56, P = 8.30 × 10-3) as susceptibility loci for MNV development in CSC. The functional enrichment analysis revealed significant enrichment of 8 pathways (GO:0051561, GO:0036444, GO:0008282, GO:1990246, GO:0015272, GO:0030955, GO:0031420, and GO:0005242) related to ion transport. CONCLUSIONS: ARMS2, CFH, COL4A3, and B3GALTL were identified as susceptibility genes for MNV development in CSC. These 4 genes are known as susceptibility genes for AMD, whereas COL4A3 and B3GALTL were previously reported to be polypoidal choroidal vasculopathy (PCV)-specific susceptibility genes. Our findings revealed the shared genetic susceptibility between PCV and MNV secondary to CSC.

A single-nucleotide-polymorphism in the 5'-flanking region of MSX1 gene as a predictive marker candidate for platinum-based therapy of esophageal carcinoma.

BACKGROUND: Platinum derivatives are important treatment options for patients with esophageal carcinoma (EC), and a predictive marker for platinum-based therapy is needed for precision medicine. PATIENTS AND METHODS: This study contained two cohorts consisting of EC patients treated using platinum-based chemoradiation therapy (CRT) as the first-line and another external cohort of nationwide clinicogenomic data from the BioBank Japan (BBJ). RESULTS: Genome-wide association study (GWAS) of therapeutic outcomes, refractory disease or not, following platinum-based CRT as first-line in 94 patients in the first cohort suggested the association of 89 SNPs using p < 0.0001. The top 10 SNPs selected from each chromosomal region by odds ratio were evaluated for progression-free survival (PFS) and overall survival (OS) hazard ratios in the first cohort, resulting in four candidates (p < 0.0025). The four selected candidates were re-evaluated in another cohort of 24 EC patients, which included patients prospectively enrolled in this study to fulfill the sample size statistically suggested by the results of the first cohort, and of the four, only rs3815544 was replicated (p < 0.0125). Furthermore, this candidate genotype of rs3815544 proceeded to the re-evaluation study in an external cohort consisting of EC patients treated with platinum derivatives and/or by radiation therapy as the first-line treatment in BBJ, which confirmed that the alternative allele (G) of rs3815544 was statistically associated with non-response (SD or PD) to platinum-based therapy in EC patients (odds ratio = 1.801, p = 0.048). The methylation QTL database as well as online clinicogenomic databases suggested that the region including rs3815544 may regulate MSX1 expression through CpG methylation, and this down-regulation was statistically associated with poor prognosis after platinum-based therapies for EC. CONCLUSION: rs3815544 is a novel candidate predictive marker for platinum-based EC therapy.

Antiemetic effects of baclofen in a shrew model of postoperative nausea and vomiting: Whole-transcriptome analysis in the nucleus of the solitary tract.

AIMS: The molecular genetic mechanisms underlying postoperative nausea and vomiting (PONV) in the brain have not been fully elucidated. This study aimed to determine the changes in whole transcriptome in the nucleus of the solitary tract (NTS) in an animal model of PONV, to screen a drug candidate and to elucidate the molecular genetic mechanisms of PONV development. METHODS: Twenty-one female musk shrews were assigned into three groups: the Surgery group (shrew PONV model, n = 9), the Sham group (n = 6), and the Naïve group (n = 6). In behavioral studies, the main outcome was the number of emetic episodes. In genetic experiments, changes in the transcriptome in the NTS were measured. In a separate study, 12 shrews were used to verify the candidate mechanism underlying PONV. RESULTS: A median of six emetic episodes occurred in both the Sham and Surgery groups. Whole-transcriptome analysis indicated the inhibition of the GABAB receptor-mediated signaling pathway in the PONV model. Baclofen (GABAB receptor agonist) administration eliminated emetic behaviors in the shrew PONV model. CONCLUSIONS: Our findings suggest that the GABAB receptor-mediated signaling pathway is involved in emesis and that baclofen may be a novel therapeutic or prophylactic agent for PONV.

Crohn's Disease and Early Exposure to Thiopurines are Independent Risk Factors for Mosaic Chromosomal Alterations in Patients with Inflammatory Bowel Diseases.

BACKGROUND AND AIMS: Mosaic chromosomal alterations [mCAs] increase the risk for haematopoietic malignancies and may be risk factors for several other diseases. Inflammatory bowel diseases [IBDs], including Crohn's disease [CD] and ulcerative colitis [UC], are associated with mCAs, and patients may be at risk for haematopoietic malignancy development and/or modification of IBD phenotypes. In the present study, we screened patients with IBD for the presence of mCAs and explored the possible pathophysiological and genetic risk factors for mCAs. METHODS: We analysed mCAs in peripheral blood from 3339 patients with IBD and investigated the clinical and genetic risk factors for mCAs. RESULTS: CD and exposure to thiopurines before the age of 20 years were identified as novel independent risk factors for mCAs [odds ratio = 2.15 and 5.68, p = 1.17e-2 and 1.60e-3, respectively]. In contrast, there were no significant associations of disease duration, anti-tumour necrosis factor alpha antibodies, or other clinical factors with mCAs. Gene ontology enrichment analysis revealed that genes specifically located in the mCAs in patients with CD were significantly associated with factors related to mucosal immune responses. A genome-wide association study revealed that ERBIN, CD96, and AC068672.2 were significantly associated with mCAs in patients with CD [p = 1.56e-8, 1.65e-8, and 4.92e-8, respectively]. CONCLUSIONS: The difference in mCAs between patients with CD and UC supports the higher incidence of haematopoietic malignancies in CD. Caution should be exercised when using thiopurines in young patients with IBD, particularly CD, in light of possible chromosomal alterations.