Search

Recent advances in human genome research have provided evidence for genotype-phenotype associations, pathogenicity, and clinical actionability of variants and genomic risk prediction of disease. However, the return of individual genomic results to healthy individuals is fraught with ethical and practical complexity.Individual genomic results were returned to BRCA1/2 pathogenic variant (PV) carriers of the Tohoku Medical Megabank cohort study participants with an information on hereditary breast and ovarian cancer syndrome (HBOC). One hundred and eighty participants, including 9 BRCA1/2 PV carriers, were asked about their willingness to receive individual genomic results, without revealing the gene name and related disorders, prior to the study. Of the 142 participants who responded, 103 showed willingness to know their genomic information. Each of the six BRCA1/2 PV carriers who consented to participate in the study received information about HBOC in person and underwent validation testing with blood resampling.All participants were in their 60s or 70s; of the four females and two males, two had a history of breast cancer and five had a family history of HBOC-related cancers. All participants appreciated the information, without remarkable negative psychological impact of the return, and intended to undergo clinical risk surveillance. Five participants were accompanied by family members while receiving the results, and three first-degree female relatives wished to undergo genomic testing at the hospital.Our results suggest that returning actionable genomic information to participants in a population-based genome cohort study is beneficial for preventing or providing early-stage intervention for associated diseases.

Disorders of sex development (DSD) comprises a congenital condition in which chromosomal, gonadal, or anatomical sex development is atypical. In this study, we screened for pathogenic variants in 32 genes associated with DSDs and central causes of hypogonadism (CHG) in a whole-genome reference panel including 8380 Japanese individuals constructed by Tohoku Medical Megabank Organization. Candidate pathogenic (P) or likely pathogenic (LP) variants were extracted from the ClinVar, InterVar, and Human Gene Mutation databases. Ninety-one candidate pathological variants were found in 25 genes; 28 novel candidate variants were identified. Nearly 1 in 40 (either ClinVar or InterVar P or LP) to 157 (both ClinVar and InterVar P or LP) individuals were found to be carriers of recessive DSD and CHG alleles. In these data, genes implicated in gonadal dysfunction did not show loss-of-function variants, with a relatively high tendency of intolerance for haploinsufficiency based on pLI and Episcore, both of which can be used for estimating haploinsufficiency. We report the types and frequencies of causative variants for DSD and CHG in the general Japanese population. This study furthers our understanding of the genetic causes and helps to refine genetic counseling of DSD and CHG.

Pharmacogenomic (PGx) testing results provide valuable information on drug selection and appropriate dosing, maximization of efficacy, and minimization of adverse effects. Although the number of large-scale, next-generation-sequencing-based PGx studies has recently increased, little is known about the risks and benefits of returning PGx results to ostensibly healthy individuals in research settings.Single-nucleotide variants of three actionable PGx genes, namely,Whereas the majority of participants were unfamiliar with the term PGx, and none had undergone PGx testing before the study, more than 80% of the participants felt that they could acquire basic PGx knowledge sufficient to understand their genomic results and were satisfied with their potential benefit and use in future prescriptions. On the other hand, some felt that the PGx concepts or terminology was difficult to fully understand and suggested that in-person return of the results was desirable.These results collectively suggest possible benefits of returning preemptive PGx information to ostensibly healthy cohort participants in a research setting.

Certain large genome cohort studies attempt to return the individual genomic results to the participants; however, the implementation process and psychosocial impacts remain largely unknown. The Tohoku Medical Megabank Project has conducted large genome cohort studies of general residents. To implement the disclosure of individual genomic results, we extracted the potential challenges and obstacles. Major challenges include the determination of genes/disorders based on the current medical system in Japan, the storage of results, prevention of misunderstanding, and collaboration of medical professionals. To overcome these challenges, we plan to conduct multilayer pilot studies, which deal with different disorders/genes. We finally chose familial hypercholesterolemia (FH) as a target disease for the first pilot study. Of the 665 eligible candidates, 33.5% were interested in the pilot study and provided consent after an educational "genetics workshop" on the basic genetics and medical facts of FH. The genetics professionals disclosed the results to the participants. All positive participants were referred to medical care, and a serial questionnaire revealed no significant psychosocial distress after the disclosure. Return of genomic results to research participants was implemented using a well-prepared protocol. To further elucidate the impact of different disorders, we will perform multilayer pilot studies with different disorders, including actionable pharmacogenomics and hereditary tumor syndromes.

