Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding.

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Clinical exome sequencing of 1000 families with complex immune phenotypes: Toward comprehensive genomic evaluations.

BACKGROUND: Prospective genetic evaluation of patients at this referral research hospital presents clinical research challenges. OBJECTIVES: This study sought not only a single-gene explanation for participants' immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune phenotype and other heritable comorbidities relevant to their presentation and health. METHODS: This study developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity. RESULTS: Probands were 50.8% female, 71.5% were ≥18 years, and had diverse immune presentations. Overall, 327 of 1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9 of 22, 40.9%). One-quarter of the molecular diagnoses (92 of 361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251 of 361 of these molecular diagnoses (69.5%) could translate into ≥1 management option. CONCLUSIONS: This program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale.

Blood DNA methylation and breast cancer risk: a meta-analysis of four prospective cohort studies.

BACKGROUND: Environmental and genetic factors play an important role in the etiology of breast cancer. Several small blood-based DNA methylation studies have reported risk associations with methylation at individual CpGs and average methylation levels; however, these findings require validation in larger prospective cohort studies. To investigate the role of blood DNA methylation on breast cancer risk, we conducted a meta-analysis of four prospective cohort studies, including a total of 1663 incident cases and 1885 controls, the largest study of blood DNA methylation and breast cancer risk to date. METHODS: We assessed associations with methylation at 365,145 CpGs present in the HumanMethylation450 (HM450K) Beadchip, after excluding CpGs that did not pass quality controls in all studies. Each of the four cohorts estimated odds ratios (ORs) and 95% confidence intervals (CI) for the association between each individual CpG and breast cancer risk. In addition, each study assessed the association between average methylation measures and breast cancer risk, adjusted and unadjusted for cell-type composition. Study-specific ORs were combined using fixed-effect meta-analysis with inverse variance weights. Stratified analyses were conducted by age at diagnosis (< 50, ≥ 50), estrogen receptor (ER) status (+/-), and time since blood collection (< 5, 5-10, > 10 years). The false discovery rate (q value) was used to account for multiple testing. RESULTS: The average age at blood draw ranged from 52.2 to 62.2 years across the four cohorts. Median follow-up time ranged from 6.6 to 8.4 years. The methylation measured at individual CpGs was not associated with breast cancer risk (q value > 0.59). In addition, higher average methylation level was not associated with risk of breast cancer (OR = 0.94, 95% CI = 0.85, 1.05; P = 0.26; P for study heterogeneity = 0.86). We found no evidence of modification of this association by age at diagnosis (P = 0.17), ER status (P = 0.88), time since blood collection (P = 0.98), or CpG location (P = 0.98). CONCLUSIONS: Our data indicate that DNA methylation measured in the blood prior to breast cancer diagnosis in predominantly postmenopausal women is unlikely to be associated with substantial breast cancer risk on the HM450K array. Larger studies or with greater methylation coverage are needed to determine if associations exist between blood DNA methylation and breast cancer risk.

Biallelic BRCA2 Mutations Shape the Somatic Mutational Landscape of Aggressive Prostate Tumors.

To identify clinically important molecular subtypes of prostate cancer (PCa), we characterized the somatic landscape of aggressive tumors via deep, whole-genome sequencing. In our discovery set of ten tumor/normal subject pairs with Gleason scores of 8-10 at diagnosis, coordinated analysis of germline and somatic variants, including single-nucleotide variants, indels, and structural variants, revealed biallelic BRCA2 disruptions in a subset of samples. Compared to the other samples, the PCa BRCA2-deficient tumors exhibited a complex and highly specific mutation signature, featuring a 2.88-fold increased somatic mutation rate, depletion of context-specific C>T substitutions, and an enrichment for deletions, especially those longer than 10 bp. We next performed a BRCA2 deficiency-targeted reanalysis of 150 metastatic PCa tumors, and each of the 18 BRCA2-mutated samples recapitulated the BRCA2 deficiency-associated mutation signature, underscoring the potent influence of these lesions on somatic mutagenesis and tumor evolution. Among all 21 individuals with BRCA2-deficient tumors, only about half carried deleterious germline alleles. Importantly, the somatic mutation signature in tumors with one germline and one somatic risk allele was indistinguishable from those with purely somatic mutations. Our observations clearly demonstrate that BRCA2-disrupted tumors represent a unique and clinically relevant molecular subtype of aggressive PCa, highlighting both the promise and utility of this mutation signature as a prognostic and treatment-selection biomarker. Further, any test designed to leverage BRCA2 status as a biomarker for PCa must consider both germline and somatic mutations and all types of deleterious mutations.

