TNRC18 engages H3K9me3 to mediate silencing of endogenous retrotransposons.

Trimethylation of histone H3 lysine 9 (H3K9me3) is crucial for the regulation of gene repression and heterochromatin formation, cell-fate determination and organismal development1. H3K9me3 also provides an essential mechanism for silencing transposable elements1-4. However, previous studies have shown that canonical H3K9me3 readers (for example, HP1 (refs. 5-9) and MPP8 (refs. 10-12)) have limited roles in silencing endogenous retroviruses (ERVs), one of the main transposable element classes in the mammalian genome13. Here we report that trinucleotide-repeat-containing 18 (TNRC18), a poorly understood chromatin regulator, recognizes H3K9me3 to mediate the silencing of ERV class I (ERV1) elements such as LTR12 (ref. 14). Biochemical, biophysical and structural studies identified the carboxy-terminal bromo-adjacent homology (BAH) domain of TNRC18 (TNRC18(BAH)) as an H3K9me3-specific reader. Moreover, the amino-terminal segment of TNRC18 is a platform for the direct recruitment of co-repressors such as HDAC-Sin3-NCoR complexes, thus enforcing optimal repression of the H3K9me3-demarcated ERVs. Point mutagenesis that disrupts the TNRC18(BAH)-mediated H3K9me3 engagement caused neonatal death in mice and, in multiple mammalian cell models, led to derepressed expression of ERVs, which affected the landscape of cis-regulatory elements and, therefore, gene-expression programmes. Collectively, we describe a new H3K9me3-sensing and regulatory pathway that operates to epigenetically silence evolutionarily young ERVs and exert substantial effects on host genome integrity, transcriptomic regulation, immunity and development.

Genomic landscape of non-small cell lung cancer in smokers and never-smokers.

We report the results of whole-genome and transcriptome sequencing of tumor and adjacent normal tissue samples from 17 patients with non-small cell lung carcinoma (NSCLC). We identified 3,726 point mutations and more than 90 indels in the coding sequence, with an average mutation frequency more than 10-fold higher in smokers than in never-smokers. Novel alterations in genes involved in chromatin modification and DNA repair pathways were identified, along with DACH1, CFTR, RELN, ABCB5, and HGF. Deep digital sequencing revealed diverse clonality patterns in both never-smokers and smokers. All validated EFGR and KRAS mutations were present in the founder clones, suggesting possible roles in cancer initiation. Analysis revealed 14 fusions, including ROS1 and ALK, as well as novel metabolic enzymes. Cell-cycle and JAK-STAT pathways are significantly altered in lung cancer, along with perturbations in 54 genes that are potentially targetable with currently available drugs.

Mapping and characterization of structural variation in 17,795 human genomes.

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Association of structural variation with cardiometabolic traits in Finns.

The contribution of genome structural variation (SV) to quantitative traits associated with cardiometabolic diseases remains largely unknown. Here, we present the results of a study examining genetic association between SVs and cardiometabolic traits in the Finnish population. We used sensitive methods to identify and genotype 129,166 high-confidence SVs from deep whole-genome sequencing (WGS) data of 4,848 individuals. We tested the 64,572 common and low-frequency SVs for association with 116 quantitative traits and tested candidate associations using exome sequencing and array genotype data from an additional 15,205 individuals. We discovered 31 genome-wide significant associations at 15 loci, including 2 loci at which SVs have strong phenotypic effects: (1) a deletion of the ALB promoter that is greatly enriched in the Finnish population and causes decreased serum albumin level in carriers (p = 1.47 × 10-54) and is also associated with increased levels of total cholesterol (p = 1.22 × 10-28) and 14 additional cholesterol-related traits, and (2) a multi-allelic copy number variant (CNV) at PDPR that is strongly associated with pyruvate (p = 4.81 × 10-21) and alanine (p = 6.14 × 10-12) levels and resides within a structurally complex genomic region that has accumulated many rearrangements over evolutionary time. We also confirmed six previously reported associations, including five led by stronger signals in single nucleotide variants (SNVs) and one linking recurrent HP gene deletion and cholesterol levels (p = 6.24 × 10-10), which was also found to be strongly associated with increased glycoprotein level (p = 3.53 × 10-35). Our study confirms that integrating SVs in trait-mapping studies will expand our knowledge of genetic factors underlying disease risk.

