Maternity Log study: a longitudinal lifelog monitoring and multiomics analysis for the early prediction of complicated pregnancy

PURPOSE

A prospective cohort study for pregnant women, the Maternity Log study, was designed to construct a time-course high-resolution reference catalogue of bioinformatic data in pregnancy and explore the associations between genomic and environmental factors and the onset of pregnancy complications, such as hypertensive disorders of pregnancy, gestational diabetes mellitus and preterm labour, using continuous lifestyle monitoring combined with multiomics data on the genome, transcriptome, proteome, metabolome and microbiome.

PARTICIPANTS

Pregnant women were recruited at the timing of first routine antenatal visits at Tohoku University Hospital, Sendai, Japan, between September 2015 and November 2016. Of the eligible women who were invited, 65.4% agreed to participate, and a total of 302 women were enrolled. The inclusion criteria were age ≥20 years and the ability to access the internet using a smartphone in the Japanese language.

FINDINGS TO DATE

Study participants uploaded daily general health information including quality of sleep, condition of bowel movements and the presence of nausea, pain and uterine contractions. Participants also collected physiological data, such as body weight, blood pressure, heart rate and body temperature, using multiple home healthcare devices. The mean upload rate for each lifelog item was ranging from 67.4% (fetal movement) to 85.3% (physical activity), and the total number of data points was over 6 million. Biospecimens, including maternal plasma, serum, urine, saliva, dental plaque and cord blood, were collected for multiomics analysis.

FUTURE PLANS

Lifelog and multiomics data will be used to construct a time-course high-resolution reference catalogue of pregnancy. The reference catalogue will allow us to discover relationships among multidimensional phenotypes and novel risk markers in pregnancy for the future personalised early prediction of pregnancy complications.

POGLUT1, the putative effector gene driven by rs2293370 in primary biliary cholangitis susceptibility locus chromosome 3q13.33

Primary biliary cholangitis (PBC) is a chronic and cholestatic autoimmune liver disease caused by the destruction of intrahepatic small bile ducts. Our previous genome-wide association study (GWAS) identified six susceptibility loci for PBC. Here, in order to further elucidate the genetic architecture of PBC, a GWAS was performed on an additional independent sample set, then a genome-wide meta-analysis with our previous GWAS was performed based on a whole-genome single nucleotide polymorphism (SNP) imputation analysis of a total of 4,045 Japanese individuals (2,060 cases and 1,985 healthy controls). A susceptibility locus on chromosome 3q13.33 (including ARHGAP31, TMEM39A, POGLUT1, TIMMDC1, and CD80) was previously identified both in the European and Chinese populations and was replicated in the Japanese population (OR = 0.7241, P = 3.5 × 10^-9). Subsequent in silico and in vitro functional analyses identified rs2293370, previously reported as the top-hit SNP in this locus in the European population, as the primary functional SNP. Moreover, e-QTL analysis indicated that the effector gene of rs2293370 was Protein O-Glucosyltransferase 1 (POGLUT1) (P = 3.4 × 10^-8). This is the first study to demonstrate that POGLUT1 and not CD80 is the effector gene regulated by the primary functional SNP rs2293370, and that increased expression of POGLUT1 might be involved in the pathogenesis of PBC.

NFKB1 and MANBA Confer Disease Susceptibility to Primary Biliary Cholangitis via Independent Putative Primary Functional Variants

BACKGROUND & AIMS

Primary biliary cholangitis (PBC) is a chronic and cholestatic liver disease that eventually leads to cirrhosis and hepatic failure. We recently identified several susceptibility genes included NFKB1 and MANBA for PBC in the Japanese population by genome-wide association study. However, the primary functional variants in the NFKB1/MANBA region and the molecular mechanism for conferring disease susceptibility to PBC have not yet been clarified.

METHODS

We performed high-density association mapping based on a single-nucleotide polymorphism (SNP) imputation analysis, using data from a whole-genome sequence reference panel of 1070 Japanese individuals and the previous genome-wide association study (1389 PBC patients, 1508 healthy controls). Among SNPs (P < 5.0 × 10^-7) in the NFKB1/MANBA region, putative primary functional variants and the molecular mechanism for conferring disease susceptibility to PBC were identified by in silico/in vitro functional analysis.

RESULTS

Among the SNPs in the NFKB1/MANBA region, rs17032850 and rs227361, which changed the binding of transcription factors lymphoid enhancer-binding factor 1 (LEF-1) and retinoid X receptor α (RXRα), respectively, were identified as putative primary functional variants that regulate gene expression. In addition, expression-quantitative trait locus data and gene editing using a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system supported the potential role of rs17032850 and rs227361 in regulating NFKB1 and MANBA expression, respectively.

CONCLUSIONS

We identified independent putative primary functional variants in NFKB1/MANBA and showed the distinct molecular mechanism by which each putative primary functional variant conferred susceptibility to PBC. Our approach was useful to dissect the pathogenesis not only of PBC, but also other digestive diseases in which NFKB1/MANBA has been reported as a susceptibility locus.

Clinical and genetic risk factors for decreased bone mineral density in Japanese patients with inflammatory bowel disease

BACKGROUND AND AIM

Patients with inflammatory bowel disease (IBD) are at a high risk of low bone mineral density (BMD). Reportedly, clinical and genetic factors cause low BMD in Caucasians; however, studies in non-Caucasian populations remain scarce.

METHODS

Clinical risk factors for low BMD were investigated in 266 Japanese patients with IBD, and a genome-wide association analysis (GWAS) was performed using linear regression with associated clinical factors as covariates. Genotyping was performed using a population-optimized genotyping array (Japonica array^® ). After quality control, the genotype data of 4 384 682 single-nucleotide polymorphisms (SNPs) from 254 patients with IBD were used for GWAS.

RESULTS

Body mass index, age, and disease duration were independently associated with the BMD of the femoral neck (P = 1.41E - 13, 1.04E - 5, and 1.58E - 3, respectively), and body mass index and sex were associated with the BMD of the lumbar spine (P = 6.90E - 10 and 6.84E - 3, respectively). In GWAS, 118 and 42 candidate SNPs of the femoral neck and lumbar spine, respectively, were identified. Among 118, 111 candidate SNPs of the femoral neck were located within the SLC22A23 gene, which is a known IBD susceptibility gene (minimum P = 1.42E - 07). Among 42, 18 candidate SNPs of the lumbar spine were located within the MECOM gene, which is associated with osteopenia (minimum P = 5.86E - 07). Interestingly, none of the known loci showed a significant association with BMD.

CONCLUSIONS

Although clinical risk factors for low BMD in IBD were similar to those in the general population, genetic risk factors were rather different.

Identification of a two-SNP PLA2R1 Haplotype and HLA-DRB1 Alleles as Primary Risk Associations in Idiopathic Membranous Nephropathy

The associations of single nucleotide polymorphisms (SNPs) in PLA2R1 and HLA-DQA1, as well as HLA-DRB1*15:01-DQB1*06:02 haplotype with idiopathic membranous nephropathy (IMN) is well known. However, the primary associations of these loci still need to be determined. We used Japanese-specific SNP genotyping array and imputation using 2,048 sequenced Japanese samples to fine-map PLA2R1 region in 98 patients and 413 controls. The most significant SNPs were replicated in a separate sample set of 130 patients and 288 controls. A two-SNP haplotype of intronic and missense SNPs showed the strongest association. The intronic SNP is strongly associated with PLA2R1 expression in the Genotype-Tissue Expression (GTEx) database, and the missense SNP is predicted to alter peptide binding with HLA-DRB1*15:01 by the Immune Epitope Database (IEDB). In HLA region, we performed relative predispositional effect (RPE) tests and identified additional risk alleles in both HLA-DRB1 and HLA-DQB1. We collapsed the risk alleles in each of HLA-DRB1 and HLA-DQB1 into single risk alleles. Reciprocal conditioning of these collapsed risk alleles showed more residual significance for HLA-DRB1 collapsed risk than HLA-DQB1 collapsed risk. These results indicate that changes in the expression levels of structurally different PLA2R protein confer risk for IMN in the presence of risk HLA-DRB1 alleles.

