Exploring the genetic diversity of the Japanese population: Insights from a large-scale whole genome sequencing analysis.

The Japanese archipelago is a terminal location for human migration, and the contemporary Japanese people represent a unique population whose genomic diversity has been shaped by multiple migrations from Eurasia. We analyzed the genomic characteristics that define the genetic makeup of the modern Japanese population from a population genetics perspective from the genomic data of 9,287 samples obtained by high-coverage whole-genome sequencing (WGS) by the National Center Biobank Network. The dataset comprised populations from the Ryukyu Islands and other parts of the Japanese archipelago (Hondo). The Hondo population underwent two episodes of population decline during the Jomon period, corresponding to the Late Neolithic, and the Edo period, corresponding to the Early Modern era, while the Ryukyu population experienced a population decline during the shell midden period of the Late Neolithic in this region. Haplotype analysis suggested increased allele frequencies for genes related to alcohol and fatty acid metabolism, which were reported as loci that had experienced positive natural selection. Two genes related to alcohol metabolism were found to be 12,500 years out of phase with the time when they began to increase in the allele frequency; this finding indicates that the genomic diversity of Japanese people has been shaped by events closely related to agriculture and food production.

A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis.

BACKGROUND AND AIMS: Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-nonspecific) and nonshared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC. APPROACH AND RESULTS: Protein tyrosine phosphatase nonreceptor type 2 ( PTPN2) was identified as a novel PBC susceptibility gene locus through GWAS and subsequent genome-wide meta-analysis involving 2181 cases and 2699 controls from the Japanese population (GWAS-lead variant: rs8098858, p = 2.6 × 10 -8 ). In silico and in vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells and plasmacytoid dendritic cells. Infiltration of PTPN2-positive T-cells and plasmacytoid dendritic cells was confirmed in the portal area of the PBC liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758. CONCLUSIONS: PTPN2 , a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC through an insufficient negative feedback loop caused by the risk allele of rs2292758 in IFN-γ signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

WDR45 variants as a major cause for a clinically variable intellectual disability syndrome from early infancy in females.

Pathogenic variants of WD repeat domain 45 (WDR45) cause neurodegeneration with brain iron accumulation 5 (NBIA5), which is characterised by progressive neurological regression and brain iron accumulation in adulthood. Early diagnosis of NBIA5 patients is difficult because they often show only a non-specific developmental delay in childhood, but it is essential for lifelong medical management. We investigated 32 females with developmental delays for coding variants of WDR45 using Sanger sequencing. Whole-genome sequencing (WGS) and X chromosome inactivation (XCI) analysis were also performed. We identified two disease-causing variants, one of which was a novel stop-loss variant, c.1051delG p.(Val351CysfsTer60), in a female with severe developmental delay from early infancy with epileptic spasms. The XCI analysis (which we originally developed) suggested a random pattern in white blood cells. WGS did not reveal any other pathogenic variants, including those in two iron transporter genes. Together with our previous findings in the WGS study, WDR45 variants accounted for 12% (6/51) of the females with developmental delay, suggesting that WDR45 is a major gene in females with developmental delay. Pathogenic variants of WDR45 result in various phenotypes that do not necessarily correlate with variant types or XCI skewing patterns.

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.

Association study of GBA1 variants with MSA based on comprehensive sequence analysis -Pitfalls in short-read sequence analysis depending on the human reference genome.

Multiple system atrophy (MSA) is a neurodegenerative disorder characterized by various combinations of autonomic failure, parkinsonism, and cerebellar ataxia. To elucidate variants associated with MSA, we have been conducting short-read-based whole-genome sequence analysis. In the process of the association studies, we initially focused on GBA1, a previously proposed susceptibility gene for MSA, to evaluate whether GBA1 variants can be efficiently identified despite its extraordinarily high homology with its pseudogene, GBA1LP. To accomplish this, we conducted a short-read whole-genome sequence analysis with alignment to GRCh38 as well as Sanger sequence analysis and compared the results. We identified five variants with inconsistencies between the two pipelines, of which three variants (p.L483P, p.A495P-p.V499V, p.L483_M489delinsW) were the results of misalignment due to minor alleles in GBA1P1 registered in GRCh38. The miscalling events in these variants were resolved by alignment to GRCh37 as the reference genome, where the major alleles are registered. In addition, a structural variant was not properly identified either by short-read or by Sanger sequence analyses. Having accomplished correct variant calling, we identified three variants pathogenic for Gaucher disease (p.S310G, p.L483P, and p.L483_M489delinsW). Of these variants, the allele frequency of p.L483P (0.003) in the MSA cases was higher than that (0.0011) in controls. The meta-analysis incorporating a previous report demonstrated a significant association of p.L483P with MSA with an odds ratio of 2.85 (95% CI; 1.05 - 7.76, p = 0.0400).