Correction: Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases.

[This corrects the article DOI: 10.1371/journal.pmed.1002174.].

Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases.

BACKGROUND: Metastasis is the main cause of cancer patient deaths and remains a poorly characterized process. It is still unclear when in tumor progression the ability to metastasize arises and whether this ability is inherent to the primary tumor or is acquired well after primary tumor formation. Next-generation sequencing and analytical methods to define clonal heterogeneity provide a means for identifying genetic events and the temporal relationships between these events in the primary and metastatic tumors within an individual. METHODS AND FINDINGS: We performed DNA whole genome and mRNA sequencing on two primary tumors, each with either four or five distinct tissue site-specific metastases, from two individuals with triple-negative/basal-like breast cancers. As evidenced by their case histories, each patient had an aggressive disease course with abbreviated survival. In each patient, the overall gene expression signatures, DNA copy number patterns, and somatic mutation patterns were highly similar across each primary tumor and its associated metastases. Almost every mutation found in the primary was found in a metastasis (for the two patients, 52/54 and 75/75). Many of these mutations were found in every tumor (11/54 and 65/75, respectively). In addition, each metastasis had fewer metastatic-specific events and shared at least 50% of its somatic mutation repertoire with the primary tumor, and all samples from each patient grouped together by gene expression clustering analysis. TP53 was the only mutated gene in common between both patients and was present in every tumor in this study. Strikingly, each metastasis resulted from multiclonal seeding instead of from a single cell of origin, and few of the new mutations, present only in the metastases, were expressed in mRNAs. Because of the clinical differences between these two patients and the small sample size of our study, the generalizability of these findings will need to be further examined in larger cohorts of patients. CONCLUSIONS: Our findings suggest that multiclonal seeding may be common amongst basal-like breast cancers. In these two patients, mutations and DNA copy number changes in the primary tumors appear to have had a biologic impact on metastatic potential, whereas mutations arising in the metastases were much more likely to be passengers.

Rapid idiosyncratic mechanisms of clinical resistance to KRAS G12C inhibition.

BACKGROUNDThe KRAS proto-oncogene is among the most frequently mutated genes in cancer, yet for 40 years it remained an elusive therapeutic target. Recently, allosteric inhibitors that covalently bind to KRAS G12C mutations have been approved for use in lung adenocarcinomas. Although responses are observed, they are often short-lived, thus making in-depth characterization of the mechanisms of resistance of paramount importance.METHODSHere, we present a rapid-autopsy case of a patient who had a KRASG12C-mutant lung adenocarcinoma who initially responded to a KRAS G12C inhibitor but then rapidly developed resistance. Using deep-RNA and whole-exome sequencing comparing pretreatment, posttreatment, and matched normal tissues, we uncover numerous mechanisms of resistance to direct KRAS inhibition.RESULTSIn addition to decreased KRAS G12C-mutant allele frequency in refractory tumors, we also found reactivation of the MAPK pathway despite no new mutations in KRAS or its downstream mediators. Tumor cell-intrinsic and non-cell autonomous mechanisms included increased complement activation, coagulation, and tumor angiogenesis, and several lines of evidence of immunologic evasion.CONCLUSIONTogether, our findings reveal numerous mechanisms of resistance to current KRAS G12C inhibitors through enrichment of clonal populations, KRAS-independent downstream signaling, and diverse remodeling of the tumor microenvironment.FUNDINGRichard and Fran Duley, Jimmy and Kay Mann, the NIH, and the North Carolina Biotechnology Center.

Genomic profiling of murine mammary tumors identifies potential personalized drug targets for p53-deficient mammary cancers.

