RefEx, a reference gene expression dataset as a web tool for the functional analysis of genes

Publication, funding, and experimental data in support of Human Reference Atlas construction and usage

Experts from 18 consortia are collaborating on the Human Reference Atlas (HRA) which aims to map the 37 trillion cells in the healthy human body. Information relevant for HRA construction and usage is held by experts, published in scholarly papers, and captured in experimental data. However, these data sources use different metadata schemas and cannot be cross-searched efficiently. This paper documents the compilation of a dataset, named HRAlit, that links the 136 HRA v1.4 digital objects (31 organs with 4,279 anatomical structures, 1,210 cell types, 2,089 biomarkers) to 583,117 experts; 7,103,180 publications; 896,680 funded projects, and 1,816 experimental datasets. The resulting HRAlit has 22 tables with 20,939,937 records including 6 junction tables with 13,170,651 relationships. The HRAlit can be mined to identify leading experts, major papers, funding trends, or alignment with existing ontologies in support of systematic HRA construction and usage.

Nonmuscular Troponin-I is required for gastrulation in sea urchin embryos

BACKGROUND

Gastrulation is one of the most important events in our lives (Barresi and Gilbert, 2020, Developmental Biology, 12th ed.). The molecular mechanisms of gastrulation in multicellular organisms are not yet fully understood, since many molecular, physical, and chemical factors are involved in the event.

RESULTS

Here, we found that one of muscle components, Troponin-I (TnI), is expressed in future gut cells, which are not muscular cells at all, and regulates gastrulation in embryos of a sea urchin, Hemicentrotus pulcherrimus. When we block the function of TnI, the invagination was inhibited in spite that the gut-cell specifier gene is normally expressed. In addition, blocking myosin activity also induced incomplete gastrulation.

CONCLUSION

These results strongly suggested that TnI regulates nonmuscular actin-myosin interactions during sea urchin gastrulation. So far, Troponin system is treated as specific only for muscle components, especially for striated muscle, but our data clearly show that TnI is involved in nonmuscular event. It is also reported that recent sensitive gene expression analysis revealed that Troponin genes are expressed in nonmuscular tissues in mammals (Ono et al., Sci Data, 2017;4:170105). These evidences propose the new evolutionary and functional scenario of the involvement of Troponin system in nonmuscular cell behaviors using actin-myosin system in bilaterians including human being.

Inhibition of 7-dehydrocholesterol reductase prevents hepatic ferroptosis under an active state of sterol synthesis

Recent evidence indicates ferroptosis is implicated in the pathophysiology of various liver diseases; however, the organ-specific regulation mechanism is poorly understood. Here, we demonstrate 7-dehydrocholesterol reductase (DHCR7), the terminal enzyme of cholesterol biosynthesis, as a regulator of ferroptosis in hepatocytes. Genetic and pharmacological inhibition (with AY9944) of DHCR7 suppress ferroptosis in human hepatocellular carcinoma Huh-7 cells. DHCR7 inhibition increases its substrate, 7-dehydrocholesterol (7-DHC). Furthermore, exogenous 7-DHC supplementation using hydroxypropyl β-cyclodextrin suppresses ferroptosis. A 7-DHC-derived oxysterol metabolite, 3β,5α-dihydroxycholest-7-en-6-one (DHCEO), is increased by the ferroptosis-inducer RSL-3 in DHCR7-deficient cells, suggesting that the ferroptosis-suppressive effect of DHCR7 inhibition is associated with the oxidation of 7-DHC. Electron spin resonance analysis reveals that 7-DHC functions as a radical trapping agent, thus protecting cells from ferroptosis. We further show that AY9944 inhibits hepatic ischemia-reperfusion injury, and genetic ablation of Dhcr7 prevents acetaminophen-induced acute liver failure in mice. These findings provide new insights into the regulatory mechanism of liver ferroptosis and suggest a potential therapeutic option for ferroptosis-related liver diseases.

Upregulated Genes in Atrial Fibrillation Blood and the Left Atrium

INTRODUCTION

Atrial fibrillation (AF) is a common arrhythmia associated with aging. Many known risk factors are associated with AF, but many senior individuals do not develop AF despite having multiple risk factors. This finding suggests that other factors may be involved in AF onset. This study aimed to identify upregulated genes in the peripheral blood and left atrium of patients with AF. These genes may serve as potential biomarkers to predict AF onset risk and its complications.

METHODS

Gene expression data were analyzed from blood (n = 3) and left atrial samples (n = 15) of patients with AF and sinus rhythm. We evaluated the significant genes identified using p value analysis of weighted average difference to confirm their rankings. We created figures for the genes using GeneMANIA and performed a functional analysis using Cytoscape3.10.1. Hub and bottleneck genes were identified based on degree and betweenness centrality. We used reference expression (RefEx) to confirm the organs in which the extracted genes were expressed. Heatmaps and Gene ontology term evaluation were performed to further elucidate the biological functions of the genes.

RESULTS

We identified 12 upregulated genes (CAST, ASAH1, MAFB, VCAN, DDIT4, FTL, HEXB, PROS1, BNIP3L, PABPC1, YBX3, and S100A6) in both the blood and left atrium of patients with AF. We analyzed the gene functions using GeneMANIA and Cytoscape. The identified genes were involved in a variety of pathways, including lysosomal function and lipid and sphingolipid catabolism. Next, we investigated whether the 12 identified genes identified were systemically expressed or had high organ specificity. Finally, RefEx was used to analyze the gene expression levels in various tissues. Four genes, FTL, ASAH1, S100A6, and PABPC1, were highly expressed in the normal heart tissue. Finally, we evaluated the expression levels of the 12 genes in the blood of patients with AF using a heatmap. Our findings suggest that the 12 genes identified in this study, especially the lysosome-related genes (FTL and ASAH1), may be involved in AF pathogenesis.

CONCLUSION

Lysosome-related genes may be important to understand the AF pathophysiology and to develop AF-related future studies.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

Clomipramine inhibits dynamin GTPase activity by L-α-phosphatidyl-L-serine stimulation

