Pan-cancer mutational signature analysis of 111,711 targeted sequenced tumors using SATS

Tumor mutational signatures have the potential to inform cancer diagnosis and treatment. However, their detection in targeted sequenced tumors is hampered by sparse mutations and variability in targeted gene panels. Here we present SATS, a scalable mutational signature analyzer addressing these challenges by leveraging tumor mutational burdens from targeted gene panels. Through analyzing simulated data, pseudo-targeted sequencing data generated by down-sampling whole exome and genome data, and samples with matched whole genome sequencing and targeted sequencing, we showed that SATS can accurately detect common mutational signatures and estimate signature burdens. Applying SATS to 111,711 targeted sequenced tumors from the AACR Project GENIE, we generated a pan-cancer catalogue of mutational signatures tailored to targeted sequencing, enabling estimation of signature burdens within individual tumors. Integrating signatures with clinical data, we demonstrated SATS's clinical utility, including identifying signatures enriched in early-onset hypermutated colorectal cancers and signatures associated with cancer prognosis and immunotherapy response.

Population-specific Mutation Patterns in Breast Tumors from African American, European American, and Kenyan Patients

Women of African descent have the highest breast cancer mortality in the United States and are more likely than women from other population groups to develop an aggressive disease. It remains uncertain to what extent breast cancer in Africa is reminiscent of breast cancer in African American or European American patients. Here, we performed whole-exome sequencing of genomic DNA from 191 breast tumor and non-cancerous adjacent tissue pairs obtained from 97 African American, 69 European American, 2 Asian American, and 23 Kenyan patients. Our analysis of the sequencing data revealed an elevated tumor mutational burden in both Kenyan and African American patients, when compared with European American patients. TP53 mutations were most prevalent, particularly in African American patients, followed by PIK3CA mutations, which showed similar frequencies in European American, African American, and the Kenyan patients. Mutations targeting TBX3 were confined to European Americans and those targeting the FBXW7 tumor suppressor to African American patients whereas mutations in the ARID1A gene that are known to confer resistance to endocrine therapy were distinctively enriched among Kenyan patients. A Kyoto Encyclopedia of Genes and Genomes pathway analysis could link FBXW7 mutations to an increased mitochondrial oxidative phosphorylation capacity in tumors carrying these mutations. Finally, Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures in tumors correlated with the occurrence of driver mutations, immune cell profiles, and neighborhood deprivation with associations ranging from being mostly modest to occasionally robust. To conclude, we found mutational profiles that were different between these patient groups. The differences concentrated among genes with low mutation frequencies in breast cancer.

SIGNIFICANCE

The study describes differences in tumor mutational profiles between African American, European American, and Kenyan breast cancer patients. It also investigates how these profiles may relate to the tumor immune environment and the neighborhood environment in which the patients had residence. Finally, it describes an overrepresentation of ARID1A gene mutations in breast tumors of the Kenyan patients.

[Chinese expert guidance on overall application of lenvatinib in hepatocellular carcinoma]

Lenvatinib mesylate is an oral receptor tyrosine kinase inhibitor against targets of vascular endothelial growth factor receptors 1-3, fibroblast growth factor receptors 1-4, platelet-derived growth factor receptor α, stem cell growth factor receptor, and rearranged during transfection, et al. Lenvatinib has been approved by the National Medical Products Administration of China on September 4, 2018, for the first-line treatment of patients with unresectable hepatocellular carcinoma who have not received systematic treatment before. Up to February 2023, Lenvatinib has been listed in China for more than 4 years, accumulating a series of post-marketing clinical research evidences. Based on the clinical practice before and after the launch of lenvatinib and referring to the clinical experience of other anti-angiogenesis inhibitors, domestic multidisciplinary experts and scholars adopt the Delphi method to formulate the Chinese Expert Guidance on Overall Application of Lenvatinib in Hepatocellular Carcinoma after repeated discussions and revisions, in order to provide reference for reasonable and effective clinical application of lenvatinib for clinicians.

Estimated Prevalence, Tumor Spectrum, and Neurofibromatosis Type 1-Like Phenotype of CDKN2A-Related Melanoma-Astrocytoma Syndrome

Importance

Knowledge about the prevalence and tumor types of CDKN2A-related melanoma-astrocytoma syndrome (MAS) is limited and could improve disease recognition.

Objective

To estimate the prevalence and describe the tumor types of MAS.

Design, Setting, and Participants

This retrospective cohort study analyzed all available MAS cases from medical centers in the US (2 sites) and Europe (2 sites) and from biomedical population genomic databases (UK Biobank [United Kingdom], Geisinger MyCode [US]) between January 1, 1976, and December 31, 2020. Patients with MAS with CDKN2A germline pathogenic variants and 1 or more neural tumors were included. Data were analyzed from June 1, 2022, to January 31, 2023.

Main Outcomes and Measures

Disease prevalence and tumor frequency.

Results

Prevalence of MAS ranged from 1 in 170 503 (n = 1 case; 95% CI, 1:30 098-1:965 887) in Geisinger MyCode (n = 170 503; mean [SD] age, 58.9 [19.1] years; 60.6% women; 96.2% White) to 1 in 39 149 (n = 12 cases; 95% CI, 1:22 396-1:68 434) in UK Biobank (n = 469 789; mean [SD] age, 70.0 [8.0] years; 54.2% women; 94.8% White). Among UK Biobank patients with MAS (n = 12) identified using an unbiased genomic ascertainment approach, brain neoplasms (4 of 12, 33%; 1 glioblastoma, 1 gliosarcoma, 1 astrocytoma, 1 unspecified type) and schwannomas (3 of 12, 25%) were the most common malignant and benign neural tumors, while cutaneous melanoma (2 of 12, 17%) and head and neck squamous cell carcinoma (2 of 12, 17%) were the most common nonneural malignant neoplasms. In a separate case series of 14 patients with MAS from the US and Europe, brain neoplasms (4 of 14, 29%; 2 glioblastomas, 2 unspecified type) and malignant peripheral nerve sheath tumor (2 of 14, 14%) were the most common neural cancers, while cutaneous melanoma (4 of 14, 29%) and sarcomas (2 of 14, 14%; 1 liposarcoma, 1 unspecified type) were the most common nonneural cancers. Cutaneous neurofibromas (7 of 14, 50%) and schwannomas (2 of 14, 14%) were also common. In 1 US family, a father and son with MAS had clinical diagnoses of neurofibromatosis type 1 (NF1). Genetic testing of the son detected a pathogenic CDKN2A splicing variant (c.151-1G>C) and was negative for NF1 genetic alterations. In UK Biobank, 2 in 150 (1.3%) individuals with clinical NF1 diagnoses had likely pathogenic variants in CDKN2A, including 1 individual with no detected variants in the NF1 gene.

Conclusions and Relevance

This cohort study estimates the prevalence and describes the tumors of MAS. Additional studies are needed in genetically diverse populations to further define population prevalence and disease phenotypes.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Spectrum and Frequency of Germline FANCM Protein-Truncating Variants in 44,803 European Female Breast Cancer Cases

FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.

Association of germline variants in telomere maintenance genes (POT1, TERF2IP, ACD, and TERT) with spitzoid morphology in familial melanoma: A multi-center case series

Background

Spitzoid morphology in familial melanoma has been associated with germline variants in POT1, a telomere maintenance gene (TMG), suggesting a link between telomere biology and spitzoid differentiation.

Objective

To assess if familial melanoma cases associated with germline variants in TMG (POT1, ACD, TERF2IP, and TERT) commonly exhibit spitzoid morphology.

Methods

In this case series, melanomas were classified as having spitzoid morphology if at least 3 of 4 dermatopathologists reported this finding in ≥25% of tumor cells. Logistic regression was used to calculate odds ratios (OR) of spitzoid morphology compared to familial melanomas from unmatched noncarriers that were previously reviewed by a National Cancer Institute dermatopathologist.

Results

Spitzoid morphology was observed in 77% (23 of 30), 75% (3 of 4), 50% (2 of 4), and 50% (1 of 2) of melanomas from individuals with germline variants in POT1, TERF2IP, ACD, and TERT, respectively. Compared to noncarriers (n = 139 melanomas), POT1 carriers (OR = 225.1, 95% confidence interval: 51.7-980.5; P < .001) and individuals with TERF2IP, ACD, and TERT variants (OR = 82.4, 95% confidence interval: 21.3-494.6; P < .001) had increased odds of spitzoid morphology.

Limitations

Findings may not be generalizable to nonfamilial melanoma cases.

Conclusion

Spitzoid morphology in familial melanoma could suggest germline alteration of TMG.

Abstract 5722: The mutational signatures of 100,477 targeted sequenced tumors

Mutational signatures, the footprints of somatic mutations, are associated with the causes of cancer and have been well-studied for tumors with whole exome or genome sequencing. However, due to the low number of detected mutations, mutational signatures were insufficiently explored for many tumors sequenced by targeted panels in clinics. It impedes the clinical application of mutational signatures. Here, we present a new method, SALMON (Signature Analyzer for Low Mutation cOuNts), to identify mutational signatures in targeted sequenced tumors based on tumor mutational burden. SALMON adjusts for panel size differences and uses a large number of targeted sequenced tumors rather than a large number of mutations per tumor (as with whole exome or genome sequencing) to overcome the challenges of mutational signature analysis for targeted sequencing. Extensive simulations and pseudo-targeted sequenced data show that SALMON can accurately detect spiky or common signatures. We applied SALMON to investigate the pan-cancer patterns of mutational signatures for 100,477 targeted sequenced tumors in AACR Project GENIE, including 14,428 lung and 11,389 breast tumors. We detected well-established signatures in tumor types that have not previously been associated with these signatures, such as the smoking signature in ovarian tumors. Interestingly, analysis of thousands of tumors per cancer type from diverse populations revealed gender discrepancies, self-described race differences, subtype heterogeneity, and metastatic enrichment of mutational signatures. For instance, most sex-biased signatures are more frequently present in males for non-gender-specific cancers. Thiopurine treatment-induced signature in glioma is enriched in Black patients. And endometrioid ovarian or uterine cancers have a higher prevalence of polymerase epsilon (POLE) deficiency-related signatures than non-endometrioid ovarian or uterine cancers, respectively. Our study demonstrates the feasibility and utility of mutational signature analysis for targeted sequenced tumors, enabling precision applications of mutational signatures in the clinical setting.

Citation Format: Donghyuk Lee, Min Hua, Difei Wang, Lei Song, Tongwu Zhang, Kai Yu, Xiaohong R. Yang, Jianxin Shi, Stephen J. Chanock, Maria Teresa Landi, Bin Zhu. The mutational signatures of 100,477 targeted sequenced tumors. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5722.

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.

The impact of coding germline variants on contralateral breast cancer risk and survival

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry

BACKGROUND

Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes.

METHODS

We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.

RESULTS

In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C.

CONCLUSIONS

Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Gene Expression Profiling Identifies Two Chordoma Subtypes Associated with Distinct Molecular Mechanisms and Clinical Outcomes

PURPOSE

Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management.

EXPERIMENTAL DESIGN

We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America.

RESULTS

Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial-mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes.

CONCLUSIONS

Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

Distinct Reproductive Risk Profiles for Intrinsic-Like Breast Cancer Subtypes: Pooled Analysis of Population-Based Studies

BACKGROUND

Reproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.

METHODS

Analyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2-like, HER2-enriched-like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.

RESULTS

Compared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2-like, and HER2-enriched-like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.

CONCLUSIONS

This large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

[Shandong hilly rural natural population cohort study: method and baseline characteristics of survey subjects]

Objective: To understand the study method and the baseline characteristics of the survey subjects of Shandong hilly rural natural population cohort study, and provide reference for the research of the prevalence and risk factors of common chronic and non-communicable diseases. Methods: Baseline survey, including questionnaire survey, physical examination, biochemical index examination and blood and saliva collection, was conducted in local residents aged 20-79 years in Kongcun and Xiaozhi townships of Pingyin county, Shandong province, from 2017 to 2019. Shandong hilly rural natural population cohort was established and main baseline characteristics of the study subjects were statistically analyzed. Results: A total of 10 296 study subjects aged 54.45 years were included in the study, in whom 40.6% were males. Among the study subjects, 88.3% had education level of junior high school or below, 62.1% were famers, and 90.7% were married. Smokers accounted for 45.6% of men and 0.9% of women, and drinkers accounted for 65.8% of men and 3.0% of women, respectively. The self-reported rates of hypertension, diabetes, coronary heart disease, stroke and tumors were 19.8%, 3.2%, 2.8%, 2.7% and 1.2%, respectively. Conclusion: The Shandong hilly rural cohort natural population study provided important evidence for assessing the risk for common chronic and non-communicable diseases and disease prevention and control in hilly rural areas.

Gene-Level Associations in Patients With and Without Pathogenic Germline Variants in CDKN2A and Pancreatic Cancer

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) is a component of familial melanoma due to germline pathogenic variants (GPVs) in CDKN2A. However, it is unclear what role this gene or other genes play in its etiology.

MATERIALS AND METHODS

We analyzed 189 cancer predisposition genes using parametric rare-variant association (RVA) tests and nonparametric permutation tests to identify gene-level associations in PDAC for patients with (CDKN2A+) and without (CDKN2A-) GPV. Exome sequencing was performed on 84 patients with PDAC, 47 CDKN2A+ and 37 CDKN2A-. After variant filtering, various RVA tests and permutation tests were run separately by CDKN2A status. Genes with the strongest nominal associations were evaluated in patients with PDAC from The Cancer Genome Atlas and the UK Biobank (UKB). A secondary analysis including only GPV from UKB was also performed.

RESULTS

