Pan-cancer analysis of non-coding recurrent mutations and their possible involvement in cancer pathogenesis.

Cancer-related mutations have been mainly identified in protein-coding regions. Recent studies have demonstrated that mutations in non-coding regions of the genome could also be a risk factor for cancer. However, the non-coding regions comprise 98% of the total length of the human genome and contain a huge number of mutations, making it difficult to interpret their impacts on pathogenesis of cancer. To comprehensively identify cancer-related non-coding mutations, we focused on recurrent mutations in non-coding regions using somatic mutation data from COSMIC and whole-genome sequencing data from The Cancer Genome Atlas (TCGA). We identified 21 574 recurrent mutations in non-coding regions that were shared by at least two different samples from both COSMIC and TCGA databases. Among them, 580 candidate cancer-related non-coding recurrent mutations were identified based on epigenomic and chromatin structure datasets. One of such mutation was located in RREB1 binding site that is thought to interact with TEAD1 promoter. Our results suggest that mutations may disrupt the binding of RREB1 to the candidate enhancer region and increase TEAD1 expression levels. Our findings demonstrate that non-coding recurrent mutations and coding mutations may contribute to the pathogenesis of cancer.

Augmented oxidative stress increases 8-oxoguanine preferentially in the transcriptionally active genomic regions.

8-Oxoguanine (8-oxoG) is the most common DNA base modification in the mammalian genome, associated with oxidative stress. Here we analysed the alterations in the distribution of 8-oxoG across the entire murine genome, before and after an elevation of oxidative stress by the use of ferric nitrilotriacetate (Fe-NTA) as an oxidative stress inducer in the renal proximal tubules. We isolated DNA fragments containing 8-oxoGs with immunoprecipitation from the murine genome, and amplified them by PCR for a distribution analysis with microarray-based comparative genomic hybridisation. The distribution profiles revealed that frequencies of 8-oxoG fluctuated with a cycle of 1-10 Mb along the chromosomes and the amplitude of the fluctuation was reduced after Fe-NTA administration. The distributions of 8-oxoG along the entire genome in the control and oxidatively stressed conditions were negatively correlated with that of gene density but positively correlated with that of Lamin B1 interaction, which corresponds to lamina-associated domains. These results on the murine genome were consistent with those on the rat genome we previously reported. We further discovered a negative correlation between the distributions of 8-oxoG and transcriptional activity along the genome. Finally, a comparison of the distributions before and after Fe-NTA administration suggested that 8-oxoGs are generated in response to the augmented oxidative stress preferentially in the transcriptionally active genomic regions, where 8-oxoGs have been less accumulated in the control condition.

Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms.

In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often observed. Although these discordant mutation sites are thought to have emerged by ancestral polymorphism or gene flow, their frequency and distribution in the genome have not yet been analyzed in detail. Using the genome sequences of all protein coding genes of 25 inbred rat strains, we analyzed the frequency and genome-wide distribution of the discordant mutation sites. From the comparison of different substrains, it was found that these loci are not substrain specific, but are common among different groups of substrains, suggesting that the discordant sites might have mainly emerged through ancestral polymorphism. It was also revealed that the discordant sites are not uniformly distributed along chromosomes, but are concentrated at certain genomic loci, such as RT1, major histocompatibility complex of rats, and olfactory receptors, indicating that genes known to be highly polymorphic tend to have more discordant sites. Our results also showed that loci with a high density of discordant sites are also rich in heterozygous variants, even though these are inbred strains.

Intratumor heterogeneity of HMCN1 mutant alleles associated with poor prognosis in patients with breast cancer.

Human breast cancers comprise a complex and highly heterogeneous population of tumor cells. Intratumor heterogeneity is an underlying cause of resistance to effective therapies and disease recurrence. To explore prognostic factors based on intratumor heterogeneity, we analyzed genomic mutations in breast cancer patients registered in The Cancer Genome Atlas. We calculated the variant allele frequency (VAF) at each mutation site and evaluated the associations of VAFs with the prognosis of breast cancer. VAFs of HMCN1 correlated with the prognosis and lymph node status. Although the detailed function of HMCN1 remains unknown, it is located in extracellular matrix and the mutation in the gene might be associated with cancer cell invasion and metastasis. This finding suggests that HMCN1 is a potential metastatic factor and can be a candidate gene for targeted breast cancer therapy.

Pan-cancer analysis of intratumor heterogeneity associated with patient prognosis using multidimensional measures.