Three dynamin isoforms play critical roles in clathrin-dependent endocytosis. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters host cells via clathrin-dependent endocytosis. We previously reported that 3-(3-chloro-10,11-dihydro-5H-dibenzo[b,f]azepin-5-yl)-N,N-dimethylpropan-1-amine (clomipramine) inhibits the GTPase activity of dynamin 1, which is in mainly neuron. Therefore, we investigated whether clomipramine inhibits the activity of other dynamin isoforms in this study. We found that, similar to its inhibitory effect on dynamin 1, clomipramine inhibited the l-α-phosphatidyl-l-serine-stimulated GTPase activity of dynamin 2, which is expressed ubiquitously, and dynamin 3, which is expressed in the lung. Inhibition of GTPase activity raises the possibility that clomipramine can suppress SARS-CoV-2 entry into host cells.

A bioinformatic investigation of proteasome and autophagy expression in the central nervous system

The ubiquitin proteasome system (UPS) and autophagy lysosome pathway (ALP) are crucial in the control of protein quality. However, data regarding the relative significance of UPS and ALP in the central nervous system (CNS) are limited. In the present study, using publicly available data, we computed the quantitative expression status of UPS- and ALP-related genes and their products in the CNS as compared with that in other tissues and cells. We obtained human and mouse gene expression datasets from the reference expression dataset (RefEx) and Genevestigator (a tool for handling curated transcriptomic data from public repositories) as well as human proteomics data from the proteomics database (ProteomicsDB). The expression levels of genes and proteins in four categories-ubiquitin, proteasome, autophagy, and lysosome--in the cells and tissues were assessed. Perturbation of the gene expression by drugs was also analyzed for the four categories. Compared with that for ubiquitin, autophagy, and lysosome, gene expression for proteasome was consistently lower in the CNS of mice but was more pronounced in humans. Neural stem cells and neurons showed low proteasome gene expression as compared with embryonic stem cells. Proteomic analyses, however, did not show trends similar to those observed in the gene expression analyses. Perturbation analyses revealed that azithromycin and vitamin D3 upregulated the expression of both UPS and ALP. Gene and proteomic expression data could offer a fresh perspective on CNS pathophysiology. Our results indicate that disproportional expression of UPS and ALP might affect CNS disorders and that this imbalance might be redressed by several therapeutic candidates.

Enhancing untargeted metabolomics using metadata-based source annotation

Human untargeted metabolomics studies annotate only ~10% of molecular features. We introduce reference-data-driven analysis to match metabolomics tandem mass spectrometry (MS/MS) data against metadata-annotated source data as a pseudo-MS/MS reference library. Applying this approach to food source data, we show that it increases MS/MS spectral usage 5.1-fold over conventional structural MS/MS library matches and allows empirical assessment of dietary patterns from untargeted data.

Optimized Mouse-on-mouse Immunohistochemical Detection of Mouse ESR2 Proteins with PPZ0506 Monoclonal Antibody

Despite the physiological significance of ESR2, a lack of well-validated detection systems for ESR2 proteins has hindered progress in ESR2 research. Thus, recent identification of a specific anti-human ESR2 monoclonal antibody (PPZ0506) and its specific cross-reactivity against mouse and rat ESR2 proteins heightened momenta toward development of appropriate immunohistochemical detection systems for rodent ESR2 proteins. Building upon our previous optimization of ESR2 immunohistochemical detection in rats using PPZ0506, in this study, we further aimed to optimize mouse-on-mouse immunohistochemical detection using PPZ0506. Our assessment of several staining conditions using paraffin-embedded ovary sections revealed that intense heat-induced antigen retrieval, appropriate blocking, and appropriate antibody dilutions were necessary for optimization of mouse-on-mouse immunohistochemistry. Subsequently, we applied the optimized immunostaining method to determine expression profiles of mouse ESR2 proteins in peripheral tissues and brain subregions. Our analyses revealed more localized distribution of mouse ESR2 proteins than previously assumed. Moreover, comparison of these results with those obtained in humans and rats using PPZ0506 revealed interspecies differences in ESR2 expression. We expect that our optimized methodology for immunohistochemical staining of mouse ESR2 proteins will help researchers to solve multiple lines of controversial evidence concerning ESR2 expression.

Oncogenic Role of ADAM32 in Hepatoblastoma: A Potential Molecular Target for Therapy

Outcomes of pediatric hepatoblastoma (HBL) have improved, but refractory cases still occur. More effective and safer drugs are needed that are based on molecular mechanisms. A disintegrin and metalloproteases (ADAMs) are expressed with high frequency in various human carcinomas and play an important role in cancer progression. In this study, we analyzed expression of ADAMs in HBL with a cDNA microarray dataset and found that the expression level of ADAM32 is particularly high. To investigate the role of ADAM32 in cancer, forced expression or knockdown experiments were conducted with HepG2 and HBL primary cells. Colony formation, cell migration and invasion, and cell viability were increased in HepG2 expressing ADAM32, whereas knockdown of ADAM32 induced a decrease in these cellular functions. Quantitative RT-PCR demonstrated an association between ADAM32 expression and the expression of genes related to cancer stem cells and epithelial–mesenchymal transition (EMT), suggesting a role of ADAM32 in cancer stemness and EMT. Furthermore, knockdown of ADAM32 increased cisplatin-induced apoptosis, and this effect was attenuated by a caspase-8 inhibitor, suggesting that ADAM32 plays a role in extrinsic apoptosis signaling. We conclude that ADAM32 plays a crucial role in progression of HBL, so it might be a promising molecular target in anticancer therapy.

New Drug Development and Clinical Trial Design by Applying Genomic Information Management

Depending on the patients’ genotype, the same drug may have different efficacies or side effects. With the cost of genomic analysis decreasing and reliability of analysis methods improving, vast amount of genomic information has been made available. Several studies in pharmacology have been based on genomic information to select the optimal drug, determine the dose, predict efficacy, and prevent side effects. This paper reviews the tissue specificity and genomic information of cancer. If the tissue specificity of cancer is low, cancer is induced in various organs based on a single gene mutation. Basket trials can be performed for carcinomas with low tissue specificity, confirming the efficacy of one drug for a single gene mutation in various carcinomas. Conversely, if the tissue specificity of cancer is high, cancer is induced in only one organ based on a single gene mutation. An umbrella trial can be performed for carcinomas with a high tissue specificity. Some drugs are effective for patients with a specific genotype. A companion diagnostic strategy that prescribes a specific drug for patients selected with a specific genotype is also reviewed. Genomic information is used in pharmacometrics to identify the relationship among pharmacokinetics, pharmacodynamics, and biomarkers of disease treatment effects. Utilizing genomic information, sophisticated clinical trials can be designed that will be better suited to the patients of specific genotypes. Genomic information also provides prospects for innovative drug development. Through proper genomic information management, factors relating to drug response and effects can be determined by selecting the appropriate data for analysis and by understanding the structure of the data. Selecting pre-processing and appropriate machine-learning libraries for use as machine-learning input features is also necessary. Professional curation of the output result is also required. Personalized medicine can be realized using a genome-based customized clinical trial design.

