Exome sequencing of liver fluke-associated cholangiocarcinoma

Opisthorchis viverrini-related cholangiocarcinoma (CCA), a fatal bile duct cancer, is a major public health concern in areas endemic for this parasite. We report here whole-exome sequencing of eight O. viverrini-related tumors and matched normal tissue. We identified and validated 206 somatic mutations in 187 genes using Sanger sequencing and selected 15 genes for mutation prevalence screening in an additional 46 individuals with CCA (cases). In addition to the known cancer-related genes TP53 (mutated in 44.4% of cases), KRAS (16.7%) and SMAD4 (16.7%), we identified somatic mutations in 10 newly implicated genes in 14.8-3.7% of cases. These included inactivating mutations in MLL3 (in 14.8% of cases), ROBO2 (9.3%), RNF43 (9.3%) and PEG3 (5.6%), and activating mutations in the GNAS oncogene (9.3%). These genes have functions that can be broadly grouped into three biological classes: (i) deactivation of histone modifiers, (ii) activation of G protein signaling and (iii) loss of genome stability. This study provides insight into the mutational landscape contributing to O. viverrini-related CCA.     

Insertional mutagenesis identifies multiple networks of cooperating genes driving intestinal tumorigenesis

The evolution of colorectal cancer suggests the involvement of many genes. To identify new drivers of intestinal cancer, we performed insertional mutagenesis using the Sleeping Beauty transposon system in mice carrying germline or somatic Apc mutations. By analyzing common insertion sites (CISs) isolated from 446 tumors, we identified many hundreds of candidate cancer drivers. Comparison to human data sets suggested that 234 CIS-targeted genes are also dysregulated in human colorectal cancers. In addition, we found 183 CIS-containing genes that are candidate Wnt targets and showed that 20 CISs-containing genes are newly discovered modifiers of canonical Wnt signaling. We also identified mutations associated with a subset of tumors containing an expanded number of Paneth cells, a hallmark of deregulated Wnt signaling, and genes associated with more severe dysplasia included those encoding members of the FGF signaling cascade. Some 70 genes had co-occurrence of CIS pairs, clustering into 38 sub-networks that may regulate tumor development.     

NMNAT1 mutations cause Leber congenital amaurosis

Leber congenital amaurosis (LCA) is an infantile-onset form of inherited retinal degeneration characterized by severe vision loss. Two-thirds of LCA cases are caused by mutations in 17 known disease-associated genes (Retinal Information Network (RetNet)). Using exome sequencing we identified a homozygous missense mutation (c.25G>A, p.Val9Met) in NMNAT1 that is likely to be disease causing in two siblings of a consanguineous Pakistani kindred affected by LCA. This mutation segregated with disease in the kindred, including in three other children with LCA. NMNAT1 resides in the previously identified LCA9 locus and encodes the nuclear isoform of nicotinamide mononucleotide adenylyltransferase, a rate-limiting enzyme in nicotinamide adenine dinucleotide (NAD(+)) biosynthesis. Functional studies showed that the p.Val9Met alteration decreased NMNAT1 enzyme activity. Sequencing NMNAT1 in 284 unrelated families with LCA identified 14 rare mutations in 13 additional affected individuals. These results are the first to link an NMNAT isoform to disease in humans and indicate that NMNAT1 mutations cause LCA.     

A genome-wide scan identifies mutations in the gene encoding phosphodiesterase 11A4 (PDE11A) in individuals with adrenocortical hyperplasia

