ERBB4 Promoter Polymorphism Is Associated with Poor Distant Disease-Free Survival in High-Risk Early Breast Cancer

A number of genetic variants have been linked to increased risk of breast cancer. Little is, however, known about the prognostic significance of hereditary factors. Here, we investigated the frequency and prognostic significance of two ERBB4 promoter region variants, −782G>T (rs62626348) and −815A>T (rs62626347), in a cohort of 1010 breast cancer patients. The frequency of nine previously described somatic ERBB4 kinase domain mutations was also analyzed. Clinical material used in the study consisted of samples from the phase III, adjuvant, FinHer breast cancer trial involving 1010 women. Tumor DNA samples were genotyped for ERBB4 variants and somatic mutations using matrix-assisted laser desorption ionization/time of flight mass spectrometry. Paraffin-embedded tumor sections from all patients were immunohistochemically stained for ErbB4 expression. Association of ERBB4 genotype to distant disease-free survival (DDFS) was assessed using Kaplan-Meier and Cox regression analyses. Genotyping was successful for 91–93% of the 1010 samples. Frequencies observed for the ERBB4 variants were 2.5% and 1.3% for −782G>T and −815A>T, respectively. Variant −815A>T was significantly associated with poor survival (HR  = 2.86 [95% CI 1.15–6.67], P = 0.017). In contrast, variant −782G>T was associated with well-differentiated cancer (P = 0.019). Two (0.2%) ERBB4 kinase domain mutations were found, both of which have previously been shown to be functional and promote cancer cell growth in vitro. These data present the germ-line ERBB4 variant −815A>T as a novel prognostic marker in high-risk early breast cancer and indicate the presence of rare but potentially oncogenic somatic ERBB4 mutations in breast cancer.     

Mutation Screening of 1,237 Cancer Genes across Six Model Cell Lines of Basal-Like Breast Cancer

Basal-like breast cancer is an aggressive subtype generally characterized as poor prognosis and lacking the expression of the three most important clinical biomarkers, estrogen receptor, progesterone receptor, and HER2. Cell lines serve as useful model systems to study cancer biology in vitro and in vivo. We performed mutational profiling of six basal-like breast cancer cell lines (HCC38, HCC1143, HCC1187, HCC1395, HCC1954, and HCC1937) and their matched normal lymphocyte DNA using targeted capture and next-generation sequencing of 1,237 cancer-associated genes, including all exons, UTRs and upstream flanking regions. In total, 658 somatic variants were identified, of which 378 were non-silent (average 63 per cell line, range 37–146) and 315 were novel (not present in the Catalogue of Somatic Mutations in Cancer database; COSMIC). 125 novel mutations were confirmed by Sanger sequencing (59 exonic, 48 3’UTR and 10 5’UTR, 1 splicing), with a validation rate of 94% of high confidence variants. Of 36 mutations previously reported for these cell lines but not detected in our exome data, 36% could not be detected by Sanger sequencing. The base replacements C/G>A/T, C/G>G/C, C/G>T/A and A/T>G/C were significantly more frequent in the coding regions compared to the non-coding regions (OR 3.2, 95% CI 2.0–5.3, P<0.0001; OR 4.3, 95% CI 2.9–6.6, P<0.0001; OR 2.4, 95% CI 1.8–3.1, P<0.0001; OR 1.8, 95% CI 1.2–2.7, P = 0.024, respectively). The single nucleotide variants within the context of T[C]T/A[G]A and T[C]A/T[G]A were more frequent in the coding than in the non-coding regions (OR 3.7, 95% CI 2.2–6.1, P<0.0001; OR 3.8, 95% CI 2.0–7.2, P = 0.001, respectively). Copy number estimations were derived from the targeted regions and correlated well to Affymetrix SNP array copy number data (Pearson correlation 0.82 to 0.96 for all compared cell lines; P<0.0001). These mutation calls across 1,237 cancer-associated genes and identification of novel variants will aid in the design and interpretation of biological experiments using these six basal-like breast cancer cell lines.     

A Bayesian ensemble approach with a disease gene network predicts damaging effects of missense variants of human cancers