Targeted therapies against basal-like breast tumors, which are typically 'triple-negative breast cancers (TNBCs)', remain an important unmet clinical need. Somatic TP53 mutations are the most common genetic event in basal-like breast tumors and TNBC. To identify additional drivers and possible drug targets of this subtype, a comparative study between human and murine tumors was performed by utilizing a murine Trp53-null mammary transplant tumor model. We show that two subsets of murine Trp53-null mammary transplant tumors resemble aspects of the human basal-like subtype. DNA-microarray, whole-genome and exome-based sequencing approaches were used to interrogate the secondary genetic aberrations of these tumors, which were then compared to human basal-like tumors to identify conserved somatic genetic features. DNA copy-number variation produced the largest number of conserved candidate personalized drug targets. These candidates were filtered using a DNA-RNA Pearson correlation cut-off and a requirement that the gene was deemed essential in at least 5% of human breast cancer cell lines from an RNA-mediated interference screen database. Five potential personalized drug target genes, which were spontaneously amplified loci in both murine and human basal-like tumors, were identified: Cul4a, Lamp1, Met, Pnpla6 and Tubgcp3 As a proof of concept, inhibition of Met using crizotinib caused Met-amplified murine tumors to initially undergo complete regression. This study identifies Met as a promising drug target in a subset of murine Trp53-null tumors, thus identifying a potential shared driver with a subset of human basal-like breast cancers. Our results also highlight the importance of comparative genomic studies for discovering personalized drug targets and for providing a preclinical model for further investigations of key tumor signaling pathways.

Rare Variation in TET2 Is Associated with Clinically Relevant Prostate Carcinoma in African Americans.

BACKGROUND: Common variants have been associated with prostate cancer risk. Unfortunately, few are reproducibly linked to aggressive disease, the phenotype of greatest clinical relevance. One possible explanation is that rare genetic variants underlie a significant proportion of the risk for aggressive disease. METHOD: To identify such variants, we performed a two-stage approach using whole-exome sequencing followed by targeted sequencing of 800 genes in 652 aggressive prostate cancer patients and 752 disease-free controls in both African and European Americans. In each population, we tested rare variants for association using two gene-based aggregation tests. We established a study-wide significance threshold of 3.125 × 10-5 to correct for multiple testing. RESULTS: TET2 in African Americans was associated with aggressive disease, with 24.4% of cases harboring a rare deleterious variant compared with 9.6% of controls (FET P = 1.84 × 10-5, OR = 3.0; SKAT-O P = 2.74 × 10-5). We report 8 additional genes with suggestive evidence of association, including the DNA repair genes PARP2 and MSH6 Finally, we observed an excess of rare truncation variants in 5 genes, including the DNA repair genes MSH6, BRCA1, and BRCA2 This adds to the growing body of evidence that DNA repair pathway defects may influence susceptibility to aggressive prostate cancer. CONCLUSIONS: Our findings suggest that rare variants influence risk of clinically relevant prostate cancer and, if validated, could serve to identify men for screening, prophylaxis, and treatment. IMPACT: This study provides evidence that rare variants in TET2 may help identify African American men at increased risk for clinically relevant prostate cancer. Cancer Epidemiol Biomarkers Prev; 25(11); 1456-63. ©2016 AACR.

Inactivation of RASA1 promotes melanoma tumorigenesis via R-Ras activation.

Inactivation of Ras GTPase activating proteins (RasGAPs) can activate Ras, increasing the risk for tumor development. Utilizing a melanoma whole genome sequencing (WGS) data from 13 patients, we identified two novel, clustered somatic missense mutations (Y472H and L481F) in RASA1 (RAS p21 protein activator 1, also called p120RasGAP). We have shown that wild type RASA1, but not identified mutants, suppresses soft agar colony formation and tumor growth of BRAF mutated melanoma cell lines via its RasGAP activity toward R-Ras (related RAS viral (r-ras) oncogene homolog) isoform. Moreover, R-Ras increased and RASA1 suppressed Ral-A activation among Ras downstream effectors. In addition to mutations, loss of RASA1 expression was frequently observed in metastatic melanoma samples on melanoma tissue microarray (TMA) and a low level of RASA1 mRNA expression was associated with decreased overall survival in melanoma patients with BRAF mutations. Thus, these data support that RASA1 is inactivated by mutation or by suppressed expression in melanoma and that RASA1 plays a tumor suppressive role by inhibiting R-Ras, a previously less appreciated member of the Ras small GTPases.

