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.

The results of genetic analysis and clinical outcomes after stent deployment in adult patients with isolated peripheral pulmonary artery stenosis

Peripheral pulmonary arterial stenosis (PPAS) is commonly observed in paediatric patients with associated congenital heart disease or systemic congenital syndromes, such as tetralogy of Fallot and Williams syndrome [1]. Meanwhile, isolated PPAS without other structural heart disease or systemic congenital syndromes is primarily reported in adulthood [2]. Recent genetic analyses have revealed that four cases of isolated PPAS in adulthood were homozygous for the ring finger protein 213 (RNF213) p.Arg4810Lys variant (also referred to as G to A alteration of rs112735431 and NM_001256071.3:c.14429G>A) [3]. RNF213 p.Arg4810Lys is a genetic determinant in moyamoya disease (MMD) in East Asians [4]. This genetic variant is unique to East Asians and is reported to have a frequency of approximately 0.8% in the general Japanese population [5, 6]. However, its prevalence in adult patients with isolated PPAS remains unknown.

The RNF213 p.Arg4810Lys mutation is closely linked to the development of PPAS. Stent deployment might be a promising treatment option, but further investigation would be warranted. https://bit.ly/40t8Gmz

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.

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.

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 (log₁₀AFP) 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.

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.

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.

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).

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.

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 GABA_B receptor‐mediated signaling pathway in the PONV model. Baclofen (GABA_B receptor agonist) administration eliminated emetic behaviors in the shrew PONV model.

Conclusions

Our findings suggest that the GABA_B 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.

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 ophthalmologic examinations, including multimodal imaging, such as fundus autofluorescence, spectral-domain optical coherence tomography, and fluorescein angiography/indocyanine green angiography and/or optimal coherence tomography 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; P_meta = 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 eight 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. The aforementioned four 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.

Genetic Background of Mesalamine-induced Fever and Diarrhea in Japanese Patients with Inflammatory Bowel Disease

BACKGROUND

Some patients with inflammatory bowel disease (IBD) who were under mesalamine treatment develop adverse reactions called "mesalamine allergy," which includes high fever and worsening diarrhea. Currently, there is no method to predict mesalamine allergy. Pharmacogenomic approaches may help identify these patients. Here we analyzed the genetic background of mesalamine intolerance in the first genome-wide association study of Japanese patients with IBD.

METHODS

Two independent pharmacogenetic IBD cohorts were analyzed: the MENDEL (n = 1523; as a discovery set) and the Tohoku (n = 788; as a replication set) cohorts. Genome-wide association studies were performed in each population, followed by a meta-analysis. In addition, we constructed a polygenic risk score model and combined genetic and clinical factors to model mesalamine intolerance.

RESULTS

In the combined cohort, mesalamine-induced fever and/or diarrhea was significantly more frequent in ulcerative colitis vs Crohn's disease. The genome-wide association studies and meta-analysis identified one significant association between rs144384547 (upstream of RGS17) and mesalamine-induced fever and diarrhea (P = 7.21e-09; odds ratio = 11.2). The estimated heritability of mesalamine allergy was 25.4%, suggesting a significant correlation with the genetic background. Furthermore, a polygenic risk score model was built to predict mesalamine allergy (P = 2.95e-2). The combined genetic/clinical prediction model yielded a higher area under the curve than did the polygenic risk score or clinical model alone (area under the curve, 0.89; sensitivity, 71.4%; specificity, 90.8%).

CONCLUSIONS

Mesalamine allergy was more common in ulcerative colitis than in Crohn's disease. We identified a novel genetic association with and developed a combined clinical/genetic model for this adverse event.

Machine learning approaches to predict gestational age in normal and complicated pregnancies via urinary metabolomics analysis

The elucidation of dynamic metabolomic changes during gestation is particularly important for the development of methods to evaluate pregnancy status or achieve earlier detection of pregnancy-related complications. Some studies have constructed models to evaluate pregnancy status and predict gestational age using omics data from blood biospecimens; however, less invasive methods are desired. Here we propose a model to predict gestational age, using urinary metabolite information. In our prospective cohort study, we collected 2741 urine samples from 187 healthy pregnant women, 23 patients with hypertensive disorders of pregnancy, and 14 patients with spontaneous preterm birth. Using gas chromatography-tandem mass spectrometry, we identified 184 urinary metabolites that showed dynamic systematic changes in healthy pregnant women according to gestational age. A model to predict gestational age during normal pregnancy progression was constructed; the correlation coefficient between actual and predicted weeks of gestation was 0.86. The predicted gestational ages of cases with hypertensive disorders of pregnancy exhibited significant progression, compared with actual gestational ages. This is the first study to predict gestational age in normal and complicated pregnancies by using urinary metabolite information. Minimally invasive urinary metabolomics might facilitate changes in the prediction of gestational age in various clinical settings.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs

BACKGROUNDS & AIMS

Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening.

METHODS

We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts.

RESULTS

We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders.

CONCLUSIONS

This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders.

LAY SUMMARY

Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis

BACKGROUND & AIMS

Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease.

METHODS

We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals).