Smooth muscle protein 22α-Cre recombination in resting cardiac fibroblasts and hematopoietic precursors

The Cre-loxP system has been widely used for cell- or organ-specific gene manipulation, but it is important to precisely understand what kind of cells the recombination takes place in. Smooth muscle 22α (SM22α)-Cre mice have been utilized to alter genes in vascular smooth muscle cells (VSMCs), activated fibroblasts or cardiomyocytes (CMs). Moreover, previous reports indicated that SM22α-Cre is expressed in adipocytes, platelets or myeloid cells. However, there have been no report of whether SM22α-Cre recombination takes place in nonCMs in hearts. Thus, we used the double-fluorescent Cre reporter mouse in which GFP is expressed when recombination occurs. Immunofluorescence analysis demonstrated that recombination occurred in resting cardiac fibroblasts (CFs) or macrophages, as well as VSMCs and CMs. Flow cytometry showed that some CFs, resident macrophages, neutrophils, T cells, and B cells were positive for GFP. These results prompted us to analyze bone marrow cells, and we observed GFP-positive hematopoietic precursor cells (HPCs). Taken together, these results indicated that SM22α-Cre-mediated recombination occurs in resting CFs and hematopoietic cell lineages, including HPCs, which is a cautionary point when using SM22α-Cre mice.

Sulfation Pathways During Neurodevelopment

Sulfate is an important nutrient that modulates a diverse range of molecular and cellular functions in mammalian physiology. Over the past 2 decades, animal studies have linked numerous sulfate maintenance genes with neurological phenotypes, including seizures, impaired neurodevelopment, and behavioral abnormalities. Despite sulfation pathways being highly conserved between humans and animals, less than one third of all known sulfate maintenance genes are clinically reportable. In this review, we curated the temporal and spatial expression of 91 sulfate maintenance genes in human fetal brain from 4 to 17 weeks post conception using the online Human Developmental Biology Resource Expression. In addition, we performed a systematic search of PubMed and Embase, identifying those sulfate maintenance genes linked to atypical neurological phenotypes in humans and animals. Those findings, together with a search of the Online Mendelian Inheritance in Man database, identified a total of 18 candidate neurological dysfunction genes that are not yet considered in clinical settings. Collectively, this article provides an overview of sulfate biology genes to inform future investigations of perturbed sulfate homeostasis associated with neurological conditions.

Loss of KAP3 decreases intercellular adhesion and impairs intracellular transport of laminin in signet ring cell carcinoma of the stomach

Signet-ring cell carcinoma (SRCC) is a unique subtype of gastric cancer that is impaired for cell-cell adhesion. The pathogenesis of SRCC remains unclear. Here, we show that expression of kinesin-associated protein 3 (KAP3), a cargo adaptor subunit of the kinesin superfamily protein 3 (KIF3), a motor protein, is specifically decreased in SRCC of the stomach. CRISPR/Cas9-mediated gene knockout experiments indicated that loss of KAP3 impairs the formation of circumferential actomyosin cables by inactivating RhoA, leading to the weakening of cell-cell adhesion. Furthermore, in KAP3 knockout cells, post-Golgi transport of laminin, a key component of the basement membrane, was inhibited, resulting in impaired basement membrane formation. Together, these findings uncover a potential role for KAP3 in the pathogenesis of SRCC of the stomach.

Development of RIKEN Plant Metabolome MetaDatabase

The advancement of metabolomics in terms of techniques for measuring small molecules has enabled the rapid detection and quantification of numerous cellular metabolites. Metabolomic data provide new opportunities to gain a deeper understanding of plant metabolism that can improve the health of both plants and humans that consume them. Although major public repositories for general metabolomic data have been established, the community still has shortcomings related to data sharing, especially in terms of data reanalysis, reusability and reproducibility. To address these issues, we developed the RIKEN Plant Metabolome MetaDatabase (RIKEN PMM, http://metabobank.riken.jp/pmm/db/plantMetabolomics), which stores mass spectrometry-based (e.g. gas chromatography-MS-based) metabolite profiling data of plants together with their detailed, structured experimental metadata, including sampling and experimental procedures. Our metadata are described as Linked Open Data based on the Resource Description Framework using standardized and controlled vocabularies, such as the Metabolomics Standards Initiative Ontology, which are to be integrated with various life and biomedical science data using the World Wide Web. RIKEN PMM implements intuitive and interactive operations for plant metabolome data, including raw data (netCDF format), mass spectra (NIST MSP format) and metabolite annotations. The feature is suitable not only for biologists who are interested in metabolomic phenotypes, but also for researchers who would like to investigate life science in general through plant metabolomic approaches.

In silico Analysis of SARS-CoV-2 ORF8-Binding Proteins Reveals the Involvement of ORF8 in Acquired-Immune and Innate-Immune Systems

SARS-CoV-2 is the causative agent of a new type of coronavirus infection, COVID-19, which has rapidly spread worldwide. The overall genome sequence homology between SARS-CoV-2 and SARS-CoV is 79%. However, the homology of the ORF8 protein between these two coronaviruses is low, at ~26%. Previously, it has been suggested that infection by the ORF8-deleted variant of SARS-CoV-2 results in less severe symptoms than in the case of wild-type SARS-CoV-2. Although we found that ORF8 is involved in the proteasome autoimmunity system, the precise role of ORF8 in infection and pathology has not been fully clarified. In this study, we determined a new network of ORF8-interacting proteins by performing in silico analysis of the binding proteins against the previously described 47 ORF8-binding proteins. We used as a dataset 431 human protein candidates from Uniprot that physically interacted with 47 ORF8-binding proteins, as identified using STRING. Homology and phylogenetic profile analyses of the protein dataset were performed on 446 eukaryotic species whose genome sequences were available in KEGG OC. Based on the phylogenetic profile results, clustering analysis was performed using Ward's method. Our phylogenetic profiling showed that the interactors of the ORF8-interacting proteins were clustered into three classes that were conserved across chordates (Class 1: 152 proteins), metazoans (Class 2: 163 proteins), and eukaryotes (Class 3: 114 proteins). Following the KEGG pathway analysis, classification of cellular localization, tissue-specific expression analysis, and a literature study on each class of the phylogenetic profiling cluster tree, we predicted that the following: protein members in Class 1 could contribute to COVID-19 pathogenesis via complement and coagulation cascades and could promote sarcoidosis; the members of Class 1 and 2, together, may contribute to the downregulation of Interferon-β; and Class 3 proteins are associated with endoplasmic reticulum stress and the degradation of human leukocyte antigen.

Placental mammals acquired functional sequences in NRK for regulating the CK2-PTEN-AKT pathway and placental cell proliferation

The molecular evolution processes underlying the acquisition of the placenta in eutherian ancestors are not fully understood. Mouse NCK-interacting kinase (NIK)-related kinase (NRK) is expressed highly in the placenta and plays a role in preventing placental hyperplasia. Here, we show the molecular evolution of NRK, which confers its function for inhibiting placental cell proliferation. Comparative genome analysis identified NRK orthologs across vertebrates, which share the kinase and citron homology (CNH) domains. Evolutionary analysis revealed that NRK underwent extensive amino acid substitutions in the ancestor of placental mammals and has been since conserved. Biochemical analysis of mouse NRK revealed that the CNH domain binds to phospholipids, and a region in NRK binds to and inhibits casein kinase-2 (CK2), which we named the CK2-inhibitory region (CIR). Cell culture experiments suggest the following: 1) Mouse NRK is localized at the plasma membrane via the CNH domain, where the CIR inhibits CK2. 2) This mitigates CK2-dependent phosphorylation and inhibition of PTEN and 3) leads to the inhibition of AKT signaling and cell proliferation. Nrk deficiency increased phosphorylation levels of PTEN and AKT in mouse placenta, supporting our hypothesis. Unlike mouse NRK, chicken NRK did not bind to phospholipids and CK2, decrease phosphorylation of AKT, or inhibit cell proliferation. Both the CNH domain and CIR have evolved under purifying selection in placental mammals. Taken together, our study suggests that placental mammals acquired the phospholipid-binding CNH domain and CIR in NRK for regulating the CK2–PTEN–AKT pathway and placental cell proliferation.