Bone dysplasias are a group of rare hereditary diseases, with up to 436 disease types. Perinatal diagnosis is clinically important for adequate personalized management and counseling. There are no reports focused on pathogenic variants of bone dysplasias in the general population. In this study, we focused on autosomal recessive bone dysplasias. We identified pathogenic variants using whole-genome reference panel data from 3552 Japanese individuals. For the first time, we were able to estimate the carrier frequencies and the proportions of potential patients. For autosomal recessive bone dysplasias, we detected 198 pathogenic variants of 54 causative genes. We estimated the variant carrier frequencies and the proportions of potential patients with variants associated with four clinically important bone dysplasias: osteogenesis imperfecta (OI), hypophosphatasia (HPP), asphyxiating thoracic dysplasia (ATD), and Ellis–van Creveld syndrome (EvC). The proportions of potential patients with OI, ATD, and EvC based on pathogenic variants classified as “pathogenic” and “likely pathogenic” by InterVar were closer to the reported incidence rates in Japanese subjects. Furthermore, the proportions of potential patients with HPP variants classified as “pathogenic” and “likely pathogenic” in InterVar and “pathogenic” in ClinVar were closer to the reported incidence rates. For bone dysplasia, the findings of this study will provide a better understanding of the variant types and frequencies in the Japanese general population, and should be useful for clinical diagnosis, genetic counseling, and personalized medicine., Bone dysplasia: Estimating carrier frequency with bioinformatics A bioinformatics approach helps estimate carrier frequency of a rare inherited bone disease which causes abnormalities in skeletal shape and structure. Autosomal recessive bone dysplasias affect bone and cartilage development and result from inheriting two mutated genes, one from each parent. Junichi Sugawara, Tohoku University, Sendai, Japan, and colleagues used mutation databases and a bioinformatics tool for variant interpretation to detect 198 pathogenic variants in 54 genes associated with autosomal recessive bone dysplasia in a whole-genome reference panel of 3,552 general Japanese individuals (3.5KJPNv2). They then estimated the frequency of people in the sample carrying bone dysplasia mutations and the expected proportion in whom the disorder could manifest, which compared well with reported incidence rates in the general population. These findings could prove useful for calculating the risk of bone dysplasia in the future children of carrier parents.

Alzheimer’s disease (AD) is a common neurological disease that causes dementia in humans. Although the reports of associated pathological genes have been increasing, the molecular mechanism leading to the accumulation of amyloid-β (Aβ) in human brain is still not well understood. To identify novel genes that cause accumulation of Aβ in AD patients, we conducted an integrative analysis by combining a human genetic association study and transcriptome analysis in mouse brain. First, we examined genome-wide gene expression levels in the hippocampus, comparing them to amyloid Aβ level in mice with mixed genetic backgrounds. Next, based on a GWAS statistics obtained by a previous study with human AD subjects, we obtained gene-based statistics from the SNP-based statistics. We combined p values from the two types of analysis across orthologous gene pairs in human and mouse into one p value for each gene to evaluate AD susceptibility. As a result, we found five genes with significant p values in this integrated analysis among the 373 genes analyzed. We also examined the gene expression level of these five genes in the hippocampus of independent human AD cases and control subjects. Two genes, LBH and SHF, showed lower expression levels in AD cases than control subjects. This is consistent with the gene expression levels of both the genes in mouse which were negatively correlated with Aβ accumulation. These results, obtained from the integrative approach, suggest that LBH and SHF are associated with the AD pathogenesis. Electronic supplementary material The online version of this article (10.1007/s00439-018-1906-z) contains supplementary material, which is available to authorized users.