In RVA tests, ERCC4 and RET showed the most compelling evidence as plausible PDAC candidate genes for CDKN2A+ patients. In contrast, the findings in CDKN2A- patients provided evidence for HMBS, EPCAM, and MRE11 as potential new candidate genes and confirmed ATM, BRCA2, and PALB2 as PDAC genes, consistent with findings in The Cancer Genome Atlas and the UKB. As expected, CDKN2A- patients were more likely to harbor GPVs from the 189 genes investigated. When including only GPVs from UKB, significant associations with PDAC were seen for ATM, BRCA2, and CDKN2A.

CONCLUSION

These results suggest that variants in other genes likely play a role in PDAC in all patients and that PDAC in CDKN2A+ patients has a distinct etiology from PDAC in CDKN2A- patients.

Impact of Transcript (p16/p14ARF) Alteration on Cancer Risk in CDKN2A Germline Pathogenic Variant Carriers

BACKGROUND

Few studies have evaluated the relationship between CDKN2A germline pathogenic variants (GPV), transcript (p16, p14ARF) alteration, and cancer risk.

METHODS

Standardized incidence ratios (SIRs) comparing cancer risk to the general population were calculated for 385 CDKN2A GPV carriers from two large cohorts (259 United States and 126 Swedish individuals) using Poisson regression; statistical significance was defined as P<.002 (Bonferroni correction). Cumulative incidence is reported for melanoma and non-melanoma cancer.

RESULTS

Incidence was significantly increased for melanoma (SIR = 159.8, 95% confidence interval [CI] = 132.1 to 193.2), pancreatic cancer (SIR = 24.1, 95% CI = 14.7 to 39.4), head and neck squamous cell carcinoma (HNSCC; SIR = 16.2, 95% CI = 9.5 to 27.6), and lung cancer (SIR = 5.6, 95% CI = 3.4 to 9.1) in GPV carriers. Similar associations were observed with p16 alteration. Combined p16 and p14ARF alteration was associated with increased incidence of esophageal cancer (SIR = 16.7, 95% CI = 5.7 to 48.9) and malignant peripheral nerve sheath tumor (SIR = 113.0, 95% CI = 16.4 to 780.9), although cancer events were limited (n < 5 for each malignancy). Cumulative incidence at age 70 for melanoma and non-melanoma cancer was 68.3% (95% CI = 68.0 to 68.6) and 35.2% (95% CI = 34.9 to 35.6), respectively. 89% of smoking-related cancers (lung, HNSCC, pancreatic, esophageal) occurred in ever smokers.

CONCLUSION

These findings highlight the impact of p16 and p14ARF alteration on cancer risk. Smoking was an important risk factor for smoking-related cancers in our study.

Rare germline deleterious variants increase susceptibility for lung cancer

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

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.

Abstract 3670: Systematic literature review of risk factor associations with breast cancer subtypes in women of African, Asian, Hispanic, and European descents

Background: Studies have reported differences in associations between breast cancer risk factors and subtypes defined by hormone receptor status, particularly estrogen receptor (ER) positive versus ER negative tumors. Most studies have been conducted in women of European descent, and an expanding body of literature in other populations suggests differences by race/ethnicity. Clarifying whether associations between risk factors and disease subtypes are consistent across racial/ethnic populations has important implications for understanding disease etiology and for improving risk prediction models. To address this question, we conducted a qualitative literature review to investigate the consistency in associations between multiple breast cancer risk factors and risk of tumor subtypes in women of African, Asian, Hispanic, and European descent.

Methods: We searched PubMed for publications between January 1, 1990 and June 8, 2021 that reported associations between breast cancer risk factors and risk of subtypes of the disease. We evaluated 19 risk factors, including reproductive, anthropometric, lifestyle, and diet factors, as well as medical history and use of menopausal hormone therapy (MHT). Subtypes were defined as ER positive or negative (specifically triple negative when available). We prioritized review studies (i.e., meta-analyses, systematic reviews, and pooled analyses) and included individual studies with populations not included in the reviews. The number of publications per risk factor ranged from 3 for calcium intake to 28 for parity, with up to 7 reviews per risk factor. Most publications reported estimates for women of European descent, followed by Asian, African, and Hispanic. Evidence of subtype heterogeneity was determined by expert review.

Results: There was strong evidence of association between reproductive factors and MHT with risk of ER positive but not ER negative disease. Parity, younger age at first birth and older age at menarche were associated with lower risk of ER positive disease, while MHT use was associated with higher risk of ER positive disease. For these risk factors, we did not find evidence that risk associations or etiologic heterogeneity varied by race/ethnicity. For the other risk factors, there was limited or no evidence of heterogeneity in risk associations by subtype, which was consistent across racial/ethnic groups. Few publications, however, specifically studied women of non-European ancestry.

Conclusion: For most breast cancer risk factors, evidence is insufficient to evaluate differences in ER-specific risk associations by race/ethnicity. More studies are needed to evaluate subtype-specific breast cancer risk factors in diverse populations, which will be important for informing the development of multi-ethnic/racial BC risk prediction models that account for subtype heterogeneity.

