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

BACKGROUND AIMS

Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-non-specific) and non-shared (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 RESULTS

Protein tyrosine phosphatase non-receptor type 2 (PTPN2) was identified as a novel PBC susceptibility gene locus through a GWAS and subsequent genome-wide meta-analysis involving 2,181 cases and 2,699 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 (Tfh) and plasmacytoid dendritic cells (pDCs). Infiltration of PTPN2-positive T-cells and pDCs were 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 via an insufficient negative feedback loop caused by the PTPN2 risk allele of rs2292758 in IFN signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

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.

Dried blood spot-based host genome analysis technique targeting pathological associations with hepatitis B: Development and clinical application in the Cambodian population

AIM

Reports of patients with hepatitis B have highlighted associations between polymorphisms in the human leukocyte antigen (HLA)-DPB1, CXCL13, and CXCR5 genes and disease pathology. Owing to its potential to contribute to the development of new diagnostic and therapeutic methods, we aimed to establish a reliable host genome analysis technique that can be used in countries with inadequate infrastructure.

METHOD

We compared multiple commercially available kits for dried blood spot (DBS)-based sample collection to develop a basic DBS-based host genome analysis technique. We then collected blood samples from Cambodian patients with hepatitis  B and performed single-nucleotide polymorphism genotyping and HLA allele typing by the DBS system.

RESULT

We were able to perform single-nucleotide polymorphism genotyping and HLA allele typing with host DNA samples obtained using a combination of a HemaSpot™ filter paper-based device and a SMITEST® EX-R&D DNA extraction kit. The accuracy of genotyping using samples obtained by this method was not inferior to one using samples obtained by venipuncture. In the Cambodian population, significant associations of HLA-DPB1*04:01 with protection against chronic hepatitis B virus (HBV) infection, and HLA-DPB1*05:01 and HLA-DPB1*13:01 with susceptibility to chronic HBV infection were identified.

CONCLUSION

Based on the DBS system, we clarified the associations of HLA-DPB1 alleles with chronic HBV infection in the Cambodian population for the first time. Because the DBS is a low-cost, durable, transportable, and easy-to-handle modality, genetic analysis based on the DBS system is a feasible strategy for obtaining a deeper understanding of HBV epidemiology, especially in middle- or low-income countries.

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.

Genome-wide association study identifies a new susceptibility locus in PLA2G4C for Multiple System Atrophy

To elucidate the molecular basis of multiple system atrophy (MSA), a neurodegenerative disease, we conducted a genome-wide association study (GWAS) in a Japanese MSA case/control series followed by replication studies in Japanese, Korean, Chinese, European and North American samples. In the GWAS stage rs2303744 on chromosome 19 showed a suggestive association ( P = 6.5 × 10 ^-7 ) that was replicated in additional Japanese samples ( P = 2.9 × 10 ^-6 . OR = 1.58; 95% confidence interval, 1.30 to 1.91), and then confirmed as highly significant in a meta-analysis of East Asian population data ( P = 5.0 × 10 ^-15 . Odds ratio= 1.49; 95% CI 1.35 to 1.72). The association of rs2303744 with MSA remained significant in combined European/North American samples ( P =0.023. Odds ratio=1.14; 95% CI 1.02 to 1.28) despite allele frequencies being quite different between these populations. rs2303744 leads to an amino acid substitution in PLA2G4C that encodes the cPLA2γ lysophospholipase/transacylase. The cPLA2γ-Ile143 isoform encoded by the MSA risk allele has significantly decreased transacylase activity compared with the alternate cPLA2γ-Val143 isoform that may perturb membrane phospholipids and α-synuclein biology.

Regulation of MEK inhibitor selumetinib sensitivity by AKT phosphorylation in the novel BRAF L525R mutant

BACKGROUND

Oncogenic mutations in BRAF genes are found in approximately 5-10% of colorectal cancers. The majority of BRAF mutations are located within exons 11-15 of the catalytic kinase domains, with BRAF V600E accounting for more than 80% of the observed BRAF mutations. Sensitivity to BRAF- and mitogen-activated protein kinase (MEK) inhibitors varies depending on BRAF mutations and tumor cell types. Previously, we newly identified, BRAF L525R-mutation, in the activation segment of the kinase in colorectal cancer patient. Here, we characterized the function of the BRAF L525R mutation.

METHODS

HEK293 cells harboring a BRAF mutation (V600E or L525R) were first characterized and then treated with cetuximab, dabrafenib, and selumetinib. Cell viability was measured using WST-1 assay and the expression of proteins involved in the extracellular signal-regulated kinase (ERK) and protein kinase B (AKT) signaling pathways was evaluated using western blot analysis.

RESULTS

The MEK inhibitor selumetinib effectively inhibited cell proliferation and ERK phosphorylation in BRAF L525R cells but not in BRAF V600E cells. Further studies revealed that AKT phosphorylation was reduced by selumetinib in BRAF L525R cells but not in BRAF V600E cells or selumetinib-resistant BRAF L525R cells. Moreover, the AKT inhibitor overcame the selumetinib resistance.

CONCLUSIONS

We established a model system harboring BRAF L525R using HEK293 cells. BRAF L525R constitutively activated ERK. AKT phosphorylation caused sensitivity and resistance to selumetinib. Our results suggest that a comprehensive network analysis may provide insights to identify effective therapies.

Development of a novel microarray data analysis tool without normalization for genotyping degraded forensic DNA

Since the arrest of the Golden State Killer in the US in April 2018, forensic geneticists have been increasingly interested in the investigative genetic genealogy (IGG) method. While this method has already been in practical use as a powerful tool for criminal investigation, we have yet to know well the limitations and potential risks. In this current study, we performed an evaluation study focusing on degraded DNA using the Affymetrix Genome-Wide Human SNP Array 6.0 platform (Thermo Fisher Scientific). We revealed one of the potential problems that occur during SNP genotype determination using a microarray-based platform. Our analysis results indicated that the SNP profiles derived from degraded DNA contained many false heterozygous SNPs. In addition, it was confirmed that the total amount of probe signal intensity on microarray chips derived from degraded DNA decreased significantly. Because the conventional analysis algorithm performs normalization during genotype determination, we concluded that noise signals could be genotype-called. To address this issue, we proposed a novel microarray data analysis method without normalization (nMAP). Although the nMAP algorithm resulted in a low call rate, it substantially improved genotyping accuracy. Finally, we confirmed the usefulness of the nMAP algorithm for kinship inferences. These findings and the nMAP algorithm will make a contribution to the advance of the IGG method.

A common variant rs2054564 in ADAMST17 is associated with susceptibility to lumbar spondylosis

The molecular pathophysiology underlying lumbar spondylosis development remains unclear. To identify genetic factors associated with lumbar spondylosis, we conducted a genome-wide association study using 83 severe lumbar spondylosis cases and 182 healthy controls and identified 65 candidate disease-associated single nucleotide polymorphisms (SNPs). Replication analysis in 510 case and 911 control subjects from five independent Japanese cohorts identified rs2054564, located in intron 7 of ADAMTS17, as a disease-associated SNP with a genome-wide significance threshold (P = 1.17 × 10-11, odds ratio = 1.92). This association was significant even after adjustment of age, sex, and body mass index (P = 7.52 × 10-11). A replication study in a Korean cohort, including 123 case and 319 control subjects, also verified the significant association of this SNP with severe lumbar spondylosis. Immunohistochemistry revealed that fibrillin-1 (FBN1) and ADAMTS17 were co-expressed in the annulus fibrosus of intervertebral discs (IVDs). ADAMTS17 overexpression in MG63 cells promoted extracellular microfibrils biogenesis, suggesting the potential role of ADAMTS17 in IVD function through interaction with fibrillin fibers. Finally, we provided evidence of FBN1 involvement in IVD function by showing that lumbar IVDs in patients with Marfan syndrome, caused by heterozygous FBN1 gene mutation, were significantly more degenerated. We identified a common SNP variant, located in ADAMTS17, associated with susceptibility to lumbar spondylosis and demonstrated the potential role of the ADAMTS17-fibrillin network in IVDs in lumbar spondylosis development.

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.

Atezolizumab plus bevacizumab versus lenvatinib or sorafenib in non-viral unresectable hepatocellular carcinoma: an international propensity score matching analysis

BACKGROUND

A growing body of evidence suggests that non-viral hepatocellular carcinoma (HCC) might benefit less from immunotherapy.

MATERIALS AND METHODS

We carried out a retrospective analysis of prospectively collected data from consecutive patients with non-viral advanced HCC, treated with atezolizumab plus bevacizumab, lenvatinib, or sorafenib, in 36 centers in 4 countries (Italy, Japan, Republic of Korea, and UK). The primary endpoint was overall survival (OS) with atezolizumab plus bevacizumab versus lenvatinib. Secondary endpoints were progression-free survival (PFS) with atezolizumab plus bevacizumab versus lenvatinib, and OS and PFS with atezolizumab plus bevacizumab versus sorafenib. For the primary and secondary endpoints, we carried out the analysis on the whole population first, and then we divided the cohort into two groups: non-alcoholic fatty liver disease (NAFLD)/non-alcoholic steatohepatitis (NASH) population and non-NAFLD/NASH population.

RESULTS

One hundred and ninety patients received atezolizumab plus bevacizumab, 569 patients received lenvatinib, and 210 patients received sorafenib. In the whole population, multivariate analysis showed that treatment with lenvatinib was associated with a longer OS [hazard ratio (HR) 0.65; 95% confidence interval (CI) 0.44-0.95; P = 0.0268] and PFS (HR 0.67; 95% CI 0.51-0.86; P = 0.002) compared to atezolizumab plus bevacizumab. In the NAFLD/NASH population, multivariate analysis confirmed that lenvatinib treatment was associated with a longer OS (HR 0.46; 95% CI 0.26-0.84; P = 0.0110) and PFS (HR 0.55; 95% CI 0.38-0.82; P = 0.031) compared to atezolizumab plus bevacizumab. In the subgroup of non-NAFLD/NASH patients, no difference in OS or PFS was observed between patients treated with lenvatinib and those treated with atezolizumab plus bevacizumab. All these results were confirmed following propensity score matching analysis. By comparing patients receiving atezolizumab plus bevacizumab versus sorafenib, no statistically significant difference in survival was observed.

CONCLUSIONS

The present analysis conducted on a large number of advanced non-viral HCC patients showed for the first time that treatment with lenvatinib is associated with a significant survival benefit compared to atezolizumab plus bevacizumab, in particular in patients with NAFLD/NASH-related HCC.

The role of cellular senescence in aortic dissection

Background

Aortic dissection (AD) is a catastrophic disease that occurs suddenly. The acute mortality is high and those who survived frequently suffer from serious complications such as aneurysm formation and distal ischemia due to progressive destruction of the aortic walls. Currently, no predictor of AD onset is available nor therapeutic intervention to specifically prevent the progressive destruction in AD, because the molecular pathogenesis is largely unknown. Clinical and experimental studies highlighted the importance of inflammation in AD, although the regulatory mechanism of inflammation remains unclear. Recently, we found that cell proliferation precedes the inflammatory response in AD. Because cell proliferation causes cellular senescence that can induce inflammatory response, we hypothesized that cellular senescence participates in AD pathogenesis.

Objective

We investigated if cellular senescence contributes to AD development and progression in mouse AD model.

Methods and results

A mouse AD model was created by continuous infusion of beta-aminopropionitrile and angiotensin II (BAPN+AngII), where AD starts to develop in 3 days and occurs to most of the mice in 14 days accompanied by frequent AD rupture and death. Infusion of BAPN+AngII resulted in the induction of senescence markers Ink4a from day 3 before AD onset and persisted for the 14 days of the observational period.

Cellular senescence, as demonstrated by the expression of senescence-associated beta-galactosidase, was evident in intimal endothelial cells, medial smooth muscle cells, adventitial macrophages and fibroblasts. We examined the role of cellular senescence in AD pathogenesis by oral administration of ABT263 which is known as “senolytics” that eliminates senescent cells. ABT263 treatment reduced the expression of the senescence marker, prevented the death by AD rupture, and ameliorated the severity of AD lesion compared to the vehicle treatment. Transcriptome analysis revealed that ABT treatment suppressed the immune and inflammatory response in AD. Quantitative RT-PCR confirmed that ABT treatment prevented the induction of p21Cip1, interleukin-6, several chemokines and their receptors by 3-day infusion of BAPN+AngII.

Conclusions

These findings demonstrated that senescence of multiple cell types precedes AD development, which is likely to induce the inflammatory response. Elimination of senescent cells effectively prevented AD progression and death. Therefore, cellular senescence represents a potential predictor and a therapeutic target for AD.

Funding Acknowledgement

Type of funding sources: None.

Genome-Wide Association Study for Chronic Hepatitis B Infection in the Thai Population

To identify novel host genetic variants that predispose to hepatitis B virus (HBV) persistence, we performed the first genome-wide association study in the Thai population involving 318 cases of chronic hepatitis B and 309 healthy controls after quality control measures. We detected the genome-wide significant association of the HLA class II region (HLA-DPA1/DPB1, rs7770370, p-value = 7.71 × 10⁻¹⁰, OR = 0.49) with HBV chronicity. Subsequent HLA allele imputation revealed HLA-DPA1*01:03 (Pc = 1.21 × 10⁻⁶, OR = 0.53), HLA-DPB1*02:01 (Pc = 2.17 × 10⁻³, OR = 0.50), and HLA-DQB1*06:09 (Pc = 2.17 × 10⁻², OR = 0.07) as protective alleles, and HLA-DPA1*02:02 (Pc = 6.32 × 10⁻⁵, OR = 1.63), HLA-DPB1*05:01 (Pc = 1.13 × 10⁻⁴, OR = 1.72), HLA-DPB1*13:01 (Pc = 4.68 × 10⁻², OR = 1.60), and HLA-DQB1*03:03 (Pc = 1.11 × 10⁻³, OR = 1.84) as risk alleles for HBV persistence. We also detected suggestive associations in the PLSCR1 (rs35766154), PDLIM5 (rs62321986), SGPL1 (rs144998273), and MGST1 (rs1828682) loci. Among single-nucleotide polymorphisms in the PLSCR1 locus, rs1061307 was identified as the primary functional variant by in silico/in vitro functional analysis. In addition to replicating the association of the HLA class II region, we detected novel candidate loci that provide new insights into the pathophysiology of chronic hepatitis B.

Genetic association of IL17 and the importance of ABO blood group antigens in saliva to COVID-19

The outbreak of COVID-19 caused by infection with SARS-CoV-2 virus has become a worldwide pandemic, and the number of patients presenting with respiratory failure is rapidly increasing in Japan. An international meta-analysis has been conducted to identify genetic factors associated with the onset and severity of COVID-19, but these factors have yet to be fully clarified. Here, we carried out genomic analysis based on a genome-wide association study (GWAS) in Japanese COVID-19 patients to determine whether genetic factors reported to be associated with the onset or severity of COVID-19 in the international meta-GWAS are replicated in the Japanese population, and whether new genetic factors exist. Although no significant genome-wide association was detected in the Japanese GWAS, an integrated analysis with the international meta-GWAS identified for the first time the involvement of the IL17A/IL17F gene in the severity of COVID-19. Among nine genes reported in the international meta-GWAS as genes involved in the onset of COVID-19, the association of FOXP4-AS1, ABO, and IFNAR2 genes was replicated in the Japanese population. Moreover, combined analysis of ABO and FUT2 genotypes revealed that the presence of oral AB antigens was significantly associated with the onset of COVID-19. FOXP4-AS1 and IFNAR2 were also significantly associated in the integrated analysis of the Japanese GWAS and international meta-GWAS when compared with severe COVID-19 cases and the general population. This made it clear that these two genes were also involved in not only the onset but also the severity of COVID-19. In particular, FOXP4-AS1 was not found to be associated with the severity of COVID-19 in the international meta-GWAS, but an integrated analysis with the Japanese GWAS revealed an association with severity. Individuals with the SNP risk allele found between IL17A and IL17F had significantly lower mRNA expression levels of IL17F, suggesting that activation of the innate immune response by IL17F may play an important role in the severity of SARS-CoV-2 infection.

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.

The role of cellular senescence in aortic dissection

Background

Aortic dissection (AD) is a catastrophic disease that occurs suddenly. The acute mortality is high and those who survived frequently suffer from serious complications such as aneurysm formation and distal ischemia due to progressive destruction of the aortic walls. Currently, no reliable predictor is available for AD development and surgical intervention is the only therapeutic option to prevent the fatal events after AD development, because the pathogenesis of AD is largely unknown.

Clinical and experimental studies highlighted the importance of inflammation in AD pathogenesis, although the trigger of inflammation remains unclear. Recently, we found that cell proliferation precedes the inflammatory response in AD. Because cell proliferation triggers cellular senescence and senescent cells secrete of proinflammatory cytokines and matrix metalloproteinases, we hypothesized that cellular senescence may participate in AD pathogenesis.

Objective

We investigated if cellular senescence contributes to AD development and progression in a mouse model of AD.

Methods and results

A mouse AD model was created by continuous infusion of beta-aminopropionitrile and angiotensin II (BAPN+AngII), where AD starts to develop in 3 days and occurs to most of the mice in 14 days accompanied by frequent AD rupture and death. Infusion of BAPN+AngII resulted in the appearance of senescent cells that are positive for senescence-associated beta-galactosidase, and expression of senescence markers Arf and Ink4a in the aortic walls. Appearance of cellular senescence occurred in one day of BAPN+AngII infusion and continued throughout the observational period of 14 days. We examined the role of cellular senescence in AD pathogenesis by oral administration of ABT263 which is known as “senolytics” that eliminates senescent cells. ABT263 treatment reduced the expression of the senescence markers. In the vehicle-treated group, the mortality was 66.7% (12/18), whereas that of ABT263-treated group was 35% (14/20, P&lt;0.05 by log-rank test). The severity of AD, as assessed by the lesion length in vehicle group was33.2±3.1 mm, whereas that in ABT263 group was 24.6±1.8 mm (P&lt;0.05).

Conclusions

These findings demonstrated that cellular senescence precedes AD development, and ABT263 effectively prevented AD progression and death, indicating the involvement of cellular senescence in AD pathogenesis. Therefore, cellular senescence represents a potential predictor and a therapeutic target for AD.

Funding Acknowledgement

Type of funding sources: None.

Involvement of FAK in aortic dissection: potential role in aortic interstitial cells

Background

Aortic dissection (AD) is a fatal disease where the intimomedial layer of the aorta suddenly fail. Although it is widely accepted that hemodynamic stress on the aortic wall triggers its destruction that is further promoted by inflammatory response as exemplified by the infiltration of neutrophils and macrophages, molecular mechanism is unknown for the link of aortic wall stress, inflammation and tissue destruction. In general, mechanical stress to the tissue is converted to the cellular response through the cell adhesion molecules and the activation of focal adhesion kinase (Fak). Although it has been reported that Fak is involved in pathogenesis of aortic aneurysm by promoting migration and activation of macrophages, its role in AD is unknown. We hypothesized that Fak may be involved in AD pathogenesis.

Purpose

We investigated the involvement of Fak in AD pathogenesis, focusing on its role in inflammatory cells.

Methods and results

We created a mouse model of AD by continuous infusion of beta-aminopropionitrile, a collagen crosslink inhibitor, and angiotensin II (BAPN + Ang II). Immunostaining for activated Fak revealed that Fak was not activated in normal aorta, but was activated in the infiltrating inflammatory cells and in interstitial cells of the aortic wall after AD development. We examined the role of Fak by oral administration of PND-1186, a specific Fak inhibitor, in mouse AD model. Vehicle-treated group showed 63.6% mortality, whereas PND-1186-treated group showed 20% mortality (P&lt;0.01, n=20 for each group) in 14 days of the observational period. The aortic arch lesion, the most critical part in AD, was improved from 1.96±0.41 mm in vehicle group to 0.66±0.29 mm in PND group (P&lt;0.05). We next examined the cell type-specific role of Fak in AD by creating macrophage and granulocyte-specific deletion of Fak driven by LysM-Cre and floxed Fak system. Unexpectedly, the genetic deletion of Fak in macrophages and granulocytes had no impact on the mortality nor the severity of AD.

Conclusions

These findings proved that Fak plays a critical role in AD progression and death. Because Fak is dispensable for macrophages and granulocytes, other cell types, possibly aortic wall interstitial cells, may be regulated by Fak in AD pathogenesis. Deciphering the role of Fak would provide the fundamental understanding of AD pathogenesis.

Funding Acknowledgement

Type of funding sources: None.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians

Background

Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors).