Phosphodiesterases (PDEs) regulate cyclic nucleotide levels. Increased cyclic AMP (cAMP) signaling has been associated with PRKAR1A or GNAS mutations and leads to adrenocortical tumors and Cushing syndrome. We investigated the genetic source of Cushing syndrome in individuals with adrenocortical hyperplasia that was not caused by known defects. We performed genome-wide SNP genotyping, including the adrenocortical tumor DNA. The region with the highest probability to harbor a susceptibility gene by loss of heterozygosity (LOH) and other analyses was 2q31-2q35. We identified mutations disrupting the expression of the PDE11A isoform-4 gene (PDE11A) in three kindreds. Tumor tissues showed 2q31-2q35 LOH, decreased protein expression and high cyclic nucleotide levels and cAMP-responsive element binding protein (CREB) phosphorylation. PDE11A codes for a dual-specificity PDE that is expressed in adrenal cortex and is partially inhibited by tadalafil and other PDE inhibitors; its germline inactivation is associated with adrenocortical hyperplasia, suggesting another means by which dysregulation of cAMP signaling causes endocrine tumors.     

Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy. Here, we performed high-resolution copy-number analysis on 125 HCC tumors and whole-exome sequencing on 24 of these tumors. We identified 135 homozygous deletions and 994 somatic mutations of genes with predicted functional consequences. We found new recurrent alterations in four genes (ARID1A, RPS6KA3, NFE2L2 and IRF2) not previously described in HCC. Functional analyses showed tumor suppressor properties for IRF2, whose inactivation, exclusively found in hepatitis B virus (HBV)-related tumors, led to impaired TP53 function. In contrast, inactivation of chromatin remodelers was frequent and predominant in alcohol-related tumors. Moreover, association of mutations in specific genes (RPS6KA3-AXIN1 and NFE2L2-CTNNB1) suggested that Wnt/β-catenin signaling might cooperate in liver carcinogenesis with both oxidative stress metabolism and Ras/mitogen-activated protein kinase (MAPK) pathways. This study provides insight into the somatic mutational landscape in HCC and identifies interactions between mutations in oncogene and tumor suppressor gene mutations related to specific risk factors.     

Fatal infantile cardioencephalomyopathy with COX deficiency and mutations in SCO2, a COX assembly gene.

Mammalian cytochrome c oxidase (COX) catalyses the transfer of reducing equivalents from cytochrome c to molecular oxygen and pumps protons across the inner mitochondrial membrane. Mitochondrial DNA (mtDNA) encodes three COX subunits (I-III) and nuclear DNA (nDNA) encodes ten. In addition, ancillary proteins are required for the correct assembly and function of COX (refs 2, 3, 4, 5, 6). Although pathogenic mutations in mtDNA-encoded COX subunits have been described, no mutations in the nDNA-encoded subunits have been uncovered in any mendelian-inherited COX deficiency disorder. In yeast, two related COX assembly genes, SCO1 and SCO2 (for synthesis of cytochrome c oxidase), enable subunits I and II to be incorporated into the holoprotein. Here we have identified mutations in the human homologue, SCO2, in three unrelated infants with a newly recognized fatal cardioencephalomyopathy and COX deficiency. Immunohistochemical studies implied that the enzymatic deficiency, which was most severe in cardiac and skeletal muscle, was due to the loss of mtDNA-encoded COX subunits. The clinical phenotype caused by mutations in human SCO2 differs from that caused by mutations in SURF1, the only other known COX assembly gene associated with a human disease, Leigh syndrome.     

De novo somatic mutations in components of the PI3K-AKT3-mTOR pathway cause hemimegalencephaly

De novo somatic mutations in focal areas are well documented in diseases such as neoplasia but are rarely reported in malformation of the developing brain. Hemimegalencephaly (HME) is characterized by overgrowth of either one of the two cerebral hemispheres. The molecular etiology of HME remains a mystery. The intractable epilepsy that is associated with HME can be relieved by the surgical treatment hemispherectomy, allowing sampling of diseased tissue. Exome sequencing and mass spectrometry analysis in paired brain-blood samples from individuals with HME (n = 20 cases) identified de novo somatic mutations in 30% of affected individuals in the PIK3CA, AKT3 and MTOR genes. A recurrent PIK3CA c.1633G>A mutation was found in four separate cases. Identified mutations were present in 8-40% of sequenced alleles in various brain regions and were associated with increased neuronal S6 protein phosphorylation in the brains of affected individuals, indicating aberrant activation of mammalian target of rapamycin (mTOR) signaling. Thus HME is probably a genetically mosaic disease caused by gain of function in phosphatidylinositol 3-kinase (PI3K)-AKT3-mTOR signaling.     

Mutations in the gene encoding mevalonate kinase cause hyper-IgD and periodic fever syndrome. International Hyper-IgD Study Group.