Functional characterization of 40 CYP2B6 allelic variants by assessing efavirenz 8-hydroxylation

Genetic variations within cytochrome P450 2B6 (CYP2B6) contribute to inter-individual variation in the metabolism of clinically important drugs, including cyclophosphamide, bupropion, methadone and efavirenz (EFZ). In this study, we performed an in vitro analysis of 40 CYP2B6 allelic variant proteins including seven novel variants identified in 1070 Japanese individuals. Wild-type and 39 variant proteins were heterologously expressed in 293FT cells to estimate the kinetic parameters (K_m, V_max, and CL_int) of EFZ 8-hydroxylation and 7-ethoxy-4-trifluoromethylcoumarin (7-ETC) O-deethylation activities. The concentrations of CYP2B6 variant holo-enzymes were measured by using carbon monoxide (CO)-reduced difference spectroscopy, and the wild-type and 28 variants showed a peak at 450 nm. The kinetic parameters were measured for the wild-type and 24 variant proteins. The values for the remaining 15 variants could not be determined because the enzymatic activity was not detected at the highest substrate concentration used. Compared to wild-type, six variants showed significantly decreased EFZ 8-hydroxylation CL_int values, while these values were significantly increased in another six variants, including CYP2B6.6. Although 7-ETC O-deethylation CL_int values of CYP2B6 variants did not differ significantly from that of CYP2B6.1, the CL_int ratios obtained for 7-ETC O-deethylation were highly correlated with EFZ 8-hydroxylation. Furthermore, three-dimensional structural modeling analysis was performed to elucidate the mechanism of changes in the kinetics of CYP2B6 variants. Our findings could provide evidence of the specific metabolic activities of the CYP2B6 proteins encoded by these variant alleles.

Regional genetic differences among Japanese populations and performance of genotype imputation using whole-genome reference panel of the Tohoku Medical Megabank Project

Background

Genotype imputation from single-nucleotide polymorphism (SNP) genotype data using a haplotype reference panel consisting of thousands of unrelated individuals from populations of interest can help to identify strongly associated variants in genome-wide association studies. The Tohoku Medical Megabank (TMM) project was established to support the development of precision medicine, together with the whole-genome sequencing of 1070 human genomes from individuals in the Miyagi region (Northeast Japan) and the construction of the 1070 Japanese genome reference panel (1KJPN). Here, we investigated the performance of 1KJPN for genotype imputation of Japanese samples not included in the TMM project and compared it with other population reference panels.

Results

We found that the 1KJPN population was more similar to other Japanese populations, Nagahama (south-central Japan) and Aki (Shikoku Island), than to East Asian populations in the 1000 Genomes Project other than JPT, suggesting that the large-scale collection (more than 1000) of Japanese genomes from the Miyagi region covered many of the genetic variations of Japanese in mainland Japan. Moreover, 1KJPN outperformed the phase 3 reference panel of the 1000 Genomes Project (1KGPp3) for Japanese samples, and IKJPN showed similar imputation rates for the TMM and other Japanese samples for SNPs with minor allele frequencies (MAFs) higher than 1%.

Conclusions

1KJPN covered most of the variants found in the samples from areas of the Japanese mainland outside the Miyagi region, implying 1KJPN is representative of the Japanese population’s genomes. 1KJPN and successive reference panels are useful genome reference panels for the mainland Japanese population. Importantly, the addition of whole genome sequences not included in the 1KJPN panel improved imputation efficiencies for SNPs with MAFs under 1% for samples from most regions of the Japanese archipelago.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4942-0) contains supplementary material, which is available to authorized users.

Strong Association of the HLA-DR/DQ Locus with Childhood Steroid-Sensitive Nephrotic Syndrome in the Japanese Population

Background Nephrotic syndrome is the most common cause of chronic glomerular disease in children. Most of these patients develop steroid-sensitive nephrotic syndrome (SSNS), but the loci conferring susceptibility to childhood SSNS are mainly unknown.Methods We conducted a genome-wide association study (GWAS) in the Japanese population; 224 patients with childhood SSNS and 419 adult healthy controls were genotyped using the Affymetrix Japonica Array in the discovery stage. Imputation for six HLA genes (HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1) was conducted on the basis of Japanese-specific references. We performed genotyping for HLA-DRB1/-DQB1 using a sequence-specific oligonucleotide-probing method on a Luminex platform. Whole-genome imputation was conducted using a phased reference panel of 2049 healthy Japanese individuals. Replication was performed in an independent Japanese sample set including 216 patients and 719 healthy controls. We genotyped candidate single-nucleotide polymorphisms using the DigiTag2 assay.Results The most significant association was detected in the HLA-DR/DQ region and replicated (rs4642516 [minor allele G], combined P_allelic=7.84×10^-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined P_allelic=1.72×10^-25; OR, 0.29; 95% CI, 0.23 to 0.37). HLA-DRB1*08:02 (Pc=1.82×10^-9; OR, 2.62; 95% CI, 1.94 to 3.54) and HLA-DQB1*06:04 (Pc=2.09×10^-12; OR, 0.10; 95% CI, 0.05 to 0.21) were considered primary HLA alleles associated with childhood SSNS. HLA-DRB1*08:02-DQB1*03:02 (Pc=7.01×10^-11; OR, 3.60; 95% CI, 2.46 to 5.29) was identified as the most significant genetic susceptibility factor.Conclusions The most significant association with childhood SSNS was detected in the HLA-DR/DQ region. Further HLA allele/haplotype analyses should enhance our understanding of molecular mechanisms underlying SSNS.

Omics research project on prospective cohort studies from the Tohoku Medical Megabank Project

Population-based prospective cohort studies are indispensable for modern medical research as they provide important knowledge on the influences of many kinds of genetic and environmental factors on the cause of disease. Although traditional cohort studies are mainly conducted using questionnaires and physical examinations, modern cohort studies incorporate omics and genomic approaches to obtain comprehensive physical information, including genetic information. Here, we report the design and midterm results of multi-omics analysis on population-based prospective cohort studies from the Tohoku Medical Megabank (TMM) Project. We have incorporated genomic and metabolomic studies in the TMM cohort study as both metabolome and genome analyses are suitable for high-throughput analysis of large-scale cohort samples. Moreover, an association study between the metabolome and genome show that metabolites are an important intermediate phenotype connecting genetic and lifestyle factors to physical and pathologic phenotypes. We apply our metabolome and genome analyses to large-scale cohort samples in the following studies.

Functional Characterization of 21 Allelic Variants of Dihydropyrimidine Dehydrogenase Identified in 1070 Japanese Individuals

Dihydropyrimidine dehydrogenase (DPD, EC 1.3.1.2), encoded by the DPYD gene, is the rate-limiting enzyme in the degradation pathway of endogenous pyrimidine and fluoropyrimidine drugs such as 5-fluorouracil (5-FU). DPD catalyzes the reduction of uracil, thymine, and 5-FU. In Caucasians, DPYD mutations, including DPYD*2A, DPYD*13, c.2846A>T, and c.1129-5923C>G/hapB3, are known to contribute to interindividual variations in the toxicity of 5-FU; however, none of these DPYD polymorphisms has been identified in the Asian population. Recently, 21 DPYD allelic variants, including some novel single-nucleotide variants (SNVs), were identified in 1070 healthy Japanese individuals by analyzing their whole-genome sequences (WGSs), but the functional alterations caused by these variants remain unknown. In this study, in vitro analysis was performed on 22 DPD allelic variants by transiently expressing wild-type DPD and 21 DPD variants in 293FT cells and characterizing their enzymatic activities using 5-FU as a substrate. DPD expression levels and dimeric forms were determined using immunoblotting and blue-native PAGE, respectively. Additionally, the values of three kinetic parameters-the Michaelis constant (K_m ), maximum velocity (V_max ), and intrinsic clearance (CL_int = V_max/K_m )-were determined for the reduction of 5-FU. Eleven variants exhibited significantly decreased intrinsic clearance compared with wild-type DPD. Moreover, the band patterns observed in the immunoblots of blue-native gels indicated that DPD dimerization is required for enzymatic activity in DPD. Thus, the detection of rare DPYD variants might facilitate severe adverse effect prediction of 5-FU-based chemotherapy in the Japanese population.

