Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

Associations of LRP5 and MTHFR Gene Variants with Osteoarthritis Prevalence in Elderly Women: A Japanese Cohort Survey Randomly Sampled from a Basic Resident Registry

Objective

Osteoarthritis (OA) is a common and degenerative joint disorder in the elderly. A greater importance of understanding the relationship between genetic factors and OA prevalence has emerged with population aging. We therefore investigated the associations of several bone disease-related genetic variants with the prevalence of OA and osteoporosis in Japanese elderly women from the Obuse study cohort, which was randomly sampled from a basic town resident registry.

Methods and Results

In total, 206 female participants (mean ± standard deviation age: 69.7 ± 11.0 years) who completed OA, bone mineral density, and genotype assessments were included. The number of patients diagnosed as having knee/hip OA and osteoporosis was 59 (28.6%) and 30 (14.6%), respectively. Fisher’s exact testing revealed significant relationships between the minor T allele of LDL receptor related protein 5 (LRP5) rs3736228 and the prevalence of knee/hip OA and osteoporosis. The respective odds ratios (ORs) of the TT genotype for knee/hip OA and osteoporosis were 7.28 (95% confidence interval [CI] 2.22–28.08) and 5.24 (95% CI 0.95–26.98). An additional subgroup analysis for knee OA revealed that the frequency of the common C allele of methylenetetrahydrofolate reductase (MTHFR) rs1801133 had a statistically significant protective association with the prevalence of knee OA (OR 0.58, 95% CI 0.35–0.97).

Conclusion

In sum, the present study demonstrated significant associations of LRP5 rs3736228 and MTHFR rs1801133 with knee/hip OA and osteoporosis prevalences and knee OA prevalence, respectively, in Japanese elderly women. These results will help further the understanding of OA pathogenesis and related genetic risk factors.

Exome Sequencing Reveals De Novo Variants in Congenital Scoliosis

Congenital scoliosis (CS) is a lateral curvature of the spine characterized by the presence of vertebral anomalies. Pathogenic genetic variants in the TBX6 gene are one of the causes of CS. However, since many clinically diagnosed cases of CS are without known TBX6 gene variations, this study aims to uncover new genes related to disease susceptibility of CS by exome sequencing (ES). This study employed ES in a cohort of 5 Japanese patients with CS and their healthy parents or a sister for a total of 16 samples among 5 families. Variant interpretation was performed using SIFT, PolyPhen-2, Mutation Taster, and CADD. Four de novo variants were identified by ES and confirmed by Sanger sequencing: 1 frameshift variant (SHISA3) and 3 missense variants (AGBL5, HDAC4, and PDE2A). ES also uncovered 1 homozygous variant in the MOCOS gene. All of these variants were predicted to be deleterious by SIFT, PolyPhen-2, Mutation Taster, and/or CADD. The number of de novo variants identified in this study was exactly what would be expected by chance. Additional functional studies or gathering matched patients using Gene Matcher are needed.

Observation of p.R4810K, a Polymorphism of the Mysterin Gene, the Susceptibility Gene for Moyamoya Disease, in Two Female Japanese Diabetic Patients with Familial Partial Lipodystrophy 1

Mysterin, which was recently shown to play an important role in maintaining cellular fat storage, has been identified to be the susceptibility gene for moyamoya disease (MMD). We encountered some female Japanese patients with partial lipodystrophy and MMD-like vascular lesions. This prompted us to examine whether mysterin variants may be present in these patients. We identified a mysterin variant, p.R4810K in two patients with MMD-like vascular lesions, who may fit the category of familial partial lipodystrophy (FPLD) 1. Our cases suggest the possibility that p.R4810K, in addition to atherogenic risk factors, might thus play a role in the development of atherosclerotic lesions in patients with FPLD1 and p.R4810K.

Olfactomedin-like protein OLFML1 inhibits Hippo signaling and mineralization in osteoblasts

Congenital scoliosis is a lateral curvature of the spine that is due to the presence of vertebral anomalies. Although genetic and environmental factors are involved in the pathogenesis of congenital scoliosis, the specific cause of only a small number of individuals has been identified to date. We identified a de novo missense mutation in the olfactomedin-like 1 (OLFML1) gene by whole-exome sequencing of a patient with congenital scoliosis. Then, we carried out further functional investigation in mice. An assessment of the tissue distribution of Olfml1 revealed it to be prominently expressed in developing skeletal tissues, specifically osteoblasts. Short hairpin RNA-mediated knockdown of Olfml1 in osteoblasts induced the translocation of Yes-associated protein (YAP) transcriptional coactivator from the cytoplasm to the nucleus, which accelerated the Hippo signaling pathway to promote osteoblast mineralization. In contrast, experimentally induced gain of function of Olfml1 retained YAP in the cytoplasm. There appears to exist a novel cell-autonomous mechanism by which osteoblasts avoid excess mineralization through Olfml1. Our results also indicate that mutation of OLFML1 leads to impaired osteoblast differentiation and abnormal development of bone tissue.

Identification of Individual Cancer-Specific Somatic Mutations for Neoantigen-Based Immunotherapy of Lung Cancer

INTRODUCTION

Two strategies for selecting neoantigens as targets for non-small cell lung cancer vaccines were compared: (1) an "off-the-shelf" approach starting with shared mutations extracted from global databases and (2) a personalized pipeline using whole-exome sequencing data on each patient's tumor.

METHODS

The Catalogue of Somatic Mutations in Cancer database was used to create a list of shared missense mutations occurring in more than 1% of patients. These mutations were then assessed for predicted binding affinity to HLA alleles of 15 lung cancer patients, and potential neoantigens (pNeoAgs) for each patient were selected on this basis. In the personalized approach, pNeoAgs were selected from missense mutations detected by whole-exome sequencing of the patient's own samples.

