Whole genome DNA sequencing provides an atlas of somatic mutagenesis in healthy human cells and identifies a tumor-prone cell type.

BACKGROUND: The lifelong accumulation of somatic mutations underlies age-related phenotypes and cancer. Mutagenic forces are thought to shape the genome of aging cells in a tissue-specific way. Whole genome analyses of somatic mutation patterns, based on both types and genomic distribution of variants, can shed light on specific processes active in different human tissues and their effect on the transition to cancer. RESULTS: To analyze somatic mutation patterns, we compile a comprehensive genetic atlas of somatic mutations in healthy human cells. High-confidence variants are obtained from newly generated and publicly available whole genome DNA sequencing data from single non-cancer cells, clonally expanded in vitro. To enable a well-controlled comparison of different cell types, we obtain single genome data (92% mean coverage) from multi-organ biopsies from the same donors. These data show multiple cell types that are protected from mutagens and display a stereotyped mutation profile, despite their origin from different tissues. Conversely, the same tissue harbors cells with distinct mutation profiles associated to different differentiation states. Analyses of mutation rate in the coding and non-coding portions of the genome identify a cell type bearing a unique mutation pattern characterized by mutation enrichment in active chromatin, regulatory, and transcribed regions. CONCLUSIONS: Our analysis of normal cells from healthy donors identifies a somatic mutation landscape that enhances the risk of tumor transformation in a specific cell population from the kidney proximal tubule. This unique pattern is characterized by high rate of mutation accumulation during adult life and specific targeting of expressed genes and regulatory regions.

Somatic mutation that affects transcription factor binding upstream of CD55 in the temporal cortex of a late-onset Alzheimer disease patient.

Alzheimer's disease (AD) is the most common neurodegenerative disease worldwide. Familial cases suggest genetic components; however, monogenetic causes are few, and the vast majority of incidences have unknown cause. Sequencing efforts have focused on germline mutations, but improved technology has opened up for studies on somatic mutations in affected brain tissue samples. Here we use ultra-deep sequencing on brain and blood from early-onset AD (EOAD) and late-onset AD (LOAD) patients and non-AD individuals (n = 16). In total, 2.86 Mb of genomic regions, previously associated with AD, were targeted included 28 genes and upstream and downstream regulatory regions. Tailored downstream bioinformatics filtering identified 11 somatic single nucleotide variants in the temporal cortex in AD patients and none in the controls. One variant was validated to be present at 0.4% allele frequency in temporal cortex of a LOAD patient. This variant was predicted to affect transcription factor binding sites upstream of the CD55 gene, contributing to AD pathogenesis by affecting the complement system. Our results suggest that future studies targeting larger portions of the genome for somatic mutation analysis are important to obtain an increased understanding for the molecular basis of both EOAD and LOAD.

Somatic mutagenesis in satellite cells associates with human skeletal muscle aging.

Human aging is associated with a decline in skeletal muscle (SkM) function and a reduction in the number and activity of satellite cells (SCs), the resident stem cells. To study the connection between SC aging and muscle impairment, we analyze the whole genome of single SC clones of the leg muscle vastus lateralis from healthy individuals of different ages (21-78 years). We find an accumulation rate of 13 somatic mutations per genome per year, consistent with proliferation of SCs in the healthy adult muscle. SkM-expressed genes are protected from mutations, but aging results in an increase in mutations in exons and promoters, targeting genes involved in SC activity and muscle function. In agreement with SC mutations affecting the whole tissue, we detect a missense mutation in a SC propagating to the muscle. Our results suggest somatic mutagenesis in SCs as a driving force in the age-related decline of SkM function.

SweGen: a whole-genome data resource of genetic variability in a cross-section of the Swedish population.

Here we describe the SweGen data set, a comprehensive map of genetic variation in the Swedish population. These data represent a basic resource for clinical genetics laboratories as well as for sequencing-based association studies by providing information on genetic variant frequencies in a cohort that is well matched to national patient cohorts. To select samples for this study, we first examined the genetic structure of the Swedish population using high-density SNP-array data from a nation-wide cohort of over 10 000 Swedish-born individuals included in the Swedish Twin Registry. A total of 1000 individuals, reflecting a cross-section of the population and capturing the main genetic structure, were selected for whole-genome sequencing. Analysis pipelines were developed for automated alignment, variant calling and quality control of the sequencing data. This resulted in a genome-wide collection of aggregated variant frequencies in the Swedish population that we have made available to the scientific community through the website https://swefreq.nbis.se. A total of 29.2 million single-nucleotide variants and 3.8 million indels were detected in the 1000 samples, with 9.9 million of these variants not present in current databases. Each sample contributed with an average of 7199 individual-specific variants. In addition, an average of 8645 larger structural variants (SVs) were detected per individual, and we demonstrate that the population frequencies of these SVs can be used for efficient filtering analyses. Finally, our results show that the genetic diversity within Sweden is substantial compared with the diversity among continental European populations, underscoring the relevance of establishing a local reference data set.

