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.

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.

Urine oxygenation predicts COVID-19 risk

Since February, 2023, the omicron variant has accounted for essentially all new coronavirus infections in Japan. If future infections involve mutant strains with the same level of infectivity and virulence as omicron, the government's basic policy will be to prevent the spread of infection, without compromising socioeconomic activities. Objectives include protecting pregnant women and elderly persons, and focusing on citizens requiring hospitalization and those at risk of serious illness, without imposing new social restrictions. Although the government tries to raise public awareness through education, most people affected by COVID-19 stay at home, and by the time patients become aware of the seriousness of their disease, it has often reached moderate or higher severity. In this review, we discuss why this situation persists even though the disease seems to have become milder with the shift from the delta variant to omicron. We also propose a pathophysiological method to determine the risk of severe illness. This assessment can be made at home in the early stages of COVID-19 infection, using urine analysis. Applicability of this method to drug discovery and development is also discussed.

Assessing the efficacy of target adaptive sampling long-read sequencing through hereditary cancer patient genomes

Innovations in sequencing technology have led to the discovery of novel mutations that cause inherited diseases. However, many patients with suspected genetic diseases remain undiagnosed. Long-read sequencing technologies are expected to significantly improve the diagnostic rate by overcoming the limitations of short-read sequencing. In addition, Oxford Nanopore Technologies (ONT) offers adaptive sampling and computationally driven target enrichment technology. This enables more affordable intensive analysis of target gene regions compared to standard non-selective long-read sequencing. In this study, we developed an efficient computational workflow for target adaptive sampling long-read sequencing (TAS-LRS) and evaluated it through application to 33 genomes collected from suspected hereditary cancer patients. Our workflow can identify single nucleotide variants with nearly the same accuracy as the short-read platform and elucidate complex forms of structural variations. We also newly identified several SINE-R/VNTR/Alu (SVA) elements affecting the APC gene in two patients with familial adenomatous polyposis, as well as their sites of origin. In addition, we demonstrated that off-target reads from adaptive sampling, which is typically discarded, can be effectively used to accurately genotype common single-nucleotide polymorphisms (SNPs) across the entire genome, enabling the calculation of a polygenic risk score. Furthermore, we identified allele-specific MLH1 promoter hypermethylation in a Lynch syndrome patient. In summary, our workflow with TAS-LRS can simultaneously capture monogenic risk variants including complex structural variations, polygenic background as well as epigenetic alterations, and will be an efficient platform for genetic disease research and diagnosis.

Cost-Effective Cas9-Mediated Targeted Sequencing of Spinocerebellar Ataxia Repeat Expansions

Hereditary repeat diseases are caused by an abnormal expansion of short tandem repeats in the genome. Among them, spinocerebellar ataxia (SCA) is a heterogeneous disease, and currently, 16 responsible repeats are known. Genetic diagnosis is obtained by analyzing the number of repeats through separate testing of each repeat. Although simultaneous detection of candidate repeats using current massively parallel sequencing technologies has been developed to avoid complicated multiple experiments, these methods are generally expensive. This study developed a cost-effective SCA repeat panel [Flongle SCA repeat panel sequencing (FLO-SCAp)] using Cas9-mediated targeted long-read sequencing and the smallest long-read sequencing apparatus, Flongle. This panel enabled the detection of repeat copy number changes, internal repeat sequences, and DNA methylation in seven patients with different repeat expansion diseases. The median (interquartile range) values of coverage and on-target rate were 39.5 (12 to 72) and 11.6% (7.5% to 16.5%), respectively. This approach was validated by comparing repeat copy number changes measured by FLO-SCAp and short-read whole-genome sequencing. A high correlation was observed between FLO-SCAp and short-read whole-genome sequencing when the repeat length was ≤250 bp (r = 0.98; P < 0.001). Thus, FLO-SCAp represents the most cost-effective method for conducting multiplex testing of repeats and can serve as the first-line diagnostic tool for SCA.

Identification of region-specific gene isoforms in the human brain using long-read transcriptome sequencing