New Therapeutic Targets for Hepatic Fibrosis in the Integrin Family, α8β1 and α11β1, Induced Specifically on Activated Stellate Cells

A huge effort has been devoted to developing drugs targeting integrins over 30 years, because of the primary roles of integrins in the cell-matrix milieu. Five αv-containing integrins, in the 24 family members, have been a central target of fibrosis. Currently, a small molecule against αvβ1 is undergoing a clinical trial for NASH-associated fibrosis as a rare agent aiming at fibrogenesis. Latent TGFβ activation, a distinct talent of αv-integrins, has been intriguing as a therapeutic target. None of the αv-integrin inhibitors, however, has been in the clinical market. αv-integrins commonly recognize an Arg-Gly-Asp (RGD) sequence, and thus the pharmacophore of inhibitors for the 5-integrins is based on the same RGD structure. The RGD preference of the integrins, at the same time, dilutes ligand specificity, as the 5-integrins share ligands containing RGD sequence such as fibronectin. With the inherent little specificity in both drugs and targets, "disease specificity" has become less important for the inhibitors than blocking as many αv-integrins. In fact, an almighty inhibitor for αv-integrins, pan-αv, was in a clinical trial. On the contrary, approved integrin inhibitors are all specific to target integrins, which are expressed in a cell-type specific manner: αIIbβ3 on platelets, α4β1, α4β7 and αLβ2 on leukocytes. Herein, "disease specific" integrins would serve as attractive targets. α8β1 and α11β1 are selectively expressed in hepatic stellate cells (HSCs) and distinctively induced upon culture activation. The exceptional specificity to activated HSCs reflects a rather "pathology specific" nature of these new integrins. The monoclonal antibodies against α8β1 and α11β1 in preclinical examinations may illuminate the road to the first medical agents.

Quantitative modeling of regular retinal microglia distribution

Microglia are resident immune cells in the central nervous system, showing a regular distribution. Advancing microscopy and image processing techniques have contributed to elucidating microglia’s morphology, dynamics, and distribution. However, the mechanism underlying the regular distribution of microglia remains to be elucidated. First, we quantitatively confirmed the regularity of the distribution pattern of microglial soma in the retina. Second, we formulated a mathematical model that includes factors that may influence regular distribution. Next, we experimentally quantified the model parameters (cell movement, process formation, and ATP dynamics). The resulting model simulation from the measured parameters showed that direct cell–cell contact is most important in generating regular cell spacing. Finally, we tried to specify the molecular pathway responsible for the repulsion between neighboring microglia.

rs6426881 in the 3'-UTR of PBX1 is involved in breast and gastric cancers via altering the binding potential of miR-522-3p

BACKGROUND

Breast and gastric cancers are the most important diseases that lead to cancer death and social healthcare challenge. Overexpression of PBX1, a proto-oncogene, is correlated with the progression and metastasis of various cancers. For the first time, in this study the researchers evaluated the relationship between rs6426881, affecting miR-522-3p binding to the PBX1, with breast and gastric cancers.

METHODS AND RESULTS

The Microarray analysis was performed for finding the relative expression level of PBX1 and hsa-miR-522-3p, based on high throughput experiments. The GSE54397, GSE112369, GSE10810, GSE241585.ER, GSE24185.PR, GSE68373, and GSE38167 datasets were analyzed. A case-control study was carried out in 123 Iranian suffering from breast cancer and 132 participants as control samples as well as 130 people suffering from gastric cancer and 54 people as control group members. SNP rs6426881 in the 3'-UTR of PBX1 was genotyped by the High-Resolution Melting (HRM) method. Association analysis revealed that rs6426881 is correlated with Estrogen and Progesterone receptors, grade, and stage of breast cancer. Furthermore, a significant relationship was observed between the genotypes and blood groups in gastric cancer, while the distribution of alleles was significantly related to smoking, status of the primary tumor, and metastasis (Chi-Square P < 0.05). Finally, Bioinformatics analyses suggested that rs6426881 contains binding sites for miR-522-3p in the 3'-UTR of PBX1 transcript. The finding suggested that TT genotype is associated with poor prognosis in breast and gastric cancer.

