Progress report of the Tohoku Medical Megabank Community-Based Cohort Study: Study profile of the repeated center-based survey during second period in Miyagi Prefecture.

BACKGROUND: The purpose of this study was to report the basic profile of the Miyagi Prefecture part of a repeated center-based survey during the second period (2nd period survey) of the Tohoku Medical Megabank Community-Based Cohort Study (TMM CommCohort Study), as well as the participants' characteristics based on their participation type in the baseline survey. METHODS: The 2nd period survey, conducted from June 2017 to March 2021, included participants of the TMM CommCohort Study (May 2013 to March 2016). In addition to the questionnaire, blood, urine, and physiological function tests were performed during the 2nd period survey. There were three main ways of participation in the baseline survey: Type 1, Type 1 additional, or Type 2 survey. The 2nd period survey was conducted in the same manner as the Type 2 survey, which was based on the community support center (CSC). RESULTS: In Miyagi Prefecture, 29,383 (57.7%) of 50,967 participants participated in the 2nd period survey. The participation rate among individuals who had visited the CSC was approximately 80%. Although some factors differed depending on the participation type in the baseline survey, the 2nd period survey respondents in the Type 1 and Type 2 survey groups at baseline had similar traits. CONCLUSIONS: The 2nd period survey of the TMM CommCohort Study provided detailed follow-up information. Following up on the health conditions of the participants will clarify the long-term effects of disasters and contribute to personalized prevention.

Returning individual genomic results to population-based cohort study participants with BRCA1/2 pathogenic variants.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Loss of CAPS2/Cadps2 leads to exocrine pancreatic cell injury and intracellular accumulation of secretory granules in mice.

The type 2 Ca2+-dependent activator protein for secretion (CAPS2/CADPS2) regulates dense-core vesicle trafficking and exocytosis and is involved in the regulated release of catecholamines, peptidergic hormones, and neuromodulators. CAPS2 is expressed in the pancreatic exocrine acinar cells that produce and secrete digestive enzymes. However, the functional role of CAPS2 in vesicular trafficking and/or exocytosis of non-regulatory proteins in the exocrine pancreas remains to be determined. Here, we analyzed the morpho-pathological indicators of the pancreatic exocrine pathway in Cadps2-deficient mouse models using histochemistry, biochemistry, and electron microscopy. We used whole exosome sequencing to identify CADPS2 variants in patients with chronic pancreatitis (CP). Caps2/Cadps2-knockout (KO) mice exhibited morphophysiological abnormalities in the exocrine pancreas, including excessive accumulation of secretory granules (zymogen granules) and their amylase content in the cytoplasm, deterioration of the fine intracellular membrane structures (disorganized rough endoplasmic reticulum, dilated Golgi cisternae, and the appearance of empty vesicles and autophagic-like vacuoles), as well as exocrine pancreatic cell injury, including acinar cell atrophy, increased fibrosis, and inflammatory cell infiltration. Pancreas-specific Cadps2 conditional KO mice exhibited pathological abnormalities in the exocrine pancreas similar to the global Cadps2 KO mice, indicating that these phenotypes were caused either directly or indirectly by CAPS2 deficiency in the pancreas. Furthermore, we identified a rare variant in the exon3 coding region of CADPS2 in a non-alcoholic patient with CP and showed that Cadps2-dex3 mice lacking CAPS2 exon3 exhibited symptoms similar to those exhibited by the Cadps2 KO and cKO mice. These results suggest that CAPS2 is critical for the proper functioning of the pancreatic exocrine pathway, and its deficiency is associated with a risk of pancreatic acinar cell pathology.

Pathological variants in genes associated with disorders of sex development and central causes of hypogonadism in a whole-genome reference panel of 8380 Japanese individuals.

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.

A Pilot Study for Return of Individual Pharmacogenomic Results to Population-Based Cohort Study Participants.

Introduction: 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. Methods: Single-nucleotide variants of three actionable PGx genes, namely, MT-RNR1, CYP2C19, and NUDT15, were returned to 161 participants in a population-based Tohoku Medical Megabank project. Informed consent was obtained from the participants after a seminar on the outline of this study. The results were sent by mail alongside sealed information letter intended for clinicians. As an exception, genetic counseling was performed for the MT-RNR1 m.1555A > G variant carriers by a medical geneticist, and consultation with an otolaryngologist was encouraged. Questionnaire surveys (QSs) were conducted five times to evaluate the participants' understanding of the topic, psychological impact, and attitude toward the study. Results: 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. Conclusions: These results collectively suggest possible benefits of returning preemptive PGx information to ostensibly healthy cohort participants in a research setting.