A novel RNA sequencing data analysis method for cell line authentication.

We have developed a novel analysis method that can interrogate the authenticity of biological samples used for generation of transcriptome profiles in public data repositories. The method uses RNA sequencing information to reveal mutations in expressed transcripts and subsequently confirms the identity of analysed cells by comparison with publicly available cell-specific mutational profiles. Cell lines constitute key model systems widely used within cancer research, but their identity needs to be confirmed in order to minimise the influence of cell contaminations and genetic drift on the analysis. Using both public and novel data, we demonstrate the use of RNA-sequencing data analysis for cell line authentication by examining the validity of COLO205, DLD1, HCT15, HCT116, HKE3, HT29 and RKO colorectal cancer cell lines. We successfully authenticate the studied cell lines and validate previous reports indicating that DLD1 and HCT15 are synonymous. We also show that the analysed HKE3 cells harbour an unexpected KRAS-G13D mutation and confirm that this cell line is a genuine KRAS dosage mutant, rather than a true isogenic derivative of HCT116 expressing only the wild type KRAS. This authentication method could be used to revisit the numerous cell line based RNA sequencing experiments available in public data repositories, analyse new experiments where whole genome sequencing is not available, as well as facilitate comparisons of data from different experiments, platforms and laboratories.

Linkage and whole genome sequencing identify a locus on 6q25-26 for formal thought disorder and implicate MEF2A regulation.

Formal thought disorder is a major feature of schizophrenia and other psychotic disorders. It is heritable, found in healthy relatives of patients with schizophrenia and other mental disorders but knowledge of specific genetic factors is lacking. The aim of this study was to search for biologically relevant high-risk variants. Formal thought disorder was assessed in participants in the Copenhagen Schizophrenia Linkage Study (N=236), a unique high-risk family study comprised of six large pedigrees. Microsatellite linkage analysis of formal thought disorder was performed and subsequent haplotype analysis of the implicated region using phased microsatellite and SNP genotypes. Whole genome sequencing (N=3) was used in the attempt to identify causative variants in the linkage region. Linkage analysis of formal thought disorder resulted in a single peak at chromosome 6(q26-q27) centred on marker D6S1277, with a maximum LOD score of 4.0. Phasing and fine mapping of the linkage peak identified a 5.5Mb haplotype (chr6:162242322-167753547, hg18) in 31 individuals, all belonging to the same pedigree sharing the haplotype from a common ancestor. The haplotype segregated with increased total thought disorder index score (P=4.9 × 10(-5)) and qualitatively severe forms of thought disturbances. Whole genome sequencing identified a novel nucleotide deletion (chr6:164377205 AG>A, hg18) predicted to disrupt the potential binding of the transcription factor MEF2A. The MEF2A binding site is located between two genes previously reported to associate with schizophrenia, QKI (HGNC:21100) and PDE10A (HGNC:8772). The findings are consistent with MEF2A deregulation conferring risk of formal thought disorder.

Quantitative multigene FISH on breast carcinomas identifies der(1;16)(q10;p10) as an early event in luminal A tumors.

In situ detection of genomic alterations in cancer provides information at the single cell level, making it possible to investigate genomic changes in cells in a tissue context. Such topological information is important when studying intratumor heterogeneity as well as alterations related to different steps in tumor progression. We developed a quantitative multigene fluorescence in situ hybridization (QM FISH) method to detect multiple genomic regions in single cells in complex tissues. As a "proof of principle" we applied the method to breast cancer samples to identify partners in whole arm (WA) translocations. WA gain of chromosome arm 1q and loss of chromosome arm 16q are among the most frequent genomic events in breast cancer. By designing five specific FISH probes based on breakpoint information from comparative genomic hybridization array (aCGH) profiles, we visualized chromosomal translocations in clinical samples at the single cell level. By analyzing aCGH data from 295 patients with breast carcinoma with known molecular subtype, we found concurrent WA gain of 1q and loss of 16q to be more frequent in luminal A tumors compared to other molecular subtypes. QM FISH applied to a subset of samples (n = 26) identified a derivative chromosome der(1;16)(q10;p10), a result of a centromere-close translocation between chromosome arms 1q and 16p. In addition, we observed that the distribution of cells with the translocation varied from sample to sample, some had a homogenous cell population while others displayed intratumor heterogeneity with cell-to-cell variation. Finally, for one tumor with both preinvasive and invasive components, the fraction of cells with translocation was lower and more heterogeneous in the preinvasive tumor cells compared to the cells in the invasive component.

Mutations in BRIP1 confer high risk of ovarian cancer.

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10(-14)). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.

Genomic architecture characterizes tumor progression paths and fate in breast cancer patients.

