Genome-wide discovery of somatic coding and noncoding mutations in pediatric endemic and sporadic Burkitt lymphoma.

Although generally curable with intensive chemotherapy in resource-rich settings, Burkitt lymphoma (BL) remains a deadly disease in older patients and in sub-Saharan Africa. Epstein-Barr virus (EBV) positivity is a feature in more than 90% of cases in malaria-endemic regions, and up to 30% elsewhere. However, the molecular features of BL have not been comprehensively evaluated when taking into account tumor EBV status or geographic origin. Through an integrative analysis of whole-genome and transcriptome data, we show a striking genome-wide increase in aberrant somatic hypermutation in EBV-positive tumors, supporting a link between EBV and activation-induced cytidine deaminase (AICDA) activity. In addition to identifying novel candidate BL genes such as SIN3A, USP7, and CHD8, we demonstrate that EBV-positive tumors had significantly fewer driver mutations, especially among genes with roles in apoptosis. We also found immunoglobulin variable region genes that were disproportionally used to encode clonal B-cell receptors (BCRs) in the tumors. These include IGHV4-34, known to produce autoreactive antibodies, and IGKV3-20, a feature described in other B-cell malignancies but not yet in BL. Our results suggest that tumor EBV status defines a specific BL phenotype irrespective of geographic origin, with particular molecular properties and distinct pathogenic mechanisms. The novel mutation patterns identified here imply rational use of DNA-damaging chemotherapy in some patients with BL and targeted agents such as the CDK4/6 inhibitor palbociclib in others, whereas the importance of BCR signaling in BL strengthens the potential benefit of inhibitors for PI3K, Syk, and Src family kinases among these patients.

Sources of erroneous sequences and artifact chimeric reads in next generation sequencing of genomic DNA from formalin-fixed paraffin-embedded samples.

Tissues used in pathology laboratories are typically stored in the form of formalin-fixed, paraffin-embedded (FFPE) samples. One important consideration in repurposing FFPE material for next generation sequencing (NGS) analysis is the sequencing artifacts that can arise from the significant damage to nucleic acids due to treatment with formalin, storage at room temperature and extraction. One such class of artifacts consists of chimeric reads that appear to be derived from non-contiguous portions of the genome. Here, we show that a major proportion of such chimeric reads align to both the 'Watson' and 'Crick' strands of the reference genome. We refer to these as strand-split artifact reads (SSARs). This study provides a conceptual framework for the mechanistic basis of the genesis of SSARs and other chimeric artifacts along with supporting experimental evidence, which have led to approaches to reduce the levels of such artifacts. We demonstrate that one of these approaches, involving S1 nuclease-mediated removal of single-stranded fragments and overhangs, also reduces sequence bias, base error rates, and false positive detection of copy number and single nucleotide variants. Finally, we describe an analytical approach for quantifying SSARs from NGS data.

Analysis of Ugandan cervical carcinomas identifies human papillomavirus clade-specific epigenome and transcriptome landscapes.

Cervical cancer is the most common cancer affecting sub-Saharan African women and is prevalent among HIV-positive (HIV+) individuals. No comprehensive profiling of cancer genomes, transcriptomes or epigenomes has been performed in this population thus far. We characterized 118 tumors from Ugandan patients, of whom 72 were HIV+, and performed extended mutation analysis on an additional 89 tumors. We detected human papillomavirus (HPV)-clade-specific differences in tumor DNA methylation, promoter- and enhancer-associated histone marks, gene expression and pathway dysregulation. Changes in histone modification at HPV integration events were correlated with upregulation of nearby genes and endogenous retroviruses.

Genetic variation in H2AFX contributes to risk of non-Hodgkin lymphoma.

Non-Hodgkin lymphoma (NHL) comprises a group of lymphoid tumors that have in common somatic translocations. H2AFX encodes a key histone involved in the detection of the DNA double-stranded breaks that can lead to translocations. H2afx is a dosage-dependent gene that protects against B-cell lymphomas in mice, making its human orthologue an ideal candidate gene for susceptibility to lymphoma. We did a population-based genetic association study of H2AFX variants in 487 NHL cases and 531 controls. Complete resequencing of the human H2AFX gene in 95 NHL cases was done to establish the spectrum of variation in affected individuals; this was followed by both direct and indirect tests for association at the level of individual single nucleotide polymorphisms (SNP) and as haplotypes. Homozygosity for the AA genotype of a SNP 417 bp upstream of the translational start of H2AFX is strongly associated [odds ratio (OR), 0.54; P = 0.001] with protection from NHL. We find a strong association of this SNP with the follicular lymphoma subtype of NHL (AA genotype: OR, 0.40; P = 0.004) and with mantle cell lymphoma (AA genotype: OR, 0.20; P = 0.01) that remains significant after adjustment for the false discovery rate, but not with diffuse large B-cell lymphoma. These data support the hypothesis that genetic variation in the H2AFX gene influences genetic susceptibility or resistance to some subtypes of NHL by contributing to the maintenance of genome stability.