Population-scale whole genome sequencing identifies 271 highly polymorphic short tandem repeats from Japanese population

Forensic DNA typing is widely used to identify missing persons and plays a central role in forensic profiling. DNA typing usually uses capillary electrophoresis fragment analysis of PCR amplification products to detect the length of short tandem repeat (STR) markers. Here, we analyzed whole genome data from 1,070 Japanese individuals generated using massively parallel short-read sequencing of 162 paired-end bases. We have analyzed 843,473 STR loci with two to six basepair repeat units and cataloged highly polymorphic STR loci in the Japanese population. To evaluate the performance of the cataloged STR loci, we compared 23 STR loci, widely used in forensic DNA typing, with capillary electrophoresis based STR genotyping results in the Japanese population. Seventeen loci had high correlations and high call rates. The other six loci had low call rates or low correlations due to either the limitations of short-read sequencing technology, the bioinformatics tool used, or the complexity of repeat patterns. With these analyses, we have also purified the suitable 218 STR loci with four basepair repeat units and 53 loci with five basepair repeat units both for short read sequencing and PCR based technologies, which would be candidates to the actual forensic DNA typing in Japanese population.

Identification of somatic mutations in postmortem human brains by whole genome sequencing and their implications for psychiatric disorders

AIM

Somatic mutations in the human brain are hypothesized to contribute to the functional diversity of brain cells as well as the pathophysiology of neuropsychiatric diseases. However, there are still few reports on somatic mutations in non-neoplastic human brain tissues. This study attempted to unveil the landscape of somatic mutations in the human brain.

METHODS

We explored the landscape of somatic mutations in human brain tissues derived from three individuals with no neuropsychiatric diseases by whole-genome deep sequencing at a depth of around 100. The candidate mutations underwent multi-layered filtering, and were validated by ultra-deep target amplicon sequencing at a depth of around 200 000.

RESULTS

Thirty-one somatic mutations were identified in the human brain, demonstrating the utility of whole-genome sequencing of bulk brain tissue. The mutations were enriched in neuron-expressed genes, and two-thirds of the identified somatic single nucleotide variants in the brain tissues were cytosine-to-thymine transitions, half of which were in CpG dinucleotides.

CONCLUSION

Our developed filtering and validation approaches will be useful to identify somatic mutations in the human brain. The vulnerability of neuron-expressed genes to mutational events suggests their potential relevance to neuropsychiatric diseases.

Allele-specific DNA methylation of disease susceptibility genes in Japanese patients with inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) has an unknown etiology; however, accumulating evidence suggests that IBD is a multifactorial disease influenced by a combination of genetic and environmental factors. The influence of genetic variants on DNA methylation in cis and cis effects on expression have been demonstrated. We hypothesized that IBD susceptibility single-nucleotide polymorphisms (SNPs) regulate susceptibility gene expressions in cis by regulating DNA methylation around SNPs. For this, we determined cis-regulated allele-specific DNA methylation (ASM) around IBD susceptibility genes in CD4+ effector/memory T cells (Tem) in lamina propria mononuclear cells (LPMCs) in patients with IBD and examined the association between the ASM SNP genotype and neighboring susceptibility gene expressions.

METHODS