A novel NODAL variant in a young embolic stroke patient with visceral heterotaxy.

BACKGROUND: Ischemic stroke in young adults can be caused by a variety of etiologies including the monogenic disorders. Visceral heterotaxy is a condition caused by abnormal left-right determinations during embryonic development. We aimed to determine the cause of a young ischemic stroke patient with visceral heterotaxy. CASE PRESENTATION: We performed neurological, radiological, and genetic evaluations in a 17-year-old male patient presenting ischemic stroke and visceral heterotaxy to determine the underlying cause of this rare disease combination. Brain magnetic resonance imaging (MRI) showed evidence of embolic stroke, abdominal computed tomography (CT) showed visceral heterotaxy, and echocardiogram showed cardiac anomaly with right-to-left-shunt (RLS). Whole genome sequencing (WGS) revealed a heterozygous missense variant (NM_018055.5: c.1016 T > C, p.(Met339Val)) in the NODAL gene, which is essential to the determination of the left-right body axis. CONCLUSIONS: Our study highlights the importance of evaluating genetic etiology in young ischemic stroke and the need for stroke risk management in visceral heterotaxy patients with RLS. To the best of our knowledge, we report the first genetically-confirmed case of visceral heterotaxy with young embolic stroke reported to date.

Comprehensive analysis of microsatellite polymorphisms in human populations.

Microsatellites (MS) are tandem repeats of short units, and have been used for population genetics, individual identification, and medical genetics. However, studies of MS on a whole-genome level are limited, and genotyping methods for MS have yet to be established. Here, we analyzed approximately 8.5 million MS regions using a previously developed MS caller for short reads (MIVcall method) for three large publicly available human genome sequencing data sets: the Korean Personal Genome Project, Simons Genome Diversity Project, and Human Genome Diversity Project. Our analysis identified 253,114 polymorphic MS. A comparison among different populations suggests that MS in the coding region evolved by random genetic drift and natural selection. In an analysis of genetic structures, MS clearly revealed population structures as SNPs and detected clusters that were not found by SNPs in African and Oceanian populations. Based on the MS polymorphisms, we selected MS marker candidates for individual identification. Finally, we applied our method to a deep sequenced ancient DNA sample. This study provides a comprehensive picture of MS polymorphisms and application to human population studies.

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

Demographic history of Ryukyu islanders at the southern part of the Japanese Archipelago inferred from whole-genome resequencing data.

The Ryukyu Islands are located in the southernmost part of the Japanese Archipelago and consist of several island groups. Each island group has its own history and culture, which differ from those of mainland Japan. People of the Ryukyu Islands are genetically subdivided; however, their detailed demographic history remains unclear. We report the results of a whole-genome sequencing analysis of a total of 50 Ryukyu islanders, focusing on genetic differentiation between Miyako and Okinawa islanders. We confirmed that Miyako and Okinawa islanders cluster differently in principal component analysis and ADMIXTURE analysis and that there is a population structure among Miyako islanders. The present study supports the hypothesis that population differentiation is primarily caused by genetic drift rather than by differences in the rate of migration from surrounding regions, such as the Japanese main islands or Taiwan. In addition, the genetic cline observed among Miyako and Okinawa islanders can be explained by recurrent migration beyond the bounds of these islands. Our analysis also suggested that the presence of multiple subpopulations during the Neolithic Ryukyu Jomon period is not crucial to explain the modern Ryukyu populations. However, the assumption of multiple subpopulations during the time of admixture with mainland Japanese is necessary to explain the modern Ryukyu populations. Our findings add insights that could help clarify the complex history of populations in the Ryukyu Islands.

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.

