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.

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.

Genome-wide association study of the risk of chronic kidney disease and kidney-related traits in the Japanese population: J-Kidney-Biobank.

Chronic kidney disease (CKD) is a syndrome characterized by a gradual loss of kidney function with decreased estimated glomerular filtration rate (eGFR), which may be accompanied by an increase in the urine albumin-to-creatinine ratio (UACR). Although trans-ethnic genome-wide association studies (GWASs) have been conducted for kidney-related traits, there have been few analyses in the Japanese population, especially for the UACR trait. In this study, we conducted a GWAS to identify loci related to multiple kidney-related traits in Japanese individuals. First, to detect loci associated with CKD, eGFR, and UACR, we performed separate GWASs with the following two datasets: 475 cases of CKD diagnosed at seven university hospitals and 3471 healthy subjects (dataset 1) and 3664 cases of CKD-suspected individuals with eGFR <60 ml/min/1.73 m2 or urinary protein ≥ 1+ and 5952 healthy subjects (dataset 2). Second, we performed a meta-analysis between these two datasets and detected the following associated loci: 10 loci for CKD, 9 loci for eGFR, and 22 loci for UACR. Among the loci detected, 22 have never been reported previously. Half of the significant loci for CKD were shared with those for eGFR, whereas most of the loci associated with UACR were different from those associated with CKD or eGFR. The GWAS of the Japanese population identified novel genetic components that were not previously detected. The results also suggest that the group primarily characterized by increased UACR possessed genetically different features from the group characterized by decreased eGFR.

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.

Establishment of the early prediction models of low-birth-weight reveals influential genetic and environmental factors: a prospective cohort study.

BACKGROUND: Low birth weight (LBW) is a leading cause of neonatal morbidity and mortality, and increases various disease risks across life stages. Prediction models of LBW have been developed before, but have limitations including small sample sizes, absence of genetic factors and no stratification of neonate into preterm and term birth groups. In this study, we challenged the development of early prediction models of LBW based on environmental and genetic factors in preterm and term birth groups, and clarified influential variables for LBW prediction. METHODS: We selected 22,711 neonates, their 21,581 mothers and 8,593 fathers from the Tohoku Medical Megabank Project Birth and Three-Generation cohort study. To establish early prediction models of LBW for preterm birth and term birth groups, we trained AI-based models using genetic and environmental factors of lifestyles. We then clarified influential environmental and genetic factors for predicting LBW in the term and preterm groups. RESULTS: We identified 2,327 (10.22%) LBW neonates consisting of 1,077 preterm births and 1,248 term births. Our early prediction models archived the area under curve 0.96 and 0.95 for term LBW and preterm LBW models, respectively. We revealed that environmental factors regarding eating habits and genetic features related to fetal growth were influential for predicting LBW in the term LBW model. On the other hand, we identified that genomic features related to toll-like receptor regulations and infection reactions are influential genetic factors for prediction in the preterm LBW model. CONCLUSIONS: We developed precise early prediction models of LBW based on lifestyle factors in the term birth group and genetic factors in the preterm birth group. Because of its accuracy and generalisability, our prediction model could contribute to risk assessment of LBW in the early stage of pregnancy and control LBW risk in the term birth group. Our prediction model could also contribute to precise prediction of LBW based on genetic factors in the preterm birth group. We then identified parental genetic and maternal environmental factors during pregnancy influencing LBW prediction, which are major targets for understanding the LBW to address serious burdens on newborns' health throughout life.

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.

Endogenization and excision of human herpesvirus 6 in human genomes.

Sequences homologous to human herpesvirus 6 (HHV-6) are integrated within the nuclear genome of about 1% of humans, but it is not clear how this came about. It is also uncertain whether integrated HHV-6 can reactivate into an infectious virus. HHV-6 integrates into telomeres, and this has recently been associated with polymorphisms affecting MOV10L1. MOV10L1 is located on the subtelomere of chromosome 22q (chr22q) and is required to make PIWI-interacting RNAs (piRNAs). As piRNAs block germline integration of transposons, piRNA-mediated repression of HHV-6 integration has been proposed to explain this association. In vitro, recombination of the HHV-6 genome along its terminal direct repeats (DRs) leads to excision from the telomere and viral reactivation, but the expected "solo-DR scar" has not been described in vivo. Here we screened for integrated HHV-6 in 7,485 Japanese subjects using whole-genome sequencing (WGS). Integrated HHV-6 was associated with polymorphisms on chr22q. However, in contrast to prior work, we find that the reported MOV10L1 polymorphism is physically linked to an ancient endogenous HHV-6A variant integrated into the telomere of chr22q in East Asians. Unexpectedly, an HHV-6B variant has also endogenized in chr22q; two endogenous HHV-6 variants at this locus thus account for 72% of all integrated HHV-6 in Japan. We also report human genomes carrying only one portion of the HHV-6B genome, a solo-DR, supporting in vivo excision and possible viral reactivation. Together these results explain the recently-reported association between integrated HHV-6 and MOV10L1/piRNAs, suggest potential exaptation of HHV-6 in its coevolution with human chr22q, and clarify the evolution and risk of reactivation of the only intact (non-retro)viral genome known to be present in human germlines.