Successful use of whole genome amplified DNA from multiple source types for high-density Illumina SNP microarrays.

BACKGROUND: The recommended genomic DNA input requirements for whole genome single nucleotide polymorphism microarrays can limit the scope of molecular epidemiological studies. We performed a large-scale evaluation of whole genome amplified DNA as input into high-density, whole-genome Illumina® Infinium® SNP microarray. RESULTS: Overall, 6622 DNA samples from 5970 individuals were obtained from three distinct biospecimen sources and genotyped using gDNA and/or wgaDNA inputs. When genotypes from the same individual were compared with standard, native gDNA input amount, we observed 99.94% mean concordance with wgaDNA input. CONCLUSIONS: Our results demonstrate that carefully conducted studies with wgaDNA inputs can yield high-quality genotyping results. These findings should enable investigators to consider expansion of ongoing studies using high-density SNP microarrays, currently challenged by small amounts of available DNA.

Genome-wide association study identifies the GLDC/IL33 locus associated with survival of osteosarcoma patients.

Survival rates for osteosarcoma, the most common primary bone cancer, have changed little over the past three decades and are particularly low for patients with metastatic disease. We conducted a multi-institutional genome-wide association study (GWAS) to identify germline genetic variants associated with overall survival in 632 patients with osteosarcoma, including 523 patients of European ancestry and 109 from Brazil. We conducted a time-to-event analysis and estimated hazard ratios (HR) and 95% confidence intervals (CI) using Cox proportional hazards models, with and without adjustment for metastatic disease. The results were combined across the European and Brazilian case sets using a random-effects meta-analysis. The strongest association after meta-analysis was for rs3765555 at 9p24.1, which was inversely associated with overall survival (HR = 1.76; 95% CI 1.41-2.18, p = 4.84 × 10-7 ). After imputation across this region, the combined analysis identified two SNPs that reached genome-wide significance. The strongest single association was with rs55933544 (HR = 1.9; 95% CI 1.5-2.4; p = 1.3 × 10-8 ), which localizes to the GLDC gene, adjacent to the IL33 gene and was consistent across both the European and Brazilian case sets. Using publicly available data, the risk allele was associated with lower expression of IL33 and low expression of IL33 was associated with poor survival in an independent set of patients with osteosarcoma. In conclusion, we have identified the GLDC/IL33 locus on chromosome 9p24.1 as associated with overall survival in patients with osteosarcoma. Further studies are needed to confirm this association and shed light on the biological underpinnings of this susceptibility locus.

Comparison against 186 canid whole-genome sequences reveals survival strategies of an ancient clonally transmissible canine tumor.

Canine transmissible venereal tumor (CTVT) is a parasitic cancer clone that has propagated for thousands of years via sexual transfer of malignant cells. Little is understood about the mechanisms that converted an ancient tumor into the world's oldest known continuously propagating somatic cell lineage. We created the largest existing catalog of canine genome-wide variation and compared it against two CTVT genome sequences, thereby separating alleles derived from the founder's genome from somatic mutations that must drive clonal transmissibility. We show that CTVT has undergone continuous adaptation to its transmissible allograft niche, with overlapping mutations at every step of immunosurveillance, particularly self-antigen presentation and apoptosis. We also identified chronologically early somatic mutations in oncogenesis- and immune-related genes that may represent key initiators of clonal transmissibility. Thus, we provide the first insights into the specific genomic aberrations that underlie CTVT's dogged perseverance in canids around the world.

Higher-than-expected population prevalence of potentially pathogenic germline TP53 variants in individuals unselected for cancer history.