Human cancers accumulate various mutations during development and consist of highly heterogeneous cell populations. This phenomenon is called intratumor heterogeneity (ITH). ITH is known to be involved in tumor growth, progression, invasion, and metastasis, presenting obstacles to accurate diagnoses and effective treatments. Numerous studies have explored the dynamics of ITH, including constructions of phylogenetic trees in cancer samples using multiregional ultradeep sequencing and simulations of evolution using statistical models. Although ITH is associated with prognosis, it is still challenging to use the characteristics of ITH as prognostic factors because of difficulties in quantifying ITH precisely. In this study, we analyzed the relationship between patient prognosis and the distribution of variant allele frequencies (VAFs) in cancer samples (n = 6,064) across 16 cancer types registered in The Cancer Genome Atlas. To measure VAF distributions multidimensionally, we adopted parameters that define the shape of VAF distributions and evaluated the relationships between these parameters and prognosis. In seven cancer types, we found significant relationships between prognosis and VAF distributions. Moreover, we observed that samples with a larger amount of mutations were not necessarily linked to worse prognosis. By evaluating the ITH from multidimensional viewpoints, it will be possible to provide a more accurate prediction of cancer prognosis.

Identification of Candidate Genes for Generalized Tonic-Clonic Seizures in Noda Epileptic Rat.

The Noda epileptic rat (NER) exhibits generalized tonic-clonic seizures (GTCS). A genetic linkage analysis identified two GTCS-associated loci, Ner1 on Chr 1 and Ner3 on Chr 5. The wild-type Ner1 and Ner3 alleles suppressed GTCS when combined in double-locus congenic lines, but not when present in single-locus congenic lines. Global expression analysis revealed that cholecystokinin B receptor (Cckbr) and suppressor of tumorigenicity 5 (St5), which map within Ner1, and PHD finger protein 24 (Phf24), which maps within Ner3, were significantly downregulated in NER. De novo BAC sequencing detected an insertion of an endogenous retrovirus sequence in intron 2 of the Phf24 gene in the NER genome, and PHF24 protein was almost absent in the NER brain. Phf24 encodes a Gαi-interacting protein involved in GABAB receptor signaling pathway. Based on these findings, we conclude that Cckbr, St5, and Phf24 are strong candidate genes for GTCS in NER.

A deletion in the intergenic region upstream of Ednrb causes head spot in the rat strain KFRS4/Kyo.

BACKGROUND: Head spot is one of the phenotypes identified in the KFRS4/Kyo rat strain. Although previous linkage analysis suggested that Ednrb, which is frequently involved in coat color variations in various animals, could be the gene responsible for this phenotype, no mutations have been identified in its coding region. RESULTS: To identify mutations causative of this phenotype in KFRS4/Kyo, we analyzed target capture sequencing data that we recently generated. Our target capture method has a unique feature, i.e., it covers not only exonic regions but also conserved non-coding sequences (CNSs) among vertebrates; therefore, it has the potential to detect regulatory mutations. We identified a deletion of approximately 50 kb in length approximately 50 kb upstream of Ednrb. A comparative analysis with the epigenomic data in the corresponding region in humans and mice showed that one of the CNSs might be an enhancer. Further comparison with Hi-C data, which provide information about chromosome conformation, indicated that the putative enhancer is spatially close to the promoter of Ednrb, suggesting that it acts as an enhancer of Ednrb. CONCLUSIONS: These in silico data analyses strongly suggest that the identified deletion in the intergenic region upstream of Ednrb, which might contain a melanocyte-specific enhancer, is the mutation causative of the head spot phenotype in the KFRS4/Kyo rat strain.

Application of target capture sequencing of exons and conserved non-coding sequences to 20 inbred rat strains.

We report sequence data obtained by our recently devised target capture method TargetEC applied to 20 inbred rat strains. This method encompasses not only all annotated exons but also highly conserved non-coding sequences shared among vertebrates. The total length of the target regions covers 146.8 Mb. On an average, we obtained 31.7 × depth of target coverage and identified 154,330 SNVs and 24,368 INDELs for each strain. This corresponds to 470,037 unique SNVs and 68,652 unique INDELs among the 20 strains. The sequence data can be accessed at DDBJ/EMBL/GenBank under accession number PRJDB4648, and the identified variants have been deposited at http://bioinfo.sls.kyushu-u.ac.jp/rat_target_capture/20_strains.vcf.gz.

Design and application of a target capture sequencing of exons and conserved non-coding sequences for the rat.