Response to COVID-19 vaccine is reduced in patients with inflammatory bowel disease, but improved with additional dose.

BACKGROUND: Coronavirus disease 2019 (COVID-19) vaccination is recommended for patients with inflammatory bowel disease (IBD); however, suppressed immune responses have been reported for fully vaccinated patients under immunosuppressive therapy, mainly from Western countries. We prospectively analyzed antibody titers of IBD patients in Asia induced by two-dose and additional dose of messengerRNA COVID-19 vaccine. METHODS: After measuring high-affinity antibody titers, factors associated with antibody titers were identified by multiple regression analyses using the following covariates: sex, age (≥60 or <60 years), disease type (Crohn's disease or ulcerative colitis), vaccine type (BNT162b2 or mRNA-1273), time from second/third vaccination, molecular-targeted agent (anti-tumor necrosis factor [TNF] agents, ustekinumab, vedolizumab, tofacitinib, or no molecular-targeted agents), thiopurine, steroid, and 5-aminosalicylic acid. RESULTS: Among 409 patients analyzed, mean titer was 1316.7 U/mL (SD, 1799.3); 403 (98.5%) were judged to be seropositive (≥0.8 U/mL), and 389 (95.1%) had neutralizing antibodies (≥15 U/mL). After the third vaccination, mean titer raised up to 21 123.8 U/mL (SD, 23 474.5); all 179 were seropositive, and 178 (99.4%) had neutralizing antibodies. In 248 patients with genetic data, there was no difference in mean titer after two/third doses between carriers and non-carriers of HLA-A24 associated with severe disease during COVID-19 infection. A multiple regression analyses using covariates revealed that older age, vaccine type (BNT162b2), time from second/third dose, anti-TNF agent, tofacitinib, and thiopurine were independently associated with lower antibody titers. CONCLUSIONS: Our findings further support the recommendation for COVID-19 vaccination in patients under immunosuppressive therapy, especially additional third dose for patients receiving anti-TNF agents and/or thiopurine or tofacitinib.

Paclitaxel-carboplatin plus bevacizumab therapy for advanced neuroendocrine carcinoma of the uterine cervix: A retrospective case series.

AIM: There is no conclusive data on the prognosis of patients who receive paclitaxel-carboplatin (TC) plus bevacizumab therapy for advanced neuroendocrine carcinoma (NEC) of the uterine cervix, a rare histological subtype of cervical cancer. Thus, the aim of this study was to determine the efficacy of TC chemotherapy plus bevacizumab and bevacizumab single maintenance therapy for advanced NEC of the cervix. METHODS: This was a retrospective review of patients who received TC plus bevacizumab therapy for metastatic, recurrent, or persistent NEC of the cervix at seven institutions between 2015 and 2020. Relevant data were extracted from the patients' medical records and analyzed. RESULTS: Seven patients, including six with small-cell NEC and one with large-cell NEC, were included for analysis. Three patients received bevacizumab single maintenance therapy following TC plus bevacizumab therapy, whereas four patients did not receive bevacizumab single maintenance therapy. The median overall survival and progression-free survival of the patients who received bevacizumab single maintenance therapy were longer than those of the patients who did not receive the therapy (34 months vs. 10.5 months and 19 months vs. 5 months, respectively). However, the patients who received bevacizumab single maintenance therapy had received cisplatin-based chemotherapy previously. CONCLUSIONS: On the premise that cisplatin-based chemotherapy is administered as the first-line treatment for advanced NEC of the cervix, bevacizumab single maintenance therapy following TC plus bevacizumab may be considered the second- or third-line treatment. However, the risk of adverse events, such as intestinal perforation, should be discussed with patients.

Bcl2l1 Deficiency in Osteoblasts Reduces the Trabecular Bone Due to Enhanced Osteoclastogenesis Likely through Osteoblast Apoptosis.

Bcl2l1 (Bcl-XL) belongs to the Bcl-2 family, Bcl2 and Bcl2-XL are major anti-apoptotic proteins, and the apoptosis of osteoblasts is a key event for bone homeostasis. As the functions of Bcl2l1 in osteoblasts and bone homeostasis remain unclear, we generated osteoblast-specific Bcl2l1-deficient (Bcl2l1fl/flCre) mice using 2.3-kb Col1a1 Cre. Trabecular bone volume and the trabecular number were lower in Bcl2l1fl/flCre mice of both sexes than in Bcl2l1fl/fl mice. In bone histomorphometric analysis, osteoclast parameters were increased in Bcl2l1fl/flCre mice, whereas osteoblast parameters and the bone formation rate were similar to those in Bcl2l1fl/fl mice. TUNEL-positive osteoblastic cells and serum TRAP5b levels were increased in Bcl2l1fl/flCre mice. The deletion of Bcl2l1 in osteoblasts induced Tnfsf11 expression, whereas the overexpression of Bcl-XL had no effect. In a co-culture of Bcl2l1-deficient primary osteoblasts and wild-type bone-marrow-derived monocyte/macrophage lineage cells, the numbers of multinucleated TRAP-positive cells and resorption pits increased. Furthermore, serum deprivation or the deletion of Bcl2l1 in primary osteoblasts increased apoptosis and ATP levels in the medium. Therefore, the reduction in trabecular bone in Bcl2l1fl/flCre mice may be due to enhanced bone resorption through osteoblast apoptosis and the release of ATP from apoptotic osteoblasts, and Bcl2l1 may inhibit bone resorption by preventing osteoblast apoptosis.

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

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

Point-of-Care Ultrasound for an Esophageal Foreign Body.

BACKGROUND: Press-through packs (PTPs) are widely used for storing drugs. The number of cases involving improper swallowing and esophageal foreign bodies (EFBs) has increased with the increasing use of PTPs. Studies have reported the utility of point-of-care ultrasound (POCUS) for evaluating EFBs. The application of POCUS for esophageal PTPs has not been reported. CASE REPORT: An 83-year-old woman complaining of neck pain and odynophagia that occurred after improperly swallowing 2 PTPs was admitted to the emergency department. EFBs were suspected, and POCUS revealed a hyperechoic material, suggestive of a PTP, in the cervical esophagus. Endoscopy was immediately performed, and the PTPs were successfully removed without complications. WHY SHOULD AN EMERGENCY PHYSICIAN BE AWARE OF THIS?: This study was significant because it documented the application of POCUS to detect esophageal PTPs. POCUS is a simple and noninvasive technique for identifying EFBs without radiation exposure. © 2022 Elsevier Inc.

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.

Genome-wide association study in patients with pulmonary Mycobacterium avium complex disease.

RATIONALE: Nontuberculous mycobacteria (NTM) are environmental mycobacteria that can cause a chronic progressive lung disease. Although epidemiological data indicate potential genetic predisposition, its nature remains unclear. OBJECTIVES: We aimed to identify host susceptibility loci for Mycobacterium avium complex (MAC), the most common NTM pathogen. METHODS: This genome-wide association study (GWAS) was conducted in Japanese patients with pulmonary MAC and healthy controls, followed by genotyping of candidate single-nucleotide polymorphisms (SNPs) in another Japanese cohort. For verification by Korean and European ancestry, we performed SNP genotyping. RESULTS: The GWAS discovery set included 475 pulmonary MAC cases and 417 controls. Both GWAS and replication analysis of 591 pulmonary MAC cases and 718 controls revealed the strongest association with chromosome 16p21, particularly with rs109592 (p=1.64×10-13, OR 0.54), which is in an intronic region of the calcineurin-like EF-hand protein 2 (CHP2). Expression quantitative trait loci analysis demonstrated an association with lung CHP2 expression. CHP2 was expressed in the lung tissue in pulmonary MAC disease. This SNP was associated with the nodular bronchiectasis subtype. Additionally, this SNP was significantly associated with the disease in patients of Korean (p=2.18×10-12, OR 0.54) and European (p=5.12×10-03, OR 0.63) ancestry. CONCLUSIONS: We identified rs109592 in the CHP2 locus as a susceptibility marker for pulmonary MAC disease.

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.