Li-Fraumeni syndrome (LFS) is an autosomal-dominant cancer predisposition disorder associated with pathogenic germline variants in TP53, with a high penetrance over an individual's lifetime. The actual population prevalence of pathogenic germline TP53 mutations is still unclear, most likely due to biased selection of cancer affected families. The aim of this study was to estimate the population prevalence of potentially pathogenic TP53 exonic variants in three sequencing databases, totaling 63,983 unrelated individuals. Potential pathogenicity was defined using an original algorithm combining bioinformatic prediction tools, suggested clinical significance, and functional data. We identified 34 different potentially pathogenic TP53 variants in 131 out of 63,983 individuals (0.2%). Twenty-eight (82%) of these variants fell within the DNA-binding domain of TP53, with an enrichment for specific variants that were not previously identified as LFS mutation hotspots, such as the p.R290H and p.N235S variants. Our findings reveal that the population prevalence of potentially pathogenic TP53 variants may be up to 10 times higher than previously estimated from family-based studies. These results point to the need for further studies aimed at evaluating cancer penetrance modifiers as well as the risk associated between cancer and rare TP53 variants.

A copy number variant at the KITLG locus likely confers risk for canine squamous cell carcinoma of the digit.

The domestic dog is a robust model for studying the genetics of complex disease susceptibility. The strategies used to develop and propagate modern breeds have resulted in an elevated risk for specific diseases in particular breeds. One example is that of Standard Poodles (STPOs), who have increased risk for squamous cell carcinoma of the digit (SCCD), a locally aggressive cancer that causes lytic bone lesions, sometimes with multiple toe recurrence. However, only STPOs of dark coat color are at high risk; light colored STPOs are almost entirely unaffected, suggesting that interactions between multiple pathways are necessary for oncogenesis. We performed a genome-wide association study (GWAS) on STPOs, comparing 31 SCCD cases to 34 unrelated black STPO controls. The peak SNP on canine chromosome 15 was statistically significant at the genome-wide level (P(raw) = 1.60 × 10(-7); P(genome) = 0.0066). Additional mapping resolved the region to the KIT Ligand (KITLG) locus. Comparison of STPO cases to other at-risk breeds narrowed the locus to a 144.9-Kb region. Haplotype mapping among 84 STPO cases identified a minimal region of 28.3 Kb. A copy number variant (CNV) containing predicted enhancer elements was found to be strongly associated with SCCD in STPOs (P = 1.72 × 10(-8)). Light colored STPOs carry the CNV risk alleles at the same frequency as black STPOs, but are not susceptible to SCCD. A GWAS comparing 24 black and 24 light colored STPOs highlighted only the MC1R locus as significantly different between the two datasets, suggesting that a compensatory mutation within the MC1R locus likely protects light colored STPOs from disease. Our findings highlight a role for KITLG in SCCD susceptibility, as well as demonstrate that interactions between the KITLG and MC1R loci are potentially required for SCCD oncogenesis. These findings highlight how studies of breed-limited diseases are useful for disentangling multigene disorders.

DNA methylation age in paired tumor and adjacent normal breast tissue in Chinese women with breast cancer.

BACKGROUND: Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong. METHODS: We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath's pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson's correlation (r), Kruskal-Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors. RESULTS: DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P < 2.2e-16) than in tumor tissue (Pearson r = 0.31, P = 7.8e-06). Although overall DNAm age or AA did not vary significantly by tissue within the same individual, luminal A tumors exhibited increased DNAm AA (P = 0.004) while HER2-enriched/basal-like tumors exhibited markedly lower DNAm AA (P = < .0001) compared with paired normal tissue. Consistent with the subtype association, tumor DNAm AA was positively correlated with ESR1 (Pearson r = 0.39, P = 6.3e-06) and PGR (Pearson r = 0.36, P = 2.4e-05) gene expression. In line with this, we found that increasing DNAm AA was associated with higher body mass index (P = 0.039) and earlier age at menarche (P = 0.035), factors that are related to cumulative exposure to estrogen. In contrast, variables indicating extensive genomic instability, such as TP53 somatic mutations, high tumor mutation/copy number alteration burden, and homologous repair deficiency were associated with lower DNAm AA. CONCLUSIONS: Our findings provide additional insights into the complexity of breast tissue aging that is associated with the interaction of hormonal, genomic, and epigenetic mechanisms in an East Asian population.

Chromosomal microarray analysis supplements exome sequencing to diagnose children with suspected inborn errors of immunity.

