Genomic characterization of cervical lymph node metastases in papillary thyroid carcinoma following the Chornobyl accident.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers (131I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions.

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus.

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Insights into the mechanisms and structure of breakage-fusion-bridge cycles in cervical cancer using long-read sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and our study applies long-read sequencing to this important chromosomal rearrangement type. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes, and only one BFB breakpoint showed chromothripsis. Five cell lines have a chromosomal region 11q BFB event, with YAP1-BIRC3-BIRC2 amplification. Indeed, YAP1 amplification is associated with a 10-year-earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that individuals with cervical cancer and YAP1-BIRC3-BIRC2 amplification, especially those of African ancestry, might benefit from targeted therapy. In summary, we uncovered valuable insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

Human leukocyte antigen-DQA1*04:01 and rs2040406 variants are associated with elevated risk of childhood Burkitt lymphoma.

Burkitt lymphoma (BL) is responsible for many childhood cancers in sub-Saharan Africa, where it is linked to recurrent or chronic infection by Epstein-Barr virus or Plasmodium falciparum. However, whether human leukocyte antigen (HLA) polymorphisms, which regulate immune response, are associated with BL has not been well investigated, which limits our understanding of BL etiology. Here we investigate this association among 4,645 children aged 0-15 years, 800 with BL, enrolled in Uganda, Tanzania, Kenya, and Malawi. HLA alleles are imputed with accuracy >90% for HLA class I and 85-89% for class II alleles. BL risk is elevated with HLA-DQA1*04:01 (adjusted odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.32-1.97, P = 3.71 × 10-6), with rs2040406(G) in HLA-DQA1 region (OR = 1.43, 95% CI = 1.26-1.63, P = 4.62 × 10-8), and with amino acid Gln at position 53 versus other variants in HLA-DQA1 (OR = 1.36, P = 2.06 × 10-6). The associations with HLA-DQA1*04:01 (OR = 1.29, P = 0.03) and rs2040406(G) (OR = 1.68, P = 0.019) persist in mutually adjusted models. The higher risk rs2040406(G) variant for BL is associated with decreased HLA-DQB1 expression in eQTLs in EBV transformed lymphocytes. Our results support the role of HLA variation in the etiology of BL and suggest that a promising area of research might be understanding the link between HLA variation and EBV control.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa.

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Glyphosate Use and Mosaic Loss of Chromosome Y among Male Farmers in the Agricultural Health Study.

BACKGROUND: Glyphosate is the most commonly used herbicide worldwide and has been implicated in the development of certain hematologic cancers. Although mechanistic studies in human cells and animals support the genotoxic effects of glyphosate, evidence in human populations is scarce. OBJECTIVES: We evaluated the association between lifetime occupational glyphosate use and mosaic loss of chromosome Y (mLOY) as a marker of genotoxicity among male farmers. METHODS: We analyzed blood-derived DNA from 1,606 farmers ≥50 years of age in the Biomarkers of Exposure and Effect in Agriculture study, a subcohort of the Agricultural Health Study. mLOY was detected using genotyping array intensity data in the pseudoautosomal region of the sex chromosomes. Cumulative lifetime glyphosate use was assessed using self-reported pesticide exposure histories. Using multivariable logistic regression, we estimated odds ratios (ORs) and 95% confidence intervals (CIs) for the associations between glyphosate use and any detectable mLOY (overall mLOY) or mLOY affecting ≥10% of cells (expanded mLOY). RESULTS: Overall, mLOY was detected in 21.4% of farmers, and 9.8% of all farmers had expanded mLOY. Increasing total lifetime days of glyphosate use was associated with expanded mLOY [highest vs. lowest quartile; OR=1.75 (95% CI: 1.00, 3.07), ptrend=0.03] but not with overall mLOY; the associations with expanded mLOY were most apparent among older (≥70 years of age) men [OR=2.30 (95% CI: 1.13, 4.67), ptrend=0.01], never smokers [OR=2.32 (95% CI: 1.04, 5.21), ptrend=0.04], and nonobese men [OR=2.04 (95% CI: 0.99, 4.19), ptrend=0.03]. Similar patterns of associations were observed for intensity-weighted lifetime days of glyphosate use. DISCUSSION: High lifetime glyphosate use could be associated with mLOY affecting a larger fraction of cells, suggesting glyphosate could confer genotoxic or selective effects relevant for clonal expansion. As the first study to investigate this association, our findings contribute novel evidence regarding the carcinogenic potential of glyphosate and require replication in future studies. https://doi.org/10.1289/EHP12834.

Mosaic chromosomal alterations in peripheral blood leukocytes of children in sub-Saharan Africa.

In high-income countries, mosaic chromosomal alterations in peripheral blood leukocytes are associated with an elevated risk of adverse health outcomes, including hematologic malignancies. We investigate mosaic chromosomal alterations in sub-Saharan Africa among 931 children with Burkitt lymphoma, an aggressive lymphoma commonly characterized by immunoglobulin-MYC chromosomal rearrangements, 3822 Burkitt lymphoma-free children, and 674 cancer-free men from Ghana. We find autosomal and X chromosome mosaic chromosomal alterations in 3.4% and 1.7% of Burkitt lymphoma-free children, and 8.4% and 3.7% of children with Burkitt lymphoma (P-values = 5.7×10-11 and 3.74×10-2, respectively). Autosomal mosaic chromosomal alterations are detected in 14.0% of Ghanaian men and increase with age. Mosaic chromosomal alterations in Burkitt lymphoma cases include gains on chromosomes 1q and 8, the latter spanning MYC, while mosaic chromosomal alterations in Burkitt lymphoma-free children include copy-neutral loss of heterozygosity on chromosomes 10, 14, and 16. Our results highlight mosaic chromosomal alterations in sub-Saharan African populations as a promising area of research.

Insights into the Mechanisms and Structure of Breakage-Fusion-Bridge Cycles in Cervical Cancer using Long-Read Sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and this is one of the first analyses of these events using long-read sequencing. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes and only one BFB breakpoint showed chromothripsis. Five cell lines have a Chr11q BFB event, with YAP1/BIRC2/BIRC3 gene amplification. Indeed, YAP1 amplification is associated with a 10-year earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that cervical cancer patients with YAP1/BIRC2/BIRC3-amplification, especially those of African American ancestry, might benefit from targeted therapy. In summary, we uncovered new insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

In utero exposure to diethylstilbestrol and blood DNA methylation in adult women: Results from a meta-analysis of two cohort studies.

BACKGROUND: Prenatal exposure to diethylstilbestrol (DES) is associated with several adverse health outcomes. Animal studies have shown associations between prenatal DES exposure and DNA methylation. OBJECTIVE: The aim of this study was to explore blood DNA methylation in women exposed and unexposed to DES in utero. METHODS: Sixty women (40 exposed and 20 unexposed) in the National Cancer Institute's Combined DES Cohort Study and 199 women (99 exposed and 100 unexposed women) in the Sister Study Cohort were included in this analysis. Within each study, robust linear regression models were used to assess associations between DES exposure and blood DNA methylation. Study-specific associations were combined using fixed-effect meta-analysis with inverse variance weights. Our analysis focused on CpG sites located within nine candidate genes identified in animal models. We further explored whether in utero DES exposure was associated with age acceleration. RESULTS: Blood DNA methylation levels at 10 CpG sites in six of the nine candidate genes were statistically significantly associated with prenatal DES exposure (P < 0.05) in this meta-analysis. Genes included EGF, EMB, EGFR, WNT11, FOS, and TGFB1, which are related to cell proliferation and differentiation. The most statistically significant CpG site was cg19830739 in gene EGF, and it was associated with lower methylation levels in women prenatally exposed to DES compared with those not exposed (P < 0.0001; false discovery rate<0.05). The association between prenatal DES exposure in utero and age acceleration was not statistically significant (P = 0.07 for meta-analyzed results). CONCLUSIONS: There are few opportunities to investigate the effects of prenatal DES exposure. These findings suggest that in utero DES exposure may be associated with differential blood DNA methylation levels, which could mediate the increased risk of several adverse health outcomes observed in exposed women. Our findings need further evaluation using larger data sets.

Germline genetic variants and pediatric rhabdomyosarcoma outcomes: a report from the Children's Oncology Group.

BACKGROUND: Relative to other pediatric cancers, survival for rhabdomyosarcoma (RMS) has not improved in recent decades, suggesting the need to enhance risk stratification. Therefore, we conducted a genome-wide association study for event-free survival (EFS) and overall survival (OS) to identify genetic variants associated with outcomes in individuals with RMS. METHODS: The study included 920 individuals with newly diagnosed RMS who were enrolled in Children's Oncology Group protocols. To assess the association of each single nucleotide polymorphism (SNP) with EFS and OS, we estimated hazard ratios (HRs) and 95% confidence intervals (CIs) using multivariable Cox proportional hazards models, adjusted for clinical covariates. All statistical tests were two sided. We also performed stratified analyses by histological subtype (alveolar and embryonal RMS) and carried out sensitivity analyses of statistically significant SNPs by PAX3/7-FOXO1 fusion status and genetic ancestry group. RESULTS: We identified that rs17321084 was associated with worse EFS (HR = 2.01, 95% CI = 1.59 to 2.53, P = 5.39 × 10-9) and rs10094840 was associated with worse OS (HR = 1.84, 95% CI = 1.48 to 2.27, P = 2.13 × 10-8). Using publicly available data, we found that rs17321084 lies in a binding region for transcription factors GATA2 and GATA3, and rs10094840 is associated with SPAG1 and RNF19A expression. We also identified that CTNNA3 rs2135732 (HR = 3.75, 95% CI = 2.34 to 5.99, P = 3.54 × 10-8) and MED31 rs74504320 (HR = 3.21, 95% CI = 2.12 to 4.86, P = 3.60 × 10-8) were associated with worse OS among individuals with alveolar RMS. CONCLUSIONS: We demonstrated that common germline variants are associated with EFS and OS among individuals with RMS. Additional replication and investigation of these SNP effects may further support their consideration in risk stratification protocols.

Extrachromosomal Amplification of Human Papillomavirus Episomes Is a Mechanism of Cervical Carcinogenesis.

SIGNIFICANCE: Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism for oncogenesis without integration.

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.

The genetic and evolutionary basis of gene expression variation in East Africans.

BACKGROUND: Mapping of quantitative trait loci (QTL) associated with molecular phenotypes is a powerful approach for identifying the genes and molecular mechanisms underlying human traits and diseases, though most studies have focused on individuals of European descent. While important progress has been made to study a greater diversity of human populations, many groups remain unstudied, particularly among indigenous populations within Africa. To better understand the genetics of gene regulation in East Africans, we perform expression and splicing QTL mapping in whole blood from a cohort of 162 diverse Africans from Ethiopia and Tanzania. We assess replication of these QTLs in cohorts of predominantly European ancestry and identify candidate genes under selection in human populations. RESULTS: We find the gene regulatory architecture of African and non-African populations is broadly shared, though there is a considerable amount of variation at individual loci across populations. Comparing our analyses to an equivalently sized cohort of European Americans, we find that QTL mapping in Africans improves the detection of expression QTLs and fine-mapping of causal variation. Integrating our QTL scans with signatures of natural selection, we find several genes related to immunity and metabolism that are highly differentiated between Africans and non-Africans, as well as a gene associated with pigmentation. CONCLUSION: Extending QTL mapping studies beyond European ancestry, particularly to diverse indigenous populations, is vital for a complete understanding of the genetic architecture of human traits and can reveal novel functional variation underlying human traits and disease.

GWAS Explorer: an open-source tool to explore, visualize, and access GWAS summary statistics in the PLCO Atlas.

The Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial is a prospective cohort study of nearly 155,000 U.S. volunteers aged 55-74 at enrollment in 1993-2001. We developed the PLCO Atlas Project, a large resource for multi-trait genome-wide association studies (GWAS), by genotyping participants with available DNA and genomic consent. Genotyping on high-density arrays and imputation was performed, and GWAS were conducted using a custom semi-automated pipeline. Association summary statistics were generated from a total of 110,562 participants of European, African and Asian ancestry. Application programming interfaces (APIs) and open-source software development kits (SKDs) enable exploring, visualizing and open data access through the PLCO Atlas GWAS Explorer website, promoting Findable, Accessible, Interoperable, and Re-usable (FAIR) principles. Currently the GWAS Explorer hosts association data for 90 traits and >78,000,000 genomic markers, focusing on cancer and cancer-related phenotypes. New traits will be posted as association data becomes available. The PLCO Atlas is a FAIR resource of high-quality genetic and phenotypic data with many potential reuse opportunities for cancer research and genetic epidemiology.

HPV16 E7 Nucleotide Variants Found in Cancer-Free Subjects Affect E7 Protein Expression and Transformation.

The human papillomavirus (HPV) type 16 E7 oncogene is critical to carcinogenesis and highly conserved. Previous studies identified a preponderance of non-synonymous E7 variants amongst HPV16-positive cancer-free controls compared to those with cervical cancer. To investigate the function of E7 variants, we constructed full-length HPV16 E7 genes and tested variants at positions H9R, D21N, N29S, E33K, T56I, D62N, S63F, S63P, T64M, E80K, D81N, P92L, and P92S (found only in controls); D14E, N29H cervical intraepithelial neoplasia (CIN2), and P6L, H51N, R77S (CIN3). We determined the steady-state level of cytoplasmic and nuclear HPV16 E7 protein. All variants from controls showed a reduced level of E7 protein, with 7/13 variants having lower protein levels. In contrast, 2/3 variants from the CIN3 precancer group had near-wild type E7 levels. We assayed the activity of representative variants in stably transfected NIH3T3 cells. The H9R, E33K, P92L, and P92S variants found in control subjects had lower transforming activity than D14E and N29H variants (CIN2), and the R77S (CIN3) had activity only slightly reduced from wild-type E7. In addition, R77S and WT E7 caused increased migration of NIH3T3 cells in a wound-healing assay compared with H9R, E33K, P92L, and P92S (controls) and D14E (CIN2). These data provide evidence that the E7 variants found in HPV16-positive cancer-free women are partially defective for transformation and cell migration, further demonstrating the importance of fully active E7 in cancer development.

Germline-somatic JAK2 interactions are associated with clonal expansion in myelofibrosis.

Myelofibrosis is a rare myeloproliferative neoplasm (MPN) with high risk for progression to acute myeloid leukemia. Our integrated genomic analysis of up to 933 myelofibrosis cases identifies 6 germline susceptibility loci, 4 of which overlap with previously identified MPN loci. Virtual karyotyping identifies high frequencies of mosaic chromosomal alterations (mCAs), with enrichment at myelofibrosis GWAS susceptibility loci and recurrently somatically mutated MPN genes (e.g., JAK2). We replicate prior MPN associations showing germline variation at the 9p24.1 risk haplotype confers elevated risk of acquiring JAK2V617F mutations, demonstrating with long-read sequencing that this relationship occurs in cis. We also describe recurrent 9p24.1 large mCAs that selectively retained JAK2V617F mutations. Germline variation associated with longer telomeres is associated with increased myelofibrosis risk. Myelofibrosis cases with high-frequency JAK2 mCAs have marked reductions in measured telomere length - suggesting a relationship between telomere biology and myelofibrosis clonal expansion. Our results advance understanding of the germline-somatic interaction at JAK2 and implicate mCAs involving JAK2 as strong promoters of clonal expansion of those mutated clones.

Breast Cancer Risk in Women from Ghana Carrying Rare Germline Pathogenic Mutations.

BACKGROUND: Risk estimates for women carrying germline mutations in breast cancer susceptibility genes are mainly based on studies of European ancestry women. METHODS: We investigated associations between pathogenic variants (PV) in 34 genes with breast cancer risk in 871 cases [307 estrogen receptor (ER)-positive, 321 ER-negative, and 243 ER-unknown] and 1,563 controls in the Ghana Breast Health Study (GBHS), and estimated lifetime risk for carriers. We compared results with those for European, Asian, and African American ancestry women. RESULTS: The frequency of PV in GBHS for nine breast cancer genes was 8.38% in cases and 1.22% in controls. Relative risk estimates for overall breast cancer were: (OR, 13.70; 95% confidence interval (CI), 4.03-46.51) for BRCA1, (OR, 7.02; 95% CI, 3.17-15.54) for BRCA2, (OR, 17.25; 95% CI, 2.15-138.13) for PALB2, 5 cases and no controls carried TP53 PVs, and 2.10, (0.72-6.14) for moderate-risk genes combined (ATM, BARD1, CHEK2, RAD51C, RAD52D). These estimates were similar to those previously reported in other populations and were modified by ER status. No other genes evaluated had mutations associated at P < 0.05 with overall risk. The estimated lifetime risks for mutation carriers in BRCA1, BRCA2, and PALB2 and moderate-risk genes were 18.4%, 9.8%, 22.4%, and 3.1%, respectively, markedly lower than in Western populations with higher baseline risks. CONCLUSIONS: We confirmed associations between PV and breast cancer risk in Ghanaian women and provide absolute risk estimates that could inform counseling in Ghana and other West African countries. IMPACT: These findings have direct relevance for breast cancer genetic counseling for women in West Africa.

Rare germline variants in PALB2 and BRCA2 in familial and sporadic chordoma.

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma.

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.