In neurological and neuropsychiatric diseases, different brain regions are affected, and differences in gene expression patterns could potentially explain this mechanism. However, limited studies have precisely explored gene expression in different regions of the human brain. In this study, we performed long-read RNA sequencing on three different brain regions of the same individuals: the cerebellum, hypothalamus, and temporal cortex. Despite stringent filtering criteria excluding isoforms predicted to be artifacts, over half of the isoforms expressed in multiple samples across multiple regions were found to be unregistered in the GENCODE reference. We then especially focused on genes with different major isoforms in each brain region, even with similar overall expression levels, and identified that many of such genes including GAS7 might have distinct roles in dendritic spine and neuronal formation in each region. We also found that DNA methylation might, in part, drive different isoform expressions in different regions. These findings highlight the significance of analyzing isoforms expressed in disease-relevant sites.

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.

Characteristics of patients with COVID-19 who have deteriorating chest X-ray findings within 48 h: a retrospective cohort study

The severity of chest X-ray (CXR) findings is a prognostic factor in patients with coronavirus disease 2019 (COVID-19). We investigated the clinical and genetic characteristics and prognosis of patients with worsening CXR findings during early hospitalization. We retrospectively included 1656 consecutive Japanese patients with COVID-19 recruited through the Japan COVID-19 Task Force. Rapid deterioration of CXR findings was defined as increased pulmonary infiltrates in ≥ 50% of the lung fields within 48 h of admission. Rapid deterioration of CXR findings was an independent risk factor for death, most severe illness, tracheal intubation, and intensive care unit admission. The presence of consolidation on CXR, comorbid cardiovascular and chronic obstructive pulmonary diseases, high body temperature, and increased serum aspartate aminotransferase, potassium, and C-reactive protein levels were independent risk factors for rapid deterioration of CXR findings. Risk variant at the ABO locus (rs529565-C) was associated with rapid deterioration of CXR findings in all patients. This study revealed the clinical features, genetic features, and risk factors associated with rapid deterioration of CXR findings, a poor prognostic factor in patients with COVID-19.

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.

Genome-wide association analysis reveals the associations of NPHP4, TYW1-AUTS2 and SEMA6D for Behçet's disease and HLA-B*46:01 for its intestinal involvement

BACKGROUND

Intestinal involvement in Behçet's disease (BD) is associated with poor prognosis and is more prevalent in East Asian than in Mediterranean populations. Identifying the genetic causes of intestinal BD is important for understanding the pathogenesis and for appropriate treatment of BD patients.

METHODS

We performed genome-wide association studies (GWAS) and imputation/replication genotyping of human leukocyte antigen (HLA) alleles for 1,689 Korean and Turkish patients with BD (including 379 patients with intestinal BD) and 2,327 healthy controls, followed by replication using 593 Japanese patients with BD (101 patients with intestinal BD) and 737 healthy controls. Stratified cross-phenotype analyses were performed for 1) overall BD, 2) intestinal BD, and 3) intestinal BD without association of overall BD.

RESULTS

We identified three novel genome-wide significant susceptibility loci including NPHP4 (rs74566205; P=1.36 × 10-8), TYW1-AUTS2 (rs60021986; P=1.14 × 10-9), and SEMA6D (rs4143322; P=5.54 × 10-9) for overall BD, and a new association with HLA-B*46:01 for intestinal BD (P=1.67 × 10-8) but not for BD without intestinal involvement. HLA peptide binding analysis revealed that Mycobacterial peptides, have a stronger binding affinity to HLA-B*46:01 compared to the known risk allele HLA-B*51:01.

CONCLUSIONS

HLA-B*46:01 is associated with the development of intestinal BD; NPHP4, TYW1-AUTS2, and SEMA6D are susceptibility loci for overall BD.

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

Association between ABO blood group/genotype and COVID-19 in a Japanese population

