jMorp: Japanese Multi-Omics Reference Panel update report 2023.

Modern medicine is increasingly focused on personalized medicine, and multi-omics data is crucial in understanding biological phenomena and disease mechanisms. Each ethnic group has its unique genetic background with specific genomic variations influencing disease risk and drug response. Therefore, multi-omics data from specific ethnic populations are essential for the effective implementation of personalized medicine. Various prospective cohort studies, such as the UK Biobank, All of Us and Lifelines, have been conducted worldwide. The Tohoku Medical Megabank project was initiated after the Great East Japan Earthquake in 2011. It collects biological specimens and conducts genome and omics analyses to build a basis for personalized medicine. Summary statistical data from these analyses are available in the jMorp web database (https://jmorp.megabank.tohoku.ac.jp), which provides a multidimensional approach to the diversity of the Japanese population. jMorp was launched in 2015 as a public database for plasma metabolome and proteome analyses and has been continuously updated. The current update will significantly expand the scale of the data (metabolome, genome, transcriptome, and metagenome). In addition, the user interface and backend server implementations were rewritten to improve the connectivity between the items stored in jMorp. This paper provides an overview of the new version of the jMorp.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Genetic etiology of truncus arteriosus excluding 22q11.2 deletion syndrome and identification of c.1617del, a prevalent variant in TMEM260, in the Japanese population.

Truncus Arteriosus (TA) is a congenital heart disease characterized by a single common blood vessel emerging from the right and left ventricles instead of the main pulmonary artery and aorta. TA accounts for 4% of all critical congenital heart diseases. The most common cause of TA is 22q11.2 deletion syndrome, accounting for 12-35% of all TA cases. However, no major causes of TA other than 22q11.2 deletion have been reported. We performed whole-genome sequencing of 11 Japanese patients having TA without 22q11.2 deletion. Among five patients, we identified pathogenic variants in TMEM260; the biallelic loss-of-function variants of which have recently been associated with structural heart defects and renal anomalies syndrome (SHDRA). In one patient, we identified a de novo pathogenic variant in GATA6, and in another patient, we identified a de novo probably pathogenic variant in NOTCH1. Notably, we identified a prevalent variant in TMEM260 (ENST00000261556.6), c.1617del (p.Trp539Cysfs*9), in 8/22 alleles among the 11 patients. The c.1617del variant was estimated to occur approximately 23 kiloyears ago. Based on the allele frequency of the c.1617del variant in the Japanese population (0.36%), approximately 26% of Japanese patients afflicted with TA could harbor homozygous c.1617del variants. This study highlights TMEM260, especially c.1617del, as a major genetic cause of TA in the Japanese population.

A fine-scale genetic map of the Japanese population.

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

Deep exome sequencing identifies enrichment of deleterious mosaic variants in neurodevelopmental disorder genes and mitochondrial tRNA regions in bipolar disorder.

Bipolar disorder (BD) is a global medical issue, afflicting around 1% of the population with manic and depressive episodes. Despite various genetic studies, the genetic architecture and pathogenesis of BD have not been fully resolved. Besides germline variants, postzygotic mosaic variants are proposed as new candidate mechanisms contributing to BD. Here, we performed extensive deep exome sequencing (DES, ~300×) and validation experiments to investigate the roles of mosaic variants in BD with 235 BD cases (194 probands of trios and 41 single cases) and 39 controls. We found an enrichment of developmental disorder (DD) genes in the genes hit by deleterious mosaic variants in BD (P = 0.000552), including a ClinVar-registered pathogenic variant in ARID2. An enrichment of deleterious mosaic variants was also observed for autism spectrum disorder (ASD) genes (P = 0.000428). The proteins coded by the DD/ASD genes with non-synonymous mosaic variants in BD form more protein-protein interaction than expected, suggesting molecular mechanisms shared with DD/ASD but restricted to a subset of cells in BD. We also found significant enrichment of mitochondrial heteroplasmic variants, another class of mosaic variants, in mitochondrial tRNA genes in BD (P = 0.0102). Among them, recurrent m.3243 A > G variants known as causal for mitochondrial diseases were found in two unrelated BD probands with allele fractions of 5-12%, lower than in mitochondrial diseases. Despite the limitation of using peripheral tissues, our DES investigation supports the possible contribution of deleterious mosaic variants in the nuclear genome responsible for severer phenotypes, such as DD/ASD, to the risk of BD and further demonstrates that the same paradigm can be applied to the mitochondrial genome. These results, as well as the enrichment of heteroplasmic mitochondrial tRNA variants in BD, add a new piece to the understanding of the genetic architecture of BD and provide general insights into the pathological roles of mosaic variants in human diseases.

Long-term clinical observation of patients with heterozygous KIF1A variants.

KIF1A-related disorders (KRDs) encompass recessive and dominant variants with wide clinical variability. Recent genetic investigations have expanded the clinical phenotypes of heterozygous KIF1A variants. However, there have been a few long-term observational studies of patients with heterozygous KIF1A variants. A retrospective chart review of consecutive patients diagnosed with spastic paraplegia at Miyagi Children's Hospital from 2016 to 2020 identified six patients with heterozygous KIF1A variants. To understand the long-term changes in clinical symptoms, we examined these patients in terms of their characteristics, clinical symptoms, results of electrophysiological and neuroimaging studies, and genetic testing. The median follow-up period was 30 years (4-44 years). This long-term observational study showed that early developmental delay and equinus gait, or unsteady gait, are the first signs of disease onset, appearing with the commencement of independent walking. In addition, later age-related progression was observed in spastic paraplegia, and the appearance of axonal neuropathy and reduced visual acuity were characteristic features of the late disease phenotype. Brain imaging showed age-related progression of cerebellar atrophy and the appearance of hyperintensity of optic radiation on T2WI and FLAIR imaging. Long-term follow-up revealed a pattern of steady progression and a variety of clinical symptoms, including spastic paraplegia, peripheral neuropathy, reduced visual acuity, and some degree of cerebellar ataxia. Clinical variability between patients was observed to some extent, and therefore, further studies are required to determine the phenotype-genotype correlation.

Treatment of ZC4H2 Variant-Associated Spastic Paraplegia with Selective Dorsal Rhizotomy and Intensive Postoperative Rehabilitation: A Case Report.

Selective dorsal rhizotomy (SDR) has been used to treat children with spastic cerebral palsy (CP), and its beneficial effect on quality of life and ambulation has been confirmed in long-term follow-up studies. However, the role of SDR in the treatment of spasticity in patients with hereditary spastic paraplegia (HSP) and related disorders is not well-established. Here, we report the first patient with the ZC4H2 variant who underwent SDR to treat spastic paraplegia. Abnormal gait was discovered during a regular checkup at the age of 3 years and 9 months, and she was diagnosed with spastic paraplegia. She was heterozygous for the ZC4H2 variant and underwent SDR at the age of 5 years and 11 months, which alleviated the spasticity. The patient underwent inpatient postoperative rehabilitation for 4 months and continued outpatient physiotherapy after discharge. The Gross Motor Function Measure-88 score and maximum walking speed decreased transiently 1 month postoperatively, but gradually recovered, and continuously improved 6 months postoperatively. SDR and postoperative intensive rehabilitation were effective in improving motor and walking functions up to 6 months after surgery, although long-term follow-up is needed to draw conclusions.

Familial Paget's disease of bone with ocular manifestations and a novel TNFRSF11A duplication variant (72dup27).

INTRODUCTION: Paget's disease of bone (PDB) is a skeletal disorder characterized by disorganized bone remodeling due to abnormal osteoclasts. Tumor necrosis factor receptor superfamily member 11A (TNFRSF11A) gene encodes the receptor activator of nuclear factor kappa B (RANK), which has a critical role in osteoclast function. There are five types of rare PDB and related osteolytic disorders due to TNFRSF11A tandem duplication variants so far, including familial expansile osteolysis (84dup18), expansile skeletal hyperphosphatasia (84dup15), early-onset familial PDB (77dup27), juvenile PDB (87dup15), and panostotic expansile bone disease (90dup12). MATERIALS AND METHODS: We reviewed a Japanese family with PDB, and performed whole-genome sequencing to identify a causative variant. RESULTS: This family had bone symptoms, hyperphosphatasia, hearing loss, tooth loss, and ocular manifestations such as angioid streaks or early-onset glaucoma. We identified a novel duplication variant of TNFRSF11A (72dup27). Angioid streaks were recognized in Juvenile Paget's disease due to loss-of-function variants in the gene TNFRSF11B, and thought to be specific for this disease. However, the novel recognition of angioid streaks in our family raised the possibility of occurrence even in bone disorders due to TNFRSF11A duplication variants and the association of RANKL-RANK signal pathway as the pathogenesis. Glaucoma has conversely not been reported in any case of Paget's disease. It is not certain whether glaucoma is coincidental or specific for PDB with 72dup27. CONCLUSION: Our new findings might suggest a broad spectrum of phenotypes in bone disorders with TNFRSF11A duplication variants.

jMorp updates in 2020: large enhancement of multi-omics data resources on the general Japanese population.

In the Tohoku Medical Megabank project, genome and omics analyses of participants in two cohort studies were performed. A part of the data is available at the Japanese Multi Omics Reference Panel (jMorp; https://jmorp.megabank.tohoku.ac.jp) as a web-based database, as reported in our previous manuscript published in Nucleic Acid Research in 2018. At that time, jMorp mainly consisted of metabolome data; however, now genome, methylome, and transcriptome data have been integrated in addition to the enhancement of the number of samples for the metabolome data. For genomic data, jMorp provides a Japanese reference sequence obtained using de novo assembly of sequences from three Japanese individuals and allele frequencies obtained using whole-genome sequencing of 8,380 Japanese individuals. In addition, the omics data include methylome and transcriptome data from ∼300 samples and distribution of concentrations of more than 755 metabolites obtained using high-throughput nuclear magnetic resonance and high-sensitivity mass spectrometry. In summary, jMorp now provides four different kinds of omics data (genome, methylome, transcriptome, and metabolome), with a user-friendly web interface. This will be a useful scientific data resource on the general population for the discovery of disease biomarkers and personalized disease prevention and early diagnosis.

Involvement of the Restoration of Cerebral Blood Flow and Maintenance of eNOS Expression in the Prophylactic Protective Effect of the Novel Ferulic Acid Derivative FAD012 against Ischemia/Reperfusion Injuries in Rats.

Tissue plasminogen activator, aiming to restore cerebral blood flow (CBF), has been used for acute ischemic strokes in clinics; however, its narrow therapeutic time window remains a serious concern. To develop novel prophylactic drugs to alleviate cerebral ischemia/reperfusion injuries, ferulic acid derivative 012 (FAD012) was synthesized and showed comparable antioxidant properties to ferulic acid (FA) and probably possesses the potent ability to cross the blood-brain barrier. A more potent cytoprotective effect of FAD012 against H2O2-induced cytotoxicity in PC12 cells was also observed. In vivo toxicity was not observed in rats given a long-term oral administration of FAD012, indicating its good tolerability. A one-week-course oral administration of FAD012 significantly alleviated middle cerebral artery occlusion (MCAO)-induced cerebral ischemia/reperfusion injuries in rats, accompanied by the restoration of CBF and endothelial nitrogen oxide synthetase (eNOS) expression. Treatment with FAD012 significantly restored the cell viability and eNOS expression damaged by H2O2, used to mimic MCAO-triggered oxidative stress, in rat brain microvascular endothelial cells. Our findings suggested that FAD012 protected the viability of vascular endothelium and maintained eNOS expression, ultimately contributing to the restoration of CBF, and may provide a rationale for the development of FAD012 into an effective prophylactic drug for patients at high risk of stroke.

The longest reported sibling survivors of a severe form of congenital myasthenic syndrome with the ALG14 pathogenic variant.

Congenital myasthenic syndromes (CMS) is a group of diseases that causes abnormalities at the neuromuscular junction owing to genetic anomalies. The pathogenic variant in ALG14 results in a severe pathological form of CMS causing end-plate acetylcholine receptor deficiency. Here, we report the cases of two siblings with CMS associated with a novel variant in ALG14. Immediately after birth, they showed hypotonia and multiple joint contractures with low Apgar scores. Ptosis, low-set ears, and high-arched palate were noted. Deep tendon reflexes were symmetrical. They showed worsening swallowing and respiratory problems; hence, nasal feeding and tracheotomy were performed. Cranial magnetic resonance imaging scans revealed delayed myelination and cerebral atrophy. Exome sequencing indicated that the siblings had novel compound heterozygous missense variants, c.590T>G (p.Val197Gly) and c.433G>A (p.Gly145Arg), in exon 4 of ALG14. Repetitive nerve stimulation test showed an abnormal decrease in compound muscle action potential. After treatment with pyridostigmine, the time off the respirator increased. Their epileptic seizures were well controlled by anti-epileptic drugs. Their clinical course is stable even now at the ages of 5 and 2 years, making them the longest reported survivors of a severe form of CMS with the ALG14 variant thus far.

Two Siblings with Cerebellar Ataxia, Mental Retardation, and Disequilibrium Syndrome 4 and a Novel Variant of ATP8A2.

Cerebellar ataxia, mental retardation, and disequilibrium syndrome 4 (CAMRQ4) is early onset neuromotor disorder and intellectual disabilities caused by variants of ATP8A2. We report sibling cases and systematically analyze previous literature to increase our understanding of CAMRQ4. Japanese siblings presented with athetotic movements at 1 and 2 months of age. They also had ptosis, ophthalmoplegia, feeding difficulty, hypotonia, and severely delayed development. One patient had retinal degeneration and optic atrophy. Flattening of the auditory brainstem responses and areflexia developed. At the last follow-up, neither patient could sit or achieve head control, although some nonverbal communication was preserved. Whole exome sequencing revealed compound heterozygous variants of ATP8A2: NM_016529.6:c.[1741C>T];[2158C>T] p.[(Arg581*)];[(Arg720*)]. The p.(Arg581*) variant has been reported, while the variant p.(Arg720*) was novel. The symptoms did not progress in the early period of development, which makes it difficult to distinguish from dyskinetic cerebral palsy, particularly in solitary cases. However, visual and hearing impairments associated with involuntary movements and severe developmental delay may be a clue to suspect CAMRQ4.

Statistical image processing quantifies the changes in cytoplasmic texture associated with aging in Caenorhabditis elegans oocytes.

BACKGROUND: Oocyte quality decreases with aging, thereby increasing errors in fertilization, chromosome segregation, and embryonic cleavage. Oocyte appearance also changes with aging, suggesting a functional relationship between oocyte quality and appearance. However, no methods are available to objectively quantify age-associated changes in oocyte appearance. RESULTS: We show that statistical image processing of Nomarski differential interference contrast microscopy images can be used to quantify age-associated changes in oocyte appearance in the nematode Caenorhabditis elegans. Max-min value (mean difference between the maximum and minimum intensities within each moving window) quantitatively characterized the difference in oocyte cytoplasmic texture between 1- and 3-day-old adults (Day 1 and Day 3 oocytes, respectively). With an appropriate parameter set, the gray level co-occurrence matrix (GLCM)-based texture feature Correlation (COR) more sensitively characterized this difference than the Max-min Value. Manipulating the smoothness of and/or adding irregular structures to the cytoplasmic texture of Day 1 oocyte images reproduced the difference in Max-min Value but not in COR between Day 1 and Day 3 oocytes. Increasing the size of granules in synthetic images recapitulated the age-associated changes in COR. Manual measurements validated that the cytoplasmic granules in oocytes become larger with aging. CONCLUSIONS: The Max-min value and COR objectively quantify age-related changes in C. elegans oocyte in Nomarski DIC microscopy images. Our methods provide new opportunities for understanding the mechanism underlying oocyte aging.

Characterization of 6-bromoferulic acid as a novel common-use matrix for matrix-assisted laser desorption/ionization time-of-flight mass spectrometry.

RATIONALE: Ferulic acid (FA) is a standard matrix used for analyzing proteins. In this study, the ability of a halogenated FA to serve as an effective MALDI matrix was investigated. Various halogenated FAs were synthesized, and the characteristics and performance of each were compared with those of the standard matrices α-cyano-4-hydroxycinnamic acid (CHCA) and 2,5-dihydrobenzoic acid (DHBA). METHODS: The abilities of 6-bromoferulic acid (6-BFA), ferulic acid (FA), and eight other halogenated FA derivatives to ionize eight synthetic peptides were examined. Absorption measurements, MM2 structure optimizations, and proton affinity (PA) calculations were also performed for 6-BFA and FA. The suitabilities of these compounds as matrices for matrix-assisted laser desorption/ionization (MALDI) for lipids, sugar chains, polymers, cyanocobalamin, synthetic peptides, and tryptic peptides originating from two types of serum proteins were also tested. RESULTS: The 6-position of FA was found to be the best site for introducing a bromine because the generated compound allowed facile detection of cyanocobalamin and several peptides. 6-BFA exhibited good sensitivity for large peptides (3-5 kDa) and peptides containing acidic amino acids or proline. 6-BFA was also shown to be a suitable matrix for tandem mass spectrometry (MS/MS) analysis when using MALDI time-of-flight (TOF) mass spectrometry (MS) with a quadrupole ion trap (QIT) system. CONCLUSIONS: The properties of 6-BFA as a MALDI matrix differed from those of DHBA and CHCA. 6-BFA appears to be a useful matrix for de novo sequencing using MALDI-QIT-TOF-MS.

Potent HIV-1 protease inhibitors incorporating squaramide-derived P2 ligands: Design, synthesis, and biological evaluation.

We describe the design, synthesis, and biological evaluation of novel HIV-1 protease inhibitors containing a squaramide-derived scaffold as the P2 ligand in combination with a (R)-hydroxyethylamine sulfonamide isostere. Inhibitor 3h with an N-methyl-3-(R)-aminotetrahydrofuranyl squaramide P2-ligand displayed an HIV-1 protease inhibitory Ki value of 0.51 nM. An energy minimized model of 3h revealed the major molecular interactions between HIV-1 protease active site and the tetrahydrofuranyl squaramide scaffold that may be responsible for its potent activity.

Design of novel HIV-1 protease inhibitors incorporating isophthalamide-derived P2-P3 ligands: Synthesis, biological evaluation and X-ray structural studies of inhibitor-HIV-1 protease complex.

Based upon molecular insights from the X-ray structures of inhibitor-bound HIV-1 protease complexes, we have designed a series of isophthalamide-derived inhibitors incorporating substituted pyrrolidines, piperidines and thiazolidines as P2-P3 ligands for specific interactions in the S2-S3 extended site. Compound 4b has shown an enzyme Ki of 0.025nM and antiviral IC50 of 69nM. An X-ray crystal structure of inhibitor 4b-HIV-1 protease complex was determined at 1.33Å resolution. We have also determined X-ray structure of 3b-bound HIV-1 protease at 1.27Å resolution. These structures revealed important molecular insight into the inhibitor-HIV-1 protease interactions in the active site.

Systematic evaluation of methyl ester bioisosteres in the context of developing alkenyldiarylmethanes (ADAMs) as non-nucleoside reverse transcriptase inhibitors (NNRTIs) for anti-HIV-1 chemotherapy.

The alkenyldiarylmethanes (ADAMs) are a class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) targeting HIV-1. Four chemically and metabolically stabilized ADAMs incorporating N-methoxyimidoyl halide replacements of the methyl esters of the lead compound were previously reported. In this study, twenty-five new ADAMs were synthesized in order to investigate the biological consequences of installing nine different methyl ester bioisosteres at three different locations. Attempts to define a universal rank order of methyl ester bioisosteres and discover the 'best' one in terms of inhibitory activity versus HIV-1 reverse transcriptase (RT) led to the realization that the potencies are critically dependent on the surrounding structure at each location, and therefore the definition of universal rank order is impossible. This investigation produced several new non-nucleoside reverse transcriptase inhibitors in which all three of the three methyl esters of the lead compound were replaced by methyl ester bioisosteres, resulting in compounds that are more potent as HIV-1 RT inhibitors and antiviral agents than the lead compound itself and are expected to also be more metabolically stable than the lead compound.

Improving spinning disk confocal microscopy by preventing pinhole cross-talk for intravital imaging.

A recent key requirement in life sciences is the observation of biological processes in their natural in vivo context. However, imaging techniques that allow fast imaging with higher resolution in 3D thick specimens are still limited. Spinning disk confocal microscopy using a Yokogawa Confocal Scanner Unit, which offers high-speed multipoint confocal live imaging, has been found to have wide utility among cell biologists. A conventional Confocal Scanner Unit configuration, however, is not optimized for thick specimens, for which the background noise attributed to "pinhole cross-talk," which is unintended pinhole transmission of out-of-focus light, limits overall performance in focal discrimination and reduces confocal capability. Here, we improve spinning disk confocal microscopy by eliminating pinhole cross-talk. First, the amount of pinhole cross-talk is reduced by increasing the interpinhole distance. Second, the generation of out-of-focus light is prevented by two-photon excitation that achieves selective-plane illumination. We evaluate the effect of these modifications and test the applicability to the live imaging of green fluorescent protein-expressing model animals. As demonstrated by visualizing the fine details of the 3D cell shape and submicron-size cytoskeletal structures inside animals, these strategies dramatically improve higher-resolution intravital imaging.

Structure-based design, synthesis, X-ray studies, and biological evaluation of novel HIV-1 protease inhibitors containing isophthalamide-derived P2-ligands.

We describe the design, synthesis and biological evaluation of a series of novel HIV-1 protease inhibitors bearing isophthalamide derivatives as the P2-P3 ligands. We have investigated a range of acyclic and heterocyclic amides as the extended P2-P3 ligands. These inhibitors displayed good to excellent HIV-1 protease inhibitory activity. Also, a number of inhibitors showed very good antiviral activity in MT cells. Compound 5n has shown an enzyme Ki of 0.17 nM and antiviral IC50 of 14 nM. An X-ray crystal structure of inhibitor 5o-bound to HIV-1 protease was determined at 1.11Å resolution. This structure revealed important molecular insight into the inhibitor-HIV-1 protease interactions in the active site.