Purpose: Though copy number variants (CNVs) have been suggested to play a significant role in inborn errors of immunity (IEI), the precise nature of this role remains largely unexplored. We sought to determine the diagnostic contribution of CNVs using genome-wide chromosomal microarray analysis (CMA) in children with IEI. Methods: We performed exome sequencing (ES) and CMA for 332 unrelated pediatric probands referred for evaluation of IEI. The analysis included primary, secondary, and incidental findings. Results: Of the 332 probands, 134 (40.4%) received molecular diagnoses. Of these, 116/134 (86.6%) were diagnosed by ES alone. An additional 15/134 (11.2%) were diagnosed by CMA alone, including two likely de novo changes. Three (2.2%) participants had diagnostic molecular findings from both ES and CMA, including two compound heterozygotes and one participant with two distinct diagnoses. Half of the participants with CMA contribution to diagnosis had CNVs in at least one non-immune gene, highlighting the clinical complexity of these cases. Overall, CMA contributed to 18/134 diagnoses (13.4%), increasing the overall diagnostic yield by 15.5% beyond ES alone. Conclusion: Pairing ES and CMA can provide a comprehensive evaluation to clarify the complex factors that contribute to both immune and non-immune phenotypes. Such a combined approach to genetic testing helps untangle complex phenotypes, not only by clarifying the differential diagnosis, but in some cases by identifying multiple diagnoses contributing to the overall clinical presentation.

Fine scale mapping of the breast cancer 16q12 locus.

Recent genome-wide association studies have identified a breast cancer susceptibility locus on 16q12 with an unknown biological basis. We used a set of single nucleotide polymorphism (SNP) markers to generate a fine-scale map and narrowed the region of association to a 133 kb DNA segment containing the largely uncharacterized hypothetical gene LOC643714, a short intergenic region and the 5' end of TOX3. Re-sequencing this segment in European subjects identified 293 common polymorphisms, including a set of 26 highly correlated candidate causal variants. By evaluation of these SNPs in five breast cancer case-control studies involving more than 23 000 subjects from populations of European and Southeast Asian ancestry, all but 14 variants could be excluded at odds of <1:100. Most of the remaining variants lie in the intergenic region, which exhibits evolutionary conservation and open chromatin conformation, consistent with a regulatory function. African-American case-control studies exhibit a different pattern of association suggestive of an additional causative variant.

Immunogenetics associated with severe coccidioidomycosis.

Disseminated coccidioidomycosis (DCM) is caused by Coccidioides, pathogenic fungi endemic to the southwestern United States and Mexico. Illness occurs in approximately 30% of those infected, less than 1% of whom develop disseminated disease. To address why some individuals allow dissemination, we enrolled patients with DCM and performed whole-exome sequencing. In an exploratory set of 67 patients with DCM, 2 had haploinsufficient STAT3 mutations, and defects in ß-glucan sensing and response were seen in 34 of 67 cases. Damaging CLEC7A and PLCG2 variants were associated with impaired production of ß-glucan-stimulated TNF-α from PBMCs compared with healthy controls. Using ancestry-matched controls, damaging CLEC7A and PLCG2 variants were overrepresented in DCM, including CLEC7A Y238* and PLCG2 R268W. A validation cohort of 111 patients with DCM confirmed the PLCG2 R268W, CLEC7A I223S, and CLEC7A Y238* variants. Stimulation with a DECTIN-1 agonist induced DUOX1/DUOXA1-derived hydrogen peroxide [H2O2] in transfected cells. Heterozygous DUOX1 or DUOXA1 variants that impaired H2O2 production were overrepresented in discovery and validation cohorts. Patients with DCM have impaired ß-glucan sensing or response affecting TNF-α and H2O2 production. Impaired Coccidioides recognition and decreased cellular response are associated with disseminated coccidioidomycosis.

FGFR2 variants and breast cancer risk: fine-scale mapping using African American studies and analysis of chromatin conformation.

