Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Divergent neurodevelopmental profiles of very-low-birth-weight infants.

BACKGROUND: Advanced perinatal medicine has decreased the mortality rate of preterm infants. Long-term neurodevelopmental outcomes of very-low-birth-weight infants (VLBWIs) remain to be investigated. METHODS: Participants were 124 VLBWIs who had in-hospital birth from 2007 to 2015. Perinatal information, developmental or intelligence quotient (DQ/IQ), and neurological comorbidities at ages 3 and 6 years were analyzed. RESULTS: Fifty-eight (47%) VLBWIs received neurodevelopmental assessments at ages 3 and 6 years. Among them, 15 (26%) showed DQ/IQ <75 at age 6 years. From age 3 to 6 years, 21 (36%) patients showed a decrease (≤-10), while 5 (9%) showed an increase (≥+10) in DQ/IQ scores. Eight (17%) with autism spectrum disorder or attention-deficit hyperactivity disorder (ASD/ADHD) showed split courses of DQ/IQ, including two with ≤-10 and one with +31 to their scores. On the other hand, all 7 VLBWIs with cerebral palsy showed DQ ≤35 at these ages. Magnetic resonance imaging detected severe brain lesions in 7 (47%) of those with DQ <75 and 1 (18%) with ASD/ADHD. CONCLUSIONS: VLBWIs show a broad spectrum of neurodevelopmental outcomes after 6 years. These divergent profiles also indicate that different risks contribute to the development of ASD/ADHD from those of cerebral palsy and epilepsy in VLBWIs. IMPACT: Very-low-birth-weight infants (VLBWIs) show divergent neurodevelopmental outcomes from age 3 to 6 years. A deep longitudinal study depicts the dynamic change in neurodevelopmental profiles of VLBWIs from age 3 to 6 years. Perinatal brain injury is associated with developmental delay, cerebral palsy and epilepsy, but not with ASD or ADHD at age 6 years.

HNF4α is required for Tkfc promoter activation by ChREBP.

Triokinase/FMN cyclase (Tkfc) is involved in fructose metabolism and is responsible for the phosphorylation of glyceraldehyde to glyceraldehyde-3-phosphate. In this study, we show that refeeding induced hepatic expression of Tkfc in mice. Luciferase reporter gene assays using the Tkfc promoter revealed the existence of two hepatocyte nuclear factor 4α (HNF4α)-responsive elements (HNF4RE1 and HNF4RE2) and one carbohydrate-responsive element-binding protein (ChREBP)-responsive element (ChoRE1). Deletion and mutation of HNF4RE1 and HNF4RE2 or ChoRE1 abolished HNF4α and ChREBP responsiveness, respectively. HNF4α and ChREBP synergistically stimulated Tkfc promoter activity. ChoRE1 mutation attenuated but maintained HNF4α responsiveness, whereas HNF4RE1 and HNF4RE2 mutations abolished ChREBP responsiveness. Moreover, Tkfc promoter activity stimulation by ChREBP was attenuated upon HNF4α knockdown. Furthermore, Tkfc expression was decreased in livers of ChREBP-/- and liver-specific HNF4-/- (Hnf4αΔHep) mice. Altogether, our data indicate that Tkfc is a target gene of ChREBP and HNF4α, and Tkfc promoter activity stimulation by ChREBP requires HNF4α.

Utility of Positional Instillation of Contrast Cystography for Diagnosing Occult Vesicoureteral Reflux in Children: A Report of Two Cases.

Positional instillation of contrast (PIC) cystography is effective for detecting occult vesicoureteral reflux (VUR), which can not be revealed by standard voiding cystourethrography (VCUG). We experienced two cases of young female patients; one had repeated urinary tract infection with a negative VUR on standard VCUG, and the other had findings suggestive of reflux hydronephrosis and intolerance of standard VCUG. They underwent PIC cystography, and occult VUR was detected in both cases. Both were successfully treated with simultaneous endoscopic injection therapy with dextranomer/hyaluronic acid. PIC cystography is useful for detecting occult VUR in children with negative VUR findings on standard VCUG or who are unable to tolerate standard VCUG.

Statistical testing in transcriptomic-neuroimaging studies: A how-to and evaluation of methods assessing spatial and gene specificity.

