Search

Abstract We report a case of Wilson disease (WD) with dilated cardiomyopathy in which whole-genome sequencing (WGS) revealed the rare co-occurrence of two novel compound heterozygous ATP7B pathogenic variants (NM_001005918.3:c.2250del/p.N751Tfs*9 and c.3496C>T/p.L1166F) and a known FLNC pathogenic variant. Our results highlight the usefulness of WGS, even in the diagnosis of well-characterized genetic diseases such as WD.

Abstract Lemierre's syndrome is a fatal condition characterized by septic thrombophlebitis of the internal jugular vein (IJV). Early diagnosis of Lemierre's syndrome in primary care remains challenging because of its rarity and common initial symptoms like fever and sore throat. The COVID‐19 pandemic has increased the difficulty of diagnosis, along with many diagnostic delays. Herein, we report a case in which an IJV thrombus was detected using point‐of‐care ultrasound (POCUS), leading to the diagnosis of Lemierre's syndrome. During the COVID‐19 pandemic, POCUS was a simple and noninvasive technique for screening Lemierre's syndrome without radiation and low droplet exposure.

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

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

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

Purpose: Achromatopsia (ACHM) is an early-onset cone dysfunction caused by 5 genes with cone-specific functions (CNGA3, CNGB3, GNAT2, PDE6C, and PDE6H) and by ATF6, a transcription factor with ubiquitous expression. To improve the relatively low variant detection ratio in these genes in a cohort of exome-sequenced Japanese patients with inherited retinal diseases (IRD), we performed genome sequencing to detect structural variants and intronic variants in patients with ACHM. Methods: Genome sequencing of 10 ACHM pedigrees was performed after exome sequencing. Structural, non-coding, and coding variants were filtered based on segregation between the affected and unaffected in each pedigree. Variant frequency and predicted damage scores were considered in identifying pathogenic variants. Results: A homozygous deletion involving exon 18 of RPGRIP1 was detected in 5 of 10 ACHM probands, and variant inheritance from each parent was confirmed. This deletion was relatively frequent (minor allele frequency = 0.0023) in the Japanese population but was only homozygous in patients with ACHM among the 199 Japanese IRD probands analyzed by the same genome sequencing pipeline. Conclusion: The deletion involving exon 18 of RPGRIP1 is a prevalent cause of ACHM in Japanese patients and contributes to the wide spectrum of RPGRIP1-associated IRD phenotypes, from Leber congenital amaurosis to ACHM.

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

Abstract Accurate genotype imputation requires large-scale reference panel datasets. When conducting genotype imputation on the Japanese population, researchers can use such datasets under collaborative studies or controlled access conditions in public databases. We developed the NBDC-DDBJ imputation server, which securely provides users with a web user interface to execute genotype imputation on the server. Our benchmarking analysis showed that the accuracy of genotype imputation was improved by leveraging controlled access datasets to increase the number of haplotypes available for analysis compared to using publicly available reference panels such as the 1000 Genomes Project. The NBDC-DDBJ imputation server facilitates the use of controlled access datasets for accurate genotype imputation.

Abstract Background Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. Results Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P

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.

Long-read sequencing technologies are useful for the multifaceted task of characterising somatic mutations, including structural variants, in cancers. Here, the authors combine short and long read sequencing for the phasing analysis, which enables them to resolve the chromosomal backgrounds of somatic mutations in Japanese non-small cell lung cancers.

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

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

Abstract Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a severe encephalopathy preceded by viral infections with high fever. AESD is a multifactorial disease, however, few disease susceptibility genes have previously been identified. Here, we conducted a genome-wide association study (GWAS) and assessed functional variants in non-coding regions to study genetic susceptibility in AESD using 254 Japanese children with AESD and 799 adult healthy controls. We also performed a microRNA enrichment analysis using GWAS statistics to search for candidate biomarkers in AESD. The variant with the lowest p-value, rs1850440, was located in the intron of serine/threonine kinase 39 gene (STK39) on chromosome 2q24.3 (p = 2.44 × 10−7, odds ratio = 1.71). The minor allele T of rs1850440 correlated with the stronger expression of STK39 in peripheral blood. This variant possessed enhancer histone modification marks in STK39, the encoded protein of which activates the p38 mitogen-activated protein kinase (MAPK) pathway. In the replication study, the odds ratios of three SNPs, including rs1850440, showed the same direction of association with that in the discovery stage GWAS. One of the candidate microRNAs identified by the microRNA enrichment analysis was associated with inflammatory responses regulated by the MAPK pathway. This study identified STK39 as a novel susceptibility locus of AESD, found microRNAs as potential biomarkers, and implicated immune responses and the MAPK cascade in its pathogenesis.

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.