Large-scale sequencing of cancer genomes has revealed many novel mutations and inter-tumoral heterogeneity. Therefore, prioritizing variants according to their potential deleterious effects has become essential. We constructed a disease gene network and proposed a Bayesian ensemble approach that integrates diverse sources to predict the functional effects of missense variants. We analyzed 23,336 missense disease mutations and 36,232 neutral polymorphisms of 12,039 human proteins. The results showed successful improvement of prediction accuracy in both sensitivity and specificity, and we demonstrated the utility of the method by applying it to somatic mutations obtained from colorectal and breast cancer cell lines. The candidate genes with predicted deleterious mutations as well as known cancer genes were significantly enriched in many KEGG pathways related to carcinogenesis, supporting genetic homogeneity of cancer at the pathway level. The breast cancer-specific network increased the prediction accuracy for breast cancer mutations. This study provides a ranked list of deleterious mutations and candidate cancer genes and suggests that mutations affecting cancer may occur in important pathways and should be interpreted on the phenotype-related network or pathway. A disease gene network may be of value in predicting functional effects of novel disease-specific mutations.     

Different Biochemical Properties Explain Why Two Equivalent Gα Subunit Mutants Cause Unrelated Diseases

There is an increasing number of disease-associated Gα mutations identified from genome-wide sequencing campaigns or targeted efforts. Albright''s Hereditary Osteodystrophy (AHO) was the first inherited disease associated with loss-of-function mutations in a G protein (Gαs) and other studies revealed gain-of-function Gα mutations in cancer. Here we attempted to solve the apparent quandary posed by the fact that the same mutation in two different G proteins appeared associated with both AHO and cancer. We first confirmed the presence of an inherited Gαs-R265H mutation from a previously described clinical case report of AHO. This mutation is structurally analogous to Gαo-R243H, an oncogenic mutant with increased activity in vitro and in cells due to rapid nucleotide exchange. We found that, contrary to Gαo-R243H, Gαs-R265H activity is compromised due to greatly impaired nucleotide binding in vitro and in cells. We obtained equivalent results when comparing another AHO mutation in Gαs (D173N) with a counterpart cancer mutation in Gαo (D151N). Gαo-R243H binds nucleotides efficiently under steady-state conditions but releases GDP much faster than the WT protein, suggesting diminished affinity for the nucleotide. These results indicate that the same disease-linked mutation in two different G proteins affects a common biochemical feature (nucleotide affinity) but to a different grade depending on the G protein (mild decrease for Gαo and severe for Gαs). We conclude that Gαs-R265H has dramatically impaired nucleotide affinity leading to the loss-of-function in AHO whereas Gαo-R243H has a mild decrease in nucleotide affinity that causes rapid nucleotide turnover and subsequent hyperactivity in cancer.     

Novel molecular and computational methods improve the accuracy of insertion site analysis in Sleeping Beauty-induced tumors

The recent development of the Sleeping Beauty (SB) system has led to the development of novel mouse models of cancer. Unlike spontaneous models, SB causes cancer through the action of mutagenic transposons that are mobilized in the genomes of somatic cells to induce mutations in cancer genes. While previous methods have successfully identified many transposon-tagged mutations in SB-induced tumors, limitations in DNA sequencing technology have prevented a comprehensive analysis of large tumor cohorts. Here we describe a novel method for producing genetic profiles of SB-induced tumors using Illumina sequencing. This method has dramatically increased the number of transposon-induced mutations identified in each tumor sample to reveal a level of genetic complexity much greater than previously appreciated. In addition, Illumina sequencing has allowed us to more precisely determine the depth of sequencing required to obtain a reproducible signature of transposon-induced mutations within tumor samples. The use of Illumina sequencing to characterize SB-induced tumors should significantly reduce sampling error that undoubtedly occurs using previous sequencing methods. As a consequence, the improved accuracy and precision provided by this method will allow candidate cancer genes to be identified with greater confidence. Overall, this method will facilitate ongoing efforts to decipher the genetic complexity of the human cancer genome by providing more accurate comparative information from Sleeping Beauty models of cancer.     