An association between coronavirus disease 2019 (COVID-19) and the ABO blood group has been reported. However, such an association has not been studied in the Japanese population on a large scale. Little is known about the association between COVID-19 and ABO genotype. This study investigated the association between COVID-19 and ABO blood group/genotype in a large Japanese population. All Japanese patients diagnosed with COVID-19 were recruited through the Japan COVID-19 Task Force between February 2020 and October 2021. We conducted a retrospective cohort study involving 1790 Japanese COVID-19 patients whose DNA was used for a genome-wide association study. We compared the ABO blood group/genotype in a healthy population (n = 611, control) and COVID-19 patients and then analyzed their associations and clinical outcomes. Blood group A was significantly more prevalent (41.6% vs. 36.8%; P = 0.038), and group O was significantly less prevalent (26.2% vs. 30.8%; P = 0.028) in the COVID-19 group than in the control group. Moreover, genotype OO was significantly less common in the COVID-19 group. Furthermore, blood group AB was identified as an independent risk factor for most severe diseases compared with blood group O [aOR (95% CI) = 1.84 (1.00-3.37)]. In ABO genotype analysis, only genotype AB was an independent risk factor for most severe diseases compared with genotype OO. Blood group O is protective, whereas group A is associated with the risk of infection. Moreover, blood group AB is associated with the risk of the "most" severe disease.

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.

Narcolepsy type I-associated DNA methylation and gene expression changes in the human leukocyte antigen region

Narcolepsy type 1 (NT1) is caused by a loss of hypothalamic orexin-producing cells, and autoreactive CD4⁺ and CD8⁺ T cells have been suggested to play a role in the autoimmune mechanism. Although NT1 showed a strong association with human leukocyte antigen (HLA)-DQB1*06:02, the responsible antigens remain unidentified. We analyzed array-based DNA methylation and gene expression data for the HLA region in CD4⁺ and CD8⁺ T cells that were separated from the peripheral blood mononuclear cells of Japanese subjects (NT1, N = 42; control, N = 42). As the large number of SNPs in the HLA region might interfere with the affinity of the array probes, we conducted a comprehensive assessment of the reliability of each probe. The criteria were based on a previous study reporting that the presence of frequent SNPs, especially on the 3' side of the probe, makes the probe unreliable. We confirmed that 90.3% of the probes after general filtering in the HLA region do not include frequent SNPs, and are thus suitable for analysis, particularly in Japanese subjects. We then performed an association analysis, and found that several CpG sites in the HLA class II region of the patients were significantly hypomethylated in CD4⁺ and CD8⁺ T cells. This association was not detected when the effect of HLA-DQB1*06:02 was considered, suggesting that the hypomethylation was possibly derived from HLA-DQB1*06:02. Further RNA sequencing revealed reduced expression levels of HLA-DQB1 alleles other than HLA-DQB1*06:02 in the patients with NT1. Our results suggest the involvement of epigenetic and expressional changes in HLA-DQB1 in the pathogenesis of NT1.

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.

Narcolepsy risk loci outline role of T cell autoimmunity and infectious triggers in narcolepsy

Narcolepsy type 1 (NT1) is caused by a loss of hypocretin/orexin transmission. Risk factors include pandemic 2009 H1N1 influenza A infection and immunization with Pandemrix®. Here, we dissect disease mechanisms and interactions with environmental triggers in a multi-ethnic sample of 6,073 cases and 84,856 controls. We fine-mapped GWAS signals within HLA (DQ0602, DQB1*03:01 and DPB1*04:02) and discovered seven novel associations (CD207, NAB1, IKZF4-ERBB3, CTSC, DENND1B, SIRPG, PRF1). Significant signals at TRA and DQB1*06:02 loci were found in 245 vaccination-related cases, who also shared polygenic risk. T cell receptor associations in NT1 modulated TRAJ*24, TRAJ*28 and TRBV*4-2 chain-usage. Partitioned heritability and immune cell enrichment analyses found genetic signals to be driven by dendritic and helper T cells. Lastly comorbidity analysis using data from FinnGen, suggests shared effects between NT1 and other autoimmune diseases. NT1 genetic variants shape autoimmunity and response to environmental triggers, including influenza A infection and immunization with Pandemrix®.

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.

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

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

Multi-population genome-wide association study implicates immune and non-immune factors in pediatric steroid-sensitive nephrotic syndrome

