Allelic Selection of Amplicons in Glioblastoma Revealed by Combining Somatic and Germline Analysis

Cancer is a disease driven by a combination of inherited risk alleles coupled with the acquisition of somatic mutations, including amplification and deletion of genomic DNA. Potential relationships between the inherited and somatic aspects of the disease have only rarely been examined on a genome-wide level. Applying a novel integrative analysis of SNP and copy number measurements, we queried the tumor and normal-tissue genomes of 178 glioblastoma patients from the Cancer Genome Atlas project for preferentially amplified alleles, under the hypothesis that oncogenic germline variants will be selectively amplified in the tumor environment. Selected alleles are revealed by allelic imbalance in amplification across samples. This general approach is based on genetic principles and provides a method for identifying important tumor-related alleles. We find that SNP alleles that are most significantly overrepresented in amplicons tend to occur in genes involved with regulation of kinase and transferase activity, and many of these genes are known contributors to gliomagenesis. The analysis also implicates variants in synapse genes. By incorporating gene expression data, we demonstrate synergy between preferential allelic amplification and expression in DOCK4 and EGFR. Our results support the notion that combining germline and tumor genetic data can identify regions relevant to cancer biology

Accurate estimation of homologue-specific DNA concentration-ratios in cancer samples allows long-range haplotyping

Interpretation of allelic copy measurements at polymorphic markers in cancer samples presents distinctive challenges and opportunities. Due to frequent gross chromosomal alterations occurring in cancer (aneuploidy), many genomic regions are present at homologous-allele imbalance. Within such regions, the unequal contribution of alleles at heterozygous markers allows for direct phasing of the haplotype derived from each individual parent. In addition, genome-wide estimates of homologue specific copy- ratios (HSCRs) are important for interpretation of the cancer genome in terms of fixed integral copy-numbers. We describe HAPSEG, a probabilistic method to interpret bi- allelic marker data in cancer samples. HAPSEG operates by partitioning the genome into segments of distinct copy number and modeling the four distinct genotypes in each segment. We describe general methods for fitting these models to data which are suit- able for both SNP microarrays and massively parallel sequencing data. In addition, we demonstrate a specially tailored error-model for interpretation of systematic variations arising in microarray platforms. The ability to directly determine haplotypes from cancer samples represents an opportunity to expand reference panels of phased chromosomes, which may have general interest in various population genetic applications. In addition, this property may be exploited to interrogate the relationship between germline risk and cancer phenotype with greater sensitivity than is possible using unphased genotype. Finally, we exploit the statistical dependency of phased genotypes to enable the fitting of more elaborate sample-level error-model parameters, allowing more accurate estimation of HSCRs in cancer samples

INVESTIGATING INVASION IN DUCTAL CARCINOMA IN SITU WITH TOPOGRAPHICAL SINGLE CELL GENOME SEQUENCING

Synchronous Ductal Carcinoma in situ (DCIS-IDC) is an early stage breast cancer invasion in which it is possible to delineate genomic evolution during invasion because of the presence of both in situ and invasive regions within the same sample. While laser capture microdissection studies of DCIS-IDC examined the relationship between the paired in situ (DCIS) and invasive (IDC) regions, these studies were either confounded by bulk tissue or limited to a small set of genes or markers. To overcome these challenges, we developed Topographic Single Cell Sequencing (TSCS), which combines laser-catapulting with single cell DNA sequencing to measure genomic copy number profiles from single tumor cells while preserving their spatial context. We applied TSCS to sequence 1,293 single cells from 10 synchronous DCIS patients. We also applied deep-exome sequencing to the in situ, invasive and normal tissues for the DCIS-IDC patients. Previous bulk tissue studies had produced several conflicting models of tumor evolution. Our data support a multiclonal invasion model, in which genome evolution occurs within the ducts and gives rise to multiple subclones that escape the ducts into the adjacent tissues to establish the invasive carcinomas. In summary, we have developed a novel method for single cell DNA sequencing, which preserves spatial context, and applied this method to understand clonal evolution during the transition between carcinoma in situ to invasive ductal carcinoma

