Synergistic role of retinoic acid signaling and Gata3 during primitive choanae formation.

Developmental defects of primitive choanae, an anatomical path to connect the embryonic nasal and oral cavity, result in disorders called choanal atresia (CA), which are associated with many congenital diseases and require immediate clinical intervention after birth. Previous studies revealed that reduced retinoid signaling underlies the etiology of CA. In the present study, by using multiple mouse models which conditionally deleted Rdh10 and Gata3 during embryogenesis, we showed that Gata3 expression is regulated by retinoid signaling during embryonic craniofacial development and plays crucial roles for development of the primitive choanae. Interestingly, Gata3 loss of function is known to cause hypoparathyroidism, sensorineural deafness and renal disease (HDR) syndrome, which exhibits CA as one of the phenotypes in humans. Our model partially phenocopies HDR syndrome with CA, and is thus a useful tool for investigating the molecular and cellular mechanisms of HDR syndrome. We further uncovered critical synergy of Gata3 and retinoid signaling during embryonic development, which will shed light on novel molecular and cellular etiology of congenital defects in primitive choanae formation.

Plant GARDEN: a portal website for cross-searching between different types of genomic and genetic resources in a wide variety of plant species.

BACKGROUND: Plant genome information is fundamental to plant research and development. Along with the increase in the number of published plant genomes, there is a need for an efficient system to retrieve various kinds of genome-related information from many plant species across plant kingdoms. Various plant databases have been developed, but no public database covers both genomic and genetic resources over a wide range of plant species. MAIN BODY: We have developed a plant genome portal site, Plant GARDEN (Genome And Resource Database Entry: https://plantgarden.jp/en/index ), to provide diverse information related to plant genomics and genetics in divergent plant species. Elasticsearch is used as a search engine, and cross-keyword search across species is available. Web-based user interfaces (WUI) for PCs and tablet computers were independently developed to make data searches more convenient. Several types of data are stored in Plant GARDEN: reference genomes, gene sequences, PCR-based DNA markers, trait-linked DNA markers identified in genetic studies, SNPs, and in/dels on publicly available sequence read archives (SRAs). The data registered in Plant GARDEN as of March 2023 included 304 assembled genome sequences, 11,331,614 gene sequences, 419,132 DNA markers, 8,225 QTLs, and 5,934 SNP lists (gvcf files). In addition, we have re-annotated all the genes registered in Plant GARDEN by using a functional annotation tool, Hayai-Annotation, to compare the orthologous relationships among genes. CONCLUSION: The aim of Plant GARDEN is to provide plant genome information for use in the fields of plant science as well as for plant-based industries, education, and other relevant areas. Therefore, we have designed a WUI that allows a diverse range of users to access such information in an easy-to-understand manner. Plant GARDEN will eventually include a wide range of plant species for which genome sequences are assembled, and thus the number of plant species in the database will continue to expand. We anticipate that Plant GARDEN will promote the understanding of genomes and gene diversity by facilitating comparisons of the registered sequences.

Disruption of a RAC1-centred network is associated with Alzheimer's disease pathology and causes age-dependent neurodegeneration.

The molecular biological mechanisms of Alzheimer's disease (AD) involve disease-associated crosstalk through many genes and include a loss of normal as well as a gain of abnormal interactions among genes. A protein domain network (PDN) is a collection of physical bindings that occur between protein domains, and the states of the PDNs in patients with AD are likely to be perturbed compared to those in normal healthy individuals. To identify PDN changes that cause neurodegeneration, we analysed the PDNs that occur among genes co-expressed in each of three brain regions at each stage of AD. Our analysis revealed that the PDNs collapsed with the progression of AD stage and identified five hub genes, including Rac1, as key players in PDN collapse. Using publicly available as well as our own gene expression data, we confirmed that the mRNA expression level of the RAC1 gene was downregulated in the entorhinal cortex (EC) of AD brains. To test the causality of these changes in neurodegeneration, we utilized Drosophila as a genetic model and found that modest knockdown of Rac1 in neurons was sufficient to cause age-dependent behavioural deficits and neurodegeneration. Finally, we identified a microRNA, hsa-miR-101-3p, as a potential regulator of RAC1 in AD brains. As the Braak neurofibrillary tangle (NFT) stage progressed, the expression levels of hsa-miR-101-3p were increased specifically in the EC. Furthermore, overexpression of hsa-miR-101-3p in the human neuronal cell line SH-SY5Y caused RAC1 downregulation. These results highlight the utility of our integrated network approach for identifying causal changes leading to neurodegeneration in AD.

Lipid Pathology of the Corpus Callosum in Schizophrenia and the Potential Role of Abnormal Gene Regulatory Networks with Reduced Microglial Marker Expression.

Structural changes in the corpus callosum have been reported in schizophrenia; however, the underlying molecular mechanism remains unclear. As the corpus callosum is high in lipid content, we analyzed the lipid contents of the corpora callosa from 15 patients with schizophrenia and 15 age- and sex-matched controls using liquid chromatography coupled to tandem mass spectrometry and identified lipid combinations associated with schizophrenia. Real-time quantitative polymerase chain reaction analyses using extended samples (schizophrenia, n = 95; control, n = 91) showed low expression levels of lipid metabolism-related genes and their potential upstream transcription factors in schizophrenia. Subsequent pathway analysis identified a gene regulatory network where nuclear factor of activated T cells 2 (NFATC2) is placed most upstream. We also observed low gene expression levels of microglial markers, inflammatory cytokines, and colony-stimulating factor 1 receptor (CSF1R), which is known to regulate the density of microglia, in the corpus callosum in schizophrenia. The interactions between CSF1R and several genes in the presently identified gene network originating from NFATC2 have been reported. Collectively, this study provides evidence regarding lipid abnormalities in the corpora callosa of patients with schizophrenia and proposes the potential role of impaired "NFATC2-relevant gene network-microglial axis" as its underlying mechanism.

Mutation in VPS33A affects metabolism of glycosaminoglycans: a new type of mucopolysaccharidosis with severe systemic symptoms.

Mucopolysaccharidoses (MPS) are a group of genetic deficiencies of lysosomal enzymes that catabolize glycosaminoglycans (GAG). Here we describe a novel MPS-like disease caused by a specific mutation in the VPS33A gene. We identified several Yakut patients showing typical manifestations of MPS: coarse facial features, skeletal abnormalities, hepatosplenomegaly, respiratory problems, mental retardation, and excess secretion of urinary GAG. However, these patients could not be diagnosed enzymatically as MPS. They showed extremely high levels of plasma heparan sulphate (HS, one of GAG); 60 times the normal reference range and 6 times that of MPS patients. Additionally, most patients developed heart, kidney, and hematopoietic disorders, which are not typical symptoms for conventional MPS, leading to a fatal outcome between 1 and 2-years old. Using whole exome and Sanger sequencing, we identified homozygous c.1492C > T (p.Arg498Trp) mutations in the VPS33A gene of 13 patients. VPS33A is involved in endocytic and autophagic pathways, but the identified mutation did not affect either of these pathways. Lysosomal over-acidification and HS accumulation were detected in patient-derived and VPS33A-depleted cells, suggesting a novel role of this gene in lysosomal functions. We hence propose a new type of MPS that is not caused by an enzymatic deficiency.

Identification of genetic risk factors of aggressive periodontitis using genomewide association studies in association with those of chronic periodontitis.

To identify the genetic risk factors for aggressive periodontitis (AgP), it is important to understand the progression and pathogenesis of AgP. The purpose of this review was to summarize the genetic risk factors for AgP identified through a case-control genomewide association study (GWAS) and replication study. The initial studies to identify novel AgP risk factors were potentially biased because they relied on previous studies. To overcome this kind of issue, an unbiased GWAS strategy was introduced to identify genetic risk factors for various diseases. Currently, three genes glycosyltransferase 6 domain containing 1 (GLT6D1), defensin α1 and α3 (DEFA1A3), and sialic acid-binding Ig-like lectin 5 (SIGLEC5) that reach the threshold for genomewide significance have been identified as genetic risk factors for AgP through a case-control GWAS.

Barth Syndrome: Different Approaches to Diagnosis.

The diagnosis of Barth syndrome is challenging owing to the wide phenotypic spectrum with allelic heterogeneity. Here we report 3 cases of Barth syndrome with phenotypic and allelic heterogeneity that were diagnosed by different approaches, including whole exome sequencing and final confirmation by reverse-transcription polymease chain reaction.

Identification of biomarker sets for predicting the efficacy of sublingual immunotherapy against pollen-induced allergic rhinitis.

Sublingual immunotherapy (SLIT) is effective against allergic rhinitis, although a substantial proportion of individuals is refractory. Herein, we describe a predictive modality to reliably identify SLIT non-responders (NRs). We conducted a 2-year clinical study in 193 adult patients with Japanese cedar pollinosis, with biweekly administration of 2000 Japanese allergy units of cedar pollen extract as the maintenance dose. After identifying high-responder (HR) patients with improved severity scores and NR patients with unchanged or exacerbated symptoms, differences in 33 HR and 34 NR patients were evaluated in terms of peripheral blood cellular profiles by flow cytometry and serum factors by ELISA and cytokine bead array, both pre- and post-SLIT. Improved clinical responses were seen in 72% of the treated patients. Pre-therapy IL-12p70 and post-therapy IgG1 serum levels were significantly different between HR and NR patients, although these parameters alone failed to distinguish NR from HR patients. However, the analysis of serum parameters in the pre-therapy samples with the Adaptive Boosting (AdaBoost) algorithm distinguished NR patients with high probability within the training data set. Cluster analysis revealed a positive correlation between serum Th1/Th2 cytokines and other cytokines/chemokines in HR patients after SLIT. Thus, processing of pre-therapy serum parameters with AdaBoost and cluster analysis can be reliably used to develop a prediction method for HR/NR patients.

HDR-del: A tool based on Hamming distance for prioritizing pathogenic chromosomal deletions in exome sequencing.

High-density oligonucleotide arrays have widely been used to detect pathogenic chromosomal deletions. In addition to high-density oligonucleotide arrays, programs using whole-exome sequencing have become available for estimating copy-number variations using depth of coverage. Here, we propose a new statistical method, HDR-del, to prioritize pathogenic chromosomal deletions based on Hamming distance in exome sequencing. In vcf (variant call format) files generated from exome sequencing, hemizygous chromosomal deletion regions lack heterozygous variants and lead to apparent long runs of homozygosity (ROH). In our Hamming distance ratio (HDR)-del approach, we calculate the "difference" in heterozygous status between an affected individual and control individuals using the HDR over all candidate chromosomal deletion regions defined as ROH longer than 1Mbp. Using a suitable test statistic, which is expected to be large for a true pathogenic deletion region, we prioritize candidate chromosomal deletion regions based on this statistic. In our approach, we were able to considerably narrow down true pathogenic chromosomal deletion regions, which were confirmed by high-density oligonucleotide arrays in four mitochondrial disease patients. Our HDR-del approach represents an easy method for detecting chromosomal deletions.

Analysis of an ADTKD family with a novel frameshift mutation in MUC1 reveals characteristic features of mutant MUC1 protein.

BACKGROUND: Medullary cystic kidney disease Type 1 is an autosomal dominant tubulointerstitial kidney disease (ADTKD). Recently, mucin 1 (MUC1) was identified as a causal gene of medullary cystic kidney disease (ADTKD-MUC1). However, the MUC1 mutation was found to be a single cytosine insertion in a single copy of the GC-rich variable number of tandem repeats (VNTRs), which are very difficult to analyze by next-generation sequencing. To date, other mutations have not been detected in ADTKD-MUC1, and the mutant MUC1 protein has not been analyzed because of the difficulty of genetically modifying the VNTR sequence. METHODS: We conducted whole-exome analyses of an ADTKD family by next-generation sequencing. We also performed histopathological analyses of a renal biopsy from a pedigree family member. We constructed a mutant protein expression vector based on the patient genome sequence and characterized the nature of the mutant protein. RESULTS: We found a novel frameshift mutation before the VNTR in the MUC1 gene. The resulting mutant MUC1 protein had a very similar amino acid sequence and predicted 3D structure to the previously reported mutant protein. Notably, the recombinant mutant MUC1 protein was trapped in the cytoplasm and appeared to self-aggregate. The patient native mutant protein was also found in urine exosomes. CONCLUSIONS: This novel frameshift mutation in the MUC1 gene and consequent mutant protein may contribute to the future discovery of the pathophysiology of ADTKD-MUC1. The mutant MUC1 protein in urine exosomes may be used for non-DNA-related diagnosis.

Identification of core gene networks and hub genes associated with progression of non-alcoholic fatty liver disease by RNA sequencing.

AIM: Non-alcoholic fatty liver disease (NAFLD) progresses because of the interaction between numerous genes. Thus, we carried out a weighted gene coexpression network analysis to identify core gene networks and key genes associated with NAFLD progression. METHODS: We enrolled 39 patients with mild NAFLD (fibrosis stages 0-2) and 21 with advanced NAFLD (fibrosis stages 3-4). Total RNA was extracted from frozen liver biopsies, and sequenced to capture a large dynamic range of expression levels. RESULTS: A total of 1777 genes differentially expressed between mild and advanced NAFLD (q-value <0.05) clustered into four modules. One module was enriched for genes that encode cell surface or extracellular matrix proteins, and are involved in cell adhesion, proliferation, and signaling. This module formed a scale-free network containing four hub genes (PAPLN, LBH, DPYSL3, and JAG1) overexpressed in advanced NAFLD. PAPLN is a component of the extracellular matrix, LBH and DPYSL3 are reported to be tumor suppressors, and JAG1 is tumorigenic. Another module formed a random network, and was enriched for genes that accumulate in the mitochondria. These genes were downregulated in advanced NAFLD, reflecting impaired mitochondrial function. However, the other two modules did not form unambiguous networks. KEGG analysis indicated that 71 differentially expressed genes were involved in "pathways in cancer". Strikingly, expression of half of all differentially expressed genes was inversely correlated with methylation of CpG sites (q-value <0.05). Among clinical parameters, serum type IV collagen 7 s was most strongly associated with the epigenetic status in NAFLD. CONCLUSIONS: Newly identified core gene networks suggest that the NAFLD liver undergoes mitochondrial dysfunction and fibrosis, and acquires tumorigenic potential epigenetically. Our data provide novel insights into the pathology and etiology of NAFLD progression, and identify potential targets for diagnosis and treatment.

Effect of Clozapine on DNA Methylation in Peripheral Leukocytes from Patients with Treatment-Resistant Schizophrenia.

Clozapine is an atypical antipsychotic, that is established as the treatment of choice for treatment-resistant schizophrenia (SCZ). To date, no study investigating comprehensive DNA methylation changes in SCZ patients treated with chronic clozapine has been reported. The purpose of the present study is to reveal the effects of clozapine on DNA methylation in treatment-resistant SCZ. We conducted a genome-wide DNA methylation profiling in peripheral leukocytes (485,764 CpG dinucleotides) from treatment-resistant SCZ patients treated with clozapine (n = 21) in a longitudinal study. Significant changes in DNA methylation were observed at 29,134 sites after one year of treatment with clozapine, and these genes were enriched for "cell substrate adhesion" and "cell matrix adhesion" gene ontology (GO) terms. Furthermore, DNA methylation changes in the CREBBP (CREB binding protein) gene were significantly correlated with the clinical improvements. Our findings provide insights into the action of clozapine in treatment-resistant SCZ.

HDR: a statistical two-step approach successfully identifies disease genes in autosomal recessive families.

In the search for sequence variants underlying disease, commonly applied filtering steps usually result in a number of candidate variants that cannot further be narrowed down. In autosomal recessive families, disease usually occurs only in one generation so that genetic linkage analysis is unlikely to help. Because homozygous recessive mutations tend to be inherited together with flanking homozygous variants, we developed a statistical method to detect pathogenic variants in autosomal recessive families: We look for differences in patterns of homozygosity around candidate variants between patients and control individuals and expect that such differences are greater for pathogenic variants than random candidate variants. In six autosomal recessive mitochondrial disease families, in which pathogenic homozygous variants have already been identified, our approach succeeded in prioritizing pathogenic mutations. Our method is applicable to single patients from recessive families with at least a few dozen control individuals from the same population; it is easy to use and is highly effective for detecting causative mutations in autosomal recessive families.

Systematic review and meta-analysis of Japanese familial Alzheimer's disease and FTDP-17.

Mutations in APP, PSEN1 and PSEN2 as the genetic causes of familial Alzheimer's disease (FAD) have been found in various ethnic populations. A substantial number of FAD pedigrees with mutations have been reported in the Japanese population; however, it remains unclear whether the genetic and clinical features of FAD in the Japanese population differ from those in other populations. To address this issue, we conducted a systematic review and meta-analysis of Japanese FAD and frontotemporal dementia with parkinsonism linked to chromosome 17 (FTDP-17) by literature search. Using this analysis, we identified 39 different PSEN1 mutations in 140 patients, 5 APP mutations in 35 patients and 16 MAPT mutations in 84 patients. There was no PSEN2 mutation among Japanese patients. The age at onset in Japanese FAD patients with PSEN1 mutations was significantly younger than that in patients with APP mutations. Kaplan-Meier analysis revealed that patients with MAPT mutations showed a shorter survival than patients with PSEN1 or APP mutations. Patients with mutations in different genes exhibit characteristic clinical presentations, suggesting that mutations in causative genes may modify the clinical presentations. By collecting and cataloging genetic and clinical information on Japanese FAD and FTDP-17, we developed an original database designated as Japanese Familial Alzheimer's Disease Database, which is accessible at http://alzdb.bri.niigata-u.ac.jp/.

KEGG OC: a large-scale automatic construction of taxonomy-based ortholog clusters.

The identification of orthologous genes in an increasing number of fully sequenced genomes is a challenging issue in recent genome science. Here we present KEGG OC (http://www.genome.jp/tools/oc/), a novel database of ortholog clusters (OCs). The current version of KEGG OC contains 1 176 030 OCs, obtained by clustering 8 357 175 genes in 2112 complete genomes (153 eukaryotes, 1830 bacteria and 129 archaea). The OCs were constructed by applying the quasi-clique-based clustering method to all possible protein coding genes in all complete genomes, based on their amino acid sequence similarities. It is computationally efficient to calculate OCs, which enables to regularly update the contents. KEGG OC has the following two features: (i) It consists of all complete genomes of a wide variety of organisms from three domains of life, and the number of organisms is the largest among the existing databases; and (ii) It is compatible with the KEGG database by sharing the same sets of genes and identifiers, which leads to seamless integration of OCs with useful components in KEGG such as biological pathways, pathway modules, functional hierarchy, diseases and drugs. The KEGG OC resources are accessible via OC Viewer that provides an interactive visualization of OCs at different taxonomic levels.

Plant Genome DataBase Japan (PGDBj): a portal website for the integration of plant genome-related databases.

The Plant Genome DataBase Japan (PGDBj, http://pgdbj.jp/?ln=en) is a portal website that aims to integrate plant genome-related information from databases (DBs) and the literature. The PGDBj is comprised of three component DBs and a cross-search engine, which provides a seamless search over the contents of the DBs. The three DBs are as follows. (i) The Ortholog DB, providing gene cluster information based on the amino acid sequence similarity. Over 500,000 amino acid sequences of 20 Viridiplantae species were subjected to reciprocal BLAST searches and clustered. Sequences from plant genome DBs (e.g. TAIR10 and RAP-DB) were also included in the cluster with a direct link to the original DB. (ii) The Plant Resource DB, integrating the SABRE DB, which provides cDNA and genome sequence resources accumulated and maintained in the RIKEN BioResource Center and National BioResource Projects. (iii) The DNA Marker DB, providing manually or automatically curated information of DNA markers, quantitative trait loci and related linkage maps, from the literature and external DBs. As the PGDBj targets various plant species, including model plants, algae, and crops important as food, fodder and biofuel, researchers in the field of basic biology as well as a wide range of agronomic fields are encouraged to perform searches using DNA sequences, gene names, traits and phenotypes of interest. The PGDBj will return the search results from the component DBs and various types of linked external DBs.

Polygenic effects on the risk of Alzheimer's disease in the Japanese population.

BACKGROUND: Polygenic effects have been proposed to account for some disease phenotypes; these effects are calculated as a polygenic risk score (PRS). This score is correlated with Alzheimer's disease (AD)-related phenotypes, such as biomarker abnormalities and brain atrophy, and is associated with conversion from mild cognitive impairment (MCI) to AD. However, the AD PRS has been examined mainly in Europeans, and owing to differences in genetic structure and lifestyle, it is unclear whether the same relationships between the PRS and AD-related phenotypes exist in non-European populations. In this study, we calculated and evaluated the AD PRS in Japanese individuals using genome-wide association study (GWAS) statistics from Europeans. METHODS: In this study, we calculated the AD PRS in 504 Japanese participants (145 cognitively unimpaired (CU) participants, 220 participants with late mild cognitive impairment (MCI), and 139 patients with mild AD dementia) enrolled in the Japanese Alzheimer's Disease Neuroimaging Initiative (J-ADNI) project. In order to evaluate the clinical value of this score, we (1) determined the polygenic effects on AD in the J-ADNI and validated it using two independent cohorts (a Japanese neuropathology (NP) cohort (n = 565) and the North American ADNI (NA-ADNI) cohort (n = 617)), (2) examined the AD-related phenotypes associated with the PRS, and (3) tested whether the PRS helps predict the conversion of MCI to AD. RESULTS: The PRS using 131 SNPs had an effect independent of APOE. The PRS differentiated between CU participants and AD patients with an area under the curve (AUC) of 0.755 when combined with the APOE variants. Similar AUC was obtained when PRS calculated by the NP and NA-ADNI cohorts was applied. In MCI patients, the PRS was associated with cerebrospinal fluid phosphorylated-tau levels (ß estimate = 0.235, p value = 0.026). MCI with a high PRS showed a significantly increased conversion to AD in APOE ε4 noncarriers with a hazard rate of 2.22. In addition, we also developed a PRS model adjusted for LD and observed similar results. CONCLUSIONS: We showed that the AD PRS is useful in the Japanese population, whose genetic structure is different from that of the European population. These findings suggest that the polygenicity of AD is partially common across ethnic differences.

Reduced plasma glucose by asparagine synthetase knockdown in the mouse liver.

Expression of the asparagine synthetase gene is dependent on extracellular glucose concentration. This gene was knocked-down in livers of male Balb/c mice by a hydrodynamic tail vein injection of small interfering RNA (siRNA) against the gene. This knockdown resulted in a significant decrease in plasma glucose concentration. These results suggested that asparagine synthetase is a novel protein that regulates plasma glucose levels.