Pediatric steroid-sensitive nephrotic syndrome (pSSNS) is the most common childhood glomerular disease. Previous genome-wide association studies (GWAS) identified a risk locus in the HLA Class II region and three additional independent risk loci. But the genetic architecture of pSSNS, and its genetically driven pathobiology, is largely unknown. Here, we conduct a multi-population GWAS meta-analysis in 38,463 participants (2440 cases). We then conduct conditional analyses and population specific GWAS. We discover twelve significant associations-eight from the multi-population meta-analysis (four novel), two from the multi-population conditional analysis (one novel), and two additional novel loci from the European meta-analysis. Fine-mapping implicates specific amino acid haplotypes in HLA-DQA1 and HLA-DQB1 driving the HLA Class II risk locus. Non-HLA loci colocalize with eQTLs of monocytes and numerous T-cell subsets in independent datasets. Colocalization with kidney eQTLs is lacking but overlap with kidney cell open chromatin suggests an uncharacterized disease mechanism in kidney cells. A polygenic risk score (PRS) associates with earlier disease onset. Altogether, these discoveries expand our knowledge of pSSNS genetic architecture across populations and provide cell-specific insights into its molecular drivers. Evaluating these associations in additional cohorts will refine our understanding of population specificity, heterogeneity, and clinical and molecular associations.

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.

Impact of respiratory bacterial infections on mortality in Japanese patients with COVID-19: a retrospective cohort study

BACKGROUND

Although cases of respiratory bacterial infections associated with coronavirus disease 2019 (COVID-19) have often been reported, their impact on the clinical course remains unclear. Herein, we evaluated and analyzed the complication rates of bacterial infections, causative organisms, patient backgrounds, and clinical outcome in Japanese patients with COVID-19.

METHODS

We performed a retrospective cohort study that included inpatients with COVID-19 from multiple centers participating in the Japan COVID-19 Taskforce (April 2020 to May 2021) and obtained demographic, epidemiological, and microbiological results and the clinical course and analyzed the cases of COVID-19 complicated by respiratory bacterial infections.

RESULTS

Of the 1,863 patients with COVID-19 included in the analysis, 140 (7.5%) had respiratory bacterial infections. Community-acquired co-infection at COVID-19 diagnosis was uncommon (55/1,863, 3.0%) and was mainly caused by Staphylococcus aureus, Klebsiella pneumoniae and Streptococcus pneumoniae. Hospital-acquired bacterial secondary infections, mostly caused by Staphylococcus aureus, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia, were diagnosed in 86 patients (4.6%). Severity-associated comorbidities were frequently observed in hospital-acquired secondary infection cases, including hypertension, diabetes, and chronic kidney disease. The study results suggest that the neutrophil-lymphocyte ratio (> 5.28) may be useful in diagnosing complications of respiratory bacterial infections. COVID-19 patients with community-acquired or hospital-acquired secondary infections had significantly increased mortality.

CONCLUSIONS

Respiratory bacterial co-infections and secondary infections are uncommon in patients with COVID-19 but may worsen outcomes. Assessment of bacterial complications is important in hospitalized patients with COVID-19, and the study findings are meaningful for the appropriate use of antimicrobial agents and management strategies.

Use of the neutrophil-to-lymphocyte ratio and an oxygen requirement to predict disease severity in patients with COVID-19

We investigated the association between complete blood count, including neutrophil-to-lymphocyte ratio (NLR) in combination with patient characteristics, and coronavirus disease (COVID-19) outcomes to identify the best prognostic indicator. We analyzed data of patients with confirmed COVID-19 from the nationwide database of the Japan COVID-19 Task Force between February 2020 and November 2021. A composite outcome was defined as the most severe condition, including noninvasive positive-pressure ventilation, high-flow nasal cannula, invasive mechanical ventilation, extracorporeal membrane oxygenation, or death. Of 2425 patients in the analysis, 472 (19.5%) experienced a composite outcome. NLR was the best predictor of composite outcomes, with an area under the curve (AUC) of 0.81, and a sensitivity and specificity of 72.3% and 75.7%, respectively, using a cut-off value of 5.04. The combination of NLR and an oxygen requirement on admission had the highest AUC (0.88). This simple combination may help identify patients at risk of progression to severe disease.

A case of hyperphosphatemic familial tumoral calcinosis due to maternal uniparental disomy of a GALNT3 variant