RESULTS

Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively).

CONCLUSIONS

This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

rs1944919 on chromosome 11q23.1 and its effector genes COLCA1/COLCA2 confer susceptibility to primary biliary cholangitis

Primary biliary cholangitis (PBC) is a chronic, progressive cholestatic liver disease in which intrahepatic bile ducts are destroyed by an autoimmune reaction. Our previous genome-wide association study (GWAS) identified chromosome 11q23.1 as a susceptibility gene locus for PBC in the Japanese population. Here, high-density association mapping based on single nucleotide polymorphism (SNP) imputation and in silico/in vitro functional analyses identified rs1944919 as the primary functional variant. Expression-quantitative trait loci analyses showed that the PBC susceptibility allele of rs1944919 was significantly associated with increased COLCA1/COLCA2 expression levels. Additionally, the effects of rs1944919 on COLCA1/COLCA2 expression levels were confirmed using genotype knock-in versions of cell lines constructed using the CRISPR/Cas9 system and differed between rs1944919-G/G clones and -T/T clones. To our knowledge, this is the first study to demonstrate the contribution of COLCA1/COLCA2 to PBC susceptibility.

High-speed parameter search of dynamic biological pathways from time-course transcriptomic profiles using high-level Petri net

Dynamic simulation promises a deeper understanding of complex molecular mechanisms of biological pathways. How to determine the reaction kinetic parameters which govern the simulation results is still an open question in the field of systems biology. (1) Background: To execute simulation experiments, it is an essential first step to search effective values of model parameters. The complexity of biological systems and the experimental measurement technology severely limit the acquirement of accurate kinetic parameters. Previously proposed genomic data assimilation (GDA) approach enables users to handle parameter estimation using time-course information. However, it highly depends on successive time points and costs massive computational resource; (2) Methods: To address this problem, we present a new high-speed parameter search method for estimating the kinetic parameters of quantitative biological pathways using time-course transcriptomic profiles. The key idea of our method is to interactively prune the search space by introducing Probabilistic Linear-time Temporal Logic (PLTL) based model checking into GDA. (3) Results and conclusion: We demonstrated the effectiveness of our method by comparing with GDA on Mus musculus transcription circuits modelled by hybrid functional Petri net with extension. As a result, our method works faster and more accurate than GDA for both time-course datasets with dense and sparse observed values.

Practical guide for managing large-scale human genome data in research

Studies in human genetics deal with a plethora of human genome sequencing data that are generated from specimens as well as available on public domains. With the development of various bioinformatics applications, maintaining the productivity of research, managing human genome data, and analyzing downstream data is essential. This review aims to guide struggling researchers to process and analyze these large-scale genomic data to extract relevant information for improved downstream analyses. Here, we discuss worldwide human genome projects that could be integrated into any data for improved analysis. Obtaining human whole-genome sequencing data from both data stores and processes is costly; therefore, we focus on the development of data format and software that manipulate whole-genome sequencing. Once the sequencing is complete and its format and data processing tools are selected, a computational platform is required. For the platform, we describe a multi-cloud strategy that balances between cost, performance, and customizability. A good quality published research relies on data reproducibility to ensure quality results, reusability for applications to other datasets, as well as scalability for the future increase of datasets. To solve these, we describe several key technologies developed in computer science, including workflow engine. We also discuss the ethical guidelines inevitable for human genomic data analysis that differ from model organisms. Finally, the future ideal perspective of data processing and analysis is summarized.

Digenic mutations in ALDH2 and ADH5 impair formaldehyde clearance and cause a multisystem disorder, AMeD syndrome

Rs671 in the aldehyde dehydrogenase 2 gene (ALDH2) is the cause of Asian alcohol flushing response after drinking. ALDH2 detoxifies endogenous aldehydes, which are the major source of DNA damage repaired by the Fanconi anemia pathway. Here, we show that the rs671 defective allele in combination with mutations in the alcohol dehydrogenase 5 gene, which encodes formaldehyde dehydrogenase (ADH5^FDH ), causes a previously unidentified disorder, AMeD (aplastic anemia, mental retardation, and dwarfism) syndrome. Cellular studies revealed that a decrease in the formaldehyde tolerance underlies a loss of differentiation and proliferation capacity of hematopoietic stem cells. Moreover, Adh5^-/-Aldh2 ^E506K/E506K double-deficient mice recapitulated key clinical features of AMeDS, showing short life span, dwarfism, and hematopoietic failure. Collectively, our results suggest that the combined deficiency of formaldehyde clearance mechanisms leads to the complex clinical features due to overload of formaldehyde-induced DNA damage, thereby saturation of DNA repair processes.

Identification of the Novel Variants in Patients With Chronic Thromboembolic Pulmonary Hypertension