DNA methylation profiling of the human major histocompatibility complex: a pilot study for the human epigenome project.

The Human Epigenome Project aims to identify, catalogue, and interpret genome-wide DNA methylation phenomena. Occurring naturally on cytosine bases at cytosine-guanine dinucleotides, DNA methylation is intimately involved in diverse biological processes and the aetiology of many diseases. Differentially methylated cytosines give rise to distinct profiles, thought to be specific for gene activity, tissue type, and disease state. The identification of such methylation variable positions will significantly improve our understanding of genome biology and our ability to diagnose disease. Here, we report the results of the pilot study for the Human Epigenome Project entailing the methylation analysis of the human major histocompatibility complex. This study involved the development of an integrated pipeline for high-throughput methylation analysis using bisulphite DNA sequencing, discovery of methylation variable positions, epigenotyping by matrix-assisted laser desorption/ionisation mass spectrometry, and development of an integrated public database available at http://www.epigenome.org. Our analysis of DNA methylation levels within the major histocompatibility complex, including regulatory exonic and intronic regions associated with 90 genes in multiple tissues and individuals, reveals a bimodal distribution of methylation profiles (i.e., the vast majority of the analysed regions were either hypo- or hypermethylated), tissue specificity, inter-individual variation, and correlation with independent gene expression data.

Sequence variant discovery in DNA repair genes from radiosensitive and radiotolerant prostate brachytherapy patients.

PURPOSE: The presence of intrinsic radiosensitivity within prostate cancer patients may be an important factor contributing to development of radiation toxicity. We investigated whether variants in genes responsible for detecting and repairing DNA damage independently contribute to toxicity following prostate brachytherapy. EXPERIMENTAL DESIGN: Genomic DNA was extracted from blood samples of 41 prostate brachytherapy patients, 21 with high and 20 with low late toxicity scores. For each patient, 242 PCR amplicons were generated containing 173 exons of eight candidate genes: ATM, BRCA1, ERCC2, H2AFX, LIG4, MDC1, MRE11A, and RAD50. These amplicons were sequenced and all sequence variants were subjected to statistical analysis to identify those associated with late radiation toxicity. RESULTS: Across 41 patients, 239 sites differed from the human genome reference sequence; 170 of these corresponded to known polymorphisms. Sixty variants, 14 of them novel, affected protein coding regions and 43 of these were missense mutations. In our patient population, the high toxicity group was enriched for individuals with at least one LIG4 coding variant (P = 0.028). One synonymous variant in MDC1, rs28986317, was associated with increased radiosensitivity (P = 0.048). A missense variant in ATM, rs1800057, associated with increased prostate cancer risk, was found exclusively in two high toxicity patients but did not reach statistical significance for association with radiosensitivity (P = 0.488). CONCLUSIONS: Our data revealed new germ-line sequence variants, indicating that existing sequence databases do not fully represent the full extent of sequence variation. Variants in three DNA repair genes were linked to increased radiosensitivity but require validation in larger populations.

Polymorphisms in DNA repair genes and epithelial ovarian cancer risk.

DNA repair gene polymorphisms and mutations are known to influence cancer risk. We studied whether polymorphisms in DNA double strand break (DSB) repair genes are associated with epithelial ovarian cancer (EOC) risk. Up to 1,600 cases and 4,241 controls from 4 separate genetic association studies from 3 countries were genotyped for 13 single nucleotide polymorphisms (SNP) in 6 genes (BRCA1, NBS1, RAD51, RAD52, XRCC2 and XRCC3) involved in homologous recombination of DNA double strand breaks. Genotype specific risks were estimated as odds ratios (OR) by unconditional logistic regression. No association was detected between EOC risk and BRCA1 Q356R, BRCA1 P871L, RAD51 g135c, RAD51 g172t, RAD52 c2259t, NBS1 L34L, NBS1 E185Q, NBS1 A399A, NBS1 P672P, XRCC2 g4324c, XRCC2 c41657t and XRCC3 T241M. The XRCC2 R188H polymorphism was associated with a modest reduction in EOC risk: OR for heterozygotes was 0.8 (95% confidence interval [CI] = 0.7-1.0) and for rare homozygotes 0.3 (0.1-0.9). The XRCC3 a4541g polymorphism, situated in the 5'UTR, and the intronic XRCC3 a17893g polymorphism were not associated with EOC risk in general, but when the serous EOC subset only was analysed, the OR for heterozygotes for a4541g was 1.0 (0.9-1.2) and for the rare homozygotes 0.5 (0.3-0.9). For the XRCC3 a17893g polymorphism, the OR for the heterozygotes and the rare homozygotes were 0.8 (0.7-0.9) and 0.9 (0.7-1.2), respectively. In our study, some polymorphisms in XRCC2 and XRCC3 genes were associated with EOC risk. Further research on the role of these genes on epithelial ovarian cancer is warranted.