Clarifying allele frequencies of disease-related genetic variants in a population is important in genomic medicine; however, such data is not yet available for the Japanese population. To estimate frequencies of actionable pathogenic variants in the Japanese population, we examined the reported pathological variants in genes recommended by the American College of Medical Genetics and Genomics (ACMG) in our reference panel of genomic variations, 2KJPN, which was created by whole-genome sequencing of 2049 individuals of the resident cohort of the Tohoku Medical Megabank Project. We searched for pathogenic variants in 2KJPN for 57 autosomal ACMG-recommended genes responsible for 26 diseases and then examined their frequencies. By referring to public databases of pathogenic variations, we identified 143 reported pathogenic variants in 2KJPN for the 57 ACMG recommended genes based on a classification system. At the individual level, 21% of the individuals were found to have at least one reported pathogenic allele. We then conducted a literature survey to review the variants and to check for evidence of pathogenicity. Our results suggest that a substantial number of people have reported pathogenic alleles for the ACMG genes, and reviewing variants is indispensable for constructing the information infrastructure of genomic medicine for the Japanese population.

The first step towards realizing personalized healthcare is to catalog the genetic variations in a population. Since the dissemination of individual-level genomic information is strictly controlled, it will be useful to construct population-level allele frequency panels with easy-to-use interfaces. In the Tohoku Medical Megabank Project, we sequenced nearly 4000 individuals from a Japanese population and constructed an allele frequency panel of 3552 individuals after removing related samples. The panel is called the 3.5KJPNv2. It was constructed by using a standard pipeline including the 1KGP and gnomAD algorithms to reduce technical biases and to allow comparisons to other populations. Our database is the first large-scale panel providing the frequencies of variants present on the X chromosome and on the mitochondria in the Japanese population. All the data are available on our original database at https://jmorp.megabank.tohoku.ac.jp., Population genetics: large database of Japanese gene variations constructed A new database provides information on the frequency of genetic variations within 3552 Japanese individuals, and facilitates comparisons with other populations. The reference panel, constructed by Kengo Kinoshita of Tohoku University, Sendai, and colleagues in Japan is also the first large-scale database to provide genetic variation frequency information on the X chromosome and mitochondrial DNA in the Japanese population. The methods used to sequence the genetic data are similar to those used in other large databases, allowing comparisons with other populations. The population size and methods used to compile the database overcome limitations in previous Japanese reference panels. This and similar databases that catalog genetic variations within populations can improve efforts towards personalizing healthcare and contribute to the study of human population genetics. The database is publicly available online.

Incidence rates of Mendelian diseases vary among ethnic groups, and frequencies of variant types of causative genes also vary among human populations. In this study, we examined to what extent we can predict population frequencies of recessive disorders from genomic data, and explored better strategies for variant interpretation and classification. We used a whole-genome reference panel from 3552 general Japanese individuals constructed by the Tohoku Medical Megabank Organization (ToMMo). Focusing on 32 genes for 17 congenital metabolic disorders included in newborn screening (NBS) in Japan, we identified reported and predicted pathogenic variants through variant annotation, interpretation, and multiple ways of classifications. The estimated carrier frequencies were compared with those from the Japanese NBS data based on 1,949,987 newborns from a previous study. The estimated carrier frequency based on genomic data with a recent guideline of variant interpretation for the PAH gene, in which defects cause hyperphenylalaninemia (HPA) and phenylketonuria (PKU), provided a closer estimate to that by the observed incidence than the other methods. In contrast, the estimated carrier frequencies for SLC25A13, which causes citrin deficiency, were much higher compared with the incidence rate. The results varied greatly among the 11 NBS diseases with single responsible genes; the possible reasons for departures from the carrier frequencies by reported incidence rates were discussed. Of note, (1) the number of pathogenic variants increases by including additional lines of evidence, (2) common variants with mild effects also contribute to the actual frequency of patients, and (3) penetrance of each variant remains unclear.

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

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.

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.