Background Although chronic thromboembolic pulmonary hypertension (CTEPH) and acute pulmonary embolism (APE) share some clinical manifestations, a limited proportion of patients with CTEPH have a history of APE. Moreover, in histopathologic studies, it has been revealed that pulmonary vasculature lesions similar to pulmonary arterial hypertension existed in patients with CTEPH. Thus, it remains unknown whether these 3 disorders also share genetic backgrounds. Methods and Results Whole exome screening was performed with DNA isolated from 51 unrelated patients with CTEPH of Japanese ancestry. The frequency of genetic variants associated with pulmonary arterial hypertension or APE in patients with CTEPH was compared with those in the integrative Japanese Genome Variation Database 3.5KJPN. Whole exome screening analysis showed 17 049 nonsynonymous variants in patients with CTEPH. Although we found 6 nonsynonymous variants that are associated with APE in patients with CTEPH, there was no nonsynonymous variant associated with pulmonary arterial hypertension. Patients with CTEPH with a history of APE had nonsynonymous variants of F5, which encodes factor V. In contrast, patients with CTEPH without a history of APE had a nonsynonymous variant of THBD, which encodes thrombomodulin. Moreover, thrombin-activatable fibrinolysis inhibitor, which is one of the pathogenic proteins in CTEPH, was significantly more activated in those who had the variants of THBD compared with those without it. Conclusions These results provide the first evidence that patients with CTEPH have some variants associated with APE, regardless of the presence or absence of a history of APE. Furthermore, the variants might be different between patients with CTEPH with and without a history of APE.

Direct reprogramming of human umbilical vein- and peripheral blood-derived endothelial cells into hepatic progenitor cells

Recent advances have enabled the direct induction of human tissue-specific stem and progenitor cells from differentiated somatic cells. However, it is not known whether human hepatic progenitor cells (hHepPCs) can be generated from other cell types by direct lineage reprogramming with defined transcription factors. Here, we show that a set of three transcription factors, FOXA3, HNF1A, and HNF6, can induce human umbilical vein endothelial cells to directly acquire the properties of hHepPCs. These induced hHepPCs (hiHepPCs) propagate in long-term monolayer culture and differentiate into functional hepatocytes and cholangiocytes by forming cell aggregates and cystic epithelial spheroids, respectively, under three-dimensional culture conditions. After transplantation, hiHepPC-derived hepatocytes and cholangiocytes reconstitute damaged liver tissues and support hepatic function. The defined transcription factors also induce hiHepPCs from endothelial cells circulating in adult human peripheral blood. These expandable and bipotential hiHepPCs may be useful in the study and treatment of human liver diseases.

The conditions to induce human hepatic progenitor cells from other cell types are unclear. Here, the authors reprogram human endothelial cells to hepatic progenitor cells by expressing FOXA3, HNF1A and HNF6, capable of giving rise to hepatocytes and cholangiocytes that reconstitute damaged liver tissues on transplantation.

A multi-ethnic meta-analysis identifies novel genes, including ACSL5, associated with amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating progressive motor neuron disease that affects people of all ethnicities. Approximately 90% of ALS cases are sporadic and thought to have multifactorial pathogenesis. To understand the genetics of sporadic ALS, we conducted a genome-wide association study using 1,173 sporadic ALS cases and 8,925 controls in a Japanese population. A combined meta-analysis of our Japanese cohort with individuals of European ancestry revealed a significant association at the ACSL5 locus (top SNP p = 2.97 × 10⁻⁸). We validated the association with ACSL5 in a replication study with a Chinese population and an independent Japanese population (1941 ALS cases, 3821 controls; top SNP p = 1.82 × 10⁻⁴). In the combined meta-analysis, the intronic ACSL5 SNP rs3736947 showed the strongest association (p = 7.81 × 10⁻¹¹). Using a gene-based analysis of the full multi-ethnic dataset, we uncovered additional genes significantly associated with ALS: ERGIC1, RAPGEF5, FNBP1, and ATXN3. These results advance our understanding of the genetic basis of sporadic ALS.

Gen Sobue, Masao Nagasaki and colleagues report a genome-wide association study for amyotrophic lateral sclerosis (ALS) in a large, multi-ethnic cohort comprising Japanese, Chinese, and European ancestry populations. They find a significant association to variants within the ACSL5 gene and identify novel associations with 4 additional genes using a gene-based approach.

Long non-coding RNA MIR4300HG polymorphisms are associated with postoperative nausea and vomiting: a genome-wide association study

Background

Genetic factors such as single-nucleotide polymorphisms (SNPs) play a key role in the development of postoperative nausea and vomiting (PONV). However, previous findings are not widely applicable to different populations because of population-specific genetic variation. We developed a Japanese-specific DNA microarray for high-throughput genotyping. The aim of the current study was to identify SNPs associated with PONV on a genome-wide scale using this microarray in a sample of Japanese surgical patients.

Methods

Associations between 659,636 SNPs and the incidence of PONV 24 h after surgery in a limited sample of 24 female patients were assessed using the microarray. After imputation of genotypes at 24,330,529 SNPs, 78 SNPs were found to be associated with the incidence of PONV. We chose 4 of the 78 SNPs to focus on by in silico functional annotation. Finally, we genotyped these 4 candidate SNPs in 255 patients using real-time PCR to verify association with the incidence of PONV.

Results

The T > C variant of rs11232965 in the long non-coding RNA MIR4300HG was significantly associated with reduced incidence of PONV among genotypes and between alleles (p = 0.01 and 0.007).