Screening and association testing of common coding variation in steroid hormone receptor co-activator and co-repressor genes in relation to breast cancer risk: the Multiethnic Cohort

<p>Abstract</p> <p>Background</p> <p>Only a limited number of studies have performed comprehensive investigations of coding variation in relation to breast cancer risk. Given the established role of estrogens in breast cancer, we hypothesized that coding variation in steroid receptor coactivator and corepressor genes may alter inter-individual response to estrogen and serve as markers of breast cancer risk.</p> <p>Methods</p> <p>We sequenced the coding exons of 17 genes (<it>EP300, CCND1, NME1, NCOA1, NCOA2, NCOA3, SMARCA4, SMARCA2, CARM1, FOXA1, MPG, NCOR1, NCOR2, CALCOCO1, PRMT1, PPARBP </it>and <it>CREBBP</it>) suggested to influence transcriptional activation by steroid hormone receptors in a multiethnic panel of women with advanced breast cancer (n = 95): African Americans, Latinos, Japanese, Native Hawaiians and European Americans. Association testing of validated coding variants was conducted in a breast cancer case-control study (1,612 invasive cases and 1,961 controls) nested in the Multiethnic Cohort. We used logistic regression to estimate odds ratios for allelic effects in ethnic-pooled analyses as well as in subgroups defined by disease stage and steroid hormone receptor status. We also investigated effect modification by established breast cancer risk factors that are associated with steroid hormone exposure.</p> <p>Results</p> <p>We identified 45 coding variants with frequencies ≥ 1% in any one ethnic group (43 non-synonymous variants). We observed nominally significant positive associations with two coding variants in ethnic-pooled analyses (<it>NCOR2</it>: His52Arg, OR = 1.79; 95% CI, 1.05–3.05; <it>CALCOCO1</it>: Arg12His, OR = 2.29; 95% CI, 1.00–5.26). A small number of variants were associated with risk in disease subgroup analyses and we observed no strong evidence of effect modification by breast cancer risk factors. Based on the large number of statistical tests conducted in this study, the nominally significant associations that we observed may be due to chance, and will need to be confirmed in other studies.</p> <p>Conclusion</p> <p>Our findings suggest that common coding variation in these candidate genes do not make a substantial contribution to breast cancer risk in the general population. Cataloging and testing of coding variants in coactivator and corepressor genes should continue and may serve as a valuable resource for investigations of other hormone-related phenotypes, such as inter-individual response to hormonal therapies used for cancer treatment and prevention.</p

Analysis of alterations in the human cancer genome

Thesis (Ph. D.)--Harvard-MIT Division of Health Sciences and Technology, 2011.Cataloged from PDF version of thesis.Includes bibliographical references.Aneuploidy, an abnormal complement of chromosomes, is present in approximately 90% of human malignancies. Despite over 100 years of research, many questions remain regarding the contribution of aneuploidy to the cancer phenotype. In this thesis, we develop computational methods to infer the presence and specific patterns of aneuploidy across thousands of primary cancer tissue specimens. We then combine these inferences with clinical and genomic features of the cancer samples to refine our understanding of both the clinical implications of aneuploidy, and how it evolves in various human cancers. We identified a signature of chromosomal instability from specific genes whose expression was consistently correlated with aneuploidy in several cancer types, and which was predictive of poor clinical outcome multiple cancer types. Current genomic characterization techniques measure somatic alterations in a cancer sample in units of genomes (DNA mass). The meaning of such measurements is highly dependent on the tumors purity and its overall ploidy; they are hence complicated to interpret and compare across samples. Ideally, copy-number should be measured in copies-per-cancer-cell. Such measurements are straightforward to interpret and, for alterations that are fixed in the cancer cell population, are simple integer values. We develop two computational methods to infer tumor purity and malignant cell ploidy directly from allelic analysis of DNA. First we describe HAPSEG, a probabilistic method to interpret bi-allelic marker data in cancer samples in order to produce genome-wide estimates of homologue specific copy-ratios. Second, we describe ABSOLUTE, a method that infers purity, ploidy, and absolute copy-numbers from the estimates produced by HAPSEG. In addition, ABSOLUTE can analyze point mutations to detect subclonal heterogeneity and somatic homozygosity. We used ABSOLUTE to analyze ovarian cancer data and discovered that 54% of somatic point mutations were, in fact, subclonal. In contrast, mutations occurring in key tumor suppressor genes, TP53 and NF1 were predominantly clonal and homozygous. Analysis of absolute allelic copy-number profiles from 3,155 cancer specimens revealed that genome-doubling events are common in human cancer, and likely occur in already aneuploid cells in many cancer types. By correlating genome-doubling status with mutation data, we found that homozygous mutations in NF1 occurred predominantly in non-doubled samples. This finding suggests that genome doubling influences the pathways of tumor progression, with recessive inactivation being less common after genome doubling.by Scott L. Carter.Ph.D