Methods

In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors.

Results

The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors.

Conclusions

Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01900-y.

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.

Exploring correlations in genetic and cultural variation across language families in northeast Asia

Culture evolves in ways that are analogous to, but distinct from, genomes. Previous studies examined similarities between cultural variation and genetic variation (population history) at small scales within language families, but few studies have empirically investigated these parallels across language families using diverse cultural data. We report an analysis comparing culture and genomes from in and around northeast Asia spanning 11 language families. We extract and summarize the variation in language (grammar, phonology, lexicon), music (song structure, performance style), and genomes (genome-wide SNPs) and test for correlations. We find that grammatical structure correlates with population history (genetic history). Recent contact and shared descent fail to explain the signal, suggesting relationships that arose before the formation of current families. Our results suggest that grammar might be a cultural indicator of population history while also demonstrating differences among cultural and genetic relationships that highlight the complex nature of human history.

Genome-wide copy number variation analysis of hepatitis B infection in a Japanese population

Genome-wide association studies have been performed to identify common genetic variants associated with hepatitis B (HB). However, little is known about copy number variations (CNVs) in HB. In this study, we performed a genome-wide CNV analysis between 1830 healthy controls and 1031 patients with HB infection after quality control. Using signal calling by the Axiom Analysis Suite and CNV detection by PennCNV software, we obtained a total of 4494 CNVs across all individuals. The genes with CNVs that were found only in the HB patients were associated with the immune system, such as antigen processing. A gene-level CNV association test revealed statistically significant CNVs in the contactin 6 (CNTN6) gene. Moreover, we also performed gene-level CNV association tests in disease subgroups, including hepatocellular carcinoma patients, liver cirrhosis patients, and HBV carriers, including asymptomatic carriers and patients with HBV-derived chronic hepatitis. Our findings from germline cells suggested that patient-specific CNVs may be inherent genetic risk factors for HB.