Background Corticotropin-releasing hormone (CRH) plays an important role in the pathophysiology of irritable bowel syndrome (IBS) and regulates the stress response through two CRH receptors (R1 and R2). Previously, we reported that a CRHR1 gene polymorphism (rs110402, rs242924, and rs7209436) and haplotypes were associated with IBS. However, the association between the CRHR2 gene and IBS was not investigated. We tested the hypothesis that genetic polymorphisms and haplotypes of CRHR2 are associated with IBS pathophysiology and negative emotion in IBS patients. Methods A total of 142 IBS patients and 142 healthy controls participated in this study. Seven single nucleotide polymorphisms (SNPs) of the CRHR2 gene (rs4722999, rs3779250, rs2240403, rs2267710, rs2190242, rs2284217, and rs2284220) were genotyped. Subjects' psychological states were evaluated using the Perceived-Stress Scale, the State-Trait Anxiety Inventory, and the Self-Rating Depression Scale. Results We found that rs4722999 and rs3779250, located in intronic region, were associated with IBS in terms of genotype frequency (rs4722999: P = 0.037; rs3779250: P = 0.017) and that the distribution of the major allele was significantly different between patients and controls. There was a significant group effect (controls vs. IBS), and a CRHR2 genotype effect was observed for three psychological scores, but the interaction was not significant. We found a haplotype of four SNPs (rs4722999, rs3779250, rs2240403, and rs2267710) and two SNPs (rs2284217 and rs2284220) in strong linkage disequilibrium (D′ > 0.90). We also found that haplotypes of the CRHR2 gene were significantly different between IBS patients and controls and that they were associated with negative emotion. Conclusion Our findings support the hypothesis that genetic polymorphisms and haplotypes of CRHR2 are related to IBS. In addition, we found associations between CRHR2 genotypes and haplotypes and negative emotion in IBS patients and controls. Further studies on IBS and the CRH system are warranted.

Although the Japanese population has a rather low genetic diversity, we recently confirmed the presence of two main clusters (the Hondo and Ryukyu clusters) through principal component analysis of genome-wide single-nucleotide polymorphism (SNP) genotypes. Understanding the genetic differences between the two main clusters requires further genome-wide analyses based on a dense SNP set and comparison of haplotype frequencies. In the present study, we determined haplotypes for the Hondo cluster of the Japanese population by detecting SNP homozygotes with 388,591 autosomal SNPs from 18,379 individuals and estimated the haplotype frequencies. Haplotypes for the Ryukyu cluster were inferred by a statistical approach using the genotype data from 504 individuals. We then compared the haplotype frequencies between the Hondo and Ryukyu clusters. In most genomic regions, the haplotype frequencies in the Hondo and Ryukyu clusters were very similar. However, in addition to the human leukocyte antigen region on chromosome 6, other genomic regions (chromosomes 3, 4, 5, 7, 10 and 12) showed dissimilarities in haplotype frequency. These regions were enriched for genes involved in the immune system, cell-cell adhesion and the intracellular signaling cascade. These differentiated genomic regions between the Hondo and Ryukyu clusters are of interest because they (1) should be examined carefully in association studies and (2) likely contain genes responsible for morphological or physiological differences between the two groups.