Hyperphosphatemic familial tumoral calcinosis (HFTC) is a rare, inherited autosomal recessive disorder caused by fibroblast growth factor-23 (FGF23), N-acetylgalactosaminyltransferase 3 (GALNT3), or Klotho (KL) gene variants. Here, we report the case of a Japanese boy who presented with a mass in his left elbow at the age of three. Laboratory test results of the patient revealed normocalcemia (10.3 mg/dL) and hyperphosphatemia (8.7 mg/dL); however, despite hyperphosphatemia, serum intact FGF23 level was low, renal tubular reabsorption of phosphate (TRP) level was inappropriately increased, and 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) level was inappropriately normal. Genetic analysis revealed maternal uniparental disomy (UPD) of chromosome 2, which included a novel GALNT3 variant (c.1780-1G>C). Reverse transcription-polymerase chain reaction (RT-PCR) analysis of GALNT3 mRNA confirmed that this variant resulted in the destruction of exon 11. We resected the mass when the patient was five years old, owing to its gradual enlargement. No relapse or new pathological lesions were observed four years after tumor resection. This is the first case report of a Japanese patient with HFTC associated with a novel GALNT3 variant, as well as the first case of HFTC caused by maternal UPD of chromosome 2 that includes the GALNT3 variant.

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.

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.

Impact of accumulative smoking exposure and chronic obstructive pulmonary disease on COVID-19 outcomes: report based on findings from the Japan COVID-19 task force

OBJECTIVES

Smoking and chronic obstructive pulmonary disease (COPD) are risk factors for severe COVID-19. However, limited literature exists on the effect of COPD and smoking on COVID-19 outcomes. This study examined the impact of smoking exposure in pack-years (PY) and COPD on COVID-19 outcomes among smokers in Japan.

METHODS

The study included 1266 smokers enrolled by the Japan COVID-19 task force between February 2020 and December 2021. PY and COPD status was self-reported by patients. Patients were classified into the non-COPD (n = 1151) and COPD (n = 115) groups; the non-COPD group was further classified into <10 PY (n = 293), 10-30 PY (n = 497), and >30 PY (n = 361). The study outcome was the need for invasive mechanical ventilation (IMV).

RESULTS

The incidence of IMV increased with increasing PY and was highest in the COPD group (<10 PY = 7.8%, 10-30 PY = 12.3%, >30 PY = 15.2%, COPD = 26.1%; P <0.001). A significant association was found for IMV requirement in the >30 PY and COPD groups through univariate (odds ratio [OR]: >30 PY = 2.11, COPD = 4.14) and multivariate (OR: >30 PY = 2.38; COPD = 7.94) analyses. Increasing PY number was also associated with increased IMV requirement in patients aged <65 years.

CONCLUSION

Cumulative smoking exposure was positively associated with COVID-19 outcomes in smokers.

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.

Risk for Prison-to-Community Tuberculosis Transmission, Thailand, 2017-2020

To determine contributions of previously incarcerated persons to tuberculosis (TB) transmission in the community, we performed a healthcare facility-based cohort study of TB patients in Thailand during 2017-2020. We used whole-genome sequencing of Mycobacterium tuberculosis isolates from patients to identify genotypic clusters and assess the association between previous incarceration and TB transmission in the community. We identified 4 large genotype clusters (>10 TB patients/cluster); 28% (14/50) of the patients in those clusters were formerly incarcerated. Formerly incarcerated TB patients were more likely than nonincarcerated patients to be included in large clusters. TB patients within the large genotype clusters were geographically dispersed throughout Chiang Rai Province. Community TB transmission in the community was associated with the presence of formerly incarcerated individuals in Thailand. To reduce the risk for prison-to-community transmission, we recommend TB screening at the time of entry and exit from prisons and follow-up screening in the community.

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

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

Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10-8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

Corrigendum: Intraocular human cytomegaloviruses of ocular diseases are distinct from those of viremia and are capable of escaping from innate and adaptive immunity by exploiting HLA-E-mediated peripheral and central tolerance

[This corrects the article DOI: 10.3389/fimmu.2022.1008220.].

Experimental testing of a prototype cantilevered liquid-nitrogen-cooled silicon mirror

This report presents testing of a prototype cantilevered liquid-nitrogen-cooled silicon mirror. This mirror was designed to be the first mirror for the new soft X-ray beamlines to be built as part of the Advanced Light Source Upgrade. Test activities focused on fracture, heat transfer, modal response and distortion, and indicated that the mirror functions as intended.