CONCLUSIONS

The rs6426881 T allele at PBX1 3'-UT is significantly related to breast and gastric cancers by altering the regulatory affinity of miR-522-3p to PBX1 3'-UTR and may be suggested as a novel prognostic biomarker for the diseases.

Nanosecond Pulsed Electric Field Only Transiently Affects the Cellular and Molecular Processes of Leydig Cells

The purpose of this study was to verify whether the nanosecond pulsed electric field, not eliciting thermal effects, permanently changes the molecular processes and gene expression of Leydig TM3 cells. The cells were exposed to a moderate electric field (80 quasi-rectangular shape pulses, 60 ns pulse width, and an electric field of 14 kV/cm). The putative disturbances were recorded over 24 h. After exposure to the nanosecond pulsed electric field, a 19% increase in cell diameter, a loss of microvilli, and a 70% reduction in cell adhesion were observed. Some cells showed the nonapoptotic externalization of phosphatidylserine through the pores in the plasma membrane. The cell proportion in the subG1 phase increased by 8% at the expense of the S and G2/M phases, and the DNA was fragmented in a small proportion of the cells. The membrane mitochondrial potential and superoxide content decreased by 37% and 23%, respectively. Microarray’s transcriptome analysis demonstrated a negative transient effect on the expression of genes involved in oxidative phosphorylation, DNA repair, cell proliferation, and the overexpression of plasma membrane proteins. We conclude that nanosecond pulsed electric field affected the physiology and gene expression of TM3 cells transiently, with a noticeable heterogeneity of cellular responses.

Multi-Omic Meta-Analysis of Transcriptomes and the Bibliome Uncovers Novel Hypoxia-Inducible Genes

Hypoxia is a condition in which cells, tissues, or organisms are deprived of sufficient oxygen supply. Aerobic organisms have a hypoxic response system, represented by hypoxia-inducible factor 1-α (HIF1A), to adapt to this condition. Due to publication bias, there has been little focus on genes other than well-known signature hypoxia-inducible genes. Therefore, in this study, we performed a meta-analysis to identify novel hypoxia-inducible genes. We searched publicly available transcriptome databases to obtain hypoxia-related experimental data, retrieved the metadata, and manually curated it. We selected the genes that are differentially expressed by hypoxic stimulation, and evaluated their relevance in hypoxia by performing enrichment analyses. Next, we performed a bibliometric analysis using gene2pubmed data to examine genes that have not been well studied in relation to hypoxia. Gene2pubmed data provides information about the relationship between genes and publications. We calculated and evaluated the number of reports and similarity coefficients of each gene to HIF1A, which is a representative gene in hypoxia studies. In this data-driven study, we report that several genes that were not known to be associated with hypoxia, including the G protein-coupled receptor 146 gene, are upregulated by hypoxic stimulation.

The mitochondrial Ca2+ uptake regulator, MICU1, is involved in cold stress-induced ferroptosis

Ferroptosis has recently attracted much interest because of its relevance to human diseases such as cancer and ischemia-reperfusion injury. We have reported that prolonged severe cold stress induces lipid peroxidation-dependent ferroptosis, but the upstream mechanism remains unknown. Here, using genome-wide CRISPR screening, we found that a mitochondrial Ca2+ uptake regulator, mitochondrial calcium uptake 1 (MICU1), is required for generating lipid peroxide and subsequent ferroptosis under cold stress. Furthermore, the gatekeeping activity of MICU1 through mitochondrial calcium uniporter (MCU) is suggested to be indispensable for cold stress-induced ferroptosis. MICU1 is required for mitochondrial Ca2+ increase, hyperpolarization of the mitochondrial membrane potential (MMP), and subsequent lipid peroxidation under cold stress. Collectively, these findings suggest that the MICU1-dependent mitochondrial Ca2+ homeostasis-MMP hyperpolarization axis is involved in cold stress-induced lipid peroxidation and ferroptosis.

Iduronate-2-sulfatase fused with anti-hTfR antibody, pabinafusp alfa, for MPS-II: A phase 2 trial in Brazil

In Hunter syndrome (mucopolysaccharidosis II [MPS-II]), systemic accumulation of glycosaminoglycans (GAGs) due to a deficiency of iduronate-2-sulfatase (IDS), caused by mutations in the IDS gene, leads to multiple somatic manifestations and in patients with the severe (neuronopathic) phenotype, also to central nervous system (CNS) involvement. These symptoms cannot be effectively treated with current enzyme-replacement therapies, as they are unable to cross the blood-brain barrier (BBB). Pabinafusp alfa, a novel IDS fused with an anti-human transferrin receptor antibody, was shown to penetrate the BBB and to address neurodegeneration in preclinical studies. Subsequent phase 1/2 and 2/3 clinical studies in Japan have shown marked reduction of GAG accumulation in the cerebrospinal fluid (CSF), along with favorable clinical responses. A 26-week, open-label, randomized, parallel-group phase 2 study was conducted in Brazil to further evaluate the safety and efficacy of intravenously administered pabinafusp alfa at 1.0, 2.0, and 4.0 mg/kg/week in MPS-II patients. The safety profiles in the three dosage groups were similar. Neurodevelopmental evaluation suggested positive neurocognitive signals despite a relatively short study period. The 2.0-mg/kg group, which demonstrated marked reductions in substrate concentrations in the CSF, serum, and urine, was considered to provide the best combination regarding safety and efficacy signals.

Loss of Ftsj1 perturbs codon-specific translation efficiency in the brain and is associated with X-linked intellectual disability

FtsJ RNA 2'-O-methyltransferase 1 (FTSJ1) gene has been implicated in X-linked intellectual disability (XLID), but the molecular pathogenesis is unknown. We show that Ftsj1 is responsible for 2'-O-methylation of 11 species of cytosolic transfer RNAs (tRNAs) at the anticodon region, and these modifications are abolished in Ftsj1 knockout (KO) mice and XLID patient-derived cells. Loss of 2'-O-methylation in Ftsj1 KO mouse selectively reduced the steady-state level of tRNA^Phe in the brain, resulting in a slow decoding at Phe codons. Ribosome profiling showed that translation efficiency is significantly reduced in a subset of genes that need to be efficiently translated to support synaptic organization and functions. Ftsj1 KO mice display immature synaptic morphology and aberrant synaptic plasticity, which are associated with anxiety-like and memory deficits. The data illuminate a fundamental role of tRNA modification in the brain through regulation of translation efficiency and provide mechanistic insights into FTSJ1-related XLID.