The risk of contracting the hepatitis B virus may be linked to the number of copies of certain genes in an individual’s genome. A Japanese team led by Masataka Kikuchi, Osaka University, and Akihiro Nakaya, University of Tokyo, looked for repeated segments of the genome, known as copy number variants (CNVs), that differed between people with hepatitis B infections and those without. Studying around 3000 individuals of Japanese descent, the researchers identified several rare CNVs associated with immune function in hepatitis-affected individuals. They also found a common CNV in a gene called CNTN6 that the hepatitis B virus often uses to integrate itself into the genome of liver cells, a process that can lead to cancer. The findings underscore the importance of CNVs as inherited risk factors for hepatitis B and other viral infections.

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.

HLA-A*11:01:01:01, HLA-C*12:02:02:01-HLA-B*52:01:02:02, Age and Sex Are Associated With Severity of Japanese COVID-19 With Respiratory Failure

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing coronavirus disease 2019 (COVID-19) was announced as an outbreak by the World Health Organization (WHO) in January 2020 and as a pandemic in March 2020. The majority of infected individuals have experienced no or only mild symptoms, ranging from fully asymptomatic cases to mild pneumonic disease. However, a minority of infected individuals develop severe respiratory symptoms. The objective of this study was to identify susceptible HLA alleles and clinical markers that can be used in risk prediction model for the early identification of severe COVID-19 among hospitalized COVID-19 patients. A total of 137 patients with mild COVID-19 (mCOVID-19) and 53 patients with severe COVID-19 (sCOVID-19) were recruited from the Center Hospital of the National Center for Global Health and Medicine (NCGM), Tokyo, Japan for the period of February–August 2020. High-resolution sequencing-based typing for eight HLA genes was performed using next-generation sequencing. In the HLA association studies, HLA-A*11:01:01:01 [P_c = 0.013, OR = 2.26 (1.27–3.91)] and HLA-C*12:02:02:01-HLA-B*52:01:01:02 [P_c = 0.020, OR = 2.25 (1.24–3.92)] were found to be significantly associated with the severity of COVID-19. After multivariate analysis controlling for other confounding factors and comorbidities, HLA-A*11:01:01:01 [P = 3.34E-03, OR = 3.41 (1.50–7.73)], age at diagnosis [P = 1.29E-02, OR = 1.04 (1.01–1.07)] and sex at birth [P = 8.88E-03, OR = 2.92 (1.31–6.54)] remained significant. The area under the curve of the risk prediction model utilizing HLA-A*11:01:01:01, age at diagnosis, and sex at birth was 0.772, with sensitivity of 0.715 and specificity of 0.717. To the best of our knowledge, this is the first article that describes associations of HLA alleles with COVID-19 at the 4-field (highest) resolution level. Early identification of potential sCOVID-19 could help clinicians prioritize medical utility and significantly decrease mortality from COVID-19.