The KANNO blood group system

The KANNO blood group system (International Society of Blood Transfusion [ISBT] 037) includes one high-prevalence antigen, KANNO1, across ethnic groups. Sporadic KANNO1- cases among East and South Asians are theoretically estimated by the DNA database library. Anti-KANNO1 has been found most often among Japanese women with current or prior pregnancy. Thus far, there are no reported cases of hemolytic transfusion reaction or hemolytic disease of the fetus and newborn due to anti-KANNO1.

Impact of upper and lower respiratory symptoms on COVID-19 outcomes: a multicenter retrospective cohort study

Background

Respiratory symptoms are associated with coronavirus disease 2019 (COVID-19) outcomes. However, the impacts of upper and lower respiratory symptoms on COVID-19 outcomes in the same population have not been compared. The objective of this study was to characterize upper and lower respiratory symptoms and compare their impacts on outcomes of hospitalized COVID-19 patients.

Methods

This was a multicenter, retrospective cohort study; the database from the Japan COVID-19 Task Force was used. A total of 3314 COVID-19 patients were included in the study, and the data on respiratory symptoms were collected. The participants were classified according to their respiratory symptoms (Group 1: no respiratory symptoms, Group 2: only upper respiratory symptoms, Group 3: only lower respiratory symptoms, and Group 4: both upper and lower respiratory symptoms). The impacts of upper and lower respiratory symptoms on the clinical outcomes were compared. The primary outcome was the percentage of patients with poor clinical outcomes, including the need for oxygen supplementation via high-flow oxygen therapy, mechanical ventilation, and extracorporeal membrane oxygenation or death.

Results

Of the 3314 COVID-19 patients, 605, 1331, 1229, and 1149 were classified as Group 1, Group 2, Group 3, and Group 4, respectively. In univariate analysis, patients in Group 2 had the best clinical outcomes among all groups (odds ratio [OR]: 0.21, 95% confidence interval [CI]: 0.11–0.39), while patients in Group 3 had the worst outcomes (OR: 3.27, 95% CI: 2.43–4.40). Group 3 patients had the highest incidence of pneumonia, other complications due to secondary infections, and thrombosis during the clinical course.

Conclusions

Upper and lower respiratory tract symptoms had vastly different impacts on the clinical outcomes of COVID-19.

National Center Biobank Network