Induced hepatic stellate cell integrin, α8β1, enhances cellular contractility and TGFβ activity in liver fibrosis

No effective therapy exists for fatal fibrosis. New therapeutic targets are needed for hepatic fibrosis because the incidence keeps increasing. The activation and differentiation of fibroblasts into myofibroblasts that causes excessive matrix deposition is central to fibrosis. Here, we investigated whether (and which) integrin receptors for matrix proteins activate hepatic stellate cells (HSCs). First, integrin α‐subunits were investigated systematically for their expression over the course of HSC activation and their distribution on fibroblasts and other systemic primary cells. The most upregulated in plate culture‐activated HSCs and specifically expressed across fibroblast linages was the α8 subunit. An anti‐α8 neutralizing mAb was evaluated in three different murine fibrosis models: for cytotoxic (CCl₄ treatment), non‐alcoholic steatohepatitis‐associated and cholestatic fibrosis. In all models, pathology and fibrosis markers (hydroxyproline and α‐smooth muscle actin) were improved following the mAb injection. We also CCl₄‐treated mice with inducible Itga8−/−; these mice were protected from increased hydroxyproline levels. Furthermore, ITGA8 was upregulated in specimens from 90 patients with liver fibrosis, indicating the relevance of our findings to liver fibrosis in people. Mechanistically, inhibition or ligand engagement of HSC α8 suppressed and enhanced myofibroblast differentiation, respectively, and HSC/fibroblast α8 activated latent TGFβ. Finally, integrin α8β1 potentially fulfils the growing need for anti‐fibrotic drugs and is an integrin not to be ignored. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

FANTOM enters 20th year: expansion of transcriptomic atlases and functional annotation of non-coding RNAs

The Functional ANnoTation Of the Mammalian genome (FANTOM) Consortium has continued to provide extensive resources in the pursuit of understanding the transcriptome, and transcriptional regulation, of mammalian genomes for the last 20 years. To share these resources with the research community, the FANTOM web-interfaces and databases are being regularly updated, enhanced and expanded with new data types. In recent years, the FANTOM Consortium's efforts have been mainly focused on creating new non-coding RNA datasets and resources. The existing FANTOM5 human and mouse miRNA atlas was supplemented with rat, dog, and chicken datasets. The sixth (latest) edition of the FANTOM project was launched to assess the function of human long non-coding RNAs (lncRNAs). From its creation until 2020, FANTOM6 has contributed to the research community a large dataset generated from the knock-down of 285 lncRNAs in human dermal fibroblasts; this is followed with extensive expression profiling and cellular phenotyping. Other updates to the FANTOM resource includes the reprocessing of the miRNA and promoter atlases of human, mouse and chicken with the latest reference genome assemblies. To facilitate the use and accessibility of all above resources we further enhanced FANTOM data viewers and web interfaces. The updated FANTOM web resource is publicly available at https://fantom.gsc.riken.jp/.

A microtubule-LUZP1 association around tight junction promotes epithelial cell apical constriction

Apical constriction is critical for epithelial morphogenesis, including neural tube formation. Vertebrate apical constriction is induced by di-phosphorylated myosin light chain (ppMLC)-driven contraction of actomyosin-based circumferential rings (CRs), also known as perijunctional actomyosin rings, around apical junctional complexes (AJCs), mainly consisting of tight junctions (TJs) and adherens junctions (AJs). Here, we revealed a ppMLC-triggered system at TJ-associated CRs for vertebrate apical constriction involving microtubules, LUZP1, and myosin phosphatase. We first identified LUZP1 via unbiased screening of microtubule-associated proteins in the AJC-enriched fraction. In cultured epithelial cells, LUZP1 was found localized at TJ-, but not at AJ-, associated CRs, and LUZP1 knockout resulted in apical constriction defects with a significant reduction in ppMLC levels within CRs. A series of assays revealed that ppMLC promotes the recruitment of LUZP1 to TJ-associated CRs, where LUZP1 spatiotemporally inhibits myosin phosphatase in a microtubule-facilitated manner. Our results uncovered a hitherto unknown microtubule-LUZP1 association at TJ-associated CRs that inhibits myosin phosphatase, contributing significantly to the understanding of vertebrate apical constriction.

Disruption in the balance between apolipoprotein A-I and mast cell chymase in chronic hypersensitivity pneumonitis

Background

Apolipoprotein A‐I (apoA‐I) has an antifibrotic effect in idiopathic pulmonary fibrosis. Although pulmonary fibrosis is associated with poor prognosis of patients with hypersensitivity pneumonitis (HP), little is known regarding the role of apoA‐I in the pathogenesis of HP.

Methods

Two‐dimensional electrophoresis, immunoblotting, and enzyme‐linked immunosorbent assays were performed for the identification and quantification of apoA‐I in bronchoalveolar lavage fluid (BALF) from patients with acute and chronic HP. To investigate the degradation of apoA‐I, apoA‐I was incubated with BALF. Moreover, the role of apoA‐I in TGF‐β1‐induced epithelial–mesenchymal transition of A549 cells was examined.

Results

The concentration of apoA‐I in the BALF was significantly lower in chronic HP (n = 56) compared with acute HP (n = 31). The expression level of apoA‐I was also low in the lung tissues of chronic HP. ApoA‐I was degraded by BALF from HP patients. The number of chymase‐positive mast cells in the alveolar parenchyma was inversely correlated with apoA‐I levels in the BALF of chronic HP patients. In vitro experiment using A549 cells, untreated apoA‐I inhibited TGF‐β1‐induced epithelial–mesenchymal transition, although this trend was not observed in the chymase‐treated apoA‐I.

Conclusions

A decrease of apoA‐I was associated with the pathogenesis of chronic HP in terms of pulmonary fibrosis and mast cell chymase attenuated the protective effect of apoA‐I against pulmonary fibrosis. Furthermore, apoA‐I could be a crucial molecule associated with lung fibrogenesis of HP.