Distinct molecular subtypes of breast carcinomas have been identified, but translation into clinical use has been limited. We have developed two platform-independent algorithms to explore genomic architectural distortion using array comparative genomic hybridization data to measure (i) whole-arm gains and losses [whole-arm aberration index (WAAI)] and (ii) complex rearrangements [complex arm aberration index (CAAI)]. By applying CAAI and WAAI to data from 595 breast cancer patients, we were able to separate the cases into eight subgroups with different distributions of genomic distortion. Within each subgroup data from expression analyses, sequencing and ploidy indicated that progression occurs along separate paths into more complex genotypes. Histological grade had prognostic impact only in the luminal-related groups, whereas the complexity identified by CAAI had an overall independent prognostic power. This study emphasizes the relation among structural genomic alterations, molecular subtype, and clinical behavior and shows that objective score of genomic complexity (CAAI) is an independent prognostic marker in breast cancer.

Inferring tumor progression from genomic heterogeneity.

Cancer progression in humans is difficult to infer because we do not routinely sample patients at multiple stages of their disease. However, heterogeneous breast tumors provide a unique opportunity to study human tumor progression because they still contain evidence of early and intermediate subpopulations in the form of the phylogenetic relationships. We have developed a method we call Sector-Ploidy-Profiling (SPP) to study the clonal composition of breast tumors. SPP involves macro-dissecting tumors, flow-sorting genomic subpopulations by DNA content, and profiling genomes using comparative genomic hybridization (CGH). Breast carcinomas display two classes of genomic structural variation: (1) monogenomic and (2) polygenomic. Monogenomic tumors appear to contain a single major clonal subpopulation with a highly stable chromosome structure. Polygenomic tumors contain multiple clonal tumor subpopulations, which may occupy the same sectors, or separate anatomic locations. In polygenomic tumors, we show that heterogeneity can be ascribed to a few clonal subpopulations, rather than a series of gradual intermediates. By comparing multiple subpopulations from different anatomic locations, we have inferred pathways of cancer progression and the organization of tumor growth.

Novel patterns of genome rearrangement and their association with survival in breast cancer.

Representational Oligonucleotide Microarray Analysis (ROMA) detects genomic amplifications and deletions with boundaries defined at a resolution of approximately 50 kb. We have used this technique to examine 243 breast tumors from two separate studies for which detailed clinical data were available. The very high resolution of this technology has enabled us to identify three characteristic patterns of genomic copy number variation in diploid tumors and to measure correlations with patient survival. One of these patterns is characterized by multiple closely spaced amplicons, or "firestorms," limited to single chromosome arms. These multiple amplifications are highly correlated with aggressive disease and poor survival even when the rest of the genome is relatively quiet. Analysis of a selected subset of clinical material suggests that a simple genomic calculation, based on the number and proximity of genomic alterations, correlates with life-table estimates of the probability of overall survival in patients with primary breast cancer. Based on this sample, we generate the working hypothesis that copy number profiling might provide information useful in making clinical decisions, especially regarding the use or not of systemic therapies (hormonal therapy, chemotherapy), in the management of operable primary breast cancer with ostensibly good prognosis, for example, small, node-negative, hormone-receptor-positive diploid cases.

PROBER: oligonucleotide FISH probe design software.

UNLABELLED: PROBER is an oligonucleotide primer design software application that designs multiple primer pairs for generating PCR probes useful for fluorescence in situ hybridization (FISH). PROBER generates Tiling Oligonucleotide Probes (TOPs) by masking repetitive genomic sequences and delineating essentially unique regions that can be amplified to yield small (100-2000 bp) DNA probes that in aggregate will generate a single, strong fluorescent signal for regions as small as a single gene. TOPs are an alternative to bacterial artificial chromosomes (BACs) that are commonly used for FISH but may be unstable, unavailable, chimeric, or non-specific to small (10-100 kb) genomic regions. PROBER can be applied to any genomic locus, with the limitation that the locus must contain at least 10 kb of essentially unique blocks. To test the software, we designed a number of probes for genomic amplifications and hemizygous deletions that were initially detected by Representational Oligonucleotide Microarray Analysis of breast cancer tumors. AVAILABILITY: http://prober.cshl.edu

Large-scale copy number polymorphism in the human genome.

The extent to which large duplications and deletions contribute to human genetic variation and diversity is unknown. Here, we show that large-scale copy number polymorphisms (CNPs) (about 100 kilobases and greater) contribute substantially to genomic variation between normal humans. Representational oligonucleotide microarray analysis of 20 individuals revealed a total of 221 copy number differences representing 76 unique CNPs. On average, individuals differed by 11 CNPs, and the average length of a CNP interval was 465 kilobases. We observed copy number variation of 70 different genes within CNP intervals, including genes involved in neurological function, regulation of cell growth, regulation of metabolism, and several genes known to be associated with disease.