Koichiro Higasa

Human genetic variation database, a reference database of genetic variations in the Japanese population

Whole-genome and -exome resequencing using next-generation sequencers is a powerful approach for identifying genomic variations that are associated with diseases. However, systematic strategies for prioritizing causative variants from many candidates to explain the disease phenotype are still far from being established, because the population-specific frequency spectrum of genetic variation has not been characterized. Here, we have collected exomic genetic variation from 1208 Japanese individuals through a collaborative effort, and aggregated the data into a prevailing catalog. In total, we identified 156 622 previously unreported variants. The allele frequencies for the majority (88.8%) were lower than 0.5% in allele frequency and predicted to be functionally deleterious. In addition, we have constructed a Japanese-specific major allele reference genome by which the number of unique mapping of the short reads in our data has increased 0.045% on average. Our results illustrate the importance of constructing an ethnicity-specific reference genome for identifying rare variants. All the collected data were centralized to a newly developed database to serve as useful resources for exploring pathogenic variations. Public access to the database is available at http://www.genome.med.kyoto-u.ac.jp/SnpDB/.

Large-scale East-Asian eQTL mapping reveals novel candidate genes for LD mapping and the genomic landscape of transcriptional effects of sequence variants

Profiles of sequence variants that influence gene transcription are very important for understanding mechanisms that affect phenotypic variation and disease susceptibility. Using genotypes at 1.4 million SNPs and a comprehensive transcriptional profile of 15,454 coding genes and 6,113 lincRNA genes obtained from peripheral blood cells of 298 Japanese individuals, we mapped expression quantitative trait loci (eQTLs). We identified 3,804 cis-eQTLs (within 500 kb from target genes) and 165 trans-eQTLs (>500 kb away or on different chromosomes). Cis-eQTLs were often located in transcribed or adjacent regions of genes; among these regions, 5′ untranslated regions and 5′ flanking regions had the largest effects. Epigenetic evidence for regulatory potential accumulated in public databases explained the magnitude of the effects of our eQTLs. Cis-eQTLs were often located near the respective target genes, if not within genes. Large effect sizes were observed with eQTLs near target genes, and effect sizes were obviously attenuated as the eQTL distance from the gene increased. Using a very stringent significance threshold, we identified 165 large-effect trans-eQTLs. We used our eQTL map to assess 8,069 disease-associated SNPs identified in 1,436 genome-wide association studies (GWAS). We identified genes that might be truly causative, but GWAS might have failed to identify for 148 out of the GWAS-identified SNPs; for example, TUFM (P = 3.3E-48) was identified for inflammatory bowel disease (early onset); ZFP90 (P = 4.4E-34) for ulcerative colitis; and IDUA (P = 2.2E-11) for Parkinsons disease. We identified four genes (P<2.0E-14) that might be related to three diseases and two hematological traits; each expression is regulated by trans-eQTLs on a different chromosome than the gene.

A Genome-Wide Association Study Identified AFF1 as a Susceptibility Locus for Systemic Lupus Eyrthematosus in Japanese

Systemic lupus erythematosus (SLE) is an autoimmune disease that causes multiple organ damage. Although recent genome-wide association studies (GWAS) have contributed to discovery of SLE susceptibility genes, few studies has been performed in Asian populations. Here, we report a GWAS for SLE examining 891 SLE cases and 3,384 controls and multi-stage replication studies examining 1,387 SLE cases and 28,564 controls in Japanese subjects. Considering that expression quantitative trait loci (eQTLs) have been implicated in genetic risks for autoimmune diseases, we integrated an eQTL study into the results of the GWAS. We observed enrichments of cis-eQTL positive loci among the known SLE susceptibility loci (30.8%) compared to the genome-wide SNPs (6.9%). In addition, we identified a novel association of a variant in the AF4/FMR2 family, member 1 (AFF1) gene at 4q21 with SLE susceptibility (rs340630; P = 8.3×10−9, odds ratio = 1.21). The risk A allele of rs340630 demonstrated a cis-eQTL effect on the AFF1 transcript with enhanced expression levels (P<0.05). As AFF1 transcripts were prominently expressed in CD4+ and CD19+ peripheral blood lymphocytes, up-regulation of AFF1 may cause the abnormality in these lymphocytes, leading to disease onset.

Identification of nine novel loci associated with white blood cell subtypes in a Japanese population.