The Interplay between Natural Selection and Susceptibility to Melanoma on Allele 374F of SLC45A2 Gene in a South European Population

We aimed to study the selective pressures interacting on SLC45A2 to investigate the interplay between selection and susceptibility to disease. Thus, we enrolled 500 volunteers from a geographically limited population (Basques from the North of Spain) and by resequencing the whole coding region and intron 5 of the 34 most and the 34 least pigmented individuals according to the reflectance distribution, we observed that the polymorphism Leu374Phe (L374F, rs16891982) was statistically associated with skin color variability within this sample. In particular, allele 374F was significantly more frequent among the individuals with lighter skin. Further genotyping an independent set of 558 individuals of a geographically wider population with known ancestry in the Spanish population also revealed that the frequency of L374F was significantly correlated with the incident UV radiation intensity. Selection tests suggest that allele 374F is being positively selected in South Europeans, thus indicating that depigmentation is an adaptive process. Interestingly, by genotyping 119 melanoma samples, we show that this variant is also associated with an increased susceptibility to melanoma in our populations. The ultimate driving force for this adaptation is unknown, but it is compatible with the vitamin D hypothesis. This shows that molecular evolution analysis can be used as a useful technology to predict phenotypic and biomedical consequences in humans

Trypanosomosis: potential driver of selection in African cattle

Trypanosomosis is a serious cause of reduction in productivity of cattle in tsetse-fly infested areas. Baoule and other local Taurine cattle breeds in Burkina Faso are trypanotolerant. Zebuine cattle, which are also kept there are susceptible to trypanosomosis but bigger in body size. Farmers have continuously been intercrossing Baoule and Zebu animals to increase production and disease tolerance. The aim of this study was to compare levels of zebuine and taurine admixture in genomic regions potentially involved in trypanotolerance with background admixture of composites to identify differences in allelic frequencies of tolerant and non tolerant animals. The study was conducted on 214 animals (90 Baoule, 90 Zebu and 34 composites), genotyped with 25 microsatellites across the genome and with 155 SNPs in 23 candidate regions. Degrees of admixture of composites were analyzed for microsatellite and SNP data separately. Average Baoule admixture based on microsatellites across the genomes of the Baoule-Zebu composites was 0.31, which was smaller than the average Baoule admixture in the trypanosomosis candidate regions of 0.37 (P=0.15). Fixation index FST measured in the overall genome based on microsatellites or with SNPs from candidate regions indicates strong differentiation between breeds. Nine out of 23 regions had FST ≥ 0.20 calculated from haplotypes or individual SNPs. The levels of admixture were significantly different from background admixture, as revealed by microsatellite data, for six out of the nine regions. Five out of the six regions showed an excess of Baoule ancestry. Information about best levels of breed composition would be useful for future breeding activities, aiming at trypanotolerant animals with higher productive capacity