Multigene Panel Sequencing Reveals Cancer-Specific and Common Somatic Mutations in Colorectal Cancer Patients: An Egyptian Experience

This study aims at identifying common pathogenic somatic mutations at different stages of colorectal carcinogenesis in Egyptian patients. Our cohort included colonoscopic biopsies collected from 120 patients: 20 biopsies from patients with inflammatory bowel disease, 38 from colonic polyp patients, and 62 from patients with colorectal cancer. On top of this, the cohort included 20 biopsies from patients with non-specific mild to moderated colitis. Targeted DNA sequencing using a customized gene panel of 96 colorectal related genes running on the Ion Torrent NGS technology was used to process the samples. Our results revealed that 69% of all cases harbored at least one somatic mutation. Fifty-seven genes were found to carry 232 somatic non-synonymous variants. The most frequently pathogenic somatic mutations were localized in TP53, APC, KRAS, and PIK3CA. In total, 16 somatic mutations were detected in the CRC group and in either the IBD or CP group. In addition, our data showed that 51% of total somatic variants were CRC-specific variants. The average number of CRC-specific variants per sample is 2.4. The top genes carrying CRC-specific mutations are APC, TP53, PIK3CA, FBXW7, ATM, and SMAD4. It seems obvious that TP53 and APC genes were the most affected genes with somatic mutations in all groups. Of interest, 85% and 28% of the APC and TP53 deleterious somatic mutations were located in Exon 14 and Exon 3, respectively. Besides, 37% and 28% of the total somatic mutations identified in APC and TP53 were CRC-specific variants, respectively. Moreover, we identified that, in 29 somatic mutations in 21 genes, their association with CRC patients was unprecedented. Ten detected variants were likely to be novel: six in PIK3CA and four variants in FBXW7. The detected P53, Wnt/βcatenin, Angiogenesis, EGFR, TGF-β and Interleukin signaling pathways were the most altered pathways in 22%, 16%, 12%, 10%, 9% and 9% of the CRC patients, respectively. These results would contribute to a better understanding of the colorectal cancer and in introducing personalized therapies for Egyptian CRC patients.     

Cancer Cell Invasion in Three-dimensional Collagen Is Regulated Differentially by Gα13 Protein and Discoidin Domain Receptor 1-Par3 Protein Signaling

Cancer cells can invade in three-dimensional collagen as single cells or as a cohesive group of cells that require coordination of cell-cell junctions and the actin cytoskeleton. To examine the role of Gα₁₃, a G₁₂ family heterotrimeric G protein, in regulating cellular invasion in three-dimensional collagen, we established a novel method to track cell invasion by membrane type 1 matrix metalloproteinase-expressing cancer cells. We show that knockdown of Gα₁₃ decreased membrane type 1 matrix metalloproteinase-driven proteolytic invasion in three-dimensional collagen and enhanced E-cadherin-mediated cell-cell adhesion. E-cadherin knockdown reversed Gα₁₃ siRNA-induced cell-cell adhesion but failed to reverse the effect of Gα₁₃ siRNA on proteolytic invasion. Instead, concurrent knockdown of E-cadherin and Gα₁₃ led to an increased number of single cells rather than groups of cells. Significantly, knockdown of discoidin domain receptor 1 (DDR1), a collagen-binding protein that also co-localizes to cell-cell junctions, reversed the effects of Gα₁₃ knockdown on cell-cell adhesion and proteolytic invasion in three-dimensional collagen. Knockdown of the polarity protein Par3, which can function downstream of DDR1, also reversed the effects of Gα₁₃ knockdown on cell-cell adhesion and proteolytic invasion in three-dimensional collagen. Overall, we show that Gα₁₃ and DDR1-Par3 differentially regulate cell-cell junctions and the actin cytoskeleton to mediate invasion in three-dimensional collagen.     

Variants of the PPARD Gene and Their Clinicopathological Significance in Colorectal Cancer

Background

Peroxisome proliferator-activated receptor delta (PPARD) is nuclear hormone receptor involved in colorectal cancer (CRC) differentiation and progression. The purpose of this study was to determine prevalence and spectrum of variants in the PPARD gene in CRC, and their contribution to clinicopathological endpoints.