Variants in DNA double-strand break repair genes and breast cancer susceptibility.

We performed genetic association studies in a population-based breast cancer case-control study analysing polymorphisms in genes involved in homologous recombination (NBS1, RAD52, RAD51, XRCC2 and XRCC3) and non-homologous end-joining (KU70/80 and LIG4). These DNA double-strand break repair genes are candidates for breast cancer susceptibility. Genotype results were available for up to 2205 cases and 1826 controls. In the homologous recombination (HR) pathway, genotype frequencies differed between cases and controls for two polymorphisms in XRCC3; T241M (P=0.015) and IVS5 A>G at nt 17893 (P=0.008). Homozygous carriers of M241 were associated with an increased risk [odds ratio (OR) MM versus TT=1.3 (95% confidence interval (CI) 1.1-1.6)], while the rare allele of IVS5A>G was associated with a dominant protective effect [OR AG versus AA=0.8 (0.7-0.9)]. The association of a rare variant in XRCC2 (R188H) was marginally significant [P=0.07; OR HH versus RR=2.6 (1.0-6.7)]. In the non-homologous end-joining (NHEJ) pathway, a polymorphism in LIG4 (T>C at nt 1977) was associated with a decrease in breast cancer risk [P=0.09; OR CC versus TT=0.7 (0.4-1.0)]. No significant association was found for 12 other polymorphisms in the other genes studied. For XRCC3, we found evidence for four common haplotypes and four rarer ones that appear to have arisen by recombination. Two haplotypes, AGC and GGC, were associated with non-significant reductions in breast cancer risk, and the rare GAT haplotype was associated with a significantly increased risk. These data provide some evidence that variants in XRCC2 and LIG4 alter breast cancer risk, together with stronger evidence that variants of XRCC3 are associated with risk. If these results can be confirmed, understanding the functional basis should improve our understanding of the role of DNA repair in breast carcinogenesis.

Polymorphisms associated with circulating sex hormone levels in postmenopausal women.

BACKGROUND: Reports suggest a relationship between circulating sex hormone levels and breast cancer risk, but genetic association studies have been inconclusive. We investigated the association between levels of sex hormones and single nucleotide polymorphisms (SNPs) in genes coding for the enzymes that regulate them. METHODS: We assayed circulating levels of estradiol, testosterone, estrone, androstenedione, 17alpha-hydroxyprogesterone, and sex hormone-binding globulin (SHBG) in 1975 normal postmenopausal women. Fifteen SNPs in the CYP17, CYP19, EDH17B2, SHBG, COMT, and CYP1B1 genes were genotyped in these postmenopausal women and in a breast cancer case-control study. Associations of genotypes with breast cancer risk were evaluated in the case-control study and with hormone levels in the postmenopausal women using multiple linear regression with assay batch, body mass index, parity, peri- or postmenopausal status, and age band as covariates. RESULTS: CYP19 SNPs (rs10046 and [TCT]+/-) were associated with differences in estradiol level (P =.0006 and P =.0003, respectively) and the estradiol : testosterone ratio (P =.000001() and P =.002). SNP rs10046 explained 1.6% of the variance (r2) in the estradiol : testosterone ratio. SHBG SNPs (5' untranslated region [5'UTR] g-a and D356N) were associated with both SHBG levels (P<10(-6) and P =.005) and the estradiol : SHBG ratio (P =().000008() and P =.01). These SNPs explained 2.4% and 0.6% of the variance in SHBG levels, respectively. SNPs in the other genes were not associated with differences in any hormone levels, and none were statistically significantly associated with breast cancer risk. CONCLUSION: Genetic variation in CYP19 and SHBG contributes to variance in circulating hormone levels between postmenopausal women, but low r2 values may explain why these genes have given inconclusive results in breast cancer case-control studies.

BRCA2 Arg372Hispolymorphism and epithelial ovarian cancer risk.

The BRCA2 372 HH genotype defined by the BRCA2 N372H nonconservative amino acid substitution polymorphism was recently reported to be associated with a small increased risk of breast cancer. We investigated whether this polymorphism was associated with ovarian cancer risk by conducting British and Australian case-control comparisons in parallel, including a total sample of 1,121 ovarian cancer cases and 2,643 controls. There was no difference in genotype frequency between control groups from the 2 studies (p = 0.9). The HH genotype was associated with an increased risk of ovarian cancer in both studies, and the risk estimate for the pooled studies was 1.36 (95% CI 1.04-1.77, p = 0.03). There was also a suggestion that this risk may be greater for ovarian cancers of the serous subtype for both studies, with an OR (95% CI) of 1.66 (1.17-2.54) for the 2 studies combined (p = 0.005). The BRCA2 372 HH genotype appears to be associated with an increased risk of ovarian cancer of a similar magnitude to that reported for breast cancer.