The return of individual genomic results to research participants: design and pilot study of Tohoku Medical Megabank Project.

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.

Estimation of the carrier frequencies and proportions of potential patients by detecting causative gene variants associated with autosomal recessive bone dysplasia using a whole-genome reference panel of Japanese individuals.

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.

Novel candidates of pathogenic variants of the BRCA1 and BRCA2 genes from a dataset of 3,552 Japanese whole genomes (3.5KJPNv2).

Identification of the population frequencies of definitely pathogenic germline variants in two major hereditary breast and ovarian cancer syndrome (HBOC) genes, BRCA1/2, is essential to estimate the number of HBOC patients. In addition, the identification of moderately penetrant HBOC gene variants that contribute to increasing the risk of breast and ovarian cancers in a population is critical to establish personalized health care. A prospective cohort subjected to genome analysis can provide both sets of information. Computational scoring and prospective cohort studies may help to identify such likely pathogenic variants in the general population. We annotated the variants in the BRCA1 and BRCA2 genes from a dataset of 3,552 whole-genome sequences obtained from members of a prospective cohorts with genome data in the Tohoku Medical Megabank Project (TMM) with InterVar software. Computational impact scores (CADD_phred and Eigen_raw) and minor allele frequencies (MAFs) of pathogenic (P) and likely pathogenic (LP) variants in ClinVar were used for filtration criteria. Familial predispositions to cancers among the 35,000 TMM genome cohort participants were analyzed to verify the identified pathogenicity. Seven potentially pathogenic variants were newly identified. The sisters of carriers of these moderately deleterious variants and definite P and LP variants among members of the TMM prospective cohort showed a statistically significant preponderance for cancer onset, from the self-reported cancer history. Filtering by computational scoring and MAF is useful to identify potentially pathogenic variants in BRCA genes in the Japanese population. These results should help to follow up the carriers of variants of uncertain significance in the HBOC genes in the longitudinal prospective cohort study.

3.5KJPNv2: an allele frequency panel of 3552 Japanese individuals including the X chromosome.

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.

Estimating carrier frequencies of newborn screening disorders using a whole-genome reference panel of 3552 Japanese individuals.

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.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare.

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.

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.

Integrated analysis of human genetic association study and mouse transcriptome suggests LBH and SHF genes as novel susceptible genes for amyloid-ß accumulation in Alzheimer's disease.

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.

Evaluation of reported pathogenic variants and their frequencies in a Japanese population based on a whole-genome reference panel of 2049 individuals.

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

Selection pressure on human STR loci and its relevance in repeat expansion disease.

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.

Corticotropin-Releasing Hormone Receptor 2 Gene Variants in Irritable Bowel Syndrome.

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.

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.

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.

Estimating copy numbers of alleles from population-scale high-throughput sequencing data.

BACKGROUND: With the recent development of microarray and high-throughput sequencing (HTS) technologies, a number of studies have revealed catalogs of copy number variants (CNVs) and their association with phenotypes and complex traits. In parallel, a number of approaches to predict CNV regions and genotypes are proposed for both microarray and HTS data. However, only a few approaches focus on haplotyping of CNV loci. RESULTS: We propose a novel approach to infer copy unit alleles and their numbers in each sample simultaneously from population-scale HTS data by variational Bayesian inference on a generative probabilistic model inspired by latent Dirichlet allocation, which is a well studied model for document classification problems. In simulation studies, we evaluated concordance between inferred and true copy unit alleles for lower-, middle-, and higher-copy number dataset, in which precision and recall were ≥ 0.9 for data with mean coverage ≥ 10× per copy unit. We also applied the approach to HTS data of 1123 samples at highly variable salivary amylase gene locus and a pseudogene locus, and confirmed consistency of the estimated alleles within samples belonging to a trio of CEPH/Utah pedigree 1463 with 11 offspring. CONCLUSIONS: Our proposed approach enables detailed analysis of copy number variations, such as association study between copy unit alleles and phenotypes or biological features including human diseases.