White blood cells (WBCs) mediate immune systems and consist of various subtypes with distinct roles. Elucidation of the mechanism that regulates the counts of the WBC subtypes would provide useful insights into both the etiology of the immune system and disease pathogenesis. In this study, we report results of genome-wide association studies (GWAS) and a replication study for the counts of the 5 main WBC subtypes (neutrophils, lymphocytes, monocytes, basophils, and eosinophils) using 14,792 Japanese subjects enrolled in the BioBank Japan Project. We identified 12 significantly associated loci that satisfied the genome-wide significance threshold of P<5.0×10−8, of which 9 loci were novel (the CDK6 locus for the neutrophil count; the ITGA4, MLZE, STXBP6 loci, and the MHC region for the monocyte count; the SLC45A3-NUCKS1, GATA2, NAALAD2, ERG loci for the basophil count). We further evaluated associations in the identified loci using 15,600 subjects from Caucasian populations. These WBC subtype-related loci demonstrated a variety of patterns of pleiotropic associations within the WBC subtypes, or with total WBC count, platelet count, or red blood cell-related traits (n = 30,454), which suggests unique and common functional roles of these loci in the processes of hematopoiesis. This study should contribute to the understanding of the genetic backgrounds of the WBC subtypes and hematological traits.

Novel germline mutation in the transmembrane domain of HER2 in familial lung adenocarcinomas

We encountered a family of Japanese descent in which multiple members developed lung cancer. Using whole-exome sequencing, we identified a novel germline mutation in the transmembrane domain of the human epidermal growth factor receptor 2 (HER2) gene (G660D). A novel somatic mutation (V659E) was also detected in the transmembrane domain of HER2 in one of 253 sporadic lung adenocarcinomas. Because the transmembrane domain of HER2 is considered to be responsible for the dimerization and subsequent activation of the HER family and downstream signaling pathways, we performed functional analyses of these HER2 mutants. Mutant HER2 G660D and V659E proteins were more stable than wild-type protein. Both the G660D and V659E mutants activated Akt. In addition, they activated p38, which is thought to promote cell proliferation in lung adenocarcinoma. Our findings strongly suggest that mutations in the transmembrane domain of HER2 may be oncogenic, causing hereditary and sporadic lung adenocarcinomas.

A Comprehensive Genomic Analysis Reveals the Genetic Landscape of Mitochondrial Respiratory Chain Complex Deficiencies

Mitochondrial disorders have the highest incidence among congenital metabolic disorders characterized by biochemical respiratory chain complex deficiencies. It occurs at a rate of 1 in 5,000 births, and has phenotypic and genetic heterogeneity. Mutations in about 1,500 nuclear encoded mitochondrial proteins may cause mitochondrial dysfunction of energy production and mitochondrial disorders. More than 250 genes that cause mitochondrial disorders have been reported to date. However exact genetic diagnosis for patients still remained largely unknown. To reveal this heterogeneity, we performed comprehensive genomic analyses for 142 patients with childhood-onset mitochondrial respiratory chain complex deficiencies. The approach includes whole mtDNA and exome analyses using high-throughput sequencing, and chromosomal aberration analyses using high-density oligonucleotide arrays. We identified 37 novel mutations in known mitochondrial disease genes and 3 mitochondria-related genes (MRPS23, QRSL1, and PNPLA4) as novel causative genes. We also identified 2 genes known to cause monogenic diseases (MECP2 and TNNI3) and 3 chromosomal aberrations (6q24.3-q25.1, 17p12, and 22q11.21) as causes in this cohort. Our approaches enhance the ability to identify pathogenic gene mutations in patients with biochemically defined mitochondrial respiratory chain complex deficiencies in clinical settings. They also underscore clinical and genetic heterogeneity and will improve patient care of this complex disorder.

ERBB4 Mutations that Disrupt the Neuregulin-ErbB4 Pathway Cause Amyotrophic Lateral Sclerosis Type 19

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disorder characterized by the degeneration of motor neurons and typically results in death within 3-5 years from onset. Familial ALS (FALS) comprises 5%-10% of ALS cases, and the identification of genes associated with FALS is indispensable to elucidating the molecular pathogenesis. We identified a Japanese family affected by late-onset, autosomal-dominant ALS in which mutations in genes known to be associated with FALS were excluded. A whole- genome sequencing and parametric linkage analysis under the assumption of an autosomal-dominant mode of inheritance with incomplete penetrance revealed the mutation c.2780G>A (p. Arg927Gln) in ERBB4. An extensive mutational analysis revealed the same mutation in a Canadian individual with familial ALS and a de novo mutation, c.3823C>T (p. Arg1275Trp), in a Japanese simplex case. These amino acid substitutions involve amino acids highly conserved among species, are predicted as probably damaging, and are located within a tyrosine kinase domain (p. Arg927Gln) or a C-terminal domain (p. Arg1275Trp), both of which mediate essential functions of ErbB4 as a receptor tyrosine kinase. Functional analysis revealed that these mutations led to a reduced autophosphorylation of ErbB4 upon neuregulin-1 (NRG-1) stimulation. Clinical presentations of the individuals with mutations were characterized by the involvement of both upper and lower motor neurons, a lack of obvious cognitive dysfunction, and relatively slow progression. This study indicates that disruption of the neuregulin-ErbB4 pathway is involved in the pathogenesis of ALS and potentially paves the way for the development of innovative therapeutic strategies such using NRGs or their agonists to upregulate ErbB4 functions.