Heterozygous calcyclin-binding protein/Siah1-interacting protein (CACYBP/SIP) gene pathogenic variant linked to a dominant family with paucity of interlobular bile duct.

Paucity of interlobular bile ducts (PILBD) is a heterogeneous disorder classified into two categories, syndromic and non-syndromic bile duct paucity. Syndromic PILBD is characterized by the presence of clinical manifestations of Alagille syndrome. Non-syndromic PILBD is caused by multiple diseases, such as metabolic and genetic disorders, infectious diseases, and inflammatory and immune disorders. We evaluated a family with a dominantly inherited PILBD, who presented with cholestasis at 1-2 months of age but spontaneously improved by 1 year of age. Next-generation sequencing analysis revealed a heterozygous CACYBP/SIP p.E177Q pathogenic variant. Calcyclin-binding protein and Siah1 interacting protein (CACYBP/SIP) form a ubiquitin ligase complex and induce proteasomal degradation of non-phosphorylated ß-catenin. Immunohistochemical analysis revealed a slight decrease in CACYBP and ß-catenin levels in the liver of patients in early infancy, which almost normalized by 13 months of age. The CACYBP/SIP p.E177Q pathogenic variant may form a more active or stable ubiquitin ligase complex that enhances the degradation of ß-catenin and delays the maturation of intrahepatic bile ducts. Our findings indicate that accurate regulation of the ß-catenin concentration is essential for the development of intrahepatic bile ducts and CACYBP/SIP pathogenic variant is a novel cause of PILDB.

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.

Maternal Baseline Characteristics and Perinatal Outcomes: The Tohoku Medical Megabank Project Birth and Three-Generation Cohort Study.

BACKGROUND: The Tohoku Medical Megabank Project Birth and Three-Generation Cohort Study was launched in 2013 to evaluate the complex interactions of genetic and environmental factors in multifactorial diseases. The present study describes the maternal baseline profile and perinatal data of participating mothers and infants. METHODS: Expectant mothers living in Miyagi Prefecture were recruited from obstetric facilities or affiliated centers between 2013 and 2017. Three sets of self-administered questionnaires were collected, and the medical records were reviewed to obtain precise information about each antenatal visit and each delivery. Biospecimens, including blood, urine, umbilical cord blood, and breast milk, were collected for the study biobank. The baseline maternal sociodemographic characteristics, results of screening tests, and obstetric outcomes were analyzed according to the maternal age group. RESULTS: A total of 23,406 pregnancies involving 23,730 fetuses resulted in 23,143 live births. Younger maternal participants had a tendency toward a higher incidence of threatened abortion and threatened premature labor, while older age groups exhibited a significantly higher rate of low lying placenta, placenta previa, gestational diabetes, and hypertensive disorders of pregnancy. CONCLUSIONS: The present study clearly shows the distribution of maternal baseline characteristics and the range of perinatal outcomes according to maternal age group. This cohort study can provide strategic information for creating breakthroughs in the pathophysiology of perinatal, developmental, and noncommunicable diseases by collaborative data visiting or sharing.

GWAS with principal component analysis identifies a gene comprehensively controlling rice architecture.

Elucidation of the genetic control of rice architecture is crucial due to the global demand for high crop yields. Rice architecture is a complex trait affected by plant height, tillering, and panicle morphology. In this study, principal component analysis (PCA) on 8 typical traits related to plant architecture revealed that the first principal component (PC), PC1, provided the most information on traits that determine rice architecture. A genome-wide association study (GWAS) using PC1 as a dependent variable was used to isolate a gene encoding rice, SPINDLY (OsSPY), that activates the gibberellin (GA) signal suppression protein SLR1. The effect of GA signaling on the regulation of rice architecture was confirmed in 9 types of isogenic plant having different levels of GA responsiveness. Further population genetics analysis demonstrated that the functional allele of OsSPY associated with semidwarfism and small panicles was selected in the process of rice breeding. In summary, the use of PCA in GWAS will aid in uncovering genes involved in traits with complex characteristics.

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.

Body mass index and colorectal cancer risk: A Mendelian randomization study.

Traditional observational studies have reported a positive association between higher body mass index (BMI) and the risk of colorectal cancer (CRC). However, evidence from other approaches to pursue the causal relationship between BMI and CRC is sparse. A two-sample Mendelian randomization (MR) study was undertaken using 68 single nucleotide polymorphisms (SNPs) from the Japanese genome-wide association study (GWAS) and 654 SNPs from the GWAS catalogue for BMI as sets of instrumental variables. For the analysis of SNP-BMI associations, we undertook a meta-analysis with 36 303 participants in the Japanese Consortium of Genetic Epidemiology studies (J-CGE), comprising normal populations. For the analysis of SNP-CRC associations, we utilized 7636 CRC cases and 37 141 controls from five studies in Japan, and undertook a meta-analysis. Mendelian randomization analysis of inverse-variance weighted method indicated that a one-unit (kg/m2 ) increase in genetically predicted BMI was associated with an odds ratio of 1.13 (95% confidence interval, 1.06-1.20; P value <.001) for CRC using the set of 68 SNPs, and an odds ratio of 1.07 (1.03-1.11, 0.001) for CRC using the set of 654 SNPs. Sensitivity analyses robustly showed increased odds ratios for CRC for every one-unit increase in genetically predicted BMI. Our MR analyses strongly support the evidence that higher BMI influences the risk of CRC. Although Asians are generally leaner than Europeans and North Americans, avoiding higher BMI seems to be important for the prevention of CRC in Asian populations.