Patterns and functional implications of rare germline variants across 12 cancer types.

Large-scale cancer sequencing data enable discovery of rare germline cancer susceptibility variants. Here we systematically analyse 4,034 cases from The Cancer Genome Atlas cancer cases representing 12 cancer types. We find that the frequency of rare germline truncations in 114 cancer-susceptibility-associated genes varies widely, from 4% (acute myeloid leukaemia (AML)) to 19% (ovarian cancer), with a notably high frequency of 11% in stomach cancer. Burden testing identifies 13 cancer genes with significant enrichment of rare truncations, some associated with specific cancers (for example, RAD51C, PALB2 and MSH6 in AML, stomach and endometrial cancers, respectively). Significant, tumour-specific loss of heterozygosity occurs in nine genes (ATM, BAP1, BRCA1/2, BRIP1, FANCM, PALB2 and RAD51C/D). Moreover, our homology-directed repair assay of 68 BRCA1 rare missense variants supports the utility of allelic enrichment analysis for characterizing variants of unknown significance. The scale of this analysis and the somatic-germline integration enable the detection of rare variants that may affect individual susceptibility to tumour development, a critical step toward precision medicine.

Clonal architectures and driver mutations in metastatic melanomas.

To reveal the clonal architecture of melanoma and associated driver mutations, whole genome sequencing (WGS) and targeted extension sequencing were used to characterize 124 melanoma cases. Significantly mutated gene analysis using 13 WGS cases and 15 additional paired extension cases identified known melanoma genes such as BRAF, NRAS, and CDKN2A, as well as a novel gene EPHA3, previously implicated in other cancer types. Extension studies using tumors from another 96 patients discovered a large number of truncation mutations in tumor suppressors (TP53 and RB1), protein phosphatases (e.g., PTEN, PTPRB, PTPRD, and PTPRT), as well as chromatin remodeling genes (e.g., ASXL3, MLL2, and ARID2). Deep sequencing of mutations revealed subclones in the majority of metastatic tumors from 13 WGS cases. Validated mutations from 12 out of 13 WGS patients exhibited a predominant UV signature characterized by a high frequency of C->T transitions occurring at the 3' base of dipyrimidine sequences while one patient (MEL9) with a hypermutator phenotype lacked this signature. Strikingly, a subclonal mutation signature analysis revealed that the founding clone in MEL9 exhibited UV signature but the secondary clone did not, suggesting different mutational mechanisms for two clonal populations from the same tumor. Further analysis of four metastases from different geographic locations in 2 melanoma cases revealed phylogenetic relationships and highlighted the genetic alterations responsible for differential drug resistance among metastatic tumors. Our study suggests that clonal evaluation is crucial for understanding tumor etiology and drug resistance in melanoma.

Integrated analysis of germline and somatic variants in ovarian cancer.

We report the first large-scale exome-wide analysis of the combined germline-somatic landscape in ovarian cancer. Here we analyse germline and somatic alterations in 429 ovarian carcinoma cases and 557 controls. We identify 3,635 high confidence, rare truncation and 22,953 missense variants with predicted functional impact. We find germline truncation variants and large deletions across Fanconi pathway genes in 20% of cases. Enrichment of rare truncations is shown in BRCA1, BRCA2 and PALB2. In addition, we observe germline truncation variants in genes not previously associated with ovarian cancer susceptibility (NF1, MAP3K4, CDKN2B and MLL3). Evidence for loss of heterozygosity was found in 100 and 76% of cases with germline BRCA1 and BRCA2 truncations, respectively. Germline-somatic interaction analysis combined with extensive bioinformatics annotation identifies 222 candidate functional germline truncation and missense variants, including two pathogenic BRCA1 and 1 TP53 deleterious variants. Finally, integrated analyses of germline and somatic variants identify significantly altered pathways, including the Fanconi, MAPK and MLL pathways.