Citation Format: Amber N. Hurson, Thomas U. Ahearn, Renske Keeman, Mustapha Abubakar, Audrey Y. Jung, Pooja Middha Kapoor, Hela Koka, Xiaohong R. Yang, Jenny Chang-Claude, Elena Martínez, Rulla M. Tamimi, Melissa A. Troester, Elisa V. Bandera, Marjanka K. Schmidt, Montserrat Garcia-Closas. Systematic literature review of risk factor associations with breast cancer subtypes in women of African, Asian, Hispanic, and European descents [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3670.

Breast cancer risks associated with missense variants in breast cancer susceptibility genes

BACKGROUND

Protein truncating variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2 are associated with increased breast cancer risk, but risks associated with missense variants in these genes are uncertain.

METHODS

We analyzed data on 59,639 breast cancer cases and 53,165 controls from studies participating in the Breast Cancer Association Consortium BRIDGES project. We sampled training (80%) and validation (20%) sets to analyze rare missense variants in ATM (1146 training variants), BRCA1 (644), BRCA2 (1425), CHEK2 (325), and PALB2 (472). We evaluated breast cancer risks according to five in silico prediction-of-deleteriousness algorithms, functional protein domain, and frequency, using logistic regression models and also mixture models in which a subset of variants was assumed to be risk-associated.

RESULTS

The most predictive in silico algorithms were Helix (BRCA1, BRCA2 and CHEK2) and CADD (ATM). Increased risks appeared restricted to functional protein domains for ATM (FAT and PIK domains) and BRCA1 (RING and BRCT domains). For ATM, BRCA1, and BRCA2, data were compatible with small subsets (approximately 7%, 2%, and 0.6%, respectively) of rare missense variants giving similar risk to those of protein truncating variants in the same gene. For CHEK2, data were more consistent with a large fraction (approximately 60%) of rare missense variants giving a lower risk (OR 1.75, 95% CI (1.47-2.08)) than CHEK2 protein truncating variants. There was little evidence for an association with risk for missense variants in PALB2. The best fitting models were well calibrated in the validation set.

CONCLUSIONS

These results will inform risk prediction models and the selection of candidate variants for functional assays and could contribute to the clinical reporting of gene panel testing for breast cancer susceptibility.

SUITOR: Selecting the number of mutational signatures through cross-validation

For de novo mutational signature analysis, the critical first step is to decide how many signatures should be expected in a cancer genomics study. An incorrect number could mislead downstream analyses. Here we present SUITOR (Selecting the nUmber of mutatIonal signaTures thrOugh cRoss-validation), an unsupervised cross-validation method that requires little assumptions and no numerical approximations to select the optimal number of signatures without overfitting the data. In vitro studies and in silico simulations demonstrated that SUITOR can correctly identify signatures, some of which were missed by other widely used methods. Applied to 2,540 whole-genome sequenced tumors across 22 cancer types, SUITOR selected signatures with the smallest prediction errors and almost all signatures of breast cancer selected by SUITOR were validated in an independent breast cancer study. SUITOR is a powerful tool to select the optimal number of mutational signatures, facilitating downstream analyses with etiological or therapeutic importance.

Mutational signatures are the footprints of exogenous exposures and endogenous mutational processes on the cancer genomes. To estimate de novo mutational signatures, the first step is to decide how many signatures should be extracted in a cancer genomics study, which determines downstream analytical steps and has been insufficiently studied. We developed SUITOR, an unsupervised cross-validation method to select the optimal number of signatures without overfitting the data. We demonstrated SUITOR’s superior performance using in vitro experimental studies, in silico simulations and in vivo pan-cancer applications of 2,540 whole-genome sequenced tumors across 22 cancer types, and validated signatures of breast cancer in additional 440 breast tumors. SUITOR advances the methodological frontier of identifying de novo mutational signatures and would help discover the causes of cancer and the means of cancer prevention and treatment.

Associations of Genetic Ancestry with Terminal Duct Lobular Unit Involution among Healthy Women

Reduced age-related terminal duct lobular unit (TDLU) involution has been linked to increased breast cancer risk and triple-negative breast cancer (TNBC). Associations of TDLU involution levels with race and ethnicity remain incompletely explored. Herein, we examined associations between genetic ancestry and TDLU involution in normal breast tissue donated by 2,014 healthy women in the US. Women of African ancestry were more likely than European women to have increased TDLU counts (odds ratio [OR]trend=1.36; 95% CI = 1.07-1.74), acini counts/TDLU (OR = 1.47; 95% CI = 1.06-2.03), and median TDLU span (ORtrend=1.44; 95% CI = 1.08-1.91), indicating lower involution; whereas East Asian descendants were associated with decreased TDLU counts (ORtrend=0.52; 95% CI = 0.35-0.78) after controlling for potential confounders. These associations are consistent with the racial variations in incidence rates of TNBC in the US and suggest opportunities for future work examining whether TDLU involution may mediate the racial differences in subtype-specific breast cancer risk.

Abstract P3-01-26: Mammographic density in relation to breast cancer risk factors among Chinese women

Higher mammographic density (MD) is a known risk factor for breast cancer (BC), however, studies examining determinants of density in Chinese populations are limited. Asian women tend to have higher MD compared to European women, and yet they have a lower overall BC incidence rate. The goal of this study was to investigate the influence of established BC risk factors on quantitatively measured MD among Chinese women. The study population included 7,351 Chinese women, who had mammograms at the National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China, between 2016 and 2017 for different indications, including but not limited to screening exam, breast pain, and suspected breast tumors. Only subjects with BI-RADS (Breast Imaging-Reporting and Data System) diagnostic classification of 1 (negative, n= 2,293), 2 (benign, n= 3,519), and 3 (probably benign, n= 1,539) were included in this analysis. VolparaDensity software was used to obtain quantitative MD measures including absolute dense breast volume and percent breast density. We conducted multivariable linear regression modeling on each density measurement to test for its associations with well-established breast cancer risk factors. The mean age of this study population was 50.1 years old (range= (25, 82)), the average body mass index (BMI) was 24.0 kg/m2 (range= (14.2, 48.9)), and about 46% of women were postmenopausal. The mean absolute dense and percent dense volumes were 58.4 (SD= 32.1) and 14.8 (SD= 7.1, equivalent to Volpara Density Grade= 3), respectively. Associations between quantitative MD measures and BC risk factors did not vary significantly across the three BI-RADS diagnostic classes and we therefore combined these patients in subsequent analyses. Most risk factors we examined showed similar associations with dense volume and percent density. For example, both dense volume and percent density were inversely associated with age, parity, longer duration (&gt; 1 year) of breastfeeding, and postmenopausal status, while they were positively associated with age at menopause. In contrast, the associations with age at menarche and BMI showed opposite directions for dense volume and percent density. While increasing BMI was positively associated with dense volume, it was inversely associated with percent density. Age at menarche, however, was inversely correlated with dense volume, but positively associated with percent density. When stratifying by menopausal status, the associations for age, BMI, and parity were similar among pre- and post-menopausal women, while the associations for age at first birth and breastfeeding were stronger among postmenopausal women. In general, observed breast cancer risk factor associations with quantitative breast density metrics in this Chinese population are similar to those previously reported in Western women. We did not observe striking differences in associations of breast density with reproductive factors and BMI among Asian women in this large analysis based on quantitative density measurements. Other risk factors such as genetic susceptibility loci associated with MD should be examined in large Asian studies to further our understanding on how MD influences breast cancer risk in Asian populations.

Citation Format: Hela Koka, Yuan Tian, Deng Lu, Kai Yu, Er-Ni Li, Changyuan Guo, Bin Zhu, Jennifer L Guida, Ariane Chan, Nan Hu, Ning Lu, Gretchen Gierach, Jing Li, Xiaohong R Yang. Mammographic density in relation to breast cancer risk factors among Chinese women [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-01-26.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Comparison of somatic mutation landscapes in Chinese versus European breast cancer patients

Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Novel loss-of-function variant in DENND5A impedes melanosomal cargo transport and predisposes to familial cutaneous melanoma

PURPOSE

More than half of the familial cutaneous melanomas have unknown genetic predisposition. This study aims at characterizing a novel melanoma susceptibility gene.

METHODS

We performed exome and targeted sequencing in melanoma-prone families without any known melanoma susceptibility genes. We analyzed the expression of candidate gene DENND5A in melanoma samples in relation to pigmentation and UV signature. Functional studies were carried out using microscopic approaches and zebrafish model.

RESULTS

We identified a novel DENND5A truncating variant that segregated with melanoma in a Swedish family and 2 additional rare DENND5A variants, 1 of which segregated with the disease in an American family. We found that DENND5A is significantly enriched in pigmented melanoma tissue. Our functional studies show that loss of DENND5A function leads to decrease in melanin content in vitro and pigmentation defects in vivo. Mechanistically, harboring the truncating variant or being suppressed leads to DENND5A losing its interaction with SNX1 and its ability to transport the SNX1-associated vesicles from melanosomes. Consequently, untethered SNX1-premelanosome protein and redundant tyrosinase are redirected to lysosomal degradation by default, causing decrease in melanin content.

CONCLUSION

Our findings provide evidence of a physiological role of DENND5A in the skin context and link its variants to melanoma susceptibility.

[Effect of placenta previa attached to cesarean scar for adverse pregnant outcomes in patients with placenta accreta spectrum disorders]

Objective: To investigate the effect of placenta previa attached to cesarean scar for adverse pregnant outcomes in patients with or without placenta accreta spectrum disorders (PAS). Methods: The clinical information of patients with cesarean section history and placenta previa during the perioperative period at Peking University First Hospital from January 1st, 2015 to December 31st, 2020 were collected retrospectively. There were 53 cases without PAS and 172 cases with PAS, 153 cases with abnormally invasive placenta (containing placenta increta and placenta percreta) and 72 cases without PAS or with placenta accreta. The pregnant outcomes including rate of postpartum hemorrhage, transfusion, hysterectomy between the above groups were compared. Multivariate analysis was performed to study the factors significantly associated with PAS. Results: Pregnant women with PAS were at higher risk of adverse pregnant outcomes than those without PAS. Patients with PAS had higher incidences of hysterectomy [12.2% (21/172) vs 0(0/53); P=0.005], postpartum hemorrhage [60.5% (104/172) vs 5.7% (3/53); P<0.01] and blood transfusion [66.9% (115/172) vs 7.5% (4/53); P<0.01]. In the subgroup analysis stratified by the type of PAS, patients with abnormally invasive placenta were at higher risk of hysterectomy [13.7% (21/153) vs 0 (0/72); P<0.01], postpartum hemorrhage [66.7% (102/153) vs 6.9% (5/72); χ²=70.873, P<0.01] and blood transfusion [74.5% (114/153) vs 6.9% (5/72); χ²=90.869, P<0.01]. After multiple logistic regression, the type of creta had the positive relation with postpartum hemorrhage (OR=27.622, 95%CI:9.873~77.280; P<0.01) and blood transfusion (OR=36.912, 95%CI:13.239~102.922; P<0.01). There were no significant correlations between adverse pregnant outcomes and the type of placenta previa or the times of cesarean section (all P>0.01). Conclusions: Placenta previa attached to cesarean scar without PAS or with placenta accreta could not act as the factor of predicting adverse pregnant outcomes in clinic. Placenta previa attached to cesarean scar with placenta increta or placenta percreta could increase the risk of adverse pregnant outcomes.

Prognostic significance of tumor-infiltrating lymphocytes in premenopausal, luminal breast cancer treated with adjuvant endocrine therapy

PURPOSE

Tumor-infiltrating lymphocytes (TILs) have strong prognostic value in triple-negative and HER2-enriched breast cancer, but their prognostic role in luminal breast cancer (LBC) is less clear. Here, we assessed the overall TIL levels and CD8+ T-cells in relation to the prognosis of LBC patients from China.

METHODS AND RESULTS

A total of 596 patients with LBC who were premenopausal and treated with adjuvant endocrine therapy were included. Among them, 160 cases were evaluated for CD8 by immunohistochemical (IHC) staining. Whole-section hematoxylin and eosin and IHC staining were visually assessed for stromal TILs (sTILs), stromal CD8+ T-cells (sCD8), and intratumoral CD8+ T-cells (iCD8). Multivariable analyses were used to test the associations between TILs and disease-free survival (DFS) and overall survival (OS) with the adjustment for clinicopathologic characteristics and treatment. High sTILs (≥10%) were associated with high histologic grade (P<0.001), luminal B/HER2- (P<0.001), luminal B/HER2+ subtype (P=0.002), and high Ki67 expression (≥25%; P=0.014). Similar associations were observed for sCD8 but not for iCD8. While sTILs and sCD8 were not associated with either DFS or OS, the presence of iCD8 (≥1%) was associated with better DFS in both univariate (HR=0.51, 95% CI 0.26-0.96, P=0.042) and multivariate (HR=0.48, 95% CI 0.25-0.92, P=0.027) analyses. Similar but less significant associations were found for iCD8 and OS (adjusted HR=0.35, 95% CI 0.11-1.10, P=0.073).

CONCLUSIONS

Among Chinese premenopausal patients with LBC, iCD8 demonstrated suggestive associations with favorable outcome. In contrast, although sTILs and sCD8 were associated with more aggressive tumor features, they did not appear to be associated with clinical outcome.

Magnitude and temporal trend of acne vulgaris burden in 204 countries and territories from 1990 to 2019: A analysis from the Global Burden of Disease Study 2019

BACKGROUND

Acne vulgaris is widespread across the world. Mapping the latest magnitudes and temporal trends of acne vulgaris provides the essential foundation for targeted public policies at the national, regional, and global levels.

METHODS

Complying with the framework of the Global Burden of Disease 2019 Study, the incidence, prevalence and disability-adjusted life years (DALYs) on acne vulgaris were retrieved in 204 countries and regions from 1990 to 2019. The average annual percentage change was calculated to depict the temporal trends in age-standardized rates (ASRs) of acne vulgaris burden by region, sex, and age.

RESULTS

Globally, 117.4 [95% uncertainty interval (UI): 103.0, 133.7] million incident cases of acne vulgaris, 231.2 million (95% UI: 208.2, 255.5) prevalent cases and 5.0 (95% UI: 3.0, 7.9) million DALYs were estimated in 2019, with about 47% increase compared with 1990. Moreover, the overall ASRs of acne vulgaris increased by around 0.55% annually over the past three decades. We observed huge disparities in ASRs of acne vulgaris with its changing trends in sex, location, and age. The ASR of acne vulgaris among women was around 1.3 times that of men, but the gender difference was narrowed for the pronounced increase among men. The ASRs of acne vulgaris were higher in high-income regions, but the increasing trend was more pronounced in other regions.

CONCLUSION

The burden rate of acne vulgaris continues to increase in almost all countries. Understanding the specific characteristics of acne vulgaris burden is essential to formulate more effective and targeted interventions for controlling acne burden.

Germline ATM variants predispose to melanoma: a joint analysis across the GenoMEL and MelaNostrum consortia

PURPOSE

Ataxia-Telangiectasia Mutated (ATM) has been implicated in the risk of several cancers, but establishing a causal relationship is often challenging. Although ATM single-nucleotide polymorphisms have been linked to melanoma, few functional alleles have been identified. Therefore, ATM impact on melanoma predisposition is unclear.

METHODS

From 22 American, Australian, and European sites, we collected 2,104 familial, multiple primary (MPM), and sporadic melanoma cases who underwent ATM genotyping via panel, exome, or genome sequencing, and compared the allele frequency (AF) of selected ATM variants classified as loss-of-function (LOF) and variants of uncertain significance (VUS) between this cohort and the gnomAD non-Finnish European (NFE) data set.

RESULTS

LOF variants were more represented in our study cohort than in gnomAD NFE, both in all (AF = 0.005 and 0.002, OR = 2.6, 95% CI = 1.56-4.11, p < 0.01), and familial + MPM cases (AF = 0.0054 and 0.002, OR = 2.97, p < 0.01). Similarly, VUS were enriched in all (AF = 0.046 and 0.033, OR = 1.41, 95% CI = 1.6-5.09, p < 0.01) and familial + MPM cases (AF = 0.053 and 0.033, OR = 1.63, p < 0.01). In a case-control comparison of two centers that provided 1,446 controls, LOF and VUS were enriched in familial + MPM cases (p = 0.027, p = 0.018).

CONCLUSION

This study, describing the largest multicenter melanoma cohort investigated for ATM germline variants, supports the role of ATM as a melanoma predisposition gene, with LOF variants suggesting a moderate-risk.

Quantitative Mammographic Density Measurements and Molecular Subtypes in Chinese Women With Breast Cancer

Background

Studies investigating associations between mammographic density (MD) and breast cancer subtypes have generated mixed results. We previously showed that having extremely dense breasts was associated with the human epidermal growth factor receptor-2 (HER2)-enriched subtype in Chinese breast cancer patients.

Methods

In this study, we reevaluated the MD-subtype association in 1549 Chinese breast cancer patients, using VolparaDensity software to obtain quantitative MD measures. All statistical tests were 2-sided.

Results

Compared with women with luminal A tumors, women with luminal B/HER2- (odds ratio [OR] = 1.20, 95% confidence interval [CI] = 1.04 to 1.38; P = .01), luminal B/HER2+ (OR = 1.22, 95% CI = 1.03 to 1.46; P = .03), and HER2-enriched tumors (OR = 1.30, 95% CI = 1.06 to 1.59; P = .01) had higher fibroglandular dense volume. These associations were stronger in patients with smaller tumors (<2 cm). In contrast, the triple-negative subtype was associated with lower nondense volume (OR = 0.82, 95% CI = 0.68 to 0.99; P = .04), and the association was only seen among older women (age 50 years or older).

Conclusion

Although biological mechanisms remain to be investigated, the associations for the HER2-enriched and luminal B subtypes with increasing MD may partially explain the higher prevalence of luminal B and HER2+ breast cancers previously reported in Asian women.

Breast cancer risk factors in relation to molecular subtypes in breast cancer patients from Kenya

BACKGROUND

Few studies have investigated risk factor heterogeneity by molecular subtypes in indigenous African populations where prevalence of traditional breast cancer (BC) risk factors, genetic background, and environmental exposures show marked differences compared to European ancestry populations.

METHODS

We conducted a case-only analysis of 838 pathologically confirmed BC cases recruited from 5 groups of public, faith-based, and private institutions across Kenya between March 2012 to May 2015. Centralized pathology review and immunohistochemistry (IHC) for key markers (ER, PR, HER2, EGFR, CK5-6, and Ki67) was performed to define subtypes. Risk factor data was collected at time of diagnosis through a questionnaire. Multivariable polytomous logistic regression models were used to determine associations between BC risk factors and tumor molecular subtypes, adjusted for clinical characteristics and risk factors.

RESULTS

The median age at menarche and first pregnancy were 14 and 21 years, median number of children was 3, and breastfeeding duration was 62 months per child. Distribution of molecular subtypes for luminal A, luminal B, HER2-enriched, and triple negative (TN) breast cancers was 34.8%, 35.8%, 10.7%, and 18.6%, respectively. After adjusting for covariates, compared to patients with ER-positive tumors, ER-negative patients were more likely to have higher parity (OR = 2.03, 95% CI = (1.11, 3.72), p = 0.021, comparing ≥ 5 to ≤ 2 children). Compared to patients with luminal A tumors, luminal B patients were more likely to have lower parity (OR = 0.45, 95% CI = 0.23, 0.87, p = 0.018, comparing ≥ 5 to ≤ 2 children); HER2-enriched patients were less likely to be obese (OR = 0.36, 95% CI = 0.16, 0.81, p = 0.013) or older age at menopause (OR = 0.38, 95% CI = 0.15, 0.997, p = 0.049). Body mass index (BMI), either overall or by menopausal status, did not vary significantly by ER status. Overall, cumulative or average breastfeeding duration did not vary significantly across subtypes.

CONCLUSIONS

In Kenya, we found associations between parity-related risk factors and ER status consistent with observations in European ancestry populations, but differing associations with BMI and breastfeeding. Inclusion of diverse populations in cancer etiology studies is needed to develop population and subtype-specific risk prediction/prevention strategies.

Rare Germline Variants in Chordoma-Related Genes and Chordoma Susceptibility

Simple Summary

Chordoma is an extremely rare bone cancer that has not been fully characterized and few risk factors have been identified, highlighting the need for improving our understanding of the disease biology. Our study aims to identify chordoma susceptibility genes by investigating 265 genes involved in chordoma-related signaling pathways and other biological processes on germline DNA of 138 chordoma patients of European ancestry compared to internal control datasets and general population databases. Results were intersected with whole genome sequencing data from 80 skull-base chordoma patients of Chinese ancestry. Several rare loss-of-function and predicted deleterious missense variants were enriched in chordoma cases in both datasets, suggesting a complex model of pathways potentially involved in chordoma development and susceptibility, warranting further investigation in larger studies.

Abstract

Background: Chordoma is a rare bone cancer with an unknown etiology. TBXT is the only chordoma susceptibility gene identified to date; germline single nucleotide variants and copy number variants in TBXT have been associated with chordoma susceptibility in familial and sporadic chordoma. However, the genetic susceptibility of chordoma remains largely unknown. In this study, we investigated rare germline genetic variants in genes involved in TBXT/chordoma-related signaling pathways and other biological processes in chordoma patients from North America and China. Methods: We identified variants that were very rare in general population and internal control datasets and showed evidence for pathogenicity in 265 genes in a whole exome sequencing (WES) dataset of 138 chordoma patients of European ancestry and in a whole genome sequencing (WGS) dataset of 80 Chinese patients with skull base chordoma. Results: Rare and likely pathogenic variants were identified in 32 of 138 European ancestry patients (23%), including genes that are part of notochord development, PI3K/AKT/mTOR, Sonic Hedgehog, SWI/SNF complex and mesoderm development pathways. Rare pathogenic variants in COL2A1, EXT1, PDK1, LRP2, TBXT and TSC2, among others, were also observed in Chinese patients. Conclusion: We identified several rare loss-of-function and predicted deleterious missense variants in germline DNA from patients with chordoma, which may influence chordoma predisposition and reflect a complex susceptibility, warranting further investigation in large studies.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4

The Impact of Longitudinal Surveillance on Tumor Thickness for Melanoma-Prone Families with and without Pathogenic Germline Variants of CDKN2A and CDK4

BACKGROUND

Skin cancer screening is routinely performed for members of melanoma-prone families, but longitudinal studies evaluating the efficacy of surveillance in this high-risk population are lacking.

METHODS

We evaluated thickness for first primary melanomas diagnosed in melanoma-prone families (≥2 individuals with melanoma) enrolled in NCT00040352 (NCI familial melanoma study) from 1976 through 2014; enrolled patients received routine skin cancer screening and education about skin self-exams. We used linear and ordinal logistic regression models adjusted for gender and age with a generalized estimating equations approach to report changes in thickness and tumor (T) stage over time, comparing outcomes for NCI cases diagnosed before (pre-study) versus after study participation (prospective) and for NCI cases versus nonfamilial cases [Surveillance, Epidemiology, and End Results (SEER) 9 registries].

RESULTS

Tumor thickness was evaluated for 293 NCI (pre-study = 246; prospective = 47) patients. Compared with NCI pre-study cases, NCI prospective melanomas were thinner (0.6 vs. 1.1 mm; P < 0.001) and more likely to be T1 stage [39/47 (83%) vs. 98/246 (40%); P < 0.001]. Similar findings (P < 0.05) were observed for familial cases with and without germline CDKN2A and CDK4 mutations. Peters-Belson modeling suggested that calendar period effects of decreasing thickness in the general population (SEER 9) did not fully explain thickness trends in NCI families.

CONCLUSIONS

Participation in a longitudinal surveillance program providing skin cancer screening and education about skin self-exams was associated with thinner melanomas for members of melanoma-prone families.

IMPACT

The study findings support the clinical benefit of screening (physician and self) for this high-risk population.

Breast Cancer Risk Factors and Survival by Tumor Subtype: Pooled Analyses from the Breast Cancer Association Consortium

BACKGROUND

It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype.

METHODS

We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer-specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype.

RESULTS

There was no evidence of heterogeneous associations between risk factors and mortality by subtype (P adj > 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5-25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06-1.34)]; current versus never smoking [1.37 (1.27-1.47)], high versus low physical activity [0.43 (0.21-0.86)], age ≥30 years versus <20 years at first pregnancy [0.79 (0.72-0.86)]; >0-<5 years versus ≥10 years since last full-term birth [1.31 (1.11-1.55)]; ever versus never use of oral contraceptives [0.91 (0.87-0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen-progestin therapy [0.61 (0.54-0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02-1.21) for current versus never smoking.

CONCLUSIONS

We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype.

IMPACT

Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Combined Associations of a Polygenic Risk Score and Classical Risk Factors With Breast Cancer Risk

We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

[Pulmonary mucosa-associated lymphoid tissue lymphoma concurrent with lung squamous cell carcinoma: a case report and literature review]

Objective: To explore the clinical manifestations, imaging features, pathological features, diagnosis and treatment of pulmonary mucosal-associated lymphoid tissue(MALT)lymphoma concurrent with lung squamous cell carcinoma, and to improve the understanding of this disease. Methods: Using "Pulmonary mucosa-associated lymphoid tissue, squamous cell carcinoma" as the search term, from January 1, 1983 to August 31, 2020, a total of 3 cases were retrieved in the PubMed database. In the Wanfang database, using "Lung mucosa-associated lymphoid tissue, lung squamous cell carcinoma" as the search term, from January 1, 1990 to August 31, 2020, a total of 1 related document was retrieved. In the CNKI database, "(lung) mucosa-associated lymphoid tissue lymphoma, (lung) squamous cell carcinoma" was used as the search term, and no relevant case reports were retrieved. Results: A 64-year-old man was admitted to the hospital because of chest tightness and shortness of breath for 10 days, cough and fever for one day. Enhanced CT of the chest showed a soft tissue mass shadow in the right lower hilar area, with obstruction of the adjacent bronchus, and local mild enhancement, suggesting of right lower lung cancer. In addition, the CT scan also showed consolidated shadows in the lower lobes of both lungs, scattered nodules, multiple lymphadenopathy in the mediastinum, and a small amount of pleural effusion on the right. Under bronchoscopy, a cauliflower-like neoplasm was seen at the opening of the lower right basal section, about 7 mm×8 mm, and biopsy showed that part of the mucosal structure was destroyed, with disappearance of the squamous epithelial layer, and the nuclei were large and deeply stained, and some were distributed in nests, with poor keratinization and a small amount of necrosis, and fibrous tissue reaction. Immunostaining revealed that the tumor was positive for p40, CK5/6 and EGFR and negative forTTF-1, NapsinA, PD-L1, p53, with about 30% Ki-67 positive cells. A puncture biopsy of the right lower lobe showed that the alveolar cavity was filled with nested lymphoid cells, consisting of small lymphocytes, central cell-like cells and monocyte-like cells, with occasionally large cells. Immunostaining revealed CD₂₀⁺, CD_79a⁺, scattered CD₃+, Bcl2+, SMA vascular+, Bcl6-, CK-, CD10-, CyclinD1-, with about 3% Ki-67 positive cells. The histopathological examinations confirmed the diagnosis of mucosal-associated lymphoid tissue extranodal marginal zone lymphoma(MALT lymphoma),and lung squamous cell carcinoma. Conclusions: Pulmonary mucosa-associated lymphoid tissue lymphoma complicated with lung squamous cell carcinoma is rare and easy to be missed and misdiagnosed. Chest CT imaging shows single or multiple nodules, mass shadows or consolidation, often accompanied by air-bronchial signs in the lesion, bronchiectasis, ground glass density around the lesion, hilar and mediastinal lymphadenopathy. Occasionally, pleural effusion can be seen. Lung biopsy is the gold standard for diagnosis.

Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women

BACKGROUND

Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking.

METHODS

We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity.

RESULTS

Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants.

CONCLUSIONS

The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

Whole genome sequencing of skull-base chordoma reveals genomic alterations associated with recurrence and chordoma-specific survival

Chordoma is a rare bone tumor with an unknown etiology and high recurrence rate. Here we conduct whole genome sequencing of 80 skull-base chordomas and identify PBRM1, a SWI/SNF (SWItch/Sucrose Non-Fermentable) complex subunit gene, as a significantly mutated driver gene. Genomic alterations in PBRM1 (12.5%) and homozygous deletions of the CDKN2A/2B locus are the most prevalent events. The combination of PBRM1 alterations and the chromosome 22q deletion, which involves another SWI/SNF gene (SMARCB1), shows strong associations with poor chordoma-specific survival (Hazard ratio [HR] = 10.55, 95% confidence interval [CI] = 2.81-39.64, p = 0.001) and recurrence-free survival (HR = 4.30, 95% CI = 2.34-7.91, p = 2.77 × 10-6). Despite the low mutation rate, extensive somatic copy number alterations frequently occur, most of which are clonal and showed highly concordant profiles between paired primary and recurrence/metastasis samples, indicating their importance in chordoma initiation. In this work, our findings provide important biological and clinical insights into skull-base chordoma.

[Study on the value of serum des-γ-carboxy prothrombin in predicting hepatocellular carcinoma recurrence after liver transplantation]

Objective: To explore the clinical value of serum des-γ-carboxy prothrombin (DCP) in predicting hepatocellular carcinoma recurrence after liver transplantation. Methods: A total of 115 cases with hepatocellular carcinoma who underwent liver transplantation in Zhongshan Hospital Affiliated to Fudan University from October 2016 to December 2018 were retrospectively analyzed. Receiver operating characteristic curve analysis, Mann-Whitney U test, Kaplan-Meier method, Log-Rank test, χ2 test, univariate and multivariate Cox regression analysis and other statistical methods were used to explore the value of DCP in predicting tumor recurrence after liver transplantation and its correlation with clinicopathological characteristics. Results: The preoperative serum DCP level in recurrent population after liver transplantation was significantly higher than that in non-recurrent population (P < 0.001). The optimal cut-off value of preoperative DCP for predicting recurrence was 200mAU/ml with the use of receiver operating characteristic curve. The sensitivity, specificity, Youden's index and the receiver operating characteristic curve was 87.90%, 57.30%, 0.452, and 0.726, respectively. Survival analysis results grouped by this cut-off value showed that patients with preoperative DCP ≥200mAU/ml had a higher probability of recurrence (P < 0.001). Further, subgroup survival analysis showed that patients with preoperative DCP≥200 mAU/ ml had a higher probability of recurrence than other cases of alpha-fetoprotein negative subgroup, cumulative tumor diameter ≤ 9 cm subgroup and Milan criteria subgroup (P < 0.05). Cox regression analysis showed that preoperative DCP≥200 mAU/ ml (P = 0.017) and cumulative tumor diameter > 9 cm (P = 0.014) was an independent risk factor for recurrence after liver transplantation. χ (2) test results showed that preoperative serum DCP level was correlated with gender, serum gamma glutamyltransferase level, serum alpha fetoprotein level, cumulative tumor diameter, vascular invasion, tumor differentiation and liver cancer transplant criteria (P < 0.05). Conclusion: Preoperative serum DCP can be used as a supplement to the existing liver cancer transplant criteria to predict hepatocellular carcinoma recurrence after liver transplantation. In addition, the accurate screening of patients with low risk of HCC recurrence after liver transplantation can improve the prognosis and efficacy of liver transplant patients.

Using whole-exome sequencing and protein interaction networks to prioritize candidate genes for germline cutaneous melanoma susceptibility

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein-protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein-protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

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.

Abstract 146: DNA methylation age of paired tumor-normal breast tissue in Chinese women with breast cancer

Background: Epigenetic age, captured by DNA methylation changes, is thought to be more informative with respect to disease risk and progression, than chronological age. A few studies have associated epigenetic age acceleration (EAA), the difference between DNA methylation age (DNAm age) and chronological age, with BC risk and subtypes. However, most of these studies were conducted in Western populations. In this study, we examined EAA in an Asian population, in which BC incidence rate is lower and age at BC onset is younger compared to most Western populations.

Methods: We performed genome-wide DNA methylation profiling of 97 tumor and 89 paired distant normal tissue samples collected from BC patients in Hong Kong (HK) using an Illumina MethylationEPIC array. Two independent datasets were used to compare results: The Cancer Genome Atlas (TCGA, n=525 tumor and n=88 adjacent normal) and healthy women from the Komen tissue bank (n=59). DNAm age was calculated using Horvath's model based on 353 CpGs. The significance of EAA was tested using a simple linear regression model. We also used a multivariate regression model to test for equality of slopes (EAA rates) among different BC subtypes or datasets, considering interaction terms between age and subtype/dataset.

Results: The average age at BC diagnosis was 58 years old (range: 33-81) and the distribution of the molecular subtype based on PAM50 was 43.0%, 29.1%, and 27.9% for luminal A, luminal B, and HER2-enriched/basal tumors, respectively. As expected, DNAm age showed a stronger correlation with chronological age in normal tissue (r=0.78, p&lt;.0001) than in tumor tissue (r=0.27, p=0.0075). The average EAAs in normal and tumor samples were 9.43 and 6.02 years, respectively. Among different BC subtypes, EAA in normal tissue did not vary by subtypes. However, in tumor samples, luminal patients showed positive EAA (average 10/13 years in luminal A/B, respectively), while HER2-enriched/basal patients showed a negative EAA (average -8 years), although the rate of EAA did not vary significantly by subtype. Analyses on TCGA data produced consistent results. When comparing the rate of DNAm acceleration in normal tissue of HK, TCGA, and Komen, HK patients had a significantly different rate (βHK=0.46) compared to TCGA (βTCGA=0.65, p=0.001) and Komen (βKomen=0.80, p&lt;.0001).

Conclusion: Consistent with previous studies, we found that EAA in tumor samples varied across tumor subtypes. We also found that HK BC patients' epigenetic age accelerated at a different rate compared to predominantly white TCGA and Komen women, suggesting a potential racial biological difference. Large studies in other Asian populations are warranted to confirm our findings, which may provide biological insight into racial heterogeneity of BC, especially with regard to age at onset.

Citation Format: Hela Koka, Bin Zhu, Shelly Lap Ah Tse, Difei Wang, Maeve Kiely, Jennifer Lyn Guida, Priscilla Lee, Feng Wang, Cherry Wu, Koon Ho Tsang, Wing-cheong Chan, Sze Hong Law, Eric Karlins, Bin Zhu, Amy Hutchinson, Belynda Hicks, Xiaohong R. Yang. DNA methylation age of paired tumor-normal breast tissue in Chinese women with breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 146.

Divergent breast cancer incidence trends by hormone receptor status in the state of Sarawak, Malaysia

Recent studies from high-risk countries such as the US, Denmark and Ireland have shown rising incidence rates of hormone receptor (HR)-positive and falling rates of HR-negative breast cancers (BC). However, it remains unclear whether a similar pattern occurs in low-risk countries. Detailed clinical and risk factor data were collected from 2,977 female invasive BC patients (≥20 years) in Sarawak General Hospital, Malaysia, representing 93% of the population. The population-at-risk was obtained from the Department of Statistics Malaysia. Secular trends in age-standardized incidence rates were assessed using estimated average annual percent changes. Associations between established BC risk factors and tumor subtypes defined by HR or joint human epidermal growth factor receptor 2 (HR/HER2) status were examined by case-case comparisons using logistic regression. From 2006 to 2015, incidence rates increased for HR-positive cancers by 4.46%/year (95% CI = 2.19-6.78) and decreased for HR-negative cancers by 2.29%/year (95% CI = -4.31 to -0.24). When further stratified by HER2, the most contrasting difference in linear trends was observed between HR+/HER2- and HR-/HER2- subtypes. After controlling for potential confounders, cases with excess body weight (OR_{overweight vs. normal} = 0.82; 95% CI = 0.69-0.98; OR_{obese vs. normal} = 0.62; 95% CI = 0.48-0.80), later age at first birth (OR_{≥26 years vs. <23 years} = 0.82; 95% CI = 0.66-1.02), nulliparity (OR_{nulliparous vs. <23 years} = 0.74; 95% CI = 0.59-0.94) and never-breastfeeding (OR_{never vs. ever} = 0.73; 95% CI = 0.55-0.97) were less frequent among HR-negative cases than among HR-positive cases. Diverging incidence trends by HR expression were similar in Sarawak and Western countries, possibly reflecting changes in the prevalence of risk factors with opposing effects by tumor subtypes in low- and high-risk populations.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.