CD4+ effector/memory T cells (Tem) were isolated from LPMCs in 15 Japanese IBD patients (ten Crohn's disease [CD] and five ulcerative colitis [UC] patients). ASM analysis was performed by methylation-sensitive SNP array analysis. We defined ASM as a changing average relative allele score ([Formula: see text]) >0.1 after digestion by methylation-sensitive restriction enzymes. Among SNPs showing [Formula: see text] >0.1, we extracted the probes located on tag-SNPs of 200 IBD susceptibility loci and around IBD susceptibility genes as candidate ASM SNPs. To validate ASM, bisulfite-pyrosequencing was performed. Transcriptome analysis was examined in 11 IBD patients (seven CD and four UC patients). The relation between rs36221701 genotype and neighboring gene expressions were analyzed.

RESULTS

We extracted six candidate ASM SNPs around IBD susceptibility genes. The top of [Formula: see text] (0.23) was rs1130368 located on HLA-DQB1. ASM around rs36221701 ([Formula: see text] = 0.14) located near SMAD3 was validated using bisulfite pyrosequencing. The SMAD3 expression was significantly associated with the rs36221701 genotype (p = 0.016).

CONCLUSIONS

We confirmed the existence of cis-regulated ASM around IBD susceptibility genes and the association between ASM SNP (rs36221701) genotype and SMAD3 expression, a susceptibility gene for IBD. These results give us supporting evidence that DNA methylation mediates genetic effects on disease susceptibility.

NUDT15 codon 139 is the best pharmacogenetic marker for predicting thiopurine-induced severe adverse events in Japanese patients with inflammatory bowel disease: a multicenter study

BACKGROUND

Despite NUDT15 variants showing significant association with thiopurine-induced adverse events (AEs) in Asians, it remains unclear which variants of NUDT15 or whether additional genetic variants should be tested to predict AEs. To clarify the best pharmacogenetic test to be used clinically, we performed association studies of NUDT15 variants and haplotypes with AEs, genome-wide association study (GWAS) to discover additional variants, and ROC analysis to select the model to predict severe AEs.

METHODS

Overall, 2630 patients with inflammatory bowel disease (IBD) were enrolled and genotyped for NUDT15 codon 139; 1291 patients were treated with thiopurines. diplotypes were analyzed in 970 patients, and GWASs of AEs were performed with 1221 patients using population-optimized genotyping array and imputation.

RESULTS

We confirmed the association of NUDT15 p.Arg139Cys with leukopenia and alopecia (p = 2.20E-63, 1.32E-69, OR = 6.59, 12.1, respectively), and found a novel association with digestive symptoms (p = 6.39E-04, OR = 1.89). Time to leukopenia was significantly shorter, and when leukopenia was diagnosed, thiopurine doses were significantly lower in Arg/Cys and Cys/Cys than in Arg/Arg. In GWASs, no additional variants were found to be associated with thiopurine-induced AEs. Despite strong correlation of leukopenia frequency with estimated enzyme activities based on the diplotypes (r² = 0.926, p = 0.0087), there were no significant differences in the AUCs of diplotypes from those of codon 139 to predict severe AEs (AUC = 0.916, 0.921, for acute severe leukopenia, AUC = 0.990, 0.991, for severe alopecia, respectively).

CONCLUSIONS

Genotyping of NUDT15 codon 139 was sufficient to predict acute severe leukopenia and alopecia in Japanese patients with IBD.

Genetic analysis of Japanese primary open-angle glaucoma patients and clinical characterization of risk alleles near CDKN2B-AS1, SIX6 and GAS7

PURPOSE

To test the genetic association between Japanese patients with primary open-angle glaucoma (POAG) and the previously reported POAG susceptibility loci and to perform genotype-phenotype analysis.

METHODS

Genetic associations for 27 SNPs from 16 loci previously linked to POAG were assessed using genome-wide SNP data of the primary cohort (565 Japanese POAG patients and 1,104 controls). Reproducibility of the assessment was tested in 607 POAG cases and 455 controls (second cohort) with a targeted genotyping approach. For POAG-associated variants, a genotype-phenotype correlation study (additive, dominant, recessive model) was performed using the objective clinical data derived from 598 eyes of 598 POAG patients.

RESULTS

Among 27 SNPs from 16 loci previously linked to POAG, genotypes for total of 20 SNPs in 13 loci were available for targeted association study. Among 8 SNPs in 3 loci that showed at least nominal association (P < 5.00E-02) in the primary cohort, a representative SNP for each loci (rs2157719 for CDKN2B-AS1, rs33912345 for SIX6, and rs9913911 for GAS7) were selected. For these SNPs the association was found significant in both the second cohort analysis and meta-analysis. The genotype-phenotype analysis revealed significant correlations between CDKN2B-AS1 (rs2157719) and decreased intraocular pressure (β = -6.89 mmHg, P = 1.70E-04; dominant model) after multiple corrections. In addition, nominal correlation was observed between CDKN2B-AS1 (rs2157719) and optic nerve head blood flow (β = -0.54 and -0.67 arbitrary units (AU), P = 2.00E-02 and 1.39E-02), between SIX6 (rs33912345) and decreased total peripapillary retinal nerve fiber layer thickness (β = -2.16 and -2.82 μm, P = 4.68E-02 and 2.40E-02, additive and recessive model, respectively) and increased optic nerve head blood flow (β = 0.44 AU, P = 2.20E-02; additive model) and between GAS7 (rs9913911) and increased cup volume (β = 0.03 mm3, P = 4.60E-02) and mean cup depth (β = 0.03 mm3, P = 4.11E-02; additive model) and decreased pattern standard deviation (β = -0.87 dB, P = 2.44E-02; dominant model).

CONCLUSION

The association between SNPs near GAS7 and POAG was found in Japanese patients for the first time. Clinical characterization of the risk variants is an important step toward understanding the pathology of the disease and optimizing treatment of patients with POAG.

Genome-wide association studies identify PRKCB as a novel genetic susceptibility locus for primary biliary cholangitis in the Japanese population

A previous genome-wide association study (GWAS) performed in 963 Japanese individuals (487 primary biliary cholangitis [PBC] cases and 476 healthy controls) identified TNFSF15 (rs4979462) and POU2AF1 (rs4938534) as strong susceptibility loci for PBC. In this study, we performed GWAS in additional 1,923 Japanese individuals (894 PBC cases and 1,029 healthy controls), and combined the results with the previous data. This GWAS, together with a subsequent replication study in an independent set of 7,024 Japanese individuals (512 PBC cases and 6,512 healthy controls), identified PRKCB (rs7404928) as a novel susceptibility locus for PBC (odds ratio [OR] = 1.26, P = 4.13 × 10-9). Furthermore, a primary functional variant of PRKCB (rs35015313) was identified by genotype imputation using a phased panel of 1,070 Japanese individuals from a prospective, general population cohort study and subsequent in vitro functional analyses. These results may lead to improved understanding of the disease pathways involved in PBC, forming a basis for prevention of PBC and development of novel therapeutics.

A Histologic Categorization of Aqueous Outflow Routes in Familial Open-Angle Glaucoma and Associations With Mutations in the MYOC Gene in Japanese Patients

Purpose

This study evaluated specific relationships between pathogenic mechanisms and genetic polymorphisms in primary open-angle glaucoma (POAG). We analyzed the morphologies of trabeculectomy specimens obtained from patients with familial POAG.

Methods

We used light microscopy and transmission electron microscopy to examine specimens obtained from 17 eyes of 14 patients with familial POAG. We also conducted exome analyses of two families and used targeted Sanger sequencing to analyze samples obtained from the remaining patients.

Results

The POAG cases examined in this study were divided into two groups based on morphologic characteristics. Group A eyes (7 eyes from 5 patients) had an abnormally thick trabecular meshwork (TM), whereas group B eyes (10 eyes from 9 patients) had a TM of normal thickness. The characteristics of the outflow routes in group A eyes were remarkable and included apoptotic TM cells, abnormally thickened TM basement membranes, fused TM beams, and occluded Schlemm's canals. All group A patients harbored mutations (F369L, P370L, T377M, and T448P) in the myocilin (MYOC) gene that were not found in group B patients.

Conclusions

Although age matching of morphologic changes in the outflow routes was impossible due to the small sample size, this study suggests that abnormal TM cells may cause sequential damage in abnormally thickened TM basement membranes, TM cell apoptosis, TM beam fusion, and the occlusion of Schlemm's canals. The four detected MYOC mutations appeared to be associated with morphologic changes in the TM and the underlying pathogenesis of a subtype of familial POAG.

Security controls in an integrated Biobank to protect privacy in data sharing: rationale and study design

Background

With the goal of realizing genome-based personalized healthcare, we have developed a biobank that integrates personal health, genome, and omics data along with biospecimens donated by volunteers of 150,000. Such a large-scale of data integration involves obvious risks of privacy violation. The research use of personal genome and health information is a topic of global discussion with regard to the protection of privacy while promoting scientific advancement. The present paper reports on our plans, current attempts, and accomplishments in addressing security problems involved in data sharing to ensure donor privacy while promoting scientific advancement.

Methods

Biospecimens and data have been collected in prospective cohort studies with the comprehensive agreement. The sample size of 150,000 participants was required for multiple researches including genome-wide screening of gene by environment interactions, haplotype phasing, and parametric linkage analysis.

Results

We established the TohokuMedicalMegabank (TMM) data sharing policy: a privacy protection rule that requires physical, personnel, and technological safeguards against privacy violation regarding the use and sharing of data. The proposed policy refers to that of NCBI and that of the Sanger Institute. The proposed policy classifies shared data according to the strength of re-identification risks. Local committees organized by TMM evaluate re-identification risk and assign a security category to a dataset. Every dataset is stored in an assigned segment of a supercomputer in accordance with its security category. A security manager should be designated to handle all security problems at individual data use locations. The proposed policy requires closed networks and IP-VPN remote connections.

Conclusion

The mission of the biobank is to distribute biological resources most productively. This mission motivated us to collect biospecimens and health data and simultaneously analyze genome/omics data in-house. The biobank also has the mission of improving the quality and quantity of the contents of the biobank. This motivated us to request users to share the results of their research as feedback to the biobank. The TMM data sharing policy has tackled every security problem originating with the missions. We believe our current implementation to be the best way to protect privacy in data sharing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12911-017-0494-5) contains supplementary material, which is available to authorized users.

Genome-Wide Association Study Identifies TLL1 Variant Associated With Development of Hepatocellular Carcinoma After Eradication of Hepatitis C Virus Infection

BACKGROUND & AIMS

There is still a risk for hepatocellular carcinoma (HCC) development after eradication of hepatitis C virus (HCV) infection with antiviral agents. We investigated genetic factors associated with the development of HCC in patients with a sustained virologic response (SVR) to treatment for chronic HCV infection.

METHODS

We obtained genomic DNA from 457 patients in Japan with a SVR to interferon-based treatment for chronic HCV infection from 2007 through 2015. We conducted a genome-wide association study (GWAS), followed by a replication analysis of 79 candidate single nucleotide polymorphisms (SNPs) in an independent set of 486 patients in Japan. The study end point was HCC diagnosis or confirmation of lack of HCC (at follow-up examinations until December 2014 in the GWAS cohort, and until January 2016 in the replication cohort). We collected clinical and laboratory data from all patients. We analyzed expression levels of candidate gene variants in human hepatic stellate cells, rats with steatohepatitis caused by a choline-deficient L-amino acid-defined diet, and a mouse model of liver injury caused by administration of carbon tetrachloride. We also analyzed expression levels in liver tissues of patients with chronic HCV infection with different stages of fibrosis or tumors vs patients without HCV infection (controls).

RESULTS