Methods and Findings

Direct sequencing of the PPARD gene was performed in 303 primary tumors, in blood samples from 50 patients with ≥3 affected first-degree relatives, 50 patients with 2 affected first-degree relatives, 50 sporadic patients, 360 healthy controls, and in 6 colon cancer cell lines. Mutation analysis revealed 22 different transversions, 7 of them were novel. Three of all variants were somatic (c.548A>G, p.Y183C, c.425-9C>T, and c.628-16G>A). Two missense mutations (p.Y183C and p.R258Q) were pathogenic using in silico predictive program. Five recurrent variants were detected in/adjacent to the exons 4 (c.1-87T>C, c.1-67G>A, c.130+3G>A, and c.1-101-8C>T) and exon 7 (c.489T>C). Variant c.489C/C detected in tumors was correlated to worse differentiation (P = 0.0397).

Conclusions

We found 7 novel variants among 22 inherited or acquired PPARD variants. Somatic and/or missense variants detected in CRC patients are rare but indicate the clinical importance of the PPARD gene.     

Single Circulating-Tumor-Cell-Targeted Sequencing to Identify Somatic Variants in Liquid Biopsies in Non-Small-Cell Lung Cancer Patients

Non-small-cell lung cancer (NSCLC) accounts for most cancer-related deaths worldwide. Liquid biopsy by a blood draw to detect circulating tumor cells (CTCs) is a tool for molecular profiling of cancer using single-cell and next-generation sequencing (NGS) technologies. The aim of the study was to identify somatic variants in single CTCs isolated from NSCLC patients by targeted NGS. Thirty-one subjects (20 NSCLC patients, 11 smokers without cancer) were enrolled for blood draws (7.5 mL). CTCs were identified by immunofluorescence, individually retrieved, and DNA-extracted. Targeted NGS was performed to detect somatic variants (single-nucleotide variants (SNVs) and insertions/deletions (Indels)) across 65 oncogenes and tumor suppressor genes. Cancer-associated variants were classified using OncoKB database. NSCLC patients had significantly higher CTC counts than control smokers (p = 0.0132; Mann–Whitney test). Analyzing 23 CTCs and 13 white blood cells across seven patients revealed a total of 644 somatic variants that occurred in all CTCs within the same subject, ranging from 1 to 137 per patient. The highest number of variants detected in ≥1 CTC within a patient was 441. A total of 18/65 (27.7%) genes were highly mutated. Mutations with oncogenic impact were identified in functional domains of seven oncogenes/tumor suppressor genes (NF1, PTCH1, TP53, SMARCB1, SMAD4, KRAS, and ERBB2). Single CTC-targeted NGS detects heterogeneous and shared mutational signatures within and between NSCLC patients. CTC single-cell genomics have potential for integration in NSCLC precision oncology.     

Inherited variants in<Emphasis Type="Italic"> MYH</Emphasis> are unlikely to contribute to the risk of lung carcinoma

The base excision repair gene MYH protects against damage to DNA from reactive oxygen species, which are commonly found in cigarette smoke. Inherited mutations in MYH predispose to colorectal adenomas and carcinomas that show a characteristic pattern of somatic G:CT:A mutations in the APC gene. A similar pattern of somatic mutations in the TP53 gene is reported in smoking-related lung cancers. We therefore tested whether germline changes in MYH may also contribute to the development of lung cancer by screening for variants in 276 patients with lung carcinoma and 106 normal controls. No patients harboured truncating mutations in MYH and only a single patient was a carrier for the G382D missense mutation. We identified three common coding region (V22M, Q324H and S501F) and intronic (157+30A>G, 462+35G>A and 1435–40G>C) variants, but none were over-represented in the patient samples, indicating that MYH variants are unlikely to predispose significantly to the risk of lung cancer.     

CT45A1 acts as a new proto-oncogene to trigger tumorigenesis and cancer metastasis