RESULTS

The list of shared mutations included 22 missense mutations for adenocarcinoma and 18 for squamous cell carcinoma (SCC), resulting in a median of 10 off-the-shelf pNeoAgs for each adenocarcinoma (range 5-13) and 9 (range 5-12) for each SCC. In contrast, a median of 59 missense mutations were identified by whole-exome sequencing (range 33-899) in adenocarcinoma and 164.5 (range 26-232) in SCC. This resulted in a median of 46 pNeoAgs (range 13-659) for adenocarcinoma and 95.5 (range 10-145) for SCC in the personalized set. We found that only one or two off-the-shelf pNeoAgs were included in the set of personalized pNeoAgs-and then in only three patients, with no overlap seen in the remaining 12 patients.

CONCLUSIONS

Use of an off-the-shelf pipeline is feasible but may not be satisfactory for most patients with non-small cell lung cancer. We recommend identifying personal mutations by comprehensive genome sequencing for developing neoantigen-targeted cancer immunotherapies.

PCSK5 mutation in a patient with the VACTERL association

BACKGROUND

The VACTERL association is a typically sporadic, non-random collection of congenital anomalies that includes vertebral defects, anal atresia, cardiac defects, tracheoesophageal fistula with esophageal atresia, renal anomalies, and limb abnormalities. Although several chromosomal aberrations and gene mutations have been reported as disease-causative, these findings have been sparsely replicated to date.

CASE PRESENTATION

In the present study, whole exome sequencing of a case with the VACTERL association uncovered a novel frameshift mutation in the PCSK5 gene, which has been reported as one of the causative genes for the VACTERL association. Although this mutation appears potentially pathogenic in its functional aspects, it was also carried by the healthy father. Furthermore, a database survey revealed several other deleterious variants in the PCSK5 gene in the general population.

CONCLUSIONS

Further studies are necessary to clarify the etiological role of the PCSK5 mutation in the VACTERL association.

PCDq: human protein complex database with quality index which summarizes different levels of evidences of protein complexes predicted from h-invitational protein-protein interactions integrative dataset

Background

Proteins interact with other proteins or biomolecules in complexes to perform cellular functions. Existing protein-protein interaction (PPI) databases and protein complex databases for human proteins are not organized to provide protein complex information or facilitate the discovery of novel subunits. Data integration of PPIs focused specifically on protein complexes, subunits, and their functions. Predicted candidate complexes or subunits are also important for experimental biologists.

Description

Based on integrated PPI data and literature, we have developed a human protein complex database with a complex quality index (PCDq), which includes both known and predicted complexes and subunits. We integrated six PPI data (BIND, DIP, MINT, HPRD, IntAct, and GNP_Y2H), and predicted human protein complexes by finding densely connected regions in the PPI networks. They were curated with the literature so that missing proteins were complemented and some complexes were merged, resulting in 1,264 complexes comprising 9,268 proteins with 32,198 PPIs. The evidence level of each subunit was assigned as a categorical variable. This indicated whether it was a known subunit, and a specific function was inferable from sequence or network analysis. To summarize the categories of all the subunits in a complex, we devised a complex quality index (CQI) and assigned it to each complex. We examined the proportion of consistency of Gene Ontology (GO) terms among protein subunits of a complex. Next, we compared the expression profiles of the corresponding genes and found that many proteins in larger complexes tend to be expressed cooperatively at the transcript level. The proportion of duplicated genes in a complex was evaluated. Finally, we identified 78 hypothetical proteins that were annotated as subunits of 82 complexes, which included known complexes. Of these hypothetical proteins, after our prediction had been made, four were reported to be actual subunits of the assigned protein complexes.

Conclusions

We constructed a new protein complex database PCDq including both predicted and curated human protein complexes. CQI is a useful source of experimentally confirmed information about protein complexes and subunits. The predicted protein complexes can provide functional clues about hypothetical proteins. PCDq is freely available at http://h-invitational.jp/hinv/pcdq/.

The H-Invitational Database (H-InvDB), a comprehensive annotation resource for human genes and transcripts

Here we report the new features and improvements in our latest release of the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/), a comprehensive annotation resource for human genes and transcripts. H-InvDB, originally developed as an integrated database of the human transcriptome based on extensive annotation of large sets of full-length cDNA (FLcDNA) clones, now provides annotation for 120 558 human mRNAs extracted from the International Nucleotide Sequence Databases (INSD), in addition to 54 978 human FLcDNAs, in the latest release H-InvDB_4.6. We mapped those human transcripts onto the human genome sequences (NCBI build 36.1) and determined 34 699 human gene clusters, which could define 34 057 (98.1%) protein-coding and 642 (1.9%) non-protein-coding loci; 858 (2.5%) transcribed loci overlapped with predicted pseudogenes. For all these transcripts and genes, we provide comprehensive annotation including gene structures, gene functions, alternative splicing variants, functional non-protein-coding RNAs, functional domains, predicted sub cellular localizations, metabolic pathways, predictions of protein 3D structure, mapping of SNPs and microsatellite repeat motifs, co-localization with orphan diseases, gene expression profiles, orthologous genes, protein-protein interactions (PPI) and annotation for gene families. The current H-InvDB annotation resources consist of two main views: Transcript view and Locus view and eight sub-databases: the DiseaseInfo Viewer, H-ANGEL, the Clustering Viewer, G-integra, the TOPO Viewer, Evola, the PPI view and the Gene family/group.