Multiscale integration of gene transcriptomic and neuroimaging data is becoming a widely used approach for exploring the molecular underpinnings of large-scale brain organization in health and disease. Proper statistical evaluation of determined associations between imaging-based phenotypic and transcriptomic data is key in these explorations, in particular to establish whether observed associations exceed "chance level" of random, nonspecific effects. Recent approaches have shown the importance of statistical models that can correct for spatial autocorrelation effects in the data to avoid inflation of reported statistics. Here, we discuss the need for examination of a second category of statistical models in transcriptomic-neuroimaging analyses, namely those that can provide "gene specificity." By means of a couple of simple examples of commonly performed transcriptomic-neuroimaging analyses, we illustrate some of the potentials and challenges of transcriptomic-imaging analyses, showing that providing gene specificity on observed transcriptomic-neuroimaging effects is of high importance to avoid reports of nonspecific effects. Through means of simulations we show that the rate of reported nonspecific effects (i.e., effects that cannot be specifically linked to a specific gene or gene-set) can run as high as 60%, with only less than 5% of transcriptomic-neuroimaging associations observed through ordinary linear regression analyses showing both spatial and gene specificity. We provide a discussion, a tutorial, and an easy-to-use toolbox for the different options of null models in transcriptomic-neuroimaging analyses.

Correction: Genome-wide analysis reveals extensive genetic overlap between schizophrenia, bipolar disorder, and intelligence.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Genome-wide analysis reveals extensive genetic overlap between schizophrenia, bipolar disorder, and intelligence.

Schizophrenia (SCZ) and bipolar disorder (BD) are severe mental disorders associated with cognitive impairment, which is considered a major determinant of functional outcome. Despite this, the etiology of the cognitive impairment is poorly understood, and no satisfactory cognitive treatments exist. Increasing evidence indicates that genetic risk for SCZ may contribute to cognitive impairment, whereas the genetic relationship between BD and cognitive function remains unclear. Here, we combined large genome-wide association study data on SCZ (n = 82,315), BD (n = 51,710), and general intelligence (n = 269,867) to investigate overlap in common genetic variants using conditional false discovery rate (condFDR) analysis. We observed substantial genetic enrichment in both SCZ and BD conditional on associations with intelligence indicating polygenic overlap. Using condFDR analysis, we leveraged this enrichment to increase statistical power and identified 75 distinct genomic loci associated with both SCZ and intelligence, and 12 loci associated with both BD and intelligence at conjunctional FDR < 0.01. Among these loci, 20 are novel for SCZ, and four are novel for BD. Most SCZ risk alleles (61 of 75, 81%) were associated with poorer cognitive performance, whereas most BD risk alleles (9 of 12, 75%) were associated with better cognitive performance. A gene set analysis of the loci shared between SCZ and intelligence implicated biological processes related to neurodevelopment, synaptic integrity, and neurotransmission; the same analysis for BD was underpowered. Altogether, the study demonstrates that both SCZ and BD share genetic influences with intelligence, albeit in a different manner, providing new insights into their genetic architectures.

MIR137 schizophrenia-associated locus controls synaptic function by regulating synaptogenesis, synapse maturation and synaptic transmission.

The MIR137 locus is a replicated genetic risk factor for schizophrenia. The risk-associated allele is reported to increase miR-137 expression and miR-137 overexpression alters synaptic transmission in mouse hippocampus. We investigated the cellular mechanisms underlying these observed effects in mouse hippocampal neurons in culture. First, we correlated the risk allele to expression of the genes in the MIR137 locus in human postmortem brain. Some evidence for increased MIR137HG expression was observed, especially in hippocampus of the disease-associated genotype. Second, in mouse hippocampal neurons, we confirmed previously observed changes in synaptic transmission upon miR-137 overexpression. Evoked synaptic transmission and spontaneous release were 50% reduced. We identified defects in release probability as the underlying cause. In contrast to previous observations, no evidence was obtained for selective synaptic vesicle docking defects. Instead, ultrastructural morphometry revealed multiple effects of miR-137 overexpression on docking, active zone length and total vesicle number. Moreover, proteomic analyses of neuronal protein showed that expression of Syt1 and Cplx1, previously reported as downregulated upon miR-137 overexpression, was unaltered. Immunocytochemistry of synapses overexpressing miR-137 showed normal Synaptotagmin1 and Complexin1 protein levels. Instead, our proteomic analyses revealed altered expression of genes involved in synaptogenesis. Concomitantly, synaptogenesis assays revealed 31% reduction in synapse formation. Taken together, these data show that miR-137 regulates synaptic function by regulating synaptogenesis, synaptic ultrastructure and synapse function. These effects are plausible contributors to the increased schizophrenia risk associated with miR-137 overexpression.

Astrocyte Subtype Vulnerability in Stem Cell Models of Vanishing White Matter.