BACKGROUND: Target capture sequencing is an efficient approach to directly identify the causative mutations of genetic disorders. To apply this strategy to laboratory rats exhibiting various phenotypes, we developed a novel target capture probe set, TargetEC (target capture for exons and conserved non-coding sequences), which can identify mutations not only in exonic regions but also in conserved non-coding sequences and thus can detect regulatory mutations. RESULTS: TargetEC covers 1,078,129 regions spanning 146.8 Mb of the genome. We applied TargetEC to four inbred rat strains (WTC/Kyo, WTC-swh/Kyo, PVG/Seac, and KFRS4/Kyo) maintained by the National BioResource Project for the Rat in Japan, and successfully identified mutations associated with these phenotypes, including one mutation detected in a conserved non-coding sequence. CONCLUSIONS: The method developed in this study can be used to efficiently identify regulatory mutations, which cannot be detected using conventional exome sequencing, and will help to deepen our understanding of the relationships between regulatory mutations and associated phenotypes.

Mutations in CDCA7 and HELLS cause immunodeficiency-centromeric instability-facial anomalies syndrome.

The life-threatening Immunodeficiency, Centromeric Instability and Facial Anomalies (ICF) syndrome is a genetically heterogeneous autosomal recessive disorder. Twenty percent of patients cannot be explained by mutations in the known ICF genes DNA methyltransferase 3B or zinc-finger and BTB domain containing 24. Here we report mutations in the cell division cycle associated 7 and the helicase, lymphoid-specific genes in 10 unexplained ICF cases. Our data highlight the genetic heterogeneity of ICF syndrome; however, they provide evidence that all genes act in common or converging pathways leading to the ICF phenotype.

Genome-wide profiling of 8-oxoguanine reveals its association with spatial positioning in nucleus.

8-Oxoguanine (8-oxoG) is one of the most common DNA lesions generated by reactive oxygen species. In this study, we analysed the genome-wide distribution profile of 8-oxoG by combining immunoprecipitation by antibodies specific for the DNA fragments containing 8-oxoG with a microarray that covers rat genome. Genome-wide mapping of 8-oxoG in normal rat kidney revealed that 8-oxoG is preferentially located at gene deserts. We did not observe differences in 8-oxoG levels between groups of genes with high and low expression, possibly because of the generally low 8-oxoG levels in genic regions compared with gene deserts. The distribution of 8-oxoG and lamina-associated domains (LADs) were strongly correlated, suggesting that the spatial location of genomic DNA in the nucleus determines the susceptibility to oxidative modifications. One possible explanation for high 8-oxoG levels in LADs is that the nuclear periphery is more susceptible to the oxidative damage caused by the extra-nuclear factors. Moreover, LADs have a rather compact conformation, which may limit the recruitment of repair components to the modified bases.

Association of microRNA-34a overexpression with proliferation is cell type-dependent.

Recently Welch et al. reported that microRNA (miRNA)-34a functions as a potential tumor suppressor in neuroblastoma cells (Oncogene 26: 5017-22, 2007). Here, we conversely show that miRNA-34a supports cell proliferation in rat oxidative stress-induced renal carcinogenesis and is overexpressed in various types of human cancers. While searching for genetically unstable chromosomal areas in rat renal carcinogenesis, we found the miRNA-34 family reciprocally overexpressed in chromosomal areas with frequent allelic loss. By in situ hybridization and reverse transcription-polymerase chain reaction, cerebral neurons and Purkinje cells showed the highest expression of a major type, miRNA-34a, followed by a variety of endocrine cells and proliferating cells including germinal center lymphocytes and mouse embryonic fibroblasts and stem cells. In contrast, normal renal tubules, hepatocytes and myocardial cells showed faint expression. After 3 weeks of ferric nitrilotriacetate (Fe-NTA)-induced oxidative stress, regenerating renal proximal tubular cells showed high miRNA-34a expression. All of the Fe-NTA-induced rat renal carcinomas and an array of human cancers (151 positive cases of 177) showed high expression of miRNA-34a. Furthermore, knockdown of miRNA-34a with small interfering RNA significantly suppressed proliferation not only of renal carcinoma cells but also of HeLa and MCF7 cells. These results indicate that miRNA-34a overexpression, an acquired trait during carcinogenesis, supports cell proliferation in the majority of cancers suggesting an unexpected link in the cellular metabolism between cancer and neuronal and/or endocrine cells, which warrants further investigation.