Abstract Background Cold agglutinin disease can cause the agglutination of red blood cells and hemolytic anemia due to cold temperature. Herein, we report a case of progressive hemolytic anemia due to cold agglutinin disease during fluid resuscitation and in the absence of exposure to cold. Case presentation A 71-year-old Japanese man was admitted to the emergency department with signs of hypotension and disturbed consciousness. He was diagnosed with diabetic ketoacidosis, and treatment with fluid resuscitation and insulin infusion was initiated. Laboratory test results obtained the following day indicated hemolytic anemia. On day 5 after admission, red blood cell agglutination was detected, and the patient was diagnosed with cold agglutinin disease. Conclusions Cold agglutinin disease should be considered in the differential diagnosis of progressive hemolytic anemia during fluid resuscitation, even if the solution is at room temperature.

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

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

Abstract Stevens–Johnson syndrome (SJS) and its severe condition with extensive skin detachment and a poor prognosis, toxic epidermal necrolysis (TEN), are immunologically mediated severe cutaneous reactions of the skin and mucous membranes such as the ocular surface. Genetic variations on the HLA-A and other autosomal genes have been identified as risk factors for cold medicine-related SJS/TEN with severe ocular complications (CM-SJS/TEN with SOC). Using a whole-genome sequencing (WGS) approach, we explored other susceptible variants of CM-SJS/TEN with SOC, especially among rare variants and structural variants (SVs). WGS was performed on samples from 133 patients with CM-SJS/TEN with SOC and 418 healthy controls to obtain single nucleotide polymorphisms (SNPs) and SVs. Genome-wide association tests were performed with these variants. Our genome-wide association test reproduced the associations of the common variants of HLA-A and loci on chromosome 16q12.1. We also identified novel associations of SVs on these loci and an aggregation of rare coding variants on the TPRM8 gene. In silico gene expression analysis on the HLA-A locus revealed that the SNP (rs12202296), which was significantly associated with susceptibility to CM-SJS/TEN with SOC, was correlated to an increase in HLA-A expression levels in the whole blood (P = 2.9 × 10−17), from the GTEx database. The majority of variants that were significantly associated with CM-SJS/TEN with SOC were found in non-coding regions, indicating the regulatory role of genetic variations in the pathogenesis of CM-SJS/TEN with SOC.

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−10, OR = 0.49) with HBV chronicity. Subsequent HLA allele imputation revealed HLA-DPA1*01:03 (Pc = 1.21 × 10−6, OR = 0.53), HLA-DPB1*02:01 (Pc = 2.17 × 10−3, OR = 0.50), and HLA-DQB1*06:09 (Pc = 2.17 × 10−2, OR = 0.07) as protective alleles, and HLA-DPA1*02:02 (Pc = 6.32 × 10−5, OR = 1.63), HLA-DPB1*05:01 (Pc = 1.13 × 10−4, OR = 1.72), HLA-DPB1*13:01 (Pc = 4.68 × 10−2, OR = 1.60), and HLA-DQB1*03:03 (Pc = 1.11 × 10−3, 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.

Background: The utility of endoscopic transplantation of epithelial cell sheets to ulcer sites after endoscopic submucosal dissection (ESD) has been shown to prevent scar stenosis after ESD of early esophageal cancer. Previously, our group reported use of an endoscopic transplantation device fabricated with a 3-dimensional printer. Cell sheets are transplanted to the esophageal wound site with the following procedure: first, a cell sheet harvested from temperature-responsive culture dishes is placed on the device's deflated balloon surface and transported to the wound site with endoscopic forceps; second, by applying pressure from inflating the balloon locally at the wound site, the cell sheet is successfully transferred and adhered to the wound tissue; third, the balloon is deflated, and the device is removed. By repeating the procedure, several cell sheets can be safely transplanted to a wider ESD area. Nonetheless, possible damage to cell sheets using this procedure has not yet been assessed. Objective: Effects of endoscopic transplantation balloon inflation on cell viability and damage of normal human epidermal keratinocyte sheets resident on the device's balloon surface were evaluated by histology after sheet placement onto lumenal surfaces in the ex vivo porcine submucosal dissection esophagus model. Endoscopic transplantation of these same cell sheets with conventional methods using a polyvinylidene fluoride (PVDF) cell sheet support membrane, balloon device transfer, and also using a novel modified balloon transfer procedure was also examined. Cell sheet transfer results obtained with these three procedures were compared. Method: Normal human epidermal keratinocyte sheets were fabricated on temperature-responsive culture inserts. By temperature reduction to 20 °C, all cells were harvested as a single contiguous cell sheet. Freshly excised porcine esophagi purchased in a slaughter house were turned inside-out, and the exposed lumenal mucosa and submucosal layers were removed by Cooper scissors. This luminal surface was then utilized as a transplantation bed in ex vivo cell sheet experiments. Cell sheets were adhered to the endoscopic transfer device balloon, expanded by balloon inflation and resulting cell viability was evaluated by trypan blue exclusion test after cell sheet trypsinization and dispersion. Cell sheets were transferred onto the esophagus lumen ex vivo using forceps and the balloon device, and also using a modified balloon transfer method. The obtained results were compared with those without balloon expansion, and evaluated for sheet thickness and lumenal histology. Finally, TUNEL staining was performed to examine cell apoptosis. Result: Cell sheets thinned after cell sheet balloon expansion, but no apoptosis was observed after these procedures. Conclusion: Expanding keratinocyte cell sheets on a balloon endoscopic transfer device did not damage the cell sheets. This sheet transplantation method using the endoscopic balloon transfer device may be considered as a future standard cell sheet endoscopic transplantation procedure for repairing the esophagus.

Abstract Background Genetic factors such as single-nucleotide polymorphisms (SNPs) play a key role in the development of postoperative nausea and vomiting (PONV). However, previous findings are not widely applicable to different populations because of population-specific genetic variation. We developed a Japanese-specific DNA microarray for high-throughput genotyping. The aim of the current study was to identify SNPs associated with PONV on a genome-wide scale using this microarray in a sample of Japanese surgical patients. Methods Associations between 659,636 SNPs and the incidence of PONV 24 h after surgery in a limited sample of 24 female patients were assessed using the microarray. After imputation of genotypes at 24,330,529 SNPs, 78 SNPs were found to be associated with the incidence of PONV. We chose 4 of the 78 SNPs to focus on by in silico functional annotation. Finally, we genotyped these 4 candidate SNPs in 255 patients using real-time PCR to verify association with the incidence of PONV. Results The T > C variant of rs11232965 in the long non-coding RNA MIR4300HG was significantly associated with reduced incidence of PONV among genotypes and between alleles (p = 0.01 and 0.007). Conclusions We identified a novel SNP (rs11232965) in the long non-coding RNA MIR4300HG that is associated with PONV. The rs11232965-SNP variant (T > C) is protective against the incidence of PONV. Trial registration This study was registered at the UMIN Clinical Trials Registry (Identifier: UMIN000022903 , date of registration: June 27, 2016, retrospectively registered.

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

Over the past few decades, obesity has become a public health issue of global concern. Even though disparities exist between human populations, e.g., the higher liver fat content of the Japanese despite a lower body mass index (BMI), studies on the genetics of obesity still largely focus on populations of European descent, leading to a dearth of genetic data on non-European populations. In this context, this study aimed to establish a broad picture of the genetic attributes of the Japanese population, by examining a representative sample of 18,889 individuals participating in the Tohoku Medical Megabank Project cohort. We applied linear mixed model methods to 17 traits related to obesity and associated diseases to estimate the heritabilities explained by common genetic variants and the genetic correlations between each pair of traits. These analyses allowed us to quantify the SNP heritability of health indicators such as BMI (0.248 ± 0.032) and HDL cholesterol (0.324 ± 0.031), and to provide one of the few estimates of the SNP heritability of cystatin C in unrelated individuals (0.260 ± 0.025). We discuss potential differences between the Japanese and people of European ancestry with respect to the genetic correlations between urinary biomarkers and adiposity traits, for which large estimates were obtained. For instance, the genetic correlations between urine potassium level and the values for weight, BMI, waist circumference, and waist-to-height ratio ranged from 0.290 to 0.559, much higher than the corresponding estimates in the UK Biobank.

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 2 versus controls (P

Chalcone synthase is the first committed enzyme in the plant flavonoid biosynthesis pathway, yet shows low product specificity in vitro. Here Waki et al. show that chalcone isomerase-like proteins bind to and reduce the catalytic promiscuity of chalcone synthase, ensuring efficient flavonoid production in planta.

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

Abstract Background Kawasaki disease (KD) is a systemic vasculitis which may be associated with coronary artery aneurysms. A notable risk factor for the development of coronary artery aneurysms is resistance to intravenous immunoglobulin (IVIG) therapy, which comprises standard treatment for the acute phase of KD. The cause of IVIG resistance in KD is largely unknown; however, the contribution of genetic factors, especially variants in immune-related genes, has been suspected. Methods To explore genetic variants related to IVIG-unresponsiveness, we designated KD patients who did not respond to both first and second courses of IVIG therapy as IVIG-unresponsive patients. Using genomic DNA from 30 IVIG-unresponsive KD patients, we performed pooled genome sequencing targeting 39 immune-related cytokine receptor genes. Results The single nucleotide variant (SNV), rs563535954 (located in the IL4R locus), was concentrated in IVIG-unresponsive KD patients. Individual genotyping showed that the minor allele of rs563535954 was present in 4/33 patients with IVIG-unresponsive KD, compared with 20/1063 individuals in the Japanese genome variation database (odds ratio = 7.19, 95% confidence interval 2.43–21.47). Furthermore, the minor allele of rs563535954 was absent in 42 KD patients who responded to IVIG treatment (P = 0.0337), indicating that a low-frequency variant, rs563535954, is associated with IVIG-unresponsiveness in KD patients. Although rs563535954 is located in the 3′-untranslated region of IL4R, there was no alternation in IL4R expression associated with the mior allele of rs563535954. However, IVIG-unresponsive patients that exhibited the minor allele of rs563535954 tended to be classified into the low-risk group (based on previously reported risk scores) for prediction of IVIG-resistance. Therefore, IVIG-unresponsiveness associated with the minor allele of rs563535954 might differ from IVIG-unresponsiveness associated with previous risk factors used to evaluate IVIG-unresponsiveness in KD. Conclusion These findings suggest that the SNV rs563535954 could serve as a predictive indicator of IVIG-unresponsiveness, thereby improving the sensitivity of risk scoring systems, and may aid in prevention of coronary artery lesions in KD patients.

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. Keywords: PBC, GWAS, Disease Susceptibility Gene, e-QTL, CRISPR/Cas9

Abstract The associations of single nucleotide polymorphisms (SNPs) in PLA2R1 and HLA-DQA1, as well as HLA-DRB1*15:01-DQB1*06:02 haplotype with idiopathic membranous nephropathy (IMN) is well known. However, the primary associations of these loci still need to be determined. We used Japanese-specific SNP genotyping array and imputation using 2,048 sequenced Japanese samples to fine-map PLA2R1 region in 98 patients and 413 controls. The most significant SNPs were replicated in a separate sample set of 130 patients and 288 controls. A two-SNP haplotype of intronic and missense SNPs showed the strongest association. The intronic SNP is strongly associated with PLA2R1 expression in the Genotype-Tissue Expression (GTEx) database, and the missense SNP is predicted to alter peptide binding with HLA-DRB1*15:01 by the Immune Epitope Database (IEDB). In HLA region, we performed relative predispositional effect (RPE) tests and identified additional risk alleles in both HLA-DRB1 and HLA-DQB1. We collapsed the risk alleles in each of HLA-DRB1 and HLA-DQB1 into single risk alleles. Reciprocal conditioning of these collapsed risk alleles showed more residual significance for HLA-DRB1 collapsed risk than HLA-DQB1 collapsed risk. These results indicate that changes in the expression levels of structurally different PLA2R protein confer risk for IMN in the presence of risk HLA-DRB1 alleles.

Abstract Background Genotype imputation from single-nucleotide polymorphism (SNP) genotype data using a haplotype reference panel consisting of thousands of unrelated individuals from populations of interest can help to identify strongly associated variants in genome-wide association studies. The Tohoku Medical Megabank (TMM) project was established to support the development of precision medicine, together with the whole-genome sequencing of 1070 human genomes from individuals in the Miyagi region (Northeast Japan) and the construction of the 1070 Japanese genome reference panel (1KJPN). Here, we investigated the performance of 1KJPN for genotype imputation of Japanese samples not included in the TMM project and compared it with other population reference panels. Results We found that the 1KJPN population was more similar to other Japanese populations, Nagahama (south-central Japan) and Aki (Shikoku Island), than to East Asian populations in the 1000 Genomes Project other than JPT, suggesting that the large-scale collection (more than 1000) of Japanese genomes from the Miyagi region covered many of the genetic variations of Japanese in mainland Japan. Moreover, 1KJPN outperformed the phase 3 reference panel of the 1000 Genomes Project (1KGPp3) for Japanese samples, and IKJPN showed similar imputation rates for the TMM and other Japanese samples for SNPs with minor allele frequencies (MAFs) higher than 1%. Conclusions 1KJPN covered most of the variants found in the samples from areas of the Japanese mainland outside the Miyagi region, implying 1KJPN is representative of the Japanese population’s genomes. 1KJPN and successive reference panels are useful genome reference panels for the mainland Japanese population. Importantly, the addition of whole genome sequences not included in the 1KJPN panel improved imputation efficiencies for SNPs with MAFs under 1% for samples from most regions of the Japanese archipelago.

Abstract Numerous genome-wide association studies (GWAS) have been performed to identify susceptibility genes to various human complex diseases. However, in many cases, neither a functional variant nor a disease susceptibility gene have been clarified. Here, we show an efficient approach for identification of a functional variant in a primary biliary cholangitis (PBC)-susceptible region, chromosome 17q12-21 (ORMDL3-GSDMB-ZPBP2-IKZF3). High-density association mapping was carried out based on SNP imputation analysis by using the whole-genome sequence data from a reference panel of 1,070 Japanese individuals (1KJPN), together with genotype data from our previous GWAS (PBC patients: n = 1,389; healthy controls: n = 1,508). Among 23 single nucleotide polymorphisms (SNPs) with P

[This corrects the article DOI: 10.1371/journal.pone.0194036.].

Forensic DNA typing is widely used to identify missing persons and plays a central role in forensic profiling. DNA typing usually uses capillary electrophoresis fragment analysis of PCR amplification products to detect the length of short tandem repeat (STR) markers. Here, we analyzed whole genome data from 1,070 Japanese individuals generated using massively parallel short-read sequencing of 162 paired-end bases. We have analyzed 843,473 STR loci with two to six basepair repeat units and cataloged highly polymorphic STR loci in the Japanese population. To evaluate the performance of the cataloged STR loci, we compared 23 STR loci, widely used in forensic DNA typing, with capillary electrophoresis based STR genotyping results in the Japanese population. Seventeen loci had high correlations and high call rates. The other six loci had low call rates or low correlations due to either the limitations of short-read sequencing technology, the bioinformatics tool used, or the complexity of repeat patterns. With these analyses, we have also purified the suitable 218 STR loci with four basepair repeat units and 53 loci with five basepair repeat units both for short read sequencing and PCR based technologies, which would be candidates to the actual forensic DNA typing in Japanese population.

Important features to be included in queueing-theoretic models of the call center operation are multiple servers, impatient customers, time-varying arrival process, and operator’s after-call work (ACW). We propose a fluid approximation technique for the queueing model with these features by extending the analysis of a similar model without ACW recently developed by Liu and Whitt (2012). Our model assumes that the service for each quantum of fluid consists of a sequence of two stages, the first stage for the conversation with a customer and the second stage for the ACW. When the duration of each stage has exponential, hyperexponential or hypo-exponential distribution, we derive the time-dependent behavior of the content of fluid in each stage of service as well as that in the waiting room. Numerical examples are shown to illustrate the system performance for the cases in which the input rate and/or the number of servers vary in sinusoidal fashion as well as in adaptive ways and in stationary cases.