There are six national centers (6NCs) for advanced and specialized medicine in Japan that conduct basic and clinical research on major diseases that have a substantial impact on national health. Disease-specific bioresources and information collected by each NC are stored in a separate biobank. The National Center Biobank Network (NCBN) was established in 2011 and coordinates the biobanks and researchers of the 6NCs via an open-access database (Catalogue Database: http://www2.ncbiobank.org/Index_en ) as an efficient means of providing registered biological resources and data for use in research communities. The NCBN resources are characterized by their high-quality and rich medical information and are available for life science research and for the development of novel testing methodologies (biomarkers), new treatments, and drugs for future health care in the scope of personalized medicine through a deeper understanding of disease pathogenesis. Here, we explain the activities of the NCBN and the characteristics of the NCBN Catalogue Database.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Intraocular human cytomegaloviruses of ocular diseases are distinct from those of viremia and are capable of escaping from innate and adaptive immunity by exploiting HLA-E-mediated peripheral and central tolerance

Human cytomegalovirus (HCMV) infections develop into CMV diseases that result in various forms of manifestations in local organs. CMV-retinitis is a form of CMV disease that develops in immunocompromised hosts with CMV-viremia after viruses in the peripheral circulation have entered the eye. In the HCMV genome, extensive diversification of the UL40 gene has produced peptide sequences that modulate NK cell effector functions when loaded onto HLA-E and are subsequently recognized by the NKG2A and NKG2C receptors. Notably, some HCMV strains carry UL40 genes that encode peptide sequences identical to the signal peptide sequences of specific HLA-A and HLA-C allotypes, which enables these CMV strains to escape HLA-E-restricted CD8+T cell responses. Variations in UL40 sequences have been studied mainly in the peripheral blood of CMV-viremia cases. In this study, we sought to investigate how ocular CMV disease develops from CMV infections. CMV gene sequences were compared between the intraocular fluids and peripheral blood of 77 clinical cases. UL40 signal peptide sequences were more diverse, and multiple sequences were typically present in CMV-viremia blood compared to intraocular fluid. Significantly stronger NK cell suppression was induced by UL40-derived peptides from intraocular HCMV compared to those identified only in peripheral blood. HCMV present in intraocular fluids were limited to those carrying a UL40 peptide sequence corresponding to the leader peptide sequence of the host's HLA class I, while UL40-derived peptides from HCMV found only in the peripheral blood were disparate from any HLA class I allotype. Overall, our analyses of CMV-retinitis inferred that specific HCMV strains with UL40 signal sequences matching the host's HLA signal peptide sequences were those that crossed the blood-ocular barrier to enter the intraocular space. UL40 peptide repertoires were the same in the intraocular fluids of all ocular CMV diseases, regardless of host immune status, implying that virus type is likely to be a common determinant in ocular CMV disease development. We thus propose a mechanism for ocular CMV disease development, in which particular HCMV types in the blood exploit peripheral and central HLA-E-mediated tolerance mechanisms and, thus, escape the antivirus responses of both innate and adaptive immunity.

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.

The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force

Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection.

Genetic mechanisms influencing COVID-19 susceptibility are not well understood. Here, the authors analyzed whole blood RNA-seq data of 465 Japanese individuals with COVID-19, highlighting thousands of fine-mapped variants affecting expression and splicing of genes, as well as the presence of COVID-19 severity-interaction eQTLs.

Low carnitine palmitoyltransferase 1 activity is a risk factor for narcolepsy type 1 and other hypersomnia

Study Objectives

Narcolepsy type 1 (NT1) is associated with metabolic abnormalities but their etiology remains largely unknown. The gene for carnitine palmitoyltransferase 1B (CPT1B) and abnormally low serum acylcarnitine levels have been linked to NT1. To elucidate the details of altered fatty acid metabolism, we determined levels of individual acylcarnitines and evaluated CPT1 activity in patients with NT1 and other hypersomnia.

Methods

Blood samples from 57 NT1, 51 other hypersomnia patients, and 61 healthy controls were analyzed. The levels of 25 major individual acylcarnitines were determined and the C0/(t[C16] + t[C18]) ratio was used as a CPT1 activity marker. We further performed transcriptome analysis using independent blood samples from 42 NT1 and 42 healthy controls to study the relevance of fatty acid metabolism. NT1-specific changes in CPT1 activity and in expression of related genes were investigated.

Results

CPT1 activity was lower in patients with NT1 (p = 0.00064) and other hypersomnia (p = 0.0014) than in controls. Regression analysis revealed that CPT1 activity was an independent risk factor for NT1 (OR: 1.68; p = 0.0031) and for other hypersomnia (OR: 1.64; p = 0.0042). There was a significant interaction between obesity (BMI &lt;25, ≥25) and the SNP rs5770917 status such that nonobese NT1 patients without risk allele had better CPT1 activity (p = 0.0089). The expression levels of carnitine-acylcarnitine translocase (CACT) and CPT2 in carnitine shuttle were lower in NT1 (p = 0.000051 and p = 0.00014, respectively).

Conclusions

These results provide evidences that abnormal fatty acid metabolism is involved in the pathophysiology of NT1 and other hypersomnia.

O-301 Genome-wide association study identified meiotic variant associated with aneuploid pregnancy loss

Study question

Which single nucleotide variant (SNVs) are associated with aneuploid pregnancy loss?

Summary answer

We identified a SNV on MEIG1 gene, which are associated with meiosis/spermiogenesis.

What is known already

Recurrent pregnancy loss (RPL) refers to the loss of two or more pregnancies, with a frequency of 5%. Chromosomal abnormalities in embryos are found in 80% of first trimester miscarriages, 86% of which are aneuploid. Recently, embryonic aneuploidy was found to be the most common cause of RPL, with a frequency of 40-50%. Most trisomy miscarriages are of maternal origin, with errors occurring during meiosis of the oocyte. Chromosome segregation abnormalities in oocytes are thought to be an event associated with increasing maternal age, but in addition, maternal genetic causes are thought to contribute.

Study design, size, duration

A Genome wide association study (GWAS) was performed on a clinically well characterized cohort of 189 women with RPL whose previous aborted conceptus was ascertained to be an aneuploid embryo. Samples were mainly collected from 2007 to 2018 mainly at Nagoya City University Hospital. For control samples, we used 1157 samples from the population-based prospective cohorts that included fertile women.

Participants/materials, setting, methods

All patients underwent a systematic examination. Patients with antiphospholipid syndrome, an abnormal chromosome in either partner, or uterine anomaly were excluded. Patients whose previously miscarried POC exhibited triploidy or 45, X were excluded. DNA was isolated from stored EDTA-blood samples and genotyped by Axiom Japonica-array v2659,503 SNVs). For the GWAS, a chi-squared test was applied to a two-by-two contingency table in allele frequency model.