Hyperimmunoglobulinaemia D and periodic fever syndrome (HIDS; MIM 260920) is a rare, apparently monogenic, autosomal recessive disorder characterized by recurrent episodes of fever accompanied with lymphadenopathy, abdominal distress, joint involvement and skin lesions. All patients have high serum IgD values (>100 U/ml) and HIDS 'attacks' are associated with an intense acute phase reaction whose exact pathophysiology remains obscure. Two other hereditary febrile disorders have been described. Familial Mediterranean fever (MIM 249100) is an autosomal recessive disorder affecting mostly populations from the Mediterranean basin and is caused by mutations in the gene MEFV (refs 5,6). Familial Hibernian fever (MIM 142680), also known as autosomal dominant familial recurrent fever, is caused by missense mutations in the gene encoding type I tumour necrosis factor receptor. Here we perform a genome-wide search to map the HIDS gene. Haplotype analysis placed the gene at 12q24 between D12S330 and D12S79. We identified the gene MVK, encoding mevalonate kinase (MK, ATP:mevalonate 5-phosphotransferase; EC 2.7.1.36), as a candidate gene. We characterized 3 missense mutations, a 92-bp loss stemming from a deletion or from exon skipping, and the absence of expression of one allele. Functional analysis demonstrated diminished MK activity in fibroblasts from HIDS patients. Our data establish MVK as the gene responsible for HIDS.     

Multiple self-healing squamous epithelioma is caused by a disease-specific spectrum of mutations in TGFBR1

Multiple self-healing squamous epithelioma (MSSE), also known as Ferguson-Smith disease (FSD), is an autosomal-dominant skin cancer condition characterized by multiple squamous-carcinoma-like locally invasive skin tumors that grow rapidly for a few weeks before spontaneously regressing, leaving scars. High-throughput genomic sequencing of a conservative estimate (24.2 Mb) of the disease locus on chromosome 9 using exon array capture identified independent mutations in TGFBR1 in three unrelated families. Subsequent dideoxy sequencing of TGFBR1 identified 11 distinct monoallelic mutations in 18 affected families, firmly establishing TGFBR1 as the causative gene. The nature of the sequence variants, which include mutations in the extracellular ligand-binding domain and a series of truncating mutations in the kinase domain, indicates a clear genotype-phenotype correlation between loss-of-function TGFBR1 mutations and MSSE. This distinguishes MSSE from the Marfan syndrome-related disorders in which missense mutations in TGFBR1 lead to developmental defects with vascular involvement but no reported predisposition to cancer.     

CCDC39 is required for assembly of inner dynein arms and the dynein regulatory complex and for normal ciliary motility in humans and dogs

Primary ciliary dyskinesia (PCD) is an inherited disorder characterized by recurrent infections of the upper and lower respiratory tract, reduced fertility in males and situs inversus in about 50% of affected individuals (Kartagener syndrome). It is caused by motility defects in the respiratory cilia that are responsible for airway clearance, the flagella that propel sperm cells and the nodal monocilia that determine left-right asymmetry. Recessive mutations that cause PCD have been identified in genes encoding components of the outer dynein arms, radial spokes and cytoplasmic pre-assembly factors of axonemal dyneins, but these mutations account for only about 50% of cases of PCD. We exploited the unique properties of dog populations to positionally clone a new PCD gene, CCDC39. We found that loss-of-function mutations in the human ortholog underlie a substantial fraction of PCD cases with axonemal disorganization and abnormal ciliary beating. Functional analyses indicated that CCDC39 localizes to ciliary axonemes and is essential for assembly of inner dynein arms and the dynein regulatory complex.     