Importance of HBsAg recognition by HLA molecules as revealed by responsiveness to different hepatitis B vaccines

Hepatitis B (HB) vaccines (Heptavax-II and Bimmugen) designed based on HBV genotypes A and C are mainly used for vaccination against HB in Japan. To determine whether there are differences in the genetic background associated with vaccine responsiveness, genome-wide association studies were performed on 555 Heptavax-II and 1193 Bimmugen recipients. Further HLA imputation and detailed analysis of the association with HLA genes showed that two haplotypes, DRB1*13:02-DQB1*06:04 and DRB1*04:05-DQB1*04:01, were significantly associated in comparison with high-responders (HBsAb > 100 mIU/mL) for the two HB vaccines. In particular, HLA-DRB1*13:02-DQB1*06:04 haplotype is of great interest in the sense that it could only be detected by direct analysis of the high-responders in vaccination with Heptavax-II or Bimmugen. Compared with healthy controls, DRB1*13:02-DQB1*06:04 was significantly less frequent in high-responders when vaccinated with Heptavax-II, indicating that high antibody titers were less likely to be obtained with Heptavax-II. As Bimmugen and Heptavax-II tended to have high and low vaccine responses to DRB1*13:02, 15 residues were found in the Heptavax-II-derived antigenic peptide predicted to have the most unstable HLA-peptide binding. Further functional analysis of selected hepatitis B patients with HLA haplotypes identified in this study is expected to lead to an understanding of the mechanisms underlying liver disease.

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.

Transmission dynamics of SARS-CoV-2 on the Diamond Princess uncovered using viral genome sequence analysis

An outbreak of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) occurred aboard the Diamond Princess cruise ship between her January 20 departure and late February 2020. Here, we used phylodynamic analyses to investigate the transmission dynamics of SARS-CoV-2 during the outbreak. Using a Bayesian coalescent-based method, the estimated mean nucleotide substitution rate of 240 SARS-CoV-2 whole-genome sequences was approximately 7.13 × 10⁻⁴ substitutions per site per year. Population dynamics and the effective reproductive number (R_e) of SARS-CoV-2 infections were estimated using a Bayesian framework. The estimated origin of the outbreak was January 21, 2020. The infection spread substantially before quarantine on February 5. The R_e peaked at 6.06 on February 4 and gradually declined to 1.51, suggesting that transmission continued slowly even after quarantine. These findings highlight the high transmissibility of SARS-CoV-2 and the need for effective measures to control outbreaks in confined settings.

Serum CCL17 level becomes a predictive marker to distinguish between mild/moderate and severe/critical disease in patients with COVID-19