Identification of GLUT12/SLC2A12 as a urate transporter that regulates the blood urate level in hyperuricemia model mice

Recent genome-wide association studies have revealed some genetic loci associated with serum uric acid levels and susceptibility to gout/hyperuricemia which contain potential candidates of physiologically important urate transporters. One of these novel loci is located upstream of SGK1 and SLC2A12, suggesting that variations in these genes increase the risks of hyperuricemia and gout. We herein focused on SLC2A12 encoding a transporter, GLUT12, the physiological function of which remains unclear. As GLUT12 belongs to the same protein family as a well-recognized urate transporter GLUT9, we hypothesized that GLUT12 mediates membrane transport of urate. Therefore, we conducted functional assays and analyzed Glut12 knockout hyperuricemia model mice, generated using the CRISPR-Cas9 system. Our results revealed that GLUT12 acts as a physiological urate transporter and its dysfunction elevates the blood urate concentration. This study provides insights into the deeper understanding of the urate regulatory system in the body, which is also important for pathophysiology of gout/hyperuricemia.

Homeobox A4 suppresses vascular remodeling by repressing YAP/TEAD transcriptional activity

The Hippo signaling pathway is involved in the pathophysiology of various cardiovascular diseases. Yes-associated protein (YAP) and transcriptional enhancer activator domain (TEAD) transcriptional factors, the main transcriptional complex of the Hippo pathway, were recently identified as modulators of phenotypic switching of vascular smooth muscle cells (VSMCs). However, the intrinsic regulator of YAP/TEAD-mediated gene expressions involved in vascular pathophysiology remains to be elucidated. Here, we identified Homeobox A4 (HOXA4) as a potent repressor of YAP/TEAD transcriptional activity using lentiviral shRNA screen. Mechanistically, HOXA4 interacts with TEADs and attenuates YAP/TEAD-mediated transcription by competing with YAP for TEAD binding. We also clarified that the expression of HOXA4 is relatively abundant in the vasculature, especially in VSMCs. In vitro experiments in human VSMCs showed HOXA4 maintains the differentiation state of VSMCs via inhibition of YAP/TEAD-induced phenotypic switching. We generated Hoxa4-deficient mice and confirmed the downregulation of smooth muscle-specific contractile genes and the exacerbation of vascular remodeling after carotid artery ligation in vivo. Our results demonstrate that HOXA4 is a repressor of VSMC phenotypic switching by inhibiting YAP/TEAD-mediated transcription.

All of gene expression (AOE): An integrated index for public gene expression databases

Gene expression data have been archived as microarray and RNA-seq datasets in two public databases, Gene Expression Omnibus (GEO) and ArrayExpress (AE). In 2018, the DNA DataBank of Japan started a similar repository called the Genomic Expression Archive (GEA). These databases are useful resources for the functional interpretation of genes, but have been separately maintained and may lack RNA-seq data, while the original sequence data are available in the Sequence Read Archive (SRA). We constructed an index for those gene expression data repositories, called All Of gene Expression (AOE), to integrate publicly available gene expression data. The web interface of AOE can graphically query data in addition to the application programming interface. By collecting gene expression data from RNA-seq in the SRA, AOE also includes data not included in GEO and AE. AOE is accessible as a search tool from the GEA website and is freely available at https://aoe.dbcls.jp/.

COL5A1 rs12722 polymorphism is not associated with passive muscle stiffness and sports-related muscle injury in Japanese athletes

BACKGROUND

Poor joint flexibility has been repeatedly proposed as a risk factor for muscle injury. The C-to-T polymorphism (rs12722) in the 3'-untranslated region of the collagen type V α1 chain gene (COL5A1) is reportedly associated with joint flexibility. Flexibility of a normal joint is largely determined by passive muscle stiffness, which is influenced by intramuscular collagenous connective tissues including type V collagen. The present study aimed to test the hypothesis that the COL5A1 rs12722 polymorphism influences joint flexibility via passive muscle stiffness, and is accordingly associated with the incidence of muscle injury.

METHODS

In Study 1, we examined whether the rs12722 polymorphism is associated with joint flexibility and passive muscle stiffness in 363 healthy young adults. Joint flexibility was evaluated by passive straight-leg-raise and sit-and-reach tests, and passive muscle stiffness was measured using ultrasound shear wave elastography. In Study 2, the association of the rs12722 polymorphism with sports-related muscle injury was assessed in 1559 Japanese athletes. Muscle injury history and severity were assessed by a questionnaire. In both Study 1 and Study 2, the rs12722 C-to-T polymorphism in the COL5A1 was determined using the TaqMan SNP Genotyping Assay.

RESULTS

Study 1 revealed that the rs12722 polymorphism had no significant effect on range of motion in passive straight-leg-raise and sit-and-reach tests. Furthermore, there was no significant difference in passive muscle stiffness of the hamstring among the rs12722 genotypes. In Study 2, rs12722 genotype frequencies did not differ between the muscle injury and no muscle injury groups. Moreover, no association was observed between rs12722 polymorphism and severity of muscle injury.

CONCLUSIONS

The present study does not support the view that COL5A1 rs12722 polymorphism has a role as a risk factor for sports-related muscle injury, or that it is a determinant for passive muscle stiffness in a Japanese population.

Construction of a simple evaluation system for the intestinal absorption of an orally administered medicine using Bombyx mori larvae

Human intestinal absorption is estimated using a human colon carcinoma cell line (Caco-2) cells from human colorectal adenocarcinoma, intestinal perfusion, or a mammalian model. These current evaluation systems are limited in their ability to estimate human intestinal absorption. In addition, in vivo evaluation systems using laboratory animals such as mice and rats entail animal ethics problems, and it is difficult to screen compounds on a large scale at the drug discovery stage. Thus, we propose the use of Bombyx mori larvae for evaluation of intestinal absorption of compounds as an alternative system in this study. First, to compare the characteristics among Caco-2 cells, human intestine, and B. mori larval midgut, we analyzed their RNA-seq data, and we found 26 drug transporters common to humans and B. mori. Next, we quantitatively developed an oral administration technique in B. mori and established a method using silkworm B. mori larvae that can easily estimate the intestinal permeability of compounds. Consequently, we could determine the dose and technique for oral administration in B. mori larvae. We also developed a B. mori model to evaluate the intestinal permeability of orally administered. Our constructed evaluation system will be useful for evaluating intestinal permeability in medical drug development.