Genome-wide association studies have identified FGFR2 as a breast cancer (BC) susceptibility gene in populations of European and Asian descent, but a causative variant has not yet been conclusively identified. We hypothesized that the weaker linkage disequilibrium across this associated region in populations of African ancestry might help refine the set of candidate-causal single nucleotide polymorphisms (SNPs) previously identified by our group. Eight candidate-causal SNPs were evaluated in 1253 African American invasive BC cases and 1245 controls. A significant association with BC risk was found with SNP rs2981578 (unadjusted per-allele odds ratio = 1.20, 95% confidence interval 1.03-1.41, P(trend) = 0.02), with the odds ratio estimate similar to that reported in European and Asian subjects. To extend the fine-mapping, genotype data from the African American studies were analyzed jointly with data from European (n = 7196 cases, 7275 controls) and Asian (n = 3901 cases, 3205 controls) studies. In the combined analysis, SNP rs2981578 was the most strongly associated. Five other SNPs were too strongly correlated to be excluded at a likelihood ratio of < 1/100 relative to rs2981578. Analysis of DNase I hypersensitive sites indicated that only two of these map to highly accessible chromatin, one of which, SNP rs2981578, has previously been implicated in up-regulating FGFR2 expression. Our results demonstrate that the association of SNPs in FGFR2 with BC risk extends to women of African American ethnicity, and illustrate the utility of combining association analysis in datasets of diverse ethnic groups with functional experiments to identify disease susceptibility variants.

Age-related DNA methylation in paired normal and tumour breast tissue in Chinese breast cancer patients.

Age-related DNA methylation is a potential mechanism contributing to breast cancer development. Studies of primarily Caucasian women have identified many CpG sites of age-related methylation in non-diseased breast tissue possibly driving cancer development over time. There is a paucity of studies involving Asian women whose ages at breast cancer onset are usually younger than Caucasians. We identified the 181 most consistent age-related methylation events in non-diseased breast tissue across published studies. Age-related methylation events were measured in adjacent normal and breast tumour tissue in an exclusively Asian population at the previously identified age-related methylation sites. Age-related methylation was found in 118 probes in adjacent normal breast tissue. Methylation of 99% of these sites was increased with age and predominantly located on CpG islands in promoter regions. To ascertain biological relevance to breast cancer, we focused on the 37 sites with overall higher methylation in tumour compared to adjacent normal samples. Some sites positively related to age, including AQP5 and CORO6, inversely correlated with gene expression. Several others have known involvement in suppression of carcinogenesis including GPC5 and SST, suggesting that perturbation of epigenetic regulation at these sites due to ageing may contribute to the progression of carcinogenesis. This study highlights an age-related methylation landscape in non-tumour tissue, consistent not just across studies, but also across different populations. We present candidate age-related methylation sites warranting further investigation as potential epigenetic drivers of breast cancer. They may serve as potential targets of site-specific demethylation intervention strategies for the prevention of age-related breast cancer.

Subsequent Neoplasm Risk Associated With Rare Variants in DNA Damage Response and Clinical Radiation Sensitivity Syndrome Genes in the Childhood Cancer Survivor Study.

PURPOSE: Radiotherapy for childhood cancer is associated with elevated subsequent neoplasm (SN) risk, but the contribution of rare variants in DNA damage response and radiation sensitivity genes to SN risk is unknown. PATIENTS AND METHODS: We conducted whole-exome sequencing in a cohort of childhood cancer survivors originally diagnosed during 1970 to 1986 (mean follow-up, 32.7 years), with reconstruction of doses to body regions from radiotherapy records. We identified patients who developed SN types previously reported to be related to radiotherapy (RT-SNs; eg, basal cell carcinoma [BCC], breast cancer, meningioma, thyroid cancer, sarcoma) and matched controls (sex, childhood cancer type/diagnosis, age, SN location, radiation dose, survival). Conditional logistic regression assessed SN risk associated with potentially protein-damaging rare variants (SnpEff, ClinVar) in 476 DNA damage response or radiation sensitivity genes with exact permutation-based P values using a Bonferroni-corrected significance threshold of P < 8.06 × 10-5. RESULTS: Among 5,105 childhood cancer survivors of European descent, 1,108 (21.7%) developed at least 1 RT-SN. Out-of-field RT-SN risk, excluding BCC, was associated with homologous recombination repair (HRR) gene variants (patient cases, 23.2%; controls, 10.8%; odds ratio [OR], 2.6; 95% CI, 1.7 to 3.9; P = 4.79 × 10-5), most notably but nonsignificantly for FANCM (patient cases, 4.0%; matched controls, 0.6%; P = 9.64 × 10-5). HRR variants were not associated with likely in/near-field RT-SNs, excluding BCC (patient cases, 12.7%; matched controls, 12.9%; P = .92). Irrespective of radiation dose, risk for RT-SNs was also associated with EXO1 variants (patient cases, 1.8%; controls, 0.4%; P = 3.31 × 10-5), another gene implicated in DNA double-strand break repair. CONCLUSION: In this large-scale discovery study, we identified novel associations between RT-SN risk after childhood cancer and potentially protein-damaging rare variants in genes involved in DNA double-strand break repair, particularly HRR. With replication, these results could affect screening recommendations for childhood cancer survivors and risk-benefit assessments of treatment approaches.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma.