COVID-19, a novel coronavirus-related illness, has spread worldwide. Patients with apparently mild/moderate symptoms can suddenly develop severe pneumonia. Therefore, almost all COVID-19 patients require hospitalization, which can reduce limited medical resources in addition to overwhelming medical facilities. To identify predictive markers for the development of severe pneumonia, a comprehensive analysis of serum chemokines and cytokines was conducted using serial serum samples from COVID-19 patients. The expression profiles were analyzed along the time axis. Serum samples of common diseases were enrolled from a BioBank to confirm the usefulness of predictive markers. Five factors, IFN-λ3, IL-6, IP-10, CXCL9, and CCL17, were identified as predicting the onset of severe/critical symptoms. The factors were classified into two categories. Category A included IFN-λ3, IL-6, IP-10, and CXCL9, and their values surged and decreased rapidly before the onset of severe pneumonia. Category B included CCL17, which provided complete separation between the mild/moderate and the severe/critical groups at an early phase of SARS-CoV-2 infection. The five markers provided a high predictive value (area under the receiver operating characteristic curve (AUROC): 0.9-1.0, p < 0.001). Low expression of CCL17 was specifically observed in pre-severe COVID-19 patients compared with other common diseases, and the predictive ability of CCL17 was confirmed in validation samples of COVID-19. The factors identified could be promising prognostic markers to distinguish between mild/moderate and severe/critical patients, enabling triage at an early phase of infection, thus avoiding overwhelming medical facilities.

The role of B cells and IgG in aortic dissection

Background

Aortic dissection (AD) is one of the destructive and fatal aortic diseases, for which molecular pathogenesis is largely unknown. Recent studies have highlighted the importance of inflammatory response in AD. We and others reported that B cells and immunoglobulins participate in pathogenesis of abdominal aortic aneurysm, another form of aortic destructive disease, by promoting inflammatory response. It is not known whether and how B cells participate in AD pathogenesis.

Methods and results

Immunohistochemical staining of human AD tissue revealed that B cells were clustered together with T cells, macrophages and neutrophils at the entry site of AD with medial disruption. B cell cluster was also observed at the site of medial disruption in mouse model of AD that was induced by continuous infusion of beta-aminopropionitrile and angiotensin II (BAPN+AngII). In muMT mouse, which is deficient for B cells and immunoglobulins due to genetic deletion of immunoglobulin heavy chain, BAPN+AngII induced significantly less severe AD compared to that in wild type. Depositions of IgG and fibrinogen, one of the endogenous antigen for natural IgG, were observed after BAPN+AngII infusion before and after AD development in wild type mice. Deposition of fibrinogen was also observed in mMT mice after BAPN+AngII infusion. The rate of aortic rupture and sudden death was approximately 42% in wild type mice, while that in muMT mouse was 12% (P&lt;0.05). Administration of mouse normal polyclonal IgG to muMT mice resulted in dramatic increase in aortic rupture and sudden death, starting at day 7 of BAPN+AngII infusion, and reaching 69% of rupture rate, indicating the critical role of IgG in AD.

Conclusion

These findings demonstrated B cells and IgG are critically involved in the destructive inflammation of AD pathogenesis. Further, the deposition of fibrinogen, one of the targets of natural IgG, precedes the development of AD. Our findings may provide the conceptual foundation of the diagnostic strategy for on-going tissue destruction and for the therapeutic opportunities to intervene the progressive tissue destruction in AD.

Funding Acknowledgement

Type of funding source: Public grant(s) – National budget only. Main funding source(s): The Japan Society for the Promotion of Science

Socioeconomic and environmental status of riverside communities of Tapajós River, Brazil

Introduction

Public health care acknowledges socioeconomic factors as one of the multiple facets to promote or protect individual or collective health. Understanding and being more acquainted with the dynamics carried out in riverside communities is fundamental to investing in policies aimed at fighting diseases and illnesses that are particular to those communities.

Objective

Describing socioeconomic characteristics, healthcare infrastructure, and occupational activities of population from riverside communities along Tapajós river, Brazil.

Methodology

A descriptive study with 96 residents of communities of Tapajós river, Brazil. Participants have answered a semi-structure questionnaire and the analyses were descriptive and the variable categories were expressed as frequencies.

Results

There was a predominance of male participants (77.09%), incomplete elementary school (58.51%), involved in agricultural sector (85.10%), with a family income up to 1 minimum wage (68.81%). However, only 21.50% have declared making their living exclusively out of agriculture. Seedlings and seeds are mainly obtained through an exchange system running inside the community. In relation to infrastructure, the water comes from artesian wells (68.88%) and 51.63% claim not treating water whatsoever. All interviewees mentioned the lack of a sewage system and residues are taken to a rudimentary cesspool (76.59%). Waste produced by the community is burned out (93.61%). The main difficulties pointed out by the interviewees were: transportation (44.94%), health (32.14%), communication (21.42%) and government cooperation (15.47%).

Conclusions

The communities in this project show similar structural dynamic based on subsistence family agriculture and poor infrastructure of basic services. The data collected can be the basis for future public policies aiming at the promotion of food production autonomy, economic autonomy, and improvement of health indicators of these communities.

Key messages

Contextualizing the structural dynamic of vulnerable populations is key to plan actions aimed at tackling and dealing with social determinants involved in the health-illness process. Contextualizing the structural dynamic of vulnerable populations is fundamental to guide strategies aimed at intervening in the social and health determinants.

Orchestration of Intracellular Circuits by G Protein-Coupled Receptor 39 for Hepatitis B Virus Proliferation

Hepatitis B virus (HBV), a highly persistent pathogen causing hepatocellular carcinoma (HCC), takes full advantage of host machinery, presenting therapeutic targets. Here we aimed to identify novel druggable host cellular factors using the reporter HBV we have recently generated. In an RNAi screen of G protein-coupled receptors (GPCRs), GPCR39 (GPR39) appeared as the top hit to facilitate HBV proliferation. Lentiviral overexpression of active GPR39 proteins and an agonist enhanced HBV replication and transcriptional activities of viral promoters, inducing the expression of CCAAT/enhancer binding protein (CEBP)-β (CEBPB). Meanwhile, GPR39 was uncovered to activate the heat shock response, upregulating the expression of proviral heat shock proteins (HSPs). In addition, glioma-associated oncogene homologue signaling, a recently reported target of GPR39, was suggested to inhibit HBV replication and eventually suppress expression of CEBPB and HSPs. Thus, GPR39 provirally governed intracellular circuits simultaneously affecting the carcinopathogenetic gene functions. GPR39 and the regulated signaling networks would serve as antiviral targets, and strategies with selective inhibitors of GPR39 functions can develop host-targeted antiviral therapies preventing HCC.

Admixture and natural selection shaped genomes of an Austronesian-speaking population in the Solomon Islands