We found a strong association between the SNP rs17047200, located within the intron of the tolloid like 1 gene (TLL1) on chromosome 4, and development of HCC; there was a genome-wide level of significance when the results of the GWAS and replication study were combined (odds ratio, 2.37; P = 2.66 × 10^-8). Multivariate analysis showed rs17047200 AT/TT to be an independent risk factor for HCC (hazard ratio, 1.78; P = .008), along with male sex, older age, lower level of albumin, advanced stage of hepatic fibrosis, presence of diabetes, and higher post-treatment level of α-fetoprotein. Combining the rs17047200 genotype with other factors, we developed prediction models for HCC development in patients with mild or advanced hepatic fibrosis. Levels of TLL1 messenger RNA (mRNA) in human hepatic stellate cells increased with activation. Levels of Tll1 mRNA increased in liver tissues of rodents with hepatic fibrogenesis compared with controls. Levels of TLL1 mRNA increased in liver tissues of patients with progression of fibrosis. Gene expression levels of TLL1 short variants, including isoform 2, were higher in patients with rs17047200 AT/TT.

CONCLUSIONS

In a GWAS, we identified the association between the SNP rs17047200, within the intron of TLL1, and development of HCC in patients who achieved an SVR to treatment for chronic HCV infection. We found levels of Tll1/TLL1 mRNA to be increased in rodent models of liver injury and liver tissues of patients with fibrosis, compared with controls. We propose that this SNP might affect splicing of TLL1 mRNA, yielding short variants with high catalytic activity that accelerates hepatic fibrogenesis and carcinogenesis. Further studies are needed to determine how rs17047200 affects TLL1 mRNA levels, splicing, and translation, as well as the prevalence of this variant among other patients with HCC. Tests for the TLL1 SNP might be used to identify patients at risk for HCC after an SVR to treatment of HCV infection.

LRRK2 but not ATG16L1 is associated with Paneth cell defect in Japanese Crohn's disease patients

BACKGROUND. Morphological patterns of Paneth cells are a prognostic biomarker in Western Crohn's disease (CD) patients, and are associated with autophagy-associated ATG16L1 and NOD2 variants. We hypothesized that genetic determinants of Paneth cell phenotype in other ethnic CD cohorts are distinct but also involved in autophagy. METHODS. We performed a hypothesis-driven analysis of 56 single nucleotide polymorphisms (SNPs) associated with CD susceptibility or known to affect Paneth cell function in 110 Japanese CD patients who underwent ileal resection. We subsequently performed a genome-wide association analysis. Paneth cell phenotype was determined by defensin-5 immunofluorescence. Selected genotype-Paneth cell defect correlations were compared to a Western CD cohort (n = 164). RESULTS. The average percentage of abnormal Paneth cells in Japanese CD was similar to Western CD (P = 0.87), and abnormal Paneth cell phenotype was also associated with early recurrence (P = 0.013). In contrast to Western CD, ATG16L1 T300A was not associated with Paneth cell defect in Japanese CD (P = 0.20). Among the 56 selected SNPs, only LRRK2 M2397T showed significant association with Paneth cell defect (P = 3.62 × 10^-4), whereas in the Western CD cohort it was not (P = 0.76). Pathway analysis of LRRK2 and other candidate genes with P less than 5 × 10^-4 showed connections with known CD susceptibility genes and links to autophagy and TNF-α networks. CONCLUSIONS. We found dichotomous effects of ATG16L1 and LRRK2 on Paneth cell defect between Japanese and Western CD. Genes affecting Paneth cell phenotype in Japanese CD were also associated with autophagy. Paneth cell phenotype also predicted prognosis in Japanese CD. FUNDING. Helmsley Charitable Trust, Doris Duke Foundation (grant 2014103), Japan Society for the Promotion of Science (KAKENHI grants JP15H04805 and JP15K15284), Crohn's and Colitis Foundation grant 274415, NIH (grants 1R56DK095820, K01DK109081, and UL1 TR000448).

Genome-wide association study using the ethnicity-specific Japonica array: identification of new susceptibility loci for cold medicine-related Stevens-Johnson syndrome with severe ocular complications

A genome-wide association study (GWAS) for cold medicine-related Stevens-Johnson syndrome (CM-SJS) with severe ocular complications (SOC) was performed in a Japanese population. A recently developed ethnicity-specific array with genome-wide imputation that was based on the whole-genome sequences of 1070 unrelated Japanese individuals was used. Validation analysis with additional samples from Japanese individuals and replication analysis using samples from Korean individuals identified two new susceptibility loci on chromosomes 15 and 16. This study might suggest the usefulness of GWAS using the ethnicity-specific array and genome-wide imputation based on large-scale whole-genome sequences. Our findings contribute to the understanding of genetic predisposition to CM-SJS with SOC.

STR-realigner: a realignment method for short tandem repeat regions

Background

In the estimation of repeat numbers in a short tandem repeat (STR) region from high-throughput sequencing data, two types of strategies are mainly taken: a strategy based on counting repeat patterns included in sequence reads spanning the region and a strategy based on estimating the difference between the actual insert size and the insert size inferred from paired-end reads. The quality of sequence alignment is crucial, especially in the former approaches although usual alignment methods have difficulty in STR regions due to insertions and deletions caused by the variations of repeat numbers.

Results

We proposed a new dynamic programming based realignment method named STR-realigner that considers repeat patterns in STR regions as prior knowledge. By allowing the size change of repeat patterns with low penalty in STR regions, accurate realignment is expected. For the performance evaluation, publicly available STR variant calling tools were applied to three types of aligned reads: synthetically generated sequencing reads aligned with BWA-MEM, those realigned with STR-realigner, those realigned with ReviSTER, and those realigned with GATK IndelRealigner. From the comparison of root mean squared errors between estimated and true STR region size, the results for the dataset realigned with STR-realigner are better than those for other cases. For real data analysis, we used a real sequencing dataset from Illumina HiSeq 2000 for a parent-offspring trio. RepeatSeq and lobSTR were applied to the sequence reads for these individuals aligned with BWA-MEM, those realigned with STR-realigner, ReviSTER, and GATK IndelRealigner. STR-realigner shows the best performance in terms of consistency of the size of estimated STR regions in Mendelian inheritance. Root mean squared error values were also calculated from the comparison of these estimated results with STR region sizes obtained from high coverage PacBio sequencing data, and the results from the realigned sequencing data with STR-realigner showed the least (the best) root mean squared error value.

Conclusions

The effectiveness of the proposed realignment method for STR regions was verified from the comparison with an existing method on both simulation datasets and real whole genome sequencing dataset.

The structural origin of metabolic quantitative diversity

Relationship between structural variants of enzymes and metabolic phenotypes in human population was investigated based on the association study of metabolite quantitative traits with whole genome sequence data for 512 individuals from a population cohort. We identified five significant associations between metabolites and non-synonymous variants. Four of these non-synonymous variants are located in enzymes involved in metabolic disorders, and structural analyses of these moderate non-synonymous variants demonstrate that they are located in peripheral regions of the catalytic sites or related regulatory domains. In contrast, two individuals with larger changes of metabolite levels were also identified, and these individuals retained rare variants, which caused non-synonymous variants located near the catalytic site. These results are the first demonstrations that variant frequency, structural location, and effect for phenotype correlate with each other in human population, and imply that metabolic individuality and susceptibility for diseases may be elicited from the moderate variants and much more deleterious but rare variants.

Conditional rod photoreceptor ablation reveals Sall1 as a microglial marker and regulator of microglial morphology in the retina