Conclusions

We identified a novel SNP (rs11232965) in the long non-coding RNA MIR4300HG that is associated with PONV. The rs11232965-SNP variant (T > C) is protective against the incidence of PONV.

Trial registration

This study was registered at the UMIN Clinical Trials Registry (Identifier: UMIN000022903, date of registration: June 27, 2016, retrospectively registered.

The Dynamics of Transcriptional Activation by Hepatic Reprogramming Factors

Specific combinations of two transcription factors (Hnf4α plus Foxa1, Foxa2, or Foxa3) can induce direct conversion of mouse fibroblasts into hepatocyte-like cells. However, the molecular mechanisms underlying hepatic reprogramming are largely unknown. Here, we show that the Foxa protein family members and Hnf4α sequentially and cooperatively bind to chromatin to activate liver-specific gene expression. Although all Foxa proteins bind to and open regions of closed chromatin as pioneer factors, Foxa3 has the unique potential of transferring from the distal to proximal regions of the transcription start site of target genes, binding RNA polymerase II, and co-traversing target genes. These distinctive characteristics of Foxa3 are essential for inducing the hepatic fate in fibroblasts. Similar functional coupling of transcription factors to RNA polymerase II may occur in other contexts whereby transcriptional activation can induce cell differentiation.

Genetic Analysis of Ulcerative Colitis in Japanese Individuals Using Population-specific SNP Array

BACKGROUND

To clarify the genetic background of ulcerative colitis (UC) in the Japanese population, we conducted a genome-wide association study (GWAS) using a population-specific single nucleotide polymorphism (SNP) array.

METHODS

We performed a GWAS and replication study including 1676 UC patients and 2381 healthy controls. The probability of colectomy was compared between genotypes of rs117506082, the top hit SNP at HLA loci, by the Kaplan-Meier method. We studied serum expression of miR-622, a newly identified candidate gene, from 32 UC patients and 8 healthy controls by quantitative reverse-transcription polymerase chain reaction.

RESULTS

In the GWAS, only the HLA loci showed genome-wide significant associations with UC (rs117506082, P = 6.69E-28). Seven nominally significant regions included 2 known loci, IL23R (rs76418789, P = 6.29E-7) and IRF8 (rs16940202, P = 1.03E-6), and 5 novel loci: MIR622 (rs9560575, P = 8.23E-7), 14q31 (rs117618617, P = 1.53E-6), KAT6B (rs12260609, P = 1.81E-6), PAX3-CCDC140-SGPP2 (rs7589797, P = 2.87E-6), and KCNA2 (rs118020656, P = 4.01E-6). Combined analysis revealed that IL23R p.G149R (rs76418789, P = 9.03E-11; odds ratio [OR], 0.51) had genome-wide significant association with UC. Patients with GG genotype of rs117506082 had a significantly lower probability of total colectomy than those with the GA+AA genotype (P = 1.72E-2). Serum expression of miR-622 in patients with inactive UC tended to be higher than in healthy controls and patients with active UC (inactive UC vs healthy controls, P = 3.03E-02; inactive UC vs active UC, P = 6.44E-02).

CONCLUSIONS

IL23R p.G149R is a susceptibility locus for UC in Japanese individuals. The GG genotype of rs117506082 at HLA loci may predict a better clinical course.

Common risk variants in NPHS1 and TNFSF15 are associated with childhood steroid-sensitive nephrotic syndrome

To understand the genetics of steroid-sensitive nephrotic syndrome (SSNS), we conducted a genome-wide association study in 987 childhood SSNS patients and 3,206 healthy controls with Japanese ancestry. Beyond known associations in the HLA-DR/DQ region, common variants in NPHS1-KIRREL2 (rs56117924, P=4.94E-20, odds ratio (OR) =1.90) and TNFSF15 (rs6478109, P=2.54E-8, OR=0.72) regions achieved genome-wide significance and were replicated in Korean, South Asian and African populations. Trans-ethnic meta-analyses including Japanese, Korean, South Asian, African, European, Hispanic and Maghrebian populations confirmed the significant associations of variants in NPHS1-KIRREL2 (P_meta=6.71E-28, OR=1.88) and TNFSF15 (P_meta=5.40E-11, OR=1.33) loci. Analysis of the NPHS1 risk alleles with glomerular NPHS1 mRNA expression from the same person revealed allele specific expression with significantly lower expression of the transcript derived from the risk haplotype (Wilcox test p=9.3E-4). Because rare pathogenic variants in NPHS1 cause congenital nephrotic syndrome of the Finnish type (CNSF), the present study provides further evidence that variation along the allele frequency spectrum in the same gene can cause or contribute to both a rare monogenic disease (CNSF) and a more complex, polygenic disease (SSNS).

Prevalence of disease-causing genes in Japanese patients with BRCA1/2-wildtype hereditary breast and ovarian cancer syndrome

Panel sequencing of susceptibility genes for hereditary breast and ovarian cancer (HBOC) syndrome has uncovered numerous germline variants; however, their pathogenic relevance and ethnic diversity remain unclear. Here, we examined the prevalence of germline variants among 568 Japanese patients with BRCA1/2-wildtype HBOC syndrome and a strong family history. Pathogenic or likely pathogenic variants were identified on 12 causal genes for 37 cases (6.5%), with recurrence for 4 SNVs/indels and 1 CNV. Comparisons with non-cancer east-Asian populations and European familial breast cancer cohorts revealed significant enrichment of PALB2, BARD1, and BLM mutations. Younger onset was associated with but not predictive of these mutations. Significant somatic loss-of-function alterations were confirmed on the wildtype alleles of genes with germline mutations, including PALB2 additional somatic truncations. This study highlights Japanese-associated germline mutations among patients with BRCA1/2 wildtype HBOC syndrome and a strong family history, and provides evidence for the medical care of this high-risk population.

Genome-Wide Association Study of Lean Nonalcoholic Fatty Liver Disease Suggests Human Leukocyte Antigen as a Novel Candidate Locus

Nonalcoholic fatty liver disease (NAFLD) is supposed to manifest its metabolic phenotype in the liver, but it is common to have lean individuals diagnosed with NAFLD, known as lean NAFLD. We conducted a two-stage analysis to identify NAFLD-associated loci in Japanese patients. In stage I, 275 metabolically healthy normal-weight patients with NAFLD were compared with 1,411 non-NAFLD controls adjusted for age, sex, and alcohol consumption by a genome-wide association study (GWAS). In stage II, human leukocyte antigen (HLA) in chromosome 6 (chr6) (P = 6.73E-08), microRNA (MIR) MIR548F3 in chr7 (P = 4.25E-07), myosin light chain 2 (MYL2) in chr12 (P = 4.39E-07), and glycoprotein precursor (GPC)6 in chr13 (P = 5.43E-07), as suggested by the GWAS, were assessed by single nucleotide polymorphism (SNP) association analysis of whole NAFLD against non-NAFLD in 9,726 members of the general population. A minor allele of the secondary lead SNP in chr6, rs2076529, was significantly associated (odds ratio [OR], 1.19; 95% confidence interval [CI], 1.11-1.28; P = 2.10E-06) and the lead SNP in chr7 was weakly associated (OR 1.15; 95% CI, 1.04-1.27; P = 6.19E-03) with increased NAFLD risk. Imputation-based typing of HLA showed a significant difference in the distribution of HLA-B, HLA-DR-beta chain 1 (DRB1), and HLA-DQ-beta chain 1 (DQB1) alleles in lean NAFLD GWAS. Next-generation sequence-based typing of HLA in 5,649 members of the general population replicated the significant difference of HLA-B allele distribution and the significant increase of the HLA-B*54:01 allele in whole NAFLD. Fecal metagenomic analysis of 3,420 members of the general population showed significant dissimilarity in beta-diversity analysis of rs2076529 and HLA-B*54:01 allele carriers from noncarriers. Veillonellaceae was increased but Verrucomicrobia was decreased in rs2076529 minor allele and HLA-B*54:01 allele carriers as in NAFLD. Conclusion: HLA was identified as a novel locus associated with NAFLD susceptibility, which might be affected by the alteration of gut microbiota.

Variants That Affect Function of Calcium Channel TRPV6 Are Associated With Early-Onset Chronic Pancreatitis

BACKGROUND & AIMS

Changes in pancreatic calcium levels affect secretion and might be involved in development of chronic pancreatitis (CP). We investigated the association of CP with the transient receptor potential cation channel subfamily V member 6 gene (TRPV6), which encodes a Ca²⁺-selective ion channel, in an international cohort of patients and in mice.

METHODS

We performed whole-exome DNA sequencing from a patient with idiopathic CP and from his parents, who did not have CP. We validated our findings by sequencing DNA from 300 patients with CP (not associated with alcohol consumption) and 1070 persons from the general population in Japan (control individuals). In replication studies, we sequenced DNA from patients with early-onset CP (20 years or younger) not associated with alcohol consumption from France (n = 470) and Germany (n = 410). We expressed TRPV6 variants in HEK293 cells and measured their activity using Ca²⁺ imaging assays. CP was induced by repeated injections of cerulein in TRPV6^mut/mut mice.

RESULTS

We identified the variants c.629C>T (p.A210V) and c.970G>A (p.D324N) in TRPV6 in the index patient. Variants that affected function of the TRPV6 product were found in 13 of 300 patients (4.3%) and 1 of 1070 control individuals (0.1%) from Japan (odds ratio [OR], 48.4; 95% confidence interval [CI], 6.3-371.7; P = 2.4 × 10^-8). Twelve of 124 patients (9.7%) with early-onset CP had such variants. In the replication set from Europe, 18 patients with CP (2.0%) carried variants that affected the function of the TRPV6 product compared with 0 control individuals (P = 6.2 × 10^-8). Variants that did not affect the function of the TRPV6 product (p.I223T and p.D324N) were overrepresented in Japanese patients vs control individuals (OR, 10.9; 95% CI, 4.5-25.9; P = 7.4 × 10^-9 for p.I223T and P = .01 for p.D324N), whereas the p.L299Q was overrepresented in European patients vs control individuals (OR, 3.0; 95% CI, 1.9-4.8; P = 1.2 × 10^-5). TRPV6^mut/mut mice given cerulein developed more severe pancreatitis than control mice, as shown by increased levels of pancreatic enzymes, histologic alterations, and pancreatic fibrosis.

CONCLUSIONS

We found that patients with early-onset CP not associated with alcohol consumption carry variants in TRPV6 that affect the function of its product, perhaps by altering Ca²⁺ balance in pancreatic cells. TRPV6 regulates Ca²⁺ homeostasis and pancreatic inflammation.

Integrated GWAS and mRNA Microarray Analysis Identified IFNG and CD40L as the Central Upstream Regulators in Primary Biliary Cholangitis

Genome‐wide association studies (GWASs) in European and East Asian populations have identified more than 40 disease‐susceptibility genes in primary biliary cholangitis (PBC). The aim of this study is to computationally identify disease pathways, upstream regulators, and therapeutic targets in PBC through integrated GWAS and messenger RNA (mRNA) microarray analysis. Disease pathways and upstream regulators were analyzed with ingenuity pathway analysis in data set 1 for GWASs (1,920 patients with PBC and 1,770 controls), which included 261 annotated genes derived from 6,760 single‐nucleotide polymorphisms (P < 0.00001), and data set 2 for mRNA microarray analysis of liver biopsy specimens (36 patients with PBC and 5 normal controls), which included 1,574 genes with fold change >2 versus controls (P < 0.05). Hierarchical cluster analysis and categorization of cell type–specific genes were performed for data set 2. There were 27 genes, 10 pathways, and 149 upstream regulators that overlapped between data sets 1 and 2. All 10 pathways were immune‐related. The most significant common upstream regulators associated with PBC disease susceptibility identified were interferon‐gamma (IFNG) and CD40 ligand (CD40L). Hierarchical cluster analysis of data set 2 revealed two distinct groups of patients with PBC by disease activity. The most significant upstream regulators associated with disease activity were IFNG and CD40L. Several molecules expressed in B cells, T cells, Kupffer cells, and natural killer–like cells were identified as potential therapeutic targets in PBC with reference to a recently reported list of cell type–specific gene expression in the liver. Conclusion: Our integrated analysis using GWAS and mRNA microarray data sets predicted that IFNG and CD40L are the central upstream regulators in both disease susceptibility and activity of PBC and identified potential downstream therapeutic targets.

Exploring the Novel Susceptibility Gene Variants for Primary Open-Angle Glaucoma in East Asian Cohorts: The GLAU-GENDISK Study

Primary open-angle glaucoma (POAG) can develop even within normal ranges of intraocular pressure, and this type of glaucoma (so-called ‘normal-tension glaucoma [NTG]’) is highly prevalent in East Asia including Korea and Japan. We conducted exome chip analysis to identify low-frequency and rare variants associated with POAG from the primary cohort (309 POAG patients and 5,400 control, all Koreans). For replication, Korean (310 POAG patients and 5,612 controls) and Japanese (565 POAG patients and 1,104 controls) cohorts were further investigated by targeted genotyping. SNP rs116121322 in LRRC27 showed nominally significant association with POAG in the discovery cohort (OR = 29.85, P = 2E–06). This SNP was validated in the Korean replication cohort but only in the NTG subgroups (OR = 9.86, P = 0.007). Japanese replication cohort did not show significant association with POAG (P .00.44). However, the meta-analysis in the entire cohort revealed significant association of rs116121322 with POAG (OR_combined = 10.28, P_combined = 1.4E–07). The LRRC27 protein expression was confirmed from human trabecular meshwork cells. For gene-based testing, METTL20 showed a significant association in POAG (P_combined = 0.002) and in the subgroup of NTG (P_combined = 0.02), whereas ZNF677 were significantly associated with only in the subgroup of high-tension glaucoma (P_combined = 1.5E–06). Our findings may provide further genetic backgrounds into the pathogenesis of POAG, especially for the patients who have lower baseline intraocular pressures.

Study Profile of the Tohoku Medical Megabank Community-Based Cohort Study

Background

We established a community-based cohort study to assess the long-term impact of the Great East Japan Earthquake on disaster victims and gene-environment interactions on the incidence of major diseases, such as cancer and cardiovascular diseases.

Methods

We asked participants to join our cohort in the health check-up settings and assessment center based settings. Inclusion criteria were aged 20 years or over and living in Miyagi or Iwate Prefecture. We obtained information on lifestyle, effect of disaster, blood, and urine information (Type 1 survey), and some detailed measurements (Type 2 survey), such as carotid echography and calcaneal ultrasound bone mineral density. All participants agreed to measure genome information and to distribute their information widely.

Results

As a result, 87,865 gave their informed consent to join our study. Participation rate at health check-up site was about 70%. The participants in the Type 1 survey were more likely to have psychological distress than those in the Type 2 survey, and women were more likely to have psychological distress than men. Additionally, coastal residents were more likely to have higher degrees of psychological distress than inland residents, regardless of sex.

Conclusion

This cohort comprised a large sample size and it contains information on the natural disaster, genome information, and metabolome information. This cohort also had several detailed measurements. Using this cohort enabled us to clarify the long-term effect of the disaster and also to establish personalized prevention based on genome, metabolome, and other omics information.

Correction to: High-resolution melt analysis enables simple genotyping of complicated polymorphisms of codon 18 rendering the NUDT15 diplotype

The correct name of the last author should be.

High-resolution melt analysis enables simple genotyping of complicated polymorphisms of codon 18 rendering the NUDT15 diplotype

BACKGROUND

The genetic variants of NUDT15 have been verified to induce adverse events (AEs) of thiopurines. Codon 139 variants are frequently observed in Asians, while multiple variants are seen in codon 18 which also cause AEs including the European ancestry. The purpose of this study is to establish a technique capable of the simple genotyping of NUDT15 codon 18 and to evaluate its efficacy.

METHODS

A high-resolution melt (HRM) technique is performed to simply determine genotypes. The accuracy of HRM analysis was evaluated with DNAs from 1245 Japanese patients with inflammatory bowel diseases. Subsequently, another group of 572 patients was analyzed to verify the method. The diplotypes and the frequency of their AEs were estimated on the basis of codon 18 and 139 genotypes.

RESULTS

The HRM analysis enabled the correct identification of the three main genotypes, ref/ref, ref/ins, and ref/V18I, in 1236 of 1241 cases. All rare genotypes including ref/del were identified as the impossible-to-determine group, the proper diagnosis rate was 99.6%. In the verification test using other samples, the diagnosis rate was 99.7%. By estimating diplotypes using both codon 18 and 139 genotypes, 2.74% and 2.13% of Japanese patients with Arg/Arg and Arg/Cys of codon 139 have a lower enzymatic activity of NUDT15 and a higher risk for adverse responses than those estimated by codon 139 genotypes alone.

CONCLUSIONS

Our study showed that HRM method enables simple genotyping of complicated codon 18 variants essential to haplotype estimation of the NUDT15.

CD45+CD326+ Cells are Predictive of Poor Prognosis in Non-Small Cell Lung Cancer Patients

PURPOSE

The epithelial-to-mesenchymal transition, the major process by which some cancer cells convert from an epithelial phenotype to a mesenchymal one, has been suggested to drive chemo-resistance and/or metastasis in patients with cancer. However, only a few studies have demonstrated the presence of CD45/CD326 doubly-positive cells (CD45/CD326 DPC) in cancer. We deployed a combination of cell surface markers to elucidate the phenotypic heterogeneity in non-small cell lung cancer (NSCLC) cells and identified a new subpopulation that is doubly-positive for epithelial and non-epithelial cell-surface markers in both NSCLC cells and patients' malignant pleural effusions.

EXPERIMENTAL DESIGN

We procured a total of 39 patients' samples, solid fresh lung cancer tissues from 21 patients and malignant pleural effusion samples from 18 others, and used FACS and fluorescence microscopy to check their surface markers. We also examined the EGFR mutations in patients with known acquired EGFR mutations.

RESULTS

Our data revealed that 0.4% to 17.9% of the solid tumor tissue cells and a higher percentage of malignant pleural effusion cells harbored CD45/CD326 DPC expressing both epithelial and nonepithelial surface markers. We selected 3 EGFR mutation patients and genetically confirmed that the newly identified cell population really originated from cancer cells. We also found that higher proportions of CD45/CD326 DPC are significantly associated with poor prognosis.

CONCLUSIONS

In conclusion, varying percentages of CD45/CD326 DPC exist in both solid cancer tissue and malignant pleural effusion in patients with NSCLC. This CD45/CD326 doubly-positive subpopulation can be an important key to clinical management of patients with NSCLC.

A low-frequency IL4R locus variant in Japanese patients with intravenous immunoglobulin therapy-unresponsive Kawasaki disease

BACKGROUND

Kawasaki disease (KD) is a systemic vasculitis which may be associated with coronary artery aneurysms. A notable risk factor for the development of coronary artery aneurysms is resistance to intravenous immunoglobulin (IVIG) therapy, which comprises standard treatment for the acute phase of KD. The cause of IVIG resistance in KD is largely unknown; however, the contribution of genetic factors, especially variants in immune-related genes, has been suspected.

METHODS

To explore genetic variants related to IVIG-unresponsiveness, we designated KD patients who did not respond to both first and second courses of IVIG therapy as IVIG-unresponsive patients. Using genomic DNA from 30 IVIG-unresponsive KD patients, we performed pooled genome sequencing targeting 39 immune-related cytokine receptor genes.

RESULTS

The single nucleotide variant (SNV), rs563535954 (located in the IL4R locus), was concentrated in IVIG-unresponsive KD patients. Individual genotyping showed that the minor allele of rs563535954 was present in 4/33 patients with IVIG-unresponsive KD, compared with 20/1063 individuals in the Japanese genome variation database (odds ratio = 7.19, 95% confidence interval 2.43-21.47). Furthermore, the minor allele of rs563535954 was absent in 42 KD patients who responded to IVIG treatment (P = 0.0337), indicating that a low-frequency variant, rs563535954, is associated with IVIG-unresponsiveness in KD patients. Although rs563535954 is located in the 3'-untranslated region of IL4R, there was no alternation in IL4R expression associated with the mior allele of rs563535954. However, IVIG-unresponsive patients that exhibited the minor allele of rs563535954 tended to be classified into the low-risk group (based on previously reported risk scores) for prediction of IVIG-resistance. Therefore, IVIG-unresponsiveness associated with the minor allele of rs563535954 might differ from IVIG-unresponsiveness associated with previous risk factors used to evaluate IVIG-unresponsiveness in KD.

CONCLUSION

These findings suggest that the SNV rs563535954 could serve as a predictive indicator of IVIG-unresponsiveness, thereby improving the sensitivity of risk scoring systems, and may aid in prevention of coronary artery lesions in KD patients.

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

Researchers in Japan have assembled a Japanese reference genome, which includes sequences missing from the international reference genome, as well as others specific to East Asian populations. A team led by Masao Nagasaki and Masayuki Yamamoto sequenced a Japanese individual using a method, which produces longer sequences than previous technologies. Using this approach, they identified thousands of sequences spanning 2.5 million bases, which were absent in the international reference genome. Many of these were sequences able to move within the genome. They showed that the majority of these sequences are also present in early humans and chimpanzees, demonstrating that their absence from the current reference is due to deletions or limitations of earlier sequencing methodologies. In addition to providing a population-specific reference, these findings demonstrate the importance of continually improving the international reference genome.

A Genome-wide Association Study Identifying RAP1A as a Novel Susceptibility Gene for Crohn's Disease in Japanese Individuals

BACKGROUND AND AIMS

Genome-wide association studies [GWASs] of European populations have identified numerous susceptibility loci for Crohn's disease [CD]. Susceptibility genes differ by ethnicity, however, so GWASs specific for Asian populations are required. This study aimed to clarify the Japanese-specific genetic background for CD by a GWAS using the Japonica array [JPA] and subsequent imputation with the 1KJPN reference panel.

METHODS

Two independent Japanese case/control sets (Tohoku region [379 CD patients, 1621 controls] and Kyushu region [334 CD patients, 462 controls]) were included. GWASs were performed separately for each population, followed by a meta-analysis. Two additional replication sets [254 + 516 CD patients and 287 + 565 controls] were analysed for top hit single nucleotide polymorphisms [SNPs] from novel genomic regions.

RESULTS

Genotype data of 4 335 144 SNPs from 713 Japanese CD patients and 2083 controls were analysed. SNPs located in TNFSF15 (rs78898421, Pmeta = 2.59 × 10-26, odds ratio [OR] = 2.10), HLA-DQB1 [rs184950714, pmeta = 3.56 × 10-19, OR = 2.05], ZNF365, and 4p14 loci were significantly associated with CD in Japanese individuals. Replication analyses were performed for four novel candidate loci [p <1 × 10-6], and rs488200 located upstream of RAP1A was significantly associated with CD [pcombined = 4.36 × 10-8, OR = 1.31]. Transcriptome analysis of CD4+ effector memory T cells from lamina propria mononuclear cells of CD patients revealed a significant association of rs488200 with RAP1A expression.

CONCLUSIONS

RAP1A is a novel susceptibility locus for CD in the Japanese population.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Maternity Log study: a longitudinal lifelog monitoring and multiomics analysis for the early prediction of complicated pregnancy

PURPOSE

A prospective cohort study for pregnant women, the Maternity Log study, was designed to construct a time-course high-resolution reference catalogue of bioinformatic data in pregnancy and explore the associations between genomic and environmental factors and the onset of pregnancy complications, such as hypertensive disorders of pregnancy, gestational diabetes mellitus and preterm labour, using continuous lifestyle monitoring combined with multiomics data on the genome, transcriptome, proteome, metabolome and microbiome.

PARTICIPANTS

Pregnant women were recruited at the timing of first routine antenatal visits at Tohoku University Hospital, Sendai, Japan, between September 2015 and November 2016. Of the eligible women who were invited, 65.4% agreed to participate, and a total of 302 women were enrolled. The inclusion criteria were age ≥20 years and the ability to access the internet using a smartphone in the Japanese language.

FINDINGS TO DATE

Study participants uploaded daily general health information including quality of sleep, condition of bowel movements and the presence of nausea, pain and uterine contractions. Participants also collected physiological data, such as body weight, blood pressure, heart rate and body temperature, using multiple home healthcare devices. The mean upload rate for each lifelog item was ranging from 67.4% (fetal movement) to 85.3% (physical activity), and the total number of data points was over 6 million. Biospecimens, including maternal plasma, serum, urine, saliva, dental plaque and cord blood, were collected for multiomics analysis.

FUTURE PLANS

Lifelog and multiomics data will be used to construct a time-course high-resolution reference catalogue of pregnancy. The reference catalogue will allow us to discover relationships among multidimensional phenotypes and novel risk markers in pregnancy for the future personalised early prediction of pregnancy complications.