Artificial intelligence powered statistical genetics in biobanks.

Large-scale, sometimes nationwide, prospective genomic cohorts biobanking rich biological specimens such as blood, urine and tissues, have been established and released their vast amount of data in several countries. These genetic and epidemiological resources are expected to allow investigators to disentangle genetic and environmental components conferring common complex diseases. There are, however, two major challenges to statistical genetics for this goal: small sample size-high dimensionality and multilayered-heterogenous endophenotypes. Rather counterintuitively, biobank data generally have small sample size relative to their data dimensionality consisting of genomic variation, lifestyle questionnaire, and sometimes their interaction. This is a widely acknowledged difficulty in data analysis, so-called "p¼n problem" in statistics or "curse of dimensionality" in machine-learning field. On the other hand, we have too many measurements of individual health status, which are endophenotypes, such as health check-up data, images, psychological test scores in addition to metabolomics and proteomics data. These endophenotypes are rich but not so tractable because of their worsen dimensionality, and substantial correlation, sometimes confusing causation among them. We have tried to overcome the problems inherent to biobank data, using statistical machine-learning and deep-learning technologies.

ALOX12 mutation in a family with dominantly inherited bleeding diathesis.

The arachidonic acid (AA) cascade plays a significant role in platelet aggregation. AA released from membrane phospholipids is metabolized by cyclooxygenase (COX) pathway to thromboxane A2 (TXA2) or by 12S-lipoxygenase (ALOX12) to 12-hydroperoxyeicosatetraenoic acid (12-HPETE). In contrast to a well-known role of the COX pathway in platelet aggregation, the role of ALOX12 is not well understood. Platelets of ALOX12-deficient mice exhibit increased sensitivity for ADP-induced aggregation. However, recent evidence strongly suggests a significant role of ALOX12 in platelet aggregation and calcium signaling. 12-HPETE potentiates thrombin- and thromboxane-induced platelet aggregation, and calcium signaling. Inhibition experiments of ALOX12 demonstrated decreased platelet aggregation and calcium signaling in stimulated platelets. We studied a family with a dominantly inherited bleeding diathesis using next-generation sequencing analysis. Platelet aggregation studies revealed that the proband's platelets had defective aggregation responses to ADP, TXA2 mimetic U46619, collagen, and AA, normal affinity of TXA2 receptor for U46619, and normal induction of GTPase activity upon stimulation with U46619. However, the production of inositol 1,4,5-triphosphate (IP3) was only increased up to 30% of the control upon U46619 stimulation, suggesting a defect in phospholipase C-ß2 (PLCB2) activation downstream from TXA2 receptors. Affected family members had no mutation of PLCB2, but had a heterozygous c.1946A > G (p.Tyr649Cys) mutation of ALOX12. ALOX12 activity in platelets from the affected members was decreased to 25-35% of the control. Our data strongly suggested that a heterozygous c.1946A > G ALOX12 mutation was a disease-causing mutation; however, further experiments are required to confirm the pathogenesis of ALOX12 mutation in platelet aggregation.

A genotype imputation method for de-identified haplotype reference information by using recurrent neural network.

Genotype imputation estimates the genotypes of unobserved variants using the genotype data of other observed variants based on a collection of haplotypes for thousands of individuals, which is known as a haplotype reference panel. In general, more accurate imputation results were obtained using a larger size of haplotype reference panel. Most of the existing genotype imputation methods explicitly require the haplotype reference panel in precise form, but the accessibility of haplotype data is often limited, due to the requirement of agreements from the donors. Since de-identified information such as summary statistics or model parameters can be used publicly, imputation methods using de-identified haplotype reference information might be useful to enhance the quality of imputation results under the condition where the access of the haplotype data is limited. In this study, we proposed a novel imputation method that handles the reference panel as its model parameters by using bidirectional recurrent neural network (RNN). The model parameters are presented in the form of de-identified information from which the restoration of the genotype data at the individual-level is almost impossible. We demonstrated that the proposed method provides comparable imputation accuracy when compared with the existing imputation methods using haplotype datasets from the 1000 Genomes Project (1KGP) and the Haplotype Reference Consortium. We also considered a scenario where a subset of haplotypes is made available only in de-identified form for the haplotype reference panel. In the evaluation using the 1KGP dataset under the scenario, the imputation accuracy of the proposed method is much higher than that of the existing imputation methods. We therefore conclude that our RNN-based method is quite promising to further promote the data-sharing of sensitive genome data under the recent movement for the protection of individuals' privacy.