OBJECTIVE: Astrocytes have gained attention as important players in neurological disease. In line with their heterogeneous character, defects in specific astrocyte subtypes have been identified. Leukodystrophy vanishing white matter (VWM) shows selective vulnerability in white matter astrocytes, but the underlying mechanisms remain unclear. Induced pluripotent stem cell technology is being extensively explored in studies of pathophysiology and regenerative medicine. However, models for distinct astrocyte subtypes for VWM are lacking, thereby hampering identification of disease-specific pathways. METHODS: Here, we characterize human and mouse pluripotent stem cell-derived gray and white matter astrocyte subtypes to generate an in vitro VWM model. We examined morphology and functionality, and used coculture methods, high-content microscopy, and RNA sequencing to study VWM cultures. RESULTS: We found intrinsic vulnerability in specific astrocyte subpopulations in VWM. When comparing VWM and control cultures, white matter-like astrocytes inhibited oligodendrocyte maturation, and showed affected pathways in both human and mouse cultures, involving the immune system and extracellular matrix. Interestingly, human white matter-like astrocytes presented additional, human-specific disease mechanisms, such as neuronal and mitochondrial functioning. INTERPRETATION: Astrocyte subtype cultures revealed disease-specific pathways in VWM. Cross-validation of human- and mouse-derived protocols identified human-specific disease aspects. This study provides new insights into VWM disease mechanisms, which helps the development of in vivo regenerative applications, and we further present strategies to study astrocyte subtype vulnerability in neurological disease. ANN NEUROL 2019;86:780-792.

Chrysin reduces the activity and protein level of mature forms of sterol regulatory element-binding proteins.

Sterol regulatory element-binding proteins (SREBPs) are transcription factors that regulate the expression of genes involved in fatty acid and cholesterol biosynthetic pathways. The present study showed that the flavonoid chrysin impairs the fatty acid synthase promoter. Chrysin reduces the expression of SREBP target genes, such as fatty acid synthase, in human hepatoma Huh-7 cells and impairs de novo synthesis of fatty acids and cholesterol. Moreover, it reduces the endogenous mature, transcriptionally active forms of SREBPs, which are generated by the proteolytic processing of precursor forms. In addition, chrysin reduces the enforced expressing mature forms of SREBPs and their transcriptional activity. The ubiquitin-proteasome system is not involved in the chrysin-mediated reduction of SREBPs mature forms. These results suggest that chrysin suppresses SREBP activity, at least partially, via the degradation of SREBPs mature forms. Abbreviations: ACC1: acetyl-CoA carboxylase 1; DMEM: Dulbecco's modified Eagle's medium; FAS: fatty acid synthase; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; 25-HC: 25-hydroxycholesterol; HMGCS: HMG-CoA synthase; LDH: lactate dehydrogenase; LPDS: lipoprotein-deficient serum; PI3K: phosphatidylinositol 3-kinase; SCD1: stearoyl-CoA desaturase; SREBPs: sterol regulatory element-binding proteins.

Development of the VIGS System in the Dioecious Plant Silene latifolia.

(1) Background: Silene latifolia is a dioecious plant, whose sex is determined by XY-type sex chromosomes. Microbotryum lychnidis-dioicae is a smut fungus that infects S. latifolia plants and causes masculinization in female flowers, as if Microbotryum were acting as a sex-determining gene. Recent large-scale sequencing efforts have promised to provide candidate genes that are involved in the sex determination machinery in plants. These candidate genes are to be analyzed for functional characterization. A virus vector can be a tool for functional gene analyses; (2) Methods: To develop a viral vector system in S. latifolia plants, we selected Apple latent spherical virus (ALSV) as an appropriate virus vector that has a wide host range; (3) Results: Following the optimization of the ALSV inoculation method, S. latifolia plants were infected with ALSV at high rates in the upper leaves. In situ hybridization analysis revealed that ALSV can migrate into the flower meristems in S. latifolia plants. Successful VIGS (virus-induced gene silencing) in S. latifolia plants was demonstrated with knockdown of the phytoene desaturase gene. Finally, the developed method was applied to floral organ genes to evaluate its usability in flowers; (4) Conclusion: The developed system enables functional gene analyses in S. latifolia plants, which can unveil gene functions and networks of S. latifolia plants, such as the mechanisms of sex determination and fungal-induced masculinization.

Rapid sex identification method of spinach (Spinacia oleracea L.) in the vegetative stage using loop-mediated isothermal amplification.