Because population stratification can cause spurious associations in case-control studies, understanding the population structure is important. Here, we examined Japanese population structure by "Eigenanalysis," using the genotypes for 140,387 SNPs in 7003 Japanese individuals, along with 60 European, 60 African, and 90 East-Asian individuals, in the HapMap project. Most Japanese individuals fell into two main clusters, Hondo and Ryukyu; the Hondo cluster includes most of the individuals from the main islands in Japan, and the Ryukyu cluster includes most of the individuals from Okinawa. The SNPs with the greatest frequency differences between the Hondo and Ryukyu clusters were found in the HLA region in chromosome 6. The nonsynonymous SNPs with the greatest frequency differences between the Hondo and Ryukyu clusters were the Val/Ala polymorphism (rs3827760) in the EDAR gene, associated with hair thickness, and the Gly/Ala polymorphism (rs17822931) in the ABCC11 gene, associated with ear-wax type. Genetic differentiation was observed, even among different regions in Honshu Island, the largest island of Japan. Simulation studies showed that the inclusion of different proportions of individuals from different regions of Japan in case and control groups can lead to an inflated rate of false-positive results when the sample sizes are large.

Short Tandem Repeats (STRs) comprise repeats of one to several base pairs. Because of the high mutability due to strand slippage during DNA synthesis, rapid evolutionary change in the number of repeating units directly shapes the range of repeat-number variation according to selection pressure. However, the remaining questions include: Why are STRs causing repeat expansion diseases maintained in the human population; and why are these limited to neurodegenerative diseases? By evaluating the genome-wide selection pressure on STRs using the database we constructed, we identified two different patterns of relationship in repeat-number polymorphisms between DNA and amino-acid sequences, although both patterns are evolutionary consequences of avoiding the formation of harmful long STRs. First, a mixture of degenerate codons is represented in poly-proline (poly-P) repeats. Second, long poly-glutamine (poly-Q) repeats are favored at the protein level; however, at the DNA level, STRs encoding long poly-Qs are frequently divided by synonymous SNPs. Furthermore, significant enrichments of apoptosis and neurodevelopment were biological processes found specifically in genes encoding poly-Qs with repeat polymorphism. This suggests the existence of a specific molecular function for polymorphic and/or long poly-Q stretches. Given that the poly-Qs causing expansion diseases were longer than other poly-Qs, even in healthy subjects, our results indicate that the evolutionary benefits of long and/or polymorphic poly-Q stretches outweigh the risks of long CAG repeats predisposing to pathological hyper-expansions. Molecular pathways in neurodevelopment requiring long and polymorphic poly-Q stretches may provide a clue to understanding why poly-Q expansion diseases are limited to neurodegenerative diseases.

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, 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.

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.

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.

To clarify the relationship between the amino acid variations of the gp120 of human immunodeficiency virus type 1 (HIV-1) and the chemokine receptors that are used as the second receptor for HIV, we evaluated amino acid site variation of gp120 between the X4 strains (use CXCR4) and the R5 strains (use CCR5) from 21 sequences of subtype B. Our analysis showed that residues 306 and 322 in the V3 loop and residue 440 in the C4 region were associated with usage of the second receptor. The polymorphism at residue 440 is clearly associated with the usage of the second receptor: The amino acid at position 440 was a basic amino acid in the R5 strains, and a nonbasic and smaller amino acid in the X4 strains, while the V3 loop of the X4 strains was more basic than that of the R5 strains. This suggests that residue 440 in the C4 region, which is close to the V3 loop in the three-dimensional structure, is critical in determining which second receptor is used. Analysis of codon frequency suggests that, in almost all cases, the difference at residue 440 between basic amino acids in the R5 strains and nonbasic amino acids in the X4 strains could be due to a single nucleotide change. These findings predict that the evolutionary changes in amino acid residue 440 may be correlated with evolutionary changes in the V3 loop. One possibility is that a change in electric charge at residue 440 compensates for a change in electric charge in the V3 loop. The amino acid polymorphism at position 440 can be useful to predict the cell tropism of a strain of HIV-1 subtype B.

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. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-673) contains supplementary material, which is available to authorized users.

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.