Mutations of PKD1 in ADPKD2 cysts suggest a pathogenic effect of trans-heterozygous mutations

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in PKD1 and PKD2. The products of these genes associate to form heteromeric complexes. Several models have been proposed to explain the mechanism of cyst formation. Here we find somatic mutations of PKD2 in 71% of ADPKD2 cysts analysed. Clonal somatic mutations of PKD1 were identified in a subset of cysts that lacked PKD2 mutations.     

Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide. We sequenced and analyzed the whole genomes of 27 HCCs, 25 of which were associated with hepatitis B or C virus infections, including two sets of multicentric tumors. Although no common somatic mutations were identified in the multicentric tumor pairs, their whole-genome substitution patterns were similar, suggesting that these tumors developed from independent mutations, although their shared etiological backgrounds may have strongly influenced their somatic mutation patterns. Statistical and functional analyses yielded a list of recurrently mutated genes. Multiple chromatin regulators, including ARID1A, ARID1B, ARID2, MLL and MLL3, were mutated in ～50% of the tumors. Hepatitis B virus genome integration in the TERT locus was frequently observed in a high clonal proportion. Our whole-genome sequencing analysis of HCCs identified the influence of etiological background on somatic mutation patterns and subsequent carcinogenesis, as well as recurrent mutations in chromatin regulators in HCCs.     

Mutations in the region encoding the von Willebrand factor A domain of matrilin-3 are associated with multiple epiphyseal dysplasia

Multiple epiphyseal dysplasia (MED) is a relatively mild and clinically variable osteochondrodysplasia, primarily characterized by delayed and irregular ossification of the epiphyses and early-onset osteoarthritis. Mutations in the genes encoding cartilage oligomeric matrix protein (COMP) and type IX collagen (COL9A2 and COL9A3) have previously been shown to cause different forms of MED (refs. 4-13). These dominant forms of MED (EDM1-3) are caused by mutations in the genes encoding structural proteins of the cartilage extracellular matrix (ECM); these proteins interact with high affinity in vitro. A recessive form of MED (EDM4) has also been reported; it is caused by a mutation in the diastrophic dysplasia sulfate transporter gene (SLC26A). A genomewide screen of family with autosomal-dominant MED not linked to the EDM1-3 genes provides significant genetic evidence for a MED locus on the short arm of chromosome 2 (2p24-p23), and a search for candidate genes identified MATN3 (ref. 18), encoding matrilin-3, within the critical region. Matrilin-3 is an oligomeric protein that is present in the cartilage ECM. We have identified two different missense mutations in the exon encoding the von Willebrand factor A (vWFA) domain of matrilin-3 in two unrelated families with MED (EDM5). These are the first mutations to be identified in any of the genes encoding the matrilin family of proteins and confirm a role for matrilin-3 in the development and homeostasis of cartilage and bone.     

Mutations in ATP2A2, encoding a Ca2+ pump, cause Darier disease

Darier disease (DD) is an autosomal-dominant skin disorder characterized by loss of adhesion between epidermal cells (acantholysis) and abnormal keratinization. Recently we constructed a 2.4-Mb, P1-derived artificial chromosome contig spanning the DD candidate region on chromosome 12q23-24.1. After screening several genes that mapped to this region, we identified mutations in the ATP2A2 gene, which encodes the sarco/endoplasmic reticulum Ca2(+)-ATPase type 2 isoform (SERCA2) and is highly expressed in keratinocytes. Thirteen mutations were identified, including frameshift deletions, in-frame deletions or insertions, splice-site mutations and non-conservative missense mutations in functional domains. Our results demonstrate that mutations in ATP2A2 cause DD and disclose a role for this pump in a Ca(2+)-signalling pathway regulating cell-to-cell adhesion and differentiation of the epidermis.     

BAP1 loss defines a new class of renal cell carcinoma

The molecular pathogenesis of renal cell carcinoma (RCC) is poorly understood. Whole-genome and exome sequencing followed by innovative tumorgraft analyses (to accurately determine mutant allele ratios) identified several putative two-hit tumor suppressor genes, including BAP1. The BAP1 protein, a nuclear deubiquitinase, is inactivated in 15% of clear cell RCCs. BAP1 cofractionates with and binds to HCF-1 in tumorgrafts. Mutations disrupting the HCF-1 binding motif impair BAP1-mediated suppression of cell proliferation but not deubiquitination of monoubiquitinated histone 2A lysine 119 (H2AK119ub1). BAP1 loss sensitizes RCC cells in vitro to genotoxic stress. Notably, mutations in BAP1 and PBRM1 anticorrelate in tumors (P = 3 × 10(-5)), and combined loss of BAP1 and PBRM1 in a few RCCs was associated with rhabdoid features (q = 0.0007). BAP1 and PBRM1 regulate seemingly different gene expression programs, and BAP1 loss was associated with high tumor grade (q = 0.0005). Our results establish the foundation for an integrated pathological and molecular genetic classification of RCC, paving the way for subtype-specific treatments exploiting genetic vulnerabilities.     