Cancer/testis antigen (CTA)-45 family (CT45) belongs to a new family of genes in phylogenetics and is absent in normal tissues except for testis, but is aberrantly overexpressed in various cancer types. Whether CT45 and other CTAs act as proto-oncogenes has not been determined. Using breast cancer as a model, we found that CT45A1, a representative CT45 family member, alone had a weak tumorigenic effect. However, its neoplastic potency was greatly enhanced in the presence of growth factors. Overexpression of CT45A1 in breast cancer cells markedly upregulated various oncogenic and metastatic genes, constitutively activated ERK and CREB signaling pathways, promoted epithelial–mesenchymal transition, and increased cell stemness, tumorigenesis, invasion, and metastasis, whereas silencing CT45A1 significantly reduced cancer cell migration and invasion. We propose that CT45A1 functions as a novel proto-oncogene to trigger oncogenesis and metastasis. CT45A1 and other CT45 members are therefore excellent targets for anticancer drug discovery and targeted tumor therapy, and valuable genes in the study of a molecular phylogenetic tree.     

Rare Variants in PLXNA4 and Parkinson’s Disease

Approximately 20% of individuals with Parkinson’s disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset familial PD followed by frequency assessment in 975 PD cases and 1014 ethnically-matched controls and linkage analysis to identify potentially causal variants. Based on the predicted penetrance and the frequencies, a variant in PLXNA4 proved to be the best candidate and PLXNA4 was screened for additional variants in 862 PD cases and 940 controls, revealing an excess of rare non-synonymous coding variants in PLXNA4 in individuals with PD. Although we cannot conclude that the variant in PLXNA4 is indeed the causative variant, these findings are interesting in the light of a surfacing role of axonal guidance mechanisms in neurodegenerative disorders but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.     

Prediction and Prioritization of Rare Oncogenic Mutations in the Cancer Kinome Using Novel Features and Multiple Classifiers

Cancer is a genetic disease that develops through a series of somatic mutations, a subset of which drive cancer progression. Although cancer genome sequencing studies are beginning to reveal the mutational patterns of genes in various cancers, identifying the small subset of “causative” mutations from the large subset of “non-causative” mutations, which accumulate as a consequence of the disease, is a challenge. In this article, we present an effective machine learning approach for identifying cancer-associated mutations in human protein kinases, a class of signaling proteins known to be frequently mutated in human cancers. We evaluate the performance of 11 well known supervised learners and show that a multiple-classifier approach, which combines the performances of individual learners, significantly improves the classification of known cancer-associated mutations. We introduce several novel features related specifically to structural and functional characteristics of protein kinases and find that the level of conservation of the mutated residue at specific evolutionary depths is an important predictor of oncogenic effect. We consolidate the novel features and the multiple-classifier approach to prioritize and experimentally test a set of rare unconfirmed mutations in the epidermal growth factor receptor tyrosine kinase (EGFR). Our studies identify T725M and L861R as rare cancer-associated mutations inasmuch as these mutations increase EGFR activity in the absence of the activating EGF ligand in cell-based assays.     

Allelic Selection of Amplicons in Glioblastoma Revealed by Combining Somatic and Germline Analysis

Cancer is a disease driven by a combination of inherited risk alleles coupled with the acquisition of somatic mutations, including amplification and deletion of genomic DNA. Potential relationships between the inherited and somatic aspects of the disease have only rarely been examined on a genome-wide level. Applying a novel integrative analysis of SNP and copy number measurements, we queried the tumor and normal-tissue genomes of 178 glioblastoma patients from the Cancer Genome Atlas project for preferentially amplified alleles, under the hypothesis that oncogenic germline variants will be selectively amplified in the tumor environment. Selected alleles are revealed by allelic imbalance in amplification across samples. This general approach is based on genetic principles and provides a method for identifying important tumor-related alleles. We find that SNP alleles that are most significantly overrepresented in amplicons tend to occur in genes involved with regulation of kinase and transferase activity, and many of these genes are known contributors to gliomagenesis. The analysis also implicates variants in synapse genes. By incorporating gene expression data, we demonstrate synergy between preferential allelic amplification and expression in DOCK4 and EGFR. Our results support the notion that combining germline and tumor genetic data can identify regions relevant to cancer biology.     