Main results and the role of chance

The mean (SD) ages and number of previous miscarriages of the patients were 36.8 (4.3) and 3.09 (1.13). GWAS data revealed 5 SNVs with suggestive significance (p &lt; 9.46e-06). The SNVs that showed the most significant associations (P = 1.06E-06, OR = 1.72) was located on meiosis/spermiogenesis associated 1 (MEIG1) gene under an allelic model after Bonferroni correction considering the number of analyzed SNVs. The SNV rs7908491 was reported as a splicing QTL in the MEIG1 gene, which is a meiosis/meiosis-associated factor and is plausibly associated with chromosome aneuploidy. This is the first GWAS in patients with RPL caused by aneuploidy.

Limitations, reasons for caution

Since this study was conducted in a single center and had a small sample size, it needs to be replicated in different centers with more subjects and on an international scale. Whole genome imputation analysis will be performed to detect SNVs with more significant associations.

Wider implications of the findings

Our findings demonstrate that a specific genotype of MEIG1 gene can be a risk factor for aneuploid pregnancy loss. The establishment of clinically applicable maternal germ cell markers could identify groups for whom PGT would be more useful or provide patients with counseling that provides prognostic information about pregnancy.

Trial registration number

not applicable

Phasing analysis of lung cancer genomes using a long read sequencer

Chromosomal backgrounds of cancerous mutations still remain elusive. Here, we conduct the phasing analysis of non-small cell lung cancer specimens of 20 Japanese patients. By the combinatory use of short and long read sequencing data, we obtain long phased blocks of 834 kb in N50 length with >99% concordance rate. By analyzing the obtained phasing information, we reveal that several cancer genomes harbor regions in which mutations are unevenly distributed to either of two haplotypes. Large-scale chromosomal rearrangement events, which resemble chromothripsis events but have smaller scales, occur on only one chromosome, and these events account for the observed biased distributions. Interestingly, the events are characteristic of EGFR mutation-positive lung adenocarcinomas. Further integration of long read epigenomic and transcriptomic data reveal that haploid chromosomes are not always at equivalent transcriptomic/epigenomic conditions. Distinct chromosomal backgrounds are responsible for later cancerous aberrations in a haplotype-specific manner.

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.

Protective association of HLA-DPB1*04:01:01 with acute encephalopathy with biphasic seizures and late reduced diffusion identified by HLA imputation

Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a severe syndrome of acute encephalopathy that affects infants and young children. AESD is a polygenic disorder preceded by common viral infections with high fever. We conducted an association study of human leukocyte antigen (HLA) regions with AESD using HLA imputation. SNP genotyping was performed on 254 Japanese patients with AESD and 799 healthy controls. We conducted 3-field HLA imputation for 14 HLA genes based on Japanese-specific references using data from our previous genome-wide association study. After quality control, 208 patients and 737 controls were included in the analysis of HLA alleles. We then compared the carrier frequencies of HLA alleles and haplotypes between the patients and controls. HLA-DPB1*04:01:01 showed a significant association with AESD, exerting a protective effect against the disease (p = 0.0053, p_corrected = 0.042, odds ratio = 0.43, 95% confidence interval = 0.21-0.80). The allele frequency of HLA-DPB1*04:01:01 was lower in East Asians than in Caucasians, which may partially account for the higher incidence of AESD in the Japanese population. The present results demonstrate the importance of fine-mapping of the HLA region to investigate disease susceptibilities and elucidate the pathogenesis of AESD.