Significant antitumor effect of an antibody against TMEM180, a new colorectal cancer-specific molecule

The present state of therapy for colorectal cancer (CRC) is far from satisfactory, highlighting the need for new targets for this disease. We identified a new CRC‐specific molecule, TMEM180, a predicted 11‐pass transmembrane protein that apparently functions as a cation symporter. We developed an anti‐TMEM180 mAb and then succeeded in humanizing the mAb. Immunohistochemistry (IHC) in CRC with the mAb showed a similar positivity rate as compared with anti‐epidermal growth factor receptor mAb, and IHC with anti‐TMEM180 mAb did not show staining in major organs used in this study. Immune electron microscopy clearly indicated that TMEM180 was present on the tumor exosome. The TMEM180 promoter region contains 10 hypoxia‐responsive element consensus sequences; accordingly, SW480 cells upregulated TMEM180 under low‐oxygen conditions. Anti‐TMEM180 mAb has in vitro antibody‐dependent cell‐mediated cytotoxicity and complement‐dependent cytotoxicity activity, and SW480 CRC xenografts were eradicated by the mAb. These data indicate that TMEM180 may be a new CRC marker and that a mAb against this protein could be used as antibody‐based therapy against CRC.

Drug discovery of anticancer drugs targeting methylenetetrahydrofolate dehydrogenase 2

Many anticancer drugs have serious adverse effects; therefore, it is necessary to target features specific to cancer cells to minimize the effects on healthy cells. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) was reported to be specifically enhanced in cancer. We confirmed the validity of MTHFD2 as a drug discovery target using clinical data. In addition, we performed in silico screening to design an anticancer drug specifically targeting MTHFD2. Analysis of the clinical data indicated that MTHFD2 was enhanced in most cancers compared with normal tissues, and affected the prognosis in cancer patients. Candidate compounds for MTHFD2 inhibitors were identified using in silico drug discovery techniques, and the important interactions for MTHFD2 binding were determined. In addition, these candidate compounds decreased levels of MTHFD2 metabolites in cancer cells. The findings of the present study may help to develop anticancer drugs targeting MTHFD2, with a view to minimizing the adverse effects of anticancer drugs.

Variation in the representation of human anatomy within digital resources: Implications for data integration

Clinical informatics makes use of anatomical representation-particularly in the form of anatomical terms. But differences and ambiguities in naming anatomical structures and partitioning the body can complicate efforts to interlink anatomical resources and integrate clinical data. To better understand differences in representations of human anatomy, we compare five digital resources: a formal ontology, a terminology, and three 3D graphics applications. Because the graphics applications offer explicit representation of the boundaries and partitions of anatomical structures, they reveal the differences in modeling of anatomy that may not be apparent through text-based representations. The variations in these resources allow us to categorize differences in representations of anatomy and to highlight the importance of this topic in the context of clinical informatics.

The Healing Effect of Human Milk Fat Globule-EGF Factor 8 Protein (MFG-E8) in A Rat Model of Parkinson's Disease

We searched for drugs that alleviate the reduction of dopaminergic neurons caused by the administration of lipopolysaccharide (LPS) to the substantia nigra of the rat brain. Human milk fat globule-EGF factor 8 protein (MFG-E8) is similar to MFG-E8-S, a short isoform, of the mouse MFG-E8. However, the function of MFG-E8-S was not clear. Rats with LPS-induced Parkinson’s disease were prepared and the effects of human MFG-E8 were examined. MFG-E8 improved the significant reduction in mesencephalic dopamine neurons induced by the administration of LPS. LPS was administered to human induced pluripotent stem cell (iPSC)-derived dopaminergic neurons to induce inflammation and the effect of MFG-E8 was examined. MFG-E8 showed no toxicity toward neurons. We reanalyzed the results using public microarray data. MFG-E8 mRNA was found to be expressed in all parts of the body, particularly by adipose-derived stem cells (ADSCs). Furthermore, we investigated the culture supernatant of ADSCs using the liquid chromatography-tandem mass spectrometry (LC–MS/MS) analysis method and successfully identified the peptide of the MFG-E8 F5/8 type C domain. The results suggested that MFG-E8-S may have a preventive effect against Parkinson’s disease.

Unbiased compound screening with a reporter gene assay highlights the role of p13 in the cardiac cellular stress response

We recently showed that a 13-kDa protein (p13), the homolog protein of formation of mitochondrial complex V assembly factor 1 in yeast, acts as a potential protective factor in pancreatic islets under diabetes. Here, we aimed to identify known compounds regulating p13 mRNA expression to obtain therapeutic insight into the cellular stress response. A luciferase reporter system was developed using the putative promoter region of the human p13 gene. Overexpression of peroxisome proliferator-activated receptor gamma coactivator 1α, a master player regulating mitochondrial metabolism, increased both reporter activity and p13 expression. Following unbiased screening with 2320 known compounds in HeLa cells, 12 pharmacological agents (including 8 cardiotonics and 2 anthracyclines) that elicited >2-fold changes in p13 mRNA expression were identified. Among them, four cardiac glycosides decreased p13 expression and concomitantly elevated cellular oxidative stress. Additional database analyses showed highest p13 expression in heart, with typically decreased expression in cardiac disease. Accordingly, our results illustrate the usefulness of unbiased compound screening as a method for identifying novel functional roles of unfamiliar genes. Our findings also highlight the importance of p13 in the cellular stress response in heart.

NBDC RDF portal: a comprehensive repository for semantic data in life sciences

In the life sciences, researchers increasingly want to access multiple databases in an integrated way. However, different databases currently use different formats and vocabularies, hindering the proper integration of heterogeneous life science data. Adopting the Resource Description Framework (RDF) has the potential to address such issues by improving database interoperability, leading to advances in automatic data processing. Based on this idea, we have advised many Japanese database development groups to expose their databases in RDF. To further promote such activities, we have developed an RDF-based life science dataset repository called the National Bioscience Database Center (NBDC) RDF portal. All the datasets in this repository have been reviewed by the NBDC to ensure interoperability and queryability. As of July 2018, the service includes 21 RDF datasets, comprising over 45.5 billion triples. It provides SPARQL endpoints for all datasets, useful metadata and the ability to download RDF files. The NBDC RDF portal can be accessed at https://integbio.jp/rdf/.