Low immunogenicity of common cancer hot spot mutations resulting in false immunogenic selection signals

Cancer is driven by somatic mutations that result in a cellular fitness advantage. This selective advantage is expected to be counterbalanced by the immune system when these driver mutations simultaneously lead to the generation of neoantigens, novel peptides that are presented at the cancer cell membrane via HLA molecules from the MHC complex. The presentability of these peptides is determined by a patient’s MHC genotype and it has been suggested that this results in MHC genotype-specific restrictions of the oncogenic mutational landscape. Here, we generated a set of virtual patients, each with an identical and prototypical MHC genotype, and show that the earlier reported HLA affinity differences between observed and unobserved mutations are unrelated to MHC genotype variation. We demonstrate how these differences are secondary to high frequencies of 13 hot spot driver mutations in 6 different genes. Several oncogenic mechanisms were identified that lower the peptides’ HLA affinity, including phospho-mimicking substitutions in BRAF, destabilizing tyrosine mutations in TP53 and glycine-rich mutational contexts in the GTP-binding KRAS domain. In line with our earlier findings, our results emphasize that HLA affinity predictions are easily misinterpreted when studying immunogenic selection processes.     

Ron Knockdown and Ron Monoclonal Antibody IMC-RON8 Sensitize Pancreatic Cancer to Histone Deacetylase Inhibitors (HDACi)

Recepteur d’origine nantais (Ron) is overexpressed in a panel of pancreatic cancer cells and tissue samples from pancreatic cancer patients. Ron can be activated by its ligand macrophage stimulating protein (MSP), thereby activating oncogenic signaling pathways. Crosstalk between Ron and EGFR, c-Met, or IGF-1R may provide a mechanism underlying drug resistance. Thus, targeting Ron may represent a novel therapeutic strategy. IMC-RON8 is the first Ron monoclonal antibody (mAb) entering clinical trial for targeting Ron overexpression. Our studies show IMC-RON8 downmodulated Ron expression in pancreatic cancer cells and significantly blocked MSP-stimulated Ron activation, downstream Akt and ERK phosphorylation, and survivin mRNA expression. IMC-RON8 hindered MSP-induced cell migration and reduced cell transformation. Histone deacetylase inhibitors (HDACi) are reported to target expression of various genes through modification of nucleosome histones and non-histone proteins. Our work shows HDACi TSA and Panobinostat (PS) decreased Ron mRNA and protein expression in pancreatic cancer cells. PS also reduced downstream signaling of pAkt, survivin, and XIAP, as well as enhanced cell apoptosis. Interestingly, PS reduced colony formation in Ron knockdown cells to a greater extent than Ron scramble control cells in colony formation and soft agarose assays. IMC-RON8 could also sensitize pancreatic cancer cells to PS, as reflected by reduced colony numbers and size in combination treatment with IMC-RON8 and PS compared to single treatment alone. The co-treatment further reduced Ron expression and pAkt, and increased PARP cleavage compared to either treatment alone. This study suggests the potential for a novel combination approach which may ultimately be of value in treatment of pancreatic cancer.     

Genetic characteristics and prognostic implications of m1A regulators in pancreatic cancer

Studies have identified the methylation of N1 adenosine (m1A), an RNA modification, playing an important role in the progression of the tumorigenesis. The present study aimed to analyze the genetic characteristics and prognostic value of m1A regulators in pancreatic cancer. In the present study, data on gene mutations, single-nucleotide variants (SNVs), and copy number variation (CNV) were obtained from 363 patients with pancreatic cancer in the Cancer Genome Atlas (TCGA) database, and survival analysis was performed using the logarithmic rank test and Cox regression model. The chi-squared test was used to examine the relationship between the changes in m1A regulatory factors and clinicopathological characteristics. And we used ICGC database to verify the reliability of prognostic markers. The results show that changes in m1A-regulating genes are related to clinical stage and that the expression of some m1A-regulating genes is positively correlated with CNV. In addition, the low expression of the ‘eraser’ gene ALKBH1 is related to the poor prognosis of patients with pancreatic cancer, and its expression level has important clinical significance for patients with pancreatic adenocarcinoma (PAAD). Mechanistically, ALKBH1 may participate in the occurrence and development of pancreatic cancer through mTOR and ErbB signaling pathway. The expression of m1A-regulating genes can be used as a prognostic marker for pancreatic cancer. These findings provide valuable clues for us to understand the epigenetics of m1A in pancreatic cancer.     