People in the Solomon Islands today are considered to have derived from Asian- and Papuan-related ancestors. Papuan-related ancestors colonized Near Oceania about 47,000 years ago, and Asian-related ancestors were Austronesian (AN)-speaking population, called Lapita, who migrated from Southeast Asia about 3,500 years ago. These two ancestral populations admixed in Near Oceania before the expansion of Lapita people into Remote Oceania. To understand the impact of the admixture on the adaptation of AN-speaking Melanesians in Near Oceania, we performed the genome-wide single nucleotide polymorphism (SNP) analysis of 21 individuals from Munda, the main town of the New Georgia Islands in the western Solomon Islands. Population samples from Munda were genetically similar to other Solomon Island population samples. The analysis of genetic contribution from the two different ancestries to the Munda genome revealed significantly higher proportions of Asian- and Papuan-related ancestries in the region containing the annexin A1 (ANXA1) gene (Asian component > 82.6%) and in the human leukocyte antigen (HLA) class II region (Papuan component > 85.4%), respectively. These regions were suspected to have undergone natural selection since the time of admixture. Our results suggest that admixture had affected adaptation of AN-speaking Melanesians in the Solomon Islands.

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.

P6494Activation of focal adhesion kinase is involved in pathogenesis of aortic dissection in mice

Background

Aortic dissection (AD) is a fatal disease where the media of the aorta suddenly fail. Currently, Molecular pathogenesis of AD is unknown. Recently, we discovered that the activity of MRTF-A, a mechanosensitive transcriptional regulator, promotes AD development. The activity of MRTF-A is regulated by mechanical stress to cells, which is transduced through focal adhesion and actin dynamics. However, it is currently unknown whether the mechanotransduction mechanism is involved in AD pathogenesis.

Purpose

We investigated the role of focal adhesion kinase (FAK), a signaling molecule that transduces mechanostress from focal adhesion to actin dynamics, in AD pathogenesis.

Methods

We created a mouse model of AD with a continuous infusion of beta-aminopropionitrile (150 mg/kg/day), a collagen crosslink inhibitor, and angiotensin II (1,000 ng/kg/min) (BAPN + AngII) by an osmotic pump. This model caused about 60% death in all mice due to AD rupture within 2 weeks. In this model, we examined the severity and mortality rate of aortic dissection after 2 weeks in mice administered with PND-1186, an orally available FAK inhibitor, and in those treated with vehicle (n=20 for each group). We performed immunochemical staining, immunofluorescence staining and Western blot for activated (phosphorylated) FAK (pFAK) to evaluate the activation status of FAK in the aortic tissue. We also performed transcriptome analysis of the aortic tissue in with and without PND-1186 with BAPN + AngII stimulation before AD development.

Results

Immunochemical staining revealed that FAK was inactive in normal mouse aorta, but was strongly activated in the aortic walls after AD development. Immunofluorescence staining showed that FAK was activated mainly in smooth muscle cells after AD development. Western blot analysis also revealed that FAK was activated in 3 days after BAPN + AngII infusion before AD development, followed by transient reduction at day 7, and re-activation after AD at day 14. Significantly, administration of PND-1186 resulted in a significant reduction in the severity of AD in the aortic arch (1.96±0.41 mm in vehicle group, 0.66±0.29 mm in PND group, P<0.05). In addition, survival rate improved from 36.4% to 80.0% by administration of PND-1186 (P<0.01). In immunofluorescence staining, the PND-1186 treated group showed weaker staining of pFAK. Transcriptome analysis showed that genes for hematopoiesis and immune system were suppressed in PND-1186 treated group.

Conclusions

These findings proved that FAK plays a central role in the pathogenesis of AD probably by transmitting pathological stress to the aortic wall to cause tissue destruction. We propose that FAK is a potential therapeutic target for limiting the fatal destruction of the aortic wall of AD.

P6493Syk activation is a defense mechanism in murine model of aortic dissection

Background

Aortic dissection (AD) is a serious clinical condition with unknown etiology that frequently results in fatal outcome. Recent studies showed essential role of inflammatory response both in promoting AD and aortic aneurysm (AA). However, the difference of the molecular pathogenesis between AD and AA is unclear. Previously, we reported that Syk, a tyrosine kinase that regulates differentiation and activation of inflammatory cells, promotes AA formation in a mouse model.

Objective

In the current study, we investigated the role of Syk in AD.

Methods and results

A mouse AD model was created by continuous infusion of beta-aminopropionitrile (125 mg/kg/day) and angiotensin II (1,000 ng/min/kg) (BAPN+AngII), which caused AD in approximately 80% of mice within 14 days. Immunohistochemical staining for activated (phosphorylated) Syk (pSyk) revealed that Syk was inactive in normal mouse aorta, but was activated in the aortic walls after AD development. Double immunofluorescence staining for pSyk and smooth muscle alpha actin showed that Syk was active not only in the infiltrating inflammatory cells, but also in smooth muscle cells in AD tissue. Western blot analysis revealed that BAPN+AngII treatment caused Syk activation at 3 days before AD development, followed by transient suppression, and reactivation at 14 days after AD development. We examined the significance of Syk activation in AD by treating mice with fostamatinib, a specific Syk inhibitor, before and during BAPN+AngII infusion. Notably, fostamatinib-treated group developed more severe AD compared to the vehicle-treated group. The AD lesion length was 3.80±0.86 mm for vehicle group and 8.87±1.69 mm for fostamatinib group (P<0.05, n=12 for each group). In addition, fostamatinib significantly worsened the mortality of mice due to the rupture of the aorta from 0% to 42% (P<0.05, n=12 for each group). Transcriptome analysis revealed that fostamatinib suppressed both positive and negative regulators of immune response, defense response and inflammatory response.

Conclusions

These findings uncovered the previously unrecognized role of Syk for protecting the aortic tissue in AD pathogenesis, and suggested fundamentally different disease mechanisms of AD and AA.

Histopathological factors affecting the extraction of high quality genomic DNA from tissue sections for next-generation sequencing

To enable the widespread application of genomic medicine, the extraction of genomic DNA from thin sections of archived formalin-fixed and paraffin-embedded (FFPE) tissue blocks for next-generation sequencing (NGS) is often necessary. However, there are currently no guidelines available on which specific regions of the microtome sections to use for macrodissection with respect to the histopathological factors observed under microscopic examination. The aim of this study was to clarify the relationship between histopathological factors and DNA quality, and to standardize the macrodissection method for more efficient implementation of NGS. FFPE tissue specimens of 218 patients from the Biomarker Research for Anti-EGFR Monoclonal Antibodies by Comprehensive Cancer Genomics study were used to investigate the relationship between 15 histopathological factors and the quantitative ratio of double-stranded DNA (dsDNA) to total nucleic acids, as well as the ∆ crossing point value of each tissue specimen. Multivariate logistic regression analysis revealed that specimen storage of ≥3 years was negatively associated with dsDNA quality (P=0.0007, OR: 4.30, 95% CI: 1.85-10.04). In contrast, the presence of a mucus pool was positively associated with dsDNA quality (P=0.0308, OR: 0.23, 95% CI: 0.06-0.87). Metastatic tumors and longer specimen storage periods were significantly associated with lower ∆Cp values (P=0.0007, OR: 4.43, 95% CI: 1.87-10.49; and P=0.0003, OR: 5.51, 95% CI: 2.18-13.95, respectively). Therefore, macrodissection should not be performed on specimens exhibiting histopathological factors associated with poor DNA quality. In particular, the use of tissue blocks with a storage period of <3 years allows the extraction of genomic DNA suitable for NGS.