Neurodegeneration has been shown to induce microglial activation and the infiltration of monocyte-derived macrophages into the CNS, resulting in the coexistence of these two populations within the same lesion, though their distinct features remain elusive. To investigate the impact of rod photoreceptor degeneration on microglial activation, we generated a toxin-mediated genetic model of rod degeneration. Rod injury induced microglial proliferation and migration toward the photoreceptors. Bone marrow transplantation revealed the invasion of monocyte-derived macrophages into the retina, with microglia and the infiltrating macrophages showing distinct distribution patterns in the retina. By comparing the gene expression profiles of the activated microglia and infiltrating macrophages, we identified microglia-specific genes, including Ak1, Ctsf, Sall1, Phlda3, and Spns2. An analysis of Sall1gfp knock-in mice showed GFP expression in the microglia of developing and mature healthy retinas. DTA injury induced the expansion of Sall1gfp(+) microglia, whereas Ly6C(+) monocyte-derived macrophages were mostly Sall1gfp(-) , supporting the idea that Sall1 is exclusively expressed in microglia within the retinal phagocyte pool. We evaluated the contribution of microglia to the phagocyte pool in rd1 mutant retinas and found that Sall1gfp(+) microglia constituted the majority of phagocytes. A Sall1 deficiency did not affect microglial colonization of the retina and the cortex, but it did change their morphology from a ramified to a more amoeboid appearance. The morphological defects observed in Sall1-deficient microglia were not rescued by the presence of wild-type non-microglial cells, suggesting that Sall1 functions cell-autonomously in microglia. Taken together, our data indicate that Sall1 regulates microglial morphology during development. GLIA 2016;64:2005-2024.

The Tohoku Medical Megabank Project: Design and Mission

The Great East Japan Earthquake (GEJE) and resulting tsunami of March 11, 2011 gave rise to devastating damage on the Pacific coast of the Tohoku region. The Tohoku Medical Megabank Project (TMM), which is being conducted by Tohoku University Tohoku Medical Megabank Organization (ToMMo) and Iwate Medical University Iwate Tohoku Medical Megabank Organization (IMM), has been launched to realize creative reconstruction and to solve medical problems in the aftermath of this disaster. We started two prospective cohort studies in Miyagi and Iwate Prefectures: a population-based adult cohort study, the TMM Community-Based Cohort Study (TMM CommCohort Study), which will recruit 80 000 participants, and a birth and three-generation cohort study, the TMM Birth and Three-Generation Cohort Study (TMM BirThree Cohort Study), which will recruit 70 000 participants, including fetuses and their parents, siblings, grandparents, and extended family members. The TMM CommCohort Study will recruit participants from 2013 to 2016 and follow them for at least 5 years. The TMM BirThree Cohort Study will recruit participants from 2013 to 2017 and follow them for at least 4 years. For children, the ToMMo Child Health Study, which adopted a cross-sectional design, was also started in November 2012 in Miyagi Prefecture. An integrated biobank will be constructed based on the two prospective cohort studies, and ToMMo and IMM will investigate the chronic medical impacts of the GEJE. The integrated biobank of TMM consists of health and clinical information, biospecimens, and genome and omics data. The biobank aims to establish a firm basis for personalized healthcare and medicine, mainly for diseases aggravated by the GEJE in the two prefectures. Biospecimens and related information in the biobank will be distributed to the research community. TMM itself will also undertake genomic and omics research. The aims of the genomic studies are: 1) to construct an integrated biobank; 2) to return genomic research results to the participants of the cohort studies, which will lead to the implementation of personalized healthcare and medicine in the affected areas in the near future; and 3) to contribute the development of personalized healthcare and medicine worldwide. Through the activities of TMM, we will clarify how to approach prolonged healthcare problems in areas damaged by large-scale disasters and how useful genomic information is for disease prevention.

Variants in the UBR1 gene are not associated with chronic pancreatitis in Japan

BACKGROUND/OBJECTIVES

The UBR1 gene encodes the enzyme ubiquitin-protein ligase E3 component n-recognin 1. Loss-of-function mutations in the UBR1 gene cause Johanson-Blizzard syndrome, which involves pancreatic exocrine insufficiency. No previous studies have examined an association of UBR1 variants with pancreatitis, in part due to the large size of the gene. This study aimed to clarify whether UBR1 variants are associated with chronic pancreatitis (CP) by the application of targeted next generation sequencing.

METHODS

Exon sequences of the UBR1 gene from 389 Japanese patients with CP (188 idiopathic, 172 alcoholic, 20 hereditary, 9 familial) were captured by the HaloPlex target enrichment technology and subjected to next generation sequencing.

RESULTS

Ninety nine point two % of the coding regions of the UBR1 gene could be sequenced by ≥ 20 reads with a mean read depth of 595 and a median depth of 399. Fifteen non-synonymous variants including three novel ones [c.4514T > C (p.I1505T), c.4828C > G (p.H1610D) and c.4856A > T (p.D1619V)] and two synonymous variants were identified in the exonic regions. The frequency of any non-synonymous or synonymous variants was not different between the patients with CP and controls.

CONCLUSIONS

Variants in the UBR1 gene were not associated with CP in Japan.

A Bayesian approach for estimating allele-specific expression from RNA-Seq data with diploid genomes

Background

RNA-sequencing (RNA-Seq) has become a popular tool for transcriptome profiling in mammals. However, accurate estimation of allele-specific expression (ASE) based on alignments of reads to the reference genome is challenging, because it contains only one allele on a mosaic haploid genome. Even with the information of diploid genome sequences, precise alignment of reads to the correct allele is difficult because of the high-similarity between the corresponding allele sequences.

Results

We propose a Bayesian approach to estimate ASE from RNA-Seq data with diploid genome sequences. In the statistical framework, the haploid choice is modeled as a hidden variable and estimated simultaneously with isoform expression levels by variational Bayesian inference. Through the simulation data analysis, we demonstrate the effectiveness of the proposed approach in terms of identifying ASE compared to the existing approach. We also show that our approach enables better quantification of isoform expression levels compared to the existing methods, TIGAR2, RSEM and Cufflinks. In the real data analysis of the human reference lymphoblastoid cell line GM12878, some autosomal genes were identified as ASE genes, and skewed paternal X-chromosome inactivation in GM12878 was identified.

Conclusions

The proposed method, called ASE-TIGAR, enables accurate estimation of gene expression from RNA-Seq data in an allele-specific manner. Our results show the effectiveness of utilizing personal genomic information for accurate estimation of ASE. An implementation of our method is available at http://nagasakilab.csml.org/ase-tigar.

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.

Mitochondria play a crucial role in ATP biosynthesis and comprise proteins encoded in both the nuclear and mitochondrial genomes. Although more than 250 mitochondrial disease-causing genes have been reported, the exact genetic causes in patients remain largely unknown. Here, we aimed to provide further insights into the pathogenic mechanisms of mitochondrial disorders. We investigated the genes encoded in the nuclear and mitochondrial genomes using comprehensive genomic analysis in 142 patients with mitochondrial respiratory chain complex deficiencies. We identified 3 novel disease-causing mitochondria-related genes (MRPS23, QRSL1, and PNPLA4) as well as other disease-causing genes and novel pathogenic mutations in known mitochondrial disease-causing genes. All pathogenic mutations in this study are validated by genetic and/or functional evidence. Our findings, including the achievement of firm genetic diagnoses for 49 of 142 patients (34.5%), were higher than the general diagnosis rate of approximately 25% and demonstrated the value of comprehensive genomic analysis. Accordingly, we have shed light on the genetic heterogeneity underlying mitochondrial disorders.

iJGVD: an integrative Japanese genome variation database based on whole-genome sequencing