Alzheimer’s disease (AD) is characterized by the accumulation of amyloid-β (Aβ). The genes that govern this process, however, have remained elusive. To this end, we combined distinct mouse strains with transcriptomics to directly identify disease-relevant genes. We show that AD model mice (APP-Tg) with DBA/2 genetic backgrounds have significantly lower levels of Aβ accumulation compared with SJL and C57BL/6 mice. We then applied brain transcriptomics to reveal the genes in DBA/2 that suppress Aβ accumulation. To avoid detecting secondarily affected genes by Aβ, we used non-Tg mice in the absence of Aβ pathology and selected candidate genes differently expressed in DBA/2 mice. Additional transcriptome analysis of APP-Tg mice with mixed genetic backgrounds revealed kinesin light chain-1 (Klc1) as an Aβ modifier, indicating a role for intracellular trafficking in Aβ accumulation. Aβ levels correlated with the expression levels of Klc1 splice variant E and the genotype of Klc1 in these APP-Tg mice. In humans, the expression levels of KLC1 variant E in brain and lymphocyte were significantly higher in AD patients compared with unaffected individuals. Finally, functional analysis using neuroblastoma cells showed that overexpression or knockdown of KLC1 variant E increases or decreases the production of Aβ, respectively. The identification of KLC1 variant E suggests that the dysfunction of intracellular trafficking is a causative factor of Aβ pathology. This unique combination of distinct mouse strains and model mice with transcriptomics is expected to be useful for the study of genetic mechanisms of other complex diseases.

Haplotype phasing is essential for identifying disease-causing variants with phase-dependent interactions as well as for the coalescent-based inference of demographic history. One of approaches for estimating haplotypes is to use phase-informative reads, which span multiple heterozygous variant positions. Although the quality of estimated variants is crucial in haplotype phasing, accurate variant calling is still challenging due to errors on sequencing and read mapping. Since some of such errors can be corrected by considering haplotype phasing, simultaneous estimation of variants and haplotypes is important. Thus, we propose a statistically unified approach for variant calling and haplotype phasing named HapMonster, where haplotype phasing information is used for improving the accuracy of variant calling and the improved variant calls are used for more accurate haplotype phasing. From the comparison with other existing methods on simulation and real sequencing data, we confirm the effectiveness of HapMonster in both variant calling and haplotype phasing.

Recent development of next generation sequencing (NGS) technologies has led to the identification of structural variants (SVs) of genomic DNA existing in the human population. Several SV detection methods utilizing NGS data have been proposed. However, there are several difficulties in analysis of NGS data, particularly with regard to handling reads from duplicated loci or low-complexity sequences of the human genome. In this paper, we propose SVEM, a novel statistical method to detect SVs with a single nucleotide resolution that can utilize multi-mapped reads on breakpoints. SVEM estimates the amount of reads on breakpoints as parameters and mapping states as latent variables using the expectation maximization algorithm. This framework enables us to handle ambiguous mapping of reads without discarding information for SV detection. SVEM is applied to simulation data and real data, and it achieves better performance than existing methods in terms of precision and recall.

We analyzed 57 HIV-1 isolates from Cameroon and the Republic of Congo, with respect to the env C2V3 and/or the pol integrase regions. The results indicated that the topology of the pol tree correlated well with that of the env tree for four clusters of subtype D, F G and H, suggesting that these trees reflect the true evolution of the overall genome structures of these subtypes. However, of 22 Cameroonian isolates that were classified as subtype A based on env, 20 of them diverged in their pol sequence into two lineages that were completely different from the prototypical subtype A, tentatively designated as subtypes A1 and A2. The subtype A1 isolates (6 out of 22) were related in their env C2V3 regions with prototypical subtype A strain, but in their pol regions, they formed an independent cluster that diverged from known HIV-1 subtypes so far reported (except for subtypes I and J). The subtype A2 isolates (14 out of 22), which represent the major epidemic type of HIV-1 in Cameroon, clustered distinctly in both the env and pol trees with the recently described A/G mosaic strains from Nigeria and Djibouti. These two lineages were not spreading in the neighboring Republic of Congo.