BACKGROUND: Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor. METHODS: We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry). RESULTS: We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8). CONCLUSIONS: These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk. IMPACT: Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in Patients With Osteosarcoma.

Importance: Osteosarcoma, the most common malignant bone tumor in children and adolescents, occurs in a high number of cancer predisposition syndromes that are defined by highly penetrant germline mutations. The germline genetic susceptibility to osteosarcoma outside of familial cancer syndromes remains unclear. Objective: To investigate the germline genetic architecture of 1244 patients with osteosarcoma. Design, Setting, and Participants: Whole-exome sequencing (n = 1104) or targeted sequencing (n = 140) of the DNA of 1244 patients with osteosarcoma from 10 participating international centers or studies was conducted from April 21, 2014, to September 1, 2017. The results were compared with the DNA of 1062 individuals without cancer assembled internally from 4 participating studies who underwent comparable whole-exome sequencing and 27 173 individuals of non-Finnish European ancestry who were identified through the Exome Aggregation Consortium (ExAC) database. In the analysis, 238 high-interest cancer-susceptibility genes were assessed followed by testing of the mutational burden across 736 additional candidate genes. Principal component analyses were used to identify 732 European patients with osteosarcoma and 994 European individuals without cancer, with outliers removed for patient-control group comparisons. Patients were subsequently compared with individuals in the ExAC group. All data were analyzed from June 1, 2017, to July 1, 2019. Main Outcomes and Measures: The frequency of rare pathogenic or likely pathogenic genetic variants. Results: Among 1244 patients with osteosarcoma (mean [SD] age at diagnosis, 16 [8.9] years [range, 2-80 years]; 684 patients [55.0%] were male), an analysis restricted to individuals with European ancestry indicated a significantly higher pathogenic or likely pathogenic variant burden in 238 high-interest cancer-susceptibility genes among patients with osteosarcoma compared with the control group (732 vs 994, respectively; P = 1.3 × 10-18). A pathogenic or likely pathogenic cancer-susceptibility gene variant was identified in 281 of 1004 patients with osteosarcoma (28.0%), of which nearly three-quarters had a variant that mapped to an autosomal-dominant gene or a known osteosarcoma-associated cancer predisposition syndrome gene. The frequency of a pathogenic or likely pathogenic cancer-susceptibility gene variant was 128 of 1062 individuals (12.1%) in the control group and 2527 of 27 173 individuals (9.3%) in the ExAC group. A higher than expected frequency of pathogenic or likely pathogenic variants was observed in genes not previously linked to osteosarcoma (eg, CDKN2A, MEN1, VHL, POT1, APC, MSH2, and ATRX) and in the Li-Fraumeni syndrome-associated gene, TP53. Conclusions and Relevance: In this study, approximately one-fourth of patients with osteosarcoma unselected for family history had a highly penetrant germline mutation requiring additional follow-up analysis and possible genetic counseling with cascade testing.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data.

BACKGROUND: Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference. RESULTS: Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3). CONCLUSIONS: Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility.

Ewing sarcoma (EWS) is a pediatric cancer characterized by the EWSR1-FLI1 fusion. We performed a genome-wide association study of 733 EWS cases and 1346 unaffected individuals of European ancestry. Our study replicates previously reported susceptibility loci at 1p36.22, 10q21.3 and 15q15.1, and identifies new loci at 6p25.1, 20p11.22 and 20p11.23. Effect estimates exhibit odds ratios in excess of 1.7, which is high for cancer GWAS, and striking in light of the rarity of EWS cases in familial cancer syndromes. Expression quantitative trait locus (eQTL) analyses identify candidate genes at 6p25.1 (RREB1) and 20p11.23 (KIZ). The 20p11.22 locus is near NKX2-2, a highly overexpressed gene in EWS. Interestingly, most loci reside near GGAA repeat sequences and may disrupt binding of the EWSR1-FLI1 fusion protein. The high locus to case discovery ratio from 733 EWS cases suggests a genetic architecture in which moderate risk SNPs constitute a significant fraction of risk.