The integrative Japanese Genome Variation Database (iJGVD; http://ijgvd.megabank.tohoku.ac.jp/) provides genomic variation data detected by whole-genome sequencing (WGS) of Japanese individuals. Specifically, the database contains variants detected by WGS of 1,070 individuals who participated in a genome cohort study of the Tohoku Medical Megabank Project. In the first release, iJGVD includes >4,300,000 autosomal single nucleotide variants (SNVs) whose minor allele frequencies are >5.0%.

Rare variant discovery by deep whole-genome sequencing of 1,070 Japanese individuals

The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of <1.0%. This detailed analysis detected signatures for purifying selection on regulatory elements as well as coding regions. We also catalogue structural variants, including 3.4 million insertions and deletions, and 25,923 genic copy-number variants. The 1KJPN was effective for imputing genotypes of the Japanese population genome wide. These data demonstrate the value of high-coverage sequencing for constructing population-specific variant panels, which covers 99.0% SNVs of minor allele frequency ≥0.1%, and its value for identifying causal rare variants of complex human disease phenotypes in genetic association studies.

The Tohoku Medical Megabank Organization establishes a biobank with detailed patient health care and genome information. Here the authors analyse whole-genome sequences of 1,070 Japanese individuals, allowing them to catalogue 21 million single-nucleotide variants including 12 million novel ones.

Somatic alteration and depleted nuclear expression of BAP1 in human esophageal squamous cell carcinoma

BRCA1-associated protein 1 (BAP1) is a deubiquitinating enzyme that is involved in the regulation of cell growth. Recently, many somatic and germline mutations of BAP1 have been reported in a broad spectrum of tumors. In this study, we identified a novel somatic non-synonymous BAP1 mutation, a phenylalanine-to-isoleucine substitution at codon 170 (F170I), in 1 of 49 patients with esophageal squamous cell carcinoma (ESCC). Multiplex ligation-dependent probe amplification (MLPA) of BAP1 gene in this ESCC tumor disclosed monoallelic deletion (LOH), suggesting BAP1 alterations on both alleles in this tumor. The deubiquitinase activity and the auto-deubiquitinase activity of F170I-mutant BAP1 were markedly suppressed compared with wild-type BAP1. In addition, wild-type BAP1 mostly localizes to the nucleus, whereas the F170I mutant preferentially localized in the cytoplasm. Microarray analysis revealed that expression of the F170I mutant drastically altered gene expression profiles compared with expressed wild-type BAP1. Gene-ontology analyses indicated that the F170I mutation altered the expression of genes involved in oncogenic pathways. We found that one candidate, TCEAL7, previously reported as a putative tumor suppressor gene, was significantly induced by wild-type BAP1 as compared to F170I mutant BAP1. Furthermore, we found that the level of BAP1 expression in the nucleus was reduced in 44% of ESCC examined by immunohistochemistry (IHC). Because the nuclear localization of BAP1 is important for its tumor suppressor function, BAP1 may be functionally inactivated in a substantial portion of ESCC. Taken together, BAP1 is likely to function as a tumor suppressor in at least a part of ESCC.

Abstract 4955: Somatic alteration and depleted nuclear expression of BAP1 in human esophageal squamous cell carcinomas

We have found a case of somatic non-synonymous mutation of BAP1, substituting phenylalanine to isoleucine at codon 170 (F170I) out of 49 patients with esophageal squamous cell carcinomas (ESCs) examined by Sanger's direct sequencing. By in vivo Ubiquitination assay, deubiquitininase activity was suppressed in F170I, comparing with wild BAP1. BAP1 localization assay revealed that F170I mutant preferred to localize within cytoplasm although wild BAP1 mostly localizes within nucleus. Microarray analysis revealed that F170I has drastically altered gene expression profile. The following gene-ontology analyses indicated that expression of genes involved in oncogenic pathways has drastically altered by F170I. We have also found that BAP1 expression was depressed around 10% of ESCs examined by immunohistochemistry. It is interesting that nuclear BAP1-depleted cases have strong cytosolic expression of BAP1, instead. As nuclear localization of BAP1 is important in tumor suppressor function, BAP1 may be functionally inactivated around 10% of ESCs.

Citation Format: Takahiro Mori, Makiko Sumii, Fumiyoshi Fujishima, Kazuko Ueno, Masao Nagasaki, Chikashi Ishioka, Natsuko Chiba. Somatic alteration and depleted nuclear expression of BAP1 in human esophageal squamous cell carcinomas. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4955. doi:10.1158/1538-7445.AM2015-4955

Inter-Individual Differences in the Oral Bacteriome Are Greater than Intra-Day Fluctuations in Individuals

Given the advent of massively parallel DNA sequencing, human microbiome is analyzed comprehensively by metagenomic approaches. However, the inter- and intra-individual variability and stability of the human microbiome remain poorly characterized, particularly at the intra-day level. This issue is of crucial importance for studies examining the effects of microbiome on human health. Here, we focused on bacteriome of oral plaques, for which repeated, time-controlled sampling is feasible. Eighty-one supragingival plaque subjects were collected from healthy individuals, examining multiple sites within the mouth at three time points (forenoon, evening, and night) over the course of 3 days. Bacterial composition was estimated by 16S rRNA sequencing and species-level profiling, resulting in identification of a total of 162 known bacterial species. We found that species compositions and their relative abundances were similar within individuals, and not between sampling time or tooth type. This suggests that species-level oral bacterial composition differs significantly between individuals, although the number of subjects is limited and the intra-individual variation also occurs. The majority of detected bacterial species (98.2%; 159/162), however, did not fluctuate over the course of the day, implying a largely stable oral microbiome on an intra-day time scale. In fact, the stability of this data set enabled us to estimate potential interactions between rare bacteria, with 40 co-occurrences supported by the existing literature. In summary, the present study provides a valuable basis for studies of the human microbiome, with significant implications in terms of biological and clinical outcomes.

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data

Background

Human leucocyte antigen (HLA) genes play an important role in determining the outcome of organ transplantation and are linked to many human diseases. Because of the diversity and polymorphisms of HLA loci, HLA typing at high resolution is challenging even with whole-genome sequencing data.

Results

We have developed a computational tool, HLA-VBSeq, to estimate the most probable HLA alleles at full (8-digit) resolution from whole-genome sequence data. HLA-VBSeq simultaneously optimizes read alignments to HLA allele sequences and abundance of reads on HLA alleles by variational Bayesian inference. We show the effectiveness of the proposed method over other methods through the analysis of predicting HLA types for HLA class I (HLA-A, -B and -C) and class II (HLA-DQA1,-DQB1 and -DRB1) loci from the simulation data of various depth of coverage, and real sequencing data of human trio samples.

Conclusions

HLA-VBSeq is an efficient and accurate HLA typing method using high-throughput sequencing data without the need of primer design for HLA loci. Moreover, it does not assume any prior knowledge about HLA allele frequencies, and hence HLA-VBSeq is broadly applicable to human samples obtained from a genetically diverse population.

TIGAR2: sensitive and accurate estimation of transcript isoform expression with longer RNA-Seq reads

BACKGROUND

High-throughput RNA sequencing (RNA-Seq) enables quantification and identification of transcripts at single-base resolution. Recently, longer sequence reads become available thanks to the development of new types of sequencing technologies as well as improvements in chemical reagents for the Next Generation Sequencers. Although several computational methods have been proposed for quantifying gene expression levels from RNA-Seq data, they are not sufficiently optimized for longer reads (e.g. >250 bp).

RESULTS

We propose TIGAR2, a statistical method for quantifying transcript isoforms from fixed and variable length RNA-Seq data. Our method models substitution, deletion, and insertion errors of sequencers based on gapped-alignments of reads to the reference cDNA sequences so that sensitive read-aligners such as Bowtie2 and BWA-MEM are effectively incorporated in our pipeline. Also, a heuristic algorithm is implemented in variational Bayesian inference for faster computation. We apply TIGAR2 to both simulation data and real data of human samples and evaluate performance of transcript quantification with TIGAR2 in comparison to existing methods.

CONCLUSIONS

TIGAR2 is a sensitive and accurate tool for quantifying transcript isoform abundances from RNA-Seq data. Our method performs better than existing methods for the fixed-length reads (100 bp, 250 bp, 500 bp, and 1000 bp of both single-end and paired-end) and variable-length reads, especially for reads longer than 250 bp.

An efficient quantitation method of next-generation sequencing libraries by using MiSeq sequencer

Library quantitation is a critical step to obtain high data output in Illumina HiSeq sequencers. Here, we introduce a library quantitation method that uses the Illumina MiSeq sequencer designated as quantitative MiSeq (qMiSeq). In this procedure, 96 dual-index libraries, including control samples, are denatured, pooled in equal volume, and sequenced by MiSeq. We found that relative concentration of each library can be determined based on the observed index ratio and can be used to determine HiSeq run condition for each library. Thus, qMiSeq provides an efficient way to quantitate a large number of libraries at a time.

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

BACKGROUND

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

RESULTS

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

CONCLUSIONS

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

SUGAR: graphical user interface-based data refiner for high-throughput DNA sequencing

Background

Next-generation sequencers (NGSs) have become one of the main tools for current biology. To obtain useful insights from the NGS data, it is essential to control low-quality portions of the data affected by technical errors such as air bubbles in sequencing fluidics.

Results

We develop a software SUGAR (subtile-based GUI-assisted refiner) which can handle ultra-high-throughput data with user-friendly graphical user interface (GUI) and interactive analysis capability. The SUGAR generates high-resolution quality heatmaps of the flowcell, enabling users to find possible signals of technical errors during the sequencing. The sequencing data generated from the error-affected regions of a flowcell can be selectively removed by automated analysis or GUI-assisted operations implemented in the SUGAR. The automated data-cleaning function based on sequence read quality (Phred) scores was applied to a public whole human genome sequencing data and we proved the overall mapping quality was improved.

Conclusion

The detailed data evaluation and cleaning enabled by SUGAR would reduce technical problems in sequence read mapping, improving subsequent variant analysis that require high-quality sequence data and mapping results. Therefore, the software will be especially useful to control the quality of variant calls to the low population cells, e.g., cancers, in a sample with technical errors of sequencing procedures.

An efficient method of exploring simulation models by assimilating literature and biological observational data

Recently, several biological simulation models of, e.g., gene regulatory networks and metabolic pathways, have been constructed based on existing knowledge of biomolecular reactions, e.g., DNA-protein and protein-protein interactions. However, since these do not always contain all necessary molecules and reactions, their simulation results can be inconsistent with observational data. Therefore, improvements in such simulation models are urgently required. A previously reported method created multiple candidate simulation models by partially modifying existing models. However, this approach was computationally costly and could not handle a large number of candidates that are required to find models whose simulation results are highly consistent with the data. In order to overcome the problem, we focused on the fact that the qualitative dynamics of simulation models are highly similar if they share a certain amount of regulatory structures. This indicates that better fitting candidates tend to share the basic regulatory structure of the best fitting candidate, which can best predict the data among candidates. Thus, instead of evaluating all candidates, we propose an efficient explorative method that can selectively and sequentially evaluate candidates based on the similarity of their regulatory structures. Furthermore, in estimating the parameter values of a candidate, e.g., synthesis and degradation rates of mRNA, for the data, those of the previously evaluated candidates can be utilized. The method is applied here to the pharmacogenomic pathways for corticosteroids in rats, using time-series microarray expression data. In the performance test, we succeeded in obtaining more than 80% of consistent solutions within 15% of the computational time as compared to the comprehensive evaluation. Then, we applied this approach to 142 literature-recorded simulation models of corticosteroid-induced genes, and consequently selected 134 newly constructed better models. The method described here was found to be capable of efficiently exploring candidate simulation models and obtaining better models within a short span of time. Furthermore, the results suggest that there may be room for improvement in literature recorded pathways and that they can be systematically updated using biological observational data.

Identification of acquired mutations by whole-genome sequencing in GATA-2 deficiency evolving into myelodysplasia and acute leukemia

Heterozygous GATA-2 germline mutations are associated with overlapping clinical manifestations termed GATA-2 deficiency, characterized by immunodeficiency and predisposition to myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML). However, there is considerable clinical heterogeneity among patients, and the molecular basis for the evolution of immunodeficiency into MDS/AML remains unknown. Thus, we conducted whole-genome sequencing on a patient with a germline GATA-2 heterozygous mutation (c. 988 C > T; p. R330X), who had a history suggestive of immunodeficiency and evolved into MDS/AML. Analysis was conducted with DNA samples from leukocytes for immunodeficiency, bone marrow mononuclear cells for MDS and bone marrow-derived mesenchymal stem cells. Whereas we did not identify a candidate genomic deletion that may contribute to the evolution into MDS, a total of 280 MDS-specific nonsynonymous single nucleotide variants were identified. By narrowing down with the single nucleotide polymorphism database, the functional missense database, and NCBI information, we finally identified three candidate mutations for EZH2, HECW2 and GATA-1, which may contribute to the evolution of the disease.