Here we describe, for the first time, recombinants between two highly divergent major groups of human immunodeficiency virus type 1 (HIV-1), M and O, within a Cameroonian woman infected with three different HIV-1 strains, a group O virus, a subtype D virus, and a recently reported IBNG (A/G)-like recombinant virus. Using nested extra-long PCR amplification, we sequenced from thepolregion to theenvregion including accessory genes of the viral genome obtained from the patient’s uncultured peripheral blood mononuclear cells and examined the phylogenetic position of each gene. Compared with sequential blood samples obtained in 1995 and 1996, there were multiple segmental exchanges between three HIV-1 strains (O, D, and IBNG) and all the recombinants appeared to be derived from a common M/O ancestor. Importantly, recombination between groups M and O occurred, even though the homology between these two groups is 69, 76, 68, and 55% in thegag,pol,vif-vpr, andenvregions, respectively. Recombination between strains with such distant lineages may contribute substantially to generating new HIV-1 variants.

A stylar S-RNase is associated with gametophytic self-incompatibility in the Rosaceae, Solanaceae, and Scrophulariaceae. This S-RNase is responsible for S-allele-specific recognition in the self-incompatible reaction, but how it functions in specific discrimination is not clear. Window analysis of the numbers of synonymous (dS) and non-synonymous (dN) substitutions in rosaceous S-RNases detected four regions with an excess of dN over dS in which positive selection may operate (PS regions). The topology of the secondary structure of the S-RNases predicted by the PHD method is very similar to that of fungal RNase Rh whose tertiary structure is known. When the sequences of S-RNases are aligned with the sequence of RNase Rh based on the predicted secondary structures, the four PS regions correspond to two surface sites on the tertiary structure of RNase Rh. These findings suggest that in S-RNases the PS regions also form two sites and are candidates for the recognition sites for S-allele-specific discrimination.

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.

Motivation: Variant calling from genome-wide sequencing data is essential for the analysis of disease-causing mutations and elucidation of disease mechanisms. However, variant calling in low coverage regions is difficult due to sequence read errors and mapping errors. Hence, variant calling approaches that are robust to low coverage data are demanded. Results: We propose a new variant calling approach that considers pedigree information and haplotyping based on sequence reads spanning two or more heterozygous positions termed phase informative reads. In our approach, genotyping and haplotyping by the assignment of each read to a haplotype based on phase informative reads are simultaneously performed. Therefore, positions with low evidence for heterozygosity are rescued by phase informative reads, and such rescued positions contribute to haplotyping in a synergistic way. In addition, pedigree information supports more accurate haplotyping as well as genotyping, especially in low coverage regions. Although heterozygous positions are useful for haplotyping, homozygous positions are not informative and weaken the information from heterozygous positions, as majority of positions are homozygous. Thus, we introduce latent variables that determine zygosity at each position to filter out homozygous positions for haplotyping. In performance evaluation with a parent–offspring trio sequencing data, our approach outperforms existing approaches in accuracy on the agreement with single nucleotide polymorphism array genotyping results. Also, performance analysis considering distance between variants showed that the use of phase informative reads is effective for accurate variant calling, and further performance improvement is expected with longer sequencing data. Contact: nagasaki@megabank.tohoku.ac.jp or kojima@megabank.tohoku.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online.

Complex diseases result from contributions of multiple genes that act in concert through pathways. Here we present a method to prioritize novel candidates of disease-susceptibility genes depending on the biological similarities to the known disease-related genes. The extent of disease-susceptibility of a gene is prioritized by analyzing seven features of human genes captured in H-InvDB. Taking rheumatoid arthritis (RA) and prostate cancer (PC) as two examples, we evaluated the efficiency of our method. Highly scored genes obtained included TNFSF12 and OSM as candidate disease genes for RA and PC, respectively. Subsequent characterization of these genes based upon an extensive literature survey reinforced the validity of these highly scored genes as possible disease-susceptibility genes. Our approach, Prioritization ANalysis of Disease Association (PANDA), is an efficient and cost-effective method to narrow down a large set of genes into smaller subsets that are most likely to be involved in the disease pathogenesis.

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