Association Between High Levels of Hepatitis B Core Antibody and Seroclearance of Hepatitis B e Antigen in Individuals With Chronic Hepatitis B Virus Infection

For chronic hepatitis B patients, hepatitis B e antigen (HBeAg) seroclearance signals a transition from an immunologically active phase to an inactive carrier state with a reduction in hepatitis B virus (HBV) DNA levels and a reduced risk of hepatocellular carcinoma (HCC).¹ Predictors of HBeAg seroclearance include lower HBV DNA levels, viral genotype, the precore mutation, and higher serum alanine aminotransferase (ALT) levels.².

Genome-Wide Association Study Confirming a Strong Effect of HLA and Identifying Variants in CSAD/lnc-ITGB7-1 on Chromosome 12q13.13 Associated With Susceptibility to Fulminant Type 1 Diabetes

The first genome-wide association study of fulminant type 1 diabetes was performed in Japanese individuals. As previously reported using a candidate gene approach, a strong association was observed with multiple single nucleotide polymorphisms (SNPs) in the HLA region, and the strongest association was observed with rs9268853 in the class II DR region (P = 1.56 × 10^-23, odds ratio [OR] 3.18). In addition, rs11170445 in CSAD/lnc-ITGB7-1 on chromosome 12q13.13 showed an association at a genome-wide significance level (P = 7.58 × 10^-9, OR 1.96). Fine mapping of the region revealed that rs3782151 in CSAD/lnc-ITGB7-1 showed the lowest P value (P = 4.60 × 10^-9, OR 1.97 [95% CI 1.57-2.48]). The risk allele of rs3782151 is a cis expression quantitative trait locus for ITGB7 that significantly increases the expression of this gene. CSAD/lnc-ITGB7-1 was found to be strongly associated with susceptibility to fulminant, but not classical, autoimmune type 1 diabetes, implicating this locus in the distinct phenotype of fulminant type 1 diabetes.

POGLUT1, the putative effector gene driven by rs2293370 in primary biliary cholangitis susceptibility locus chromosome 3q13.33

Primary biliary cholangitis (PBC) is a chronic and cholestatic autoimmune liver disease caused by the destruction of intrahepatic small bile ducts. Our previous genome-wide association study (GWAS) identified six susceptibility loci for PBC. Here, in order to further elucidate the genetic architecture of PBC, a GWAS was performed on an additional independent sample set, then a genome-wide meta-analysis with our previous GWAS was performed based on a whole-genome single nucleotide polymorphism (SNP) imputation analysis of a total of 4,045 Japanese individuals (2,060 cases and 1,985 healthy controls). A susceptibility locus on chromosome 3q13.33 (including ARHGAP31, TMEM39A, POGLUT1, TIMMDC1, and CD80) was previously identified both in the European and Chinese populations and was replicated in the Japanese population (OR = 0.7241, P = 3.5 × 10^-9). Subsequent in silico and in vitro functional analyses identified rs2293370, previously reported as the top-hit SNP in this locus in the European population, as the primary functional SNP. Moreover, e-QTL analysis indicated that the effector gene of rs2293370 was Protein O-Glucosyltransferase 1 (POGLUT1) (P = 3.4 × 10^-8). This is the first study to demonstrate that POGLUT1 and not CD80 is the effector gene regulated by the primary functional SNP rs2293370, and that increased expression of POGLUT1 might be involved in the pathogenesis of PBC.

NFKB1 and MANBA Confer Disease Susceptibility to Primary Biliary Cholangitis via Independent Putative Primary Functional Variants

BACKGROUND & AIMS

Primary biliary cholangitis (PBC) is a chronic and cholestatic liver disease that eventually leads to cirrhosis and hepatic failure. We recently identified several susceptibility genes included NFKB1 and MANBA for PBC in the Japanese population by genome-wide association study. However, the primary functional variants in the NFKB1/MANBA region and the molecular mechanism for conferring disease susceptibility to PBC have not yet been clarified.

METHODS

We performed high-density association mapping based on a single-nucleotide polymorphism (SNP) imputation analysis, using data from a whole-genome sequence reference panel of 1070 Japanese individuals and the previous genome-wide association study (1389 PBC patients, 1508 healthy controls). Among SNPs (P < 5.0 × 10^-7) in the NFKB1/MANBA region, putative primary functional variants and the molecular mechanism for conferring disease susceptibility to PBC were identified by in silico/in vitro functional analysis.

RESULTS

Among the SNPs in the NFKB1/MANBA region, rs17032850 and rs227361, which changed the binding of transcription factors lymphoid enhancer-binding factor 1 (LEF-1) and retinoid X receptor α (RXRα), respectively, were identified as putative primary functional variants that regulate gene expression. In addition, expression-quantitative trait locus data and gene editing using a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system supported the potential role of rs17032850 and rs227361 in regulating NFKB1 and MANBA expression, respectively.

CONCLUSIONS

We identified independent putative primary functional variants in NFKB1/MANBA and showed the distinct molecular mechanism by which each putative primary functional variant conferred susceptibility to PBC. Our approach was useful to dissect the pathogenesis not only of PBC, but also other digestive diseases in which NFKB1/MANBA has been reported as a susceptibility locus.

Work-related psychosocial factors and metabolic syndrome onset among workers: a systematic review and meta-analysis

BACKGROUND

Work-related psychosocial factors have been associated with metabolic syndrome. However, no systematic reviews or meta-analyses have evaluated this association.

METHODS

A systematic literature search was conducted, using PubMed, Embase, PsycINFO, PsycARTICLES and the Japan Medical Abstracts Society. Eligible studies included those that examined the previously mentioned association; had a longitudinal or prospective cohort design; were conducted among workers; provided sufficient data for calculating odds ratios, relative risks or hazard ratios with 95% confidence intervals; were original articles in English or Japanese; and were published no later than 2016. Study characteristics, exposure and outcome variables and association measures of studies were extracted by the investigators independently.

RESULTS

Among 4,664 identified studies, 8 were eligible for review and meta-analysis. The pooled risk of adverse work-related stress on metabolic syndrome onset was significant and positive (RR = 1.47; 95% CI, 1.22-1.78). Sensitivity analyses limiting only the effects of job strain and shift work also indicated a significant positive relationship (RR = 1.75; 95% CI, 1.09-2.79; and RR = 1.59; 95% CI, 1.00-2.54, P = 0.049 respectively).

CONCLUSION

This study reveals a strong positive association between work-related psychosocial factors and an elevated risk of metabolic syndrome onset. The effects of job strain and shift work on metabolic syndrome appear to be significant.

A Glycomarker for Short-term Prediction of Hepatocellular Carcinoma: A Longitudinal Study With Serial Measurements

Objectives

Wisteria floribunda agglutinin-positive human Mac-2-binding protein (WFA⁺-M2BP) is a glycomarker. The present community-based long-term follow-up study repeatedly determined the serum WFA⁺-M2BP level and examined its short- and long-term associations with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC).

Methods