Comprehensive Molecular Diagnosis of a Large Cohort of Japanese Retinitis Pigmentosa and Usher Syndrome Patients by Next-Generation Sequencing

PURPOSE Retinitis pigmentosa (RP), a major cause of blindness in developed countries, has multiple causative genes; its prevalence differs by ethnicity. Usher syndrome is the most common form of syndromic RP and is accompanied by hearing impairment. Although molecular diagnosis is challenging, recent technological advances such as targeted high-throughput resequencing are efficient screening tools. METHODS We performed comprehensive molecular testing in 329 Japanese RP and Usher syndrome patients by using a custom capture panel that covered the coding exons and exon/intron boundaries of all 193 known inherited eye disease genes combined with Illumina HiSequation 2500. Candidate variants were screened using systematic data analyses, and their potential pathogenicity was assessed according to the frequency of the variants in normal populations, in silico prediction tools, and compatibility with known phenotypes or inheritance patterns. RESULTS Molecular diagnoses were made in 115/317 RP patients (36.3%) and 6/12 Usher syndrome patients (50%). We identified 104 distinct mutations, including 66 novel mutations. EYS, USH2A, and RHO were common causative genes. In particular, mutations in EYS accounted for 15.0% of the autosomal recessive/simplex RP patients or 10.7% of the entire RP cohort. Among the 189 previously reported mutations detected in the current study, 55 (29.1%) were found commonly in Japanese or other public databases and were excluded from molecular diagnoses. CONCLUSIONS By screening a large cohort of patients, this study catalogued the genetic variations involved in RP and Usher syndrome in a Japanese population and highlighted the different distribution of causative genes among populations.

CSF1R mutations identified in three families with autosomal dominantly inherited leukoencephalopathy

Genetic and phenotypic heterogeneities are considerably high in adult‐onset leukoencephalopathy, in which comprehensive mutational analyses of the candidate genes by conventional methods are too laborious. We applied exome sequencing to conduct a comprehensive mutational analysis of genes for autosomal dominant leukoencephalopathies. Genomic DNA samples from four patients of three families with autosomal dominantly inherited adult‐onset leukodystrophy were subjected to exome sequencing. On the basis of the results, 21 patients with adult‐onset sporadic leukodystrophy and one patient with pathologically proven HDLS were additionally screened for CSF1R mutations. Exome sequencing identified heterozygous CSF1R mutations (p.I794T and p.R777W) in two families. I794T has recently been reported as a causative mutation for hereditary diffuse leukoencephalopathy with spheroids (HDLS), and R777W is a novel mutation. Although mutational analysis of CSF1R in 21 sporadic cases revealed no mutations, another novel CSF1R mutation, p.C653Y, was identified in one patient with autopsy‐proven HDSL. These variants were located in the PTK domain where the causative mutations cluster. Functional prediction of the mutant CSF1R as well as cross‐species conservation of the affected amino acids supports the notion that these variants are pathogenic for HDLS. Exome sequencing is useful for a comprehensive mutational analysis of causative genes for hereditary leukoencephalopathies, and CSF1R should be considered a candidate gene for patients with autosomal dominant leukoencephalopathies.

Molecular basis of hereditary methaemoglobinaemia, types I and II: two novel mutations in the NADH-cytochrome b5 reductase gene

Hereditary methaemoglobinaemia, caused by deficiency of NADH‐cytochrome b5 reductase (b5R), has been classified into two types, an erythrocyte (type I) and a generalized (type II). We analysed the b5R gene of two Thai patients and found two novel mutations. The patient with type II was homozygous for a C‐to‐T substitution in codon 83 that changes Arg (CGA) to a stop codon (TGA), resulting in a truncated b5R without the catalytic portion. The patient with type I was homozygous for a C‐to‐T substitution in codon 178 causing replacement of Ala (GCG) with Val (GTG). To characterize effects of this missense mutation, we investigated enzymatic properties of mutant b5R (Ala 178 Val). Although the mutant enzyme showed normal catalytic activity, less stability and different spectra were observed. These results suggest that this substitution influenced enzyme stability due to the slight change of structure. In conclusion, the nonsense mutation led to type II because of malfunction of the truncated protein. On the other hand, the missense mutation caused type I, due to degradation of the unstable mutant enzyme with normal activities in patients erythrocytes, because of the lack of compensation by new protein synthesis during the long life‐span of erythrocytes.