Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Hereditary cancer variants and homologous recombination deficiency in biliary tract cancer.

BACKGROUND & AIMS: The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs. METHODS: We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status. RESULTS: Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p <0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants without a second hit or variants of other homologous recombination-related genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes. CONCLUSIONS: In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs. IMPACT AND IMPLICATIONS: We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Association of an IGHV3-66 gene variant with Kawasaki disease.

In a meta-analysis of three GWAS for susceptibility to Kawasaki disease (KD) conducted in Japan, Korea, and Taiwan and follow-up studies with a total of 11,265 subjects (3428 cases and 7837 controls), a significantly associated SNV in the immunoglobulin heavy variable gene (IGHV) cluster in 14q33.32 was identified (rs4774175; OR = 1.20, P = 6.0 × 10-9). Investigation of nonsynonymous SNVs of the IGHV cluster in 9335 Japanese subjects identified the C allele of rs6423677, located in IGHV3-66, as the most significant reproducible association (OR = 1.25, P = 6.8 × 10-10 in 3603 cases and 5731 controls). We observed highly skewed allelic usage of IGHV3-66, wherein the rs6423677 A allele was nearly abolished in the transcripts in peripheral blood mononuclear cells of both KD patients and healthy adults. Association of the high-expression allele with KD strongly indicates some active roles of B-cells or endogenous immunoglobulins in the disease pathogenesis. Considering that significant association of SNVs in the IGHV region with disease susceptibility was previously known only for rheumatic heart disease (RHD), a complication of acute rheumatic fever (ARF), these observations suggest that common B-cell related mechanisms may mediate the symptomology of KD and ARF as well as RHD.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Combined Genetic and Genealogic Studies Uncover a Large BAP1 Cancer Syndrome Kindred Tracing Back Nine Generations to a Common Ancestor from the 1700s.

We recently discovered an inherited cancer syndrome caused by BRCA1-Associated Protein 1 (BAP1) germline mutations, with high incidence of mesothelioma, uveal melanoma and other cancers and very high penetrance by age 55. To identify families with the BAP1 cancer syndrome, we screened patients with family histories of multiple mesotheliomas and melanomas and/or multiple cancers. We identified four families that shared an identical BAP1 mutation: they lived across the US and did not appear to be related. By combining family histories, molecular genetics, and genealogical approaches, we uncovered a BAP1 cancer syndrome kindred of ~80,000 descendants with a core of 106 individuals, whose members descend from a couple born in Germany in the early 1700s who immigrated to North America. Their descendants spread throughout the country with mutation carriers affected by multiple malignancies. Our data show that, once a proband is identified, extended analyses of these kindreds, using genomic and genealogical studies to identify the most recent common ancestor, allow investigators to uncover additional branches of the family that may carry BAP1 mutations. Using this knowledge, we have identified new branches of this family carrying BAP1 mutations. We have also implemented early-detection strategies that help identify cancers at early-stage, when they can be cured (melanomas) or are more susceptible to therapy (MM and other malignancies).