Whole Exome Sequencing Reveals Homozygous Mutations in RAI1, OTOF,and SLC26A4 Genes Associated with Nonsyndromic Hearing Loss in Altaian Families (South Siberia)

Hearing loss (HL) is one of the most common sensorineural disorders and several dozen genes contribute to its pathogenesis. Establishing a genetic diagnosis of HL is of great importance for clinical evaluation of deaf patients and for estimating recurrence risks for their families. Efforts to identify genes responsible for HL have been challenged by high genetic heterogeneity and different ethnic-specific prevalence of inherited deafness. Here we present the utility of whole exome sequencing (WES) for identifying candidate causal variants for previously unexplained nonsyndromic HL of seven patients from four unrelated Altaian families (the Altai Republic, South Siberia). The WES analysis revealed homozygous missense mutations in three genes associated with HL. Mutation c.2168A>G (SLC26A4) was found in one family, a novel mutation c.1111G>C (OTOF) was revealed in another family, and mutation c.5254G>A (RAI1) was found in two families. Sanger sequencing was applied for screening of identified variants in an ethnically diverse cohort of other patients with HL (n = 116) and in Altaian controls (n = 120). Identified variants were found only in patients of Altaian ethnicity (n = 93). Several lines of evidences support the association of homozygosity for discovered variants c.5254G>A (RAI1), c.1111C>G (OTOF), and c.2168A>G (SLC26A4) with HL in Altaian patients. Local prevalence of identified variants implies possible founder effect in significant number of HL cases in indigenous population of the Altai region. Notably, this is the first reported instance of patients with RAI1 missense mutation whose HL is not accompanied by specific traits typical for Smith-Magenis syndrome. Presumed association of RAI1 gene variant c.5254G>A with isolated HL needs to be proved by further experimental studies.     

Mitochondrial DNA Polymerase POLG1 Disease Mutations and Germline Variants Promote Tumorigenic Properties

Germline mutations in mitochondrial DNA polymerase gamma (POLG1) induce mitochondrial DNA (mtDNA) mutations, depletion, and decrease oxidative phosphorylation. Earlier, we identified somatic mutations in POLG1 and the contribution of these mutations in human cancer. However, a role for germline variations in POLG1 in human cancers is unknown. In this study, we examined a role for disease associated germline variants of POLG1, POLG1 gene expression, copy number variation and regulation in human cancers. We analyzed the mutations, expression and copy number variation in POLG1 in several cancer databases and validated the analyses in primary breast tumors and breast cancer cell lines. We discovered 5-aza-2''-deoxycytidine led epigenetic regulation of POLG1, mtDNA-encoded genes and increased mitochondrial respiration. We conducted comprehensive race based bioinformatics analyses of POLG1 gene in more than 33,000 European-Americans and 5,000 African-Americans. We identified a mitochondrial disease causing missense variation in polymerase domain of POLG1 protein at amino acid 1143 (E1143G) to be 25 times more prevalent in European-Americans (allele frequency 0.03777) when compared to African-American (allele frequency 0.00151) population. We identified T251I and P587L missense variations in exonuclease and linker region of POLG1 also to be more prevalent in European-Americans. Expression of these variants increased glucose consumption, decreased ATP production and increased matrigel invasion. Interestingly, conditional expression of these variants revealed that matrigel invasion properties conferred by these germline variants were reversible suggesting a role of epigenetic regulators. Indeed, we identified a set of miRNA whose expression was reversible after variant expression was turned off. Together, our studies demonstrate altered genetic and epigenetic regulation of POLG1 in human cancers and suggest a role for POLG1 germline variants in promoting tumorigenic properties.