A total of 921 participants with antibodies against HCV seropositive, but seronegative for hepatitis B surface antigen were enrolled from seven townships in Taiwan during 1991–1992. The participants were regularly followed and their serum WFA⁺-M2BP levels were measured at baseline and follow-up. HCC was ascertained through active follow-up and computerized data linkage with the National Cancer Registration System until December 31, 2013. Cox proportional hazards and logistic regression models were applied to estimate the magnitude of associations between serum WFA⁺-M2BP levels and HCC.

Results

During a median follow-up of 21.7 years, 122 new-onset HCC cases were identified. Elevated serum WFA⁺-M2BP levels were associated with an increased risk of HCC (p < 0.001). Patients with increasing changes in serum WFA⁺-M2BP levels, relative to their baseline levels, had a 4.36-fold risk of HCC. The areas under receiver operating curves (AUROCs) of WFA⁺-M2BP for predicting HCC showed that the prediction efficacy was significantly higher while closer to HCC diagnosis (p = 0.024). The AUROC was 0.91 for predicting HCC within 1 year by including the predictors of age, sex, alanine aminotransferase, alpha-fetoprotein (AFP) and WFA⁺-M2BP.

Conclusions

Serum WFA⁺-M2BP level may elevate before HCC onset and is a short-term predictor of HCC among patients infected with HCV.

Key HLA-DRB1-DQB1 haplotypes and role of the BTNL2 gene for response to a hepatitis B vaccine

Approximately 5-10% of individuals who are vaccinated with a hepatitis B (HB) vaccine designed based on the hepatitis B virus (HBV) genotype C fail to acquire protective levels of antibodies. Here, host genetic factors behind low immune response to this HB vaccine were investigated by a genome-wide association study (GWAS) and Human Leukocyte Antigen (HLA) association tests. The GWAS and HLA association tests were carried out using a total of 1,193 Japanese individuals including 107 low responders, 351 intermediate responders, and 735 high responders. Classical HLA class II alleles were statistically imputed using the genome-wide SNP typing data. The GWAS identified independent associations of HLA-DRB1-DQB1, HLA-DPB1 and BTNL2 genes with immune response to a HB vaccine designed based on the HBV genotype C. Five HLA-DRB1-DQB1 haplotypes and two DPB1 alleles showed significant associations with response to the HB vaccine in a comparison of three groups of 1,193 HB vaccinated individuals. When frequencies of DRB1-DQB1 haplotypes and DPB1 alleles were compared between low immune responders and HBV patients, significant associations were identified for three DRB1-DQB1 haplotypes, and no association was identified for any of the DPB1 alleles. In contrast, no association was identified for DRB1-DQB1 haplotypes and DPB1 alleles in a comparison between high immune responders and healthy individuals. Conclusion: The findings in this study clearly show the importance of HLA-DR-DQ (i.e., recognition of a vaccine related HB surface antigen (HBsAg) by specific DR-DQ haplotypes) and BTNL2 molecules (i.e., high immune response to HB vaccine) for response to a HB vaccine designed based on the HBV genotype C. (Hepatology 2018).

Increased expression and altered localization of cathepsin Z are associated with progression to jaundice stage in primary biliary cholangitis

Our recent genome-wide association study found that the NELFCD/CTSZ locus was significantly associated with progression of primary biliary cholangitis (PBC) to jaundice stage in the Japanese population. In this study, we investigated the role of cathepsin Z in the etiology and pathology of PBC. Serum cathepsin Z levels were measured using enzyme-linked immunosorbent assay. The expression and localization of cathepsin Z in liver specimens were analyzed by western blotting and immunohistochemistry. In PBC patients, serum cathepsin Z levels were significantly increased with disease progression. In addition, its levels were positively correlated with alanine transaminase, aspartate transaminase and total bilirubin, and were negatively correlated with platelet count and albumin. Cathepsin Z expression was markedly increased in hepatocytes at later stages of PBC, and its localization was altered from the peri-bile canaliculus to the cytoplasm, where a fraction was no longer colocalized with endosomal/lysosomal vesicles. Similar altered expression of cathepsin Z was observed in end-stage of other cholestatic liver diseases including sepsis, obstructive jaundice, and Alagille syndrome. Our results indicate that altered expression and localization of cathepsin Z in hepatocytes are characteristic features of PBC and other cholestatic liver diseases, and are implicated in the progression of PBC.

Discerning the Origins of the Negritos, First Sundaland People: Deep Divergence and Archaic Admixture

Human presence in Southeast Asia dates back to at least 40,000 years ago, when the current islands formed a continental shelf called Sundaland. In the Philippine Islands, Peninsular Malaysia, and Andaman Islands, there exist indigenous groups collectively called Negritos whose ancestry can be traced to the "First Sundaland People." To understand the relationship between these Negrito groups and their demographic histories, we generated genome-wide single nucleotide polymorphism data in the Philippine Negritos and compared them with existing data from other populations. Phylogenetic tree analyses show that Negritos are basal to other East and Southeast Asians, and that they diverged from West Eurasians at least 38,000 years ago. We also found relatively high traces of Denisovan admixture in the Philippine Negritos, but not in the Malaysian and Andamanese groups, suggesting independent introgression and/or parallel losses involving Denisovan introgressed regions. Shared genetic loci between all three Negrito groups could be related to skin pigmentation, height, facial morphology and malarial resistance. These results show the unique status of Negrito groups as descended from the First Sundaland People.

Induction of IFN-λ3 as an additional effect of nucleotide, not nucleoside, analogues: a new potential target for HBV infection

OBJECTIVE

The clinical significance of polymorphisms in the interleukin-28B gene encoding interferon (IFN)-λ3, which has antiviral effects, is known in chronic HCV but not in HBV infection. Thus, we measured IFN-λ3 levels in patients with HBV and investigated its clinical significance and association with nucleos(t)ide (NUC) analogue administration.

DESIGN

Serum IFN-λ3 level was measured in 254 patients with HBV with varying clinical conditions using our own high sensitivity method. The resulting values were compared with various clinical variables. In addition, cell lines originating from various organs were cultured with NUCs, and the production of IFN-λ3 was evaluated.

RESULTS

Higher serum IFN-λ3 levels were detected in the patients treated with nucleotide analogues (adefovir or tenofovir) compared with those treated with nucleoside analogues (lamivudine or entecavir). There were no other differences in the clinical background between the two groups. A rise in the serum IFN-λ3 levels was observed during additional administration of the nucleotide analogues. In vitro experiments showed that the nucleotide analogues directly and dose-dependently induced IFN-λ3 production only in colon cancer cells. Furthermore, the supernatant from cultured adefovir-treated colon cancer cells significantly induced IFN-stimulated genes (ISGs) and inhibited hepatitis B surface antigen (HBsAg) production in hepatoma cells, as compared with the supernatant from entecavir-treated cells.

CONCLUSIONS

We discovered that the nucleotide analogues show an additional pharmacological effect by inducing IFN-λ3 production, which further induces ISGs and results in a reduction of HBsAg production. These findings provide novel insights for HBV treatment and suggest IFN-λ3 induction as a possible target.