A de novo paradigm for mental retardation

The per-generation mutation rate in humans is high. De novo mutations may compensate for allele loss due to severely reduced fecundity in common neurodevelopmental and psychiatric diseases, explaining a major paradox in evolutionary genetic theory. Here we used a family based exome sequencing approach to test this de novo mutation hypothesis in ten individuals with unexplained mental retardation. We identified and validated unique non-synonymous de novo mutations in nine genes. Six of these, identified in six different individuals, are likely to be pathogenic based on gene function, evolutionary conservation and mutation impact. Our findings provide strong experimental support for a de novo paradigm for mental retardation. Together with de novo copy number variation, de novo point mutations of large effect could explain the majority of all mental retardation cases in the population.     

Use of human tissue to assess the oncogenic activity of melanoma-associated mutations

Multiple genetic alterations occur in melanoma, a lethal skin malignancy of increasing incidence. These include mutations that activate Ras and two of its effector cascades, Raf and phosphoinositide 3-kinase (PI3K). Induction of Ras and Raf can be caused by active N-Ras and B-Raf mutants as well as by gene amplification. Activation of PI3K pathway components occurs by PTEN loss and by AKT3 amplification. Melanomas also commonly show impairment of the p16(INK4A)-CDK4-Rb and ARF-HDM2-p53 tumor suppressor pathways. CDKN2A mutations can produce p16(INK4A) and ARF protein loss. Rb bypass can also occur through activating CDK4 mutations as well as by CDK4 amplification. In addition to ARF deletion, p53 pathway disruption can result from dominant negative TP53 mutations. TERT amplification also occurs in melanoma. The extent to which these mutations can induce human melanocytic neoplasia is unknown. Here we characterize pathways sufficient to generate human melanocytic neoplasia and show that genetically altered human tissue facilitates functional analysis of mutations observed in human tumors.     

HRPT2, encoding parafibromin,is mutated in hyperparathyroidism-jaw tumor syndrome

We report here the identification of a gene associated with the hyperparathyroidism-jaw tumor (HPT-JT) syndrome. A single locus associated with HPT-JT (HRPT2) was previously mapped to chromosomal region 1q25-q32. We refined this region to a critical interval of 12 cM by genotyping in 26 affected kindreds. Using a positional candidate approach, we identified thirteen different heterozygous, germline, inactivating mutations in a single gene in fourteen families with HPT-JT. The proposed role of HRPT2 as a tumor suppressor was supported by mutation screening in 48 parathyroid adenomas with cystic features, which identified three somatic inactivating mutations, all located in exon 1. None of these mutations were detected in normal controls, and all were predicted to cause deficient or impaired protein function. HRPT2 is a ubiquitously expressed, evolutionarily conserved gene encoding a predicted protein of 531 amino acids, for which we propose the name parafibromin. Our findings suggest that HRPT2 is a tumor-suppressor gene, the inactivation of which is directly involved in predisposition to HPT-JT and in development of some sporadic parathyroid tumors.     

Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodeling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.     

High-throughput oncogene mutation profiling in human cancer

Systematic efforts are underway to decipher the genetic changes associated with tumor initiation and progression. However, widespread clinical application of this information is hampered by an inability to identify critical genetic events across the spectrum of human tumors with adequate sensitivity and scalability. Here, we have adapted high-throughput genotyping to query 238 known oncogene mutations across 1,000 human tumor samples. This approach established robust mutation distributions spanning 17 cancer types. Of 17 oncogenes analyzed, we found 14 to be mutated at least once, and 298 (30%) samples carried at least one mutation. Moreover, we identified previously unrecognized oncogene mutations in several tumor types and observed an unexpectedly high number of co-occurring mutations. These results offer a new dimension in tumor genetics, where mutations involving multiple cancer genes may be interrogated simultaneously and in 'real time' to guide cancer classification and rational therapeutic intervention.