iSVP: an integrated structural variant calling pipeline from high-throughput sequencing data

BACKGROUND

Structural variations (SVs), such as insertions, deletions, inversions, and duplications, are a common feature in human genomes, and a number of studies have reported that such SVs are associated with human diseases. Although the progress of next generation sequencing (NGS) technologies has led to the discovery of a large number of SVs, accurate and genome-wide detection of SVs remains challenging. Thus far, various calling algorithms based on NGS data have been proposed. However, their strategies are diverse and there is no tool able to detect a full range of SVs accurately.

RESULTS

We focused on evaluating the performance of existing deletion calling algorithms for various spanning ranges from low- to high-coverage simulation data. The simulation data was generated from a whole genome sequence with artificial SVs constructed based on the distribution of variants obtained from the 1000 Genomes Project. From the simulation analysis, deletion calls of various deletion sizes were obtained with each caller, and it was found that the performance was quite different according to the type of algorithms and targeting deletion size. Based on these results, we propose an integrated structural variant calling pipeline (iSVP) that combines existing methods with a newly devised filtering and merging processes. It achieved highly accurate deletion calling with >90% precision and >90% recall on the 30× read data for a broad range of size. We applied iSVP to the whole-genome sequence data of a CEU HapMap sample, and detected a large number of deletions, including notable peaks around 300 bp and 6,000 bp, which corresponded to Alus and long interspersed nuclear elements, respectively. In addition, many of the predicted deletions were highly consistent with experimentally validated ones by other studies.

CONCLUSIONS

We present iSVP, a new deletion calling pipeline to obtain a genome-wide landscape of deletions in a highly accurate manner. From simulation and real data analysis, we show that iSVP is broadly applicable to human whole-genome sequencing data, which will elucidate relationships between SVs across genomes and associated diseases or biological functions.

Recurrent mutations in multiple components of the cohesin complex in myeloid neoplasms

Cohesin is a multimeric protein complex that is involved in the cohesion of sister chromatids, post-replicative DNA repair and transcriptional regulation. Here we report recurrent mutations and deletions involving multiple components of the cohesin complex, including STAG2, RAD21, SMC1A and SMC3, in different myeloid neoplasms. These mutations and deletions were mostly mutually exclusive and occurred in 12.1% (19/157) of acute myeloid leukemia, 8.0% (18/224) of myelodysplastic syndromes, 10.2% (9/88) of chronic myelomonocytic leukemia, 6.3% (4/64) of chronic myelogenous leukemia and 1.3% (1/77) of classical myeloproliferative neoplasms. Cohesin-mutated leukemic cells showed reduced amounts of chromatin-bound cohesin components, suggesting a substantial loss of cohesin binding sites on chromatin. The growth of leukemic cell lines harboring a mutation in RAD21 (Kasumi-1 cells) or having severely reduced expression of RAD21 and STAG2 (MOLM-13 cells) was suppressed by forced expression of wild-type RAD21 and wild-type RAD21 and STAG2, respectively. These findings suggest a role for compromised cohesin functions in myeloid leukemogenesis.