MAIN CONCLUSION: A LAMP-mediated, simple and rapid method for sex identification in spinach was developed. Nutrient compositional analysis showed a higher iron content in male than female plants. Spinach (Spinacia oleracea L.) is a dioecious plant with its sex determined by the XY system. Male and female floral organs differ morphologically, but plants do not differ in the vegetative stage before flowering. PCR with Y chromosome markers has been used to determine the sex of dioecious plants before flowering. In this study, we developed a genotype-specific loop-mediated isothermal amplification (LAMP) for sex identification of individual vegetative-stage spinach plants, using primers designed for the genomic region flanked by male-specific markers. LAMP could specifically detect spinach males. The method was further modified to omit DNA purification and use just an aliquot of crude leaf extract homogenized in water. We compared the nutrient composition of males and females, finding higher amounts of iron in the males. Our method could therefore be used for rapidly discriminating male plants in the field, which is useful for efficient hybrid breeding.

Tonsillectomy reduces recurrence of IgA nephropathy in mesangial hypercellularity type categorized by the Oxford classification.

BACKGROUND: In patients with IgA nephropathy (IgAN), recurrence after steroid pulse therapy is associated with reduced renal survival. However, the predictors of recurrence have not yet been clarified. METHODS: All patients who received 6-month steroid pulse therapy from 2004 to 2010 in our four affiliated hospitals and achieved a reduction of proteinuria to <0.4 g/day 1 year after treatment were retrospectively evaluated. The primary outcome was proteinuria ≥1.0 g/day during follow-up or additional antiproteinuric therapy. Two histological classifications were evaluated, the Oxford Classification with a split system and Japanese histological grades (HGs) with a lumped system. RESULTS: During a median follow-up of 3.4 years, 27 (26.7 %) of the 101 patients showed recurrence. Multivariate analysis showed that HG was the only significant predictor of recurrence, with HG 2+3+4 vs HG 1 having a hazard ratio of 7.38 (95 % confidence interval 1.52-133). Furthermore, in patients with mesangial hypercellularity according to the Oxford Classification, cumulative rate of recurrence-free survival was greater in patients with steroid therapy plus tonsillectomy compared with those who received steroid therapy alone (Log-rank test, P = 0.022). However, this association was not observed in patients without mesangial hypercellularity. CONCLUSIONS: HG is a novel predictor of recurrence after steroid pulse therapy in patients with IgAN. Moreover, the combination of steroid pulse therapy plus tonsillectomy may indicate a lower risk of recurrence in patients with mesangial hypercellularity, as defined by the Oxford Classification.

DPP-4 inhibitor and alpha-glucosidase inhibitor equally improve endothelial function in patients with type 2 diabetes: EDGE study.

BACKGROUND: Alpha glucosidase inhibitor (GI) attenuates postprandial hyperglycemia (PPH) and reduces the risk of cardiovascular events in patients with impaired glucose tolerance or type 2 diabetes. Dipeptidyl peptidase 4 (DPP-4) inhibitors also attenuate PPH. PPH is one of the factors leading to endothelial dysfunction which is an early event in the pathogenesis of atherosclerosis. Furthermore, DPP-4 inhibitors protect endothelial function through a GLP-1-dependent mechanism. However, the impact of these two types of drugs on endothelial dysfunction in patients with type 2 diabetes has not been fully elucidated. We compared the effects of sitagliptin, a DPP-4 inhibitor, and voglibose, an alpha GI, on endothelial function in patients with diabetes. METHODS: We conducted a randomized prospective multicenter study in 66 patients with type 2 diabetes who did not achieve the treatment goal with sulfonylurea, metformin or pioglitazone treatment; 31 patients received sitagliptin treatment and 35 patients, voglibose treatment. The flow-mediated dilatation (FMD) of the brachial artery was measured in the fasting state at baseline and after 12 weeks of treatment. The primary endpoint was a change in FMD (ΔFMD) from the baseline to the end of follow-up. The effects of sitagliptin and voglibose on FMD were assessed by ANCOVA after adjustment for the baseline FMD, age, sex, current smoking, diabetes duration and body mass index. Secondary efficacy measures included changes in HbA1c, GIP, GLP-1, C-peptide, CD34, lipid profile, oxidative stress markers, inflammatory markers and eGFR and any adverse events. RESULTS: ΔFMD was significantly improved after 12 weeks of treatment in both groups, and there was no significant difference in ΔFMD between the two groups. There were no significant differences in changes in HbA1c, GIP, GLP-1, C-peptide, lipid profile, oxidative stress marker, inflammatory marker and eGFR between the two groups. Compared with voglibose, sitagliptin significantly increased the circulating CD34, a marker of endothelial progenitor cells. Adverse events were observed in 5 patients in only the voglibose group (diarrhea 1, nausea 1, edema 2 and abdominal fullness 1). CONCLUSIONS: Sitagliptin improved endothelial dysfunction just as well as voglibose in patients with type 2 diabetes. Sitagliptin had protective effects on endothelial function without adverse events. TRIAL REGISTRATION: registered at http://www.umin.ac.jp/ctrj/ under UMIN000003951.