Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Germline pathogenic variants associated with triple-negative breast cancer in US Hispanic and Guatemalan women using hospital and community-based recruitment strategies

PURPOSE

Recruit and sequence breast cancer subjects in Guatemalan and US Hispanic populations. Identify optimum strategies to recruit Latin American and Hispanic women into genetic studies of breast cancer.

METHODS

We used targeted gene sequencing to identify pathogenic variants in 19 familial breast cancer susceptibility genes in DNA from unselected Hispanic breast cancer cases in the US and Guatemala. Recruitment across the US was achieved through community-based strategies. In addition, we obtained patients receiving cancer treatment at major hospitals in Texas and Guatemala.

RESULTS

We recruited 287 Hispanic US women, 38 (13%) from community-based and 249 (87%) from hospital-based strategies. In addition, we ascertained 801 Guatemalan women using hospital-based recruitment. In our experience, a hospital-based approach was more efficient than community-based recruitment. In this study, we sequenced 103 US and 137 Guatemalan women and found 11 and 10 pathogenic variants, respectively. The most frequently mutated genes were BRCA1, BRCA2, CHEK2, and ATM. In addition, an analysis of 287 US Hispanic patients with pathology reports showed a significantly higher percentage of triple-negative disease in patients with pathogenic variants (41% vs. 15%). Finally, an analysis of mammography usage in 801 Guatemalan patients found reduced screening in women with a lower socioeconomic status (p < 0.001).

CONCLUSION

Guatemalan and US Hispanic women have rates of hereditary breast cancer pathogenic variants similar to other populations and are more likely to have early age at diagnosis, a family history, and a more aggressive disease. Patient recruitment was higher using hospital-based versus community enrollment. This data supports genetic testing in breast cancer patients to reduce breast cancer mortality in Hispanic women.

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

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

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

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

Five-Year Follow-Up of a Slowly Resorbable Biosynthetic P4HB Mesh (Phasix) in VHWG Grade 3 Incisional Hernia Repair

Objective

To assess the 5-year recurrence rate of incisional hernia repair in Ventral Hernia Working Group (VHWG) 3 hernia with a slowly resorbable mesh.

Summary Background Data

Incisional hernia recurs frequently after initial repair. In potentially contaminated hernia, recurrences rise to 40%. Recently, the biosynthetic Phasix mesh has been developed that is resorbed in 12-18 months. Resorbable meshes might be a solution for incisional hernia repair to decrease short- and long-term (mesh) complications. However, long-term outcomes after resorption are scarce.

Methods

Patients with VHWG grade 3 incisional midline hernia, who participated in the Phasix trial (Clinilcaltrials.gov: NCT02720042) were included by means of physical examination and computed tomography (CT). Primary outcome was hernia recurrence; secondary outcomes comprised of long-term mesh complications, reoperations, and abdominal wall pain [visual analogue score (VAS): 0-10].

Results

In total, 61/84 (72.6%) patients were seen. Median follow-up time was 60.0 [interquartile range (IQR): 55-64] months. CT scan was made in 39 patients (68.4%). A recurrence rate of 15.9% (95% confidence interval: 6.9-24.8) was calculated after 5 years. Four new recurrences (6.6%) were found between 2 and 5 years. Two were asymptomatic. In total, 13/84 recurrences were found. No long-term mesh complications and/or interventions occurred. VAS scores were 0 (IQR: 0-2).

Conclusions

Hernia repair with Phasix mesh in high-risk patients (VHWG 3, body mass index >28) demonstrated a recurrence rate of 15.9%, low pain scores, no mesh-related complications or reoperations for chronic pain between the 2- and 5-year follow-up. Four new recurrences occurred, 2 were asymptomatic. The poly-4-hydroxybutyrate mesh is a safe mesh for hernia repair in VHWG 3 patients, which avoids long-term mesh complications like pain and mesh infection.

Analysis of Several Common APOBEC-type Mutations in Bladder Tumors Suggests Links to Viral Infection

FGFR3 and PIK3CA are among the most frequently mutated genes in bladder tumors. We hypothesized that recurrent mutations in these genes might be caused by common carcinogenic exposures such as smoking and other factors. We analyzed 2,816 bladder tumors with available data on FGFR3 and/or PIK3CA mutations, focusing on the most recurrent mutations detected in ≥10% of tumors. Compared to tumors with other FGFR3/PIK3CA mutations, FGFR3-Y375C was more common in tumors from smokers than never-smokers (P = 0.009), while several APOBEC-type driver mutations were enriched in never-smokers: FGFR3-S249C (P = 0.013) and PIK3CA-E542K/PIK3CA-E545K (P = 0.009). To explore possible causes of these APOBEC-type mutations, we analyzed RNA sequencing (RNA-seq) data from 798 bladder tumors and detected several viruses, with BK polyomavirus (BKPyV) being the most common. We then performed IHC staining for polyomavirus (PyV) Large T-antigen (LTAg) in an independent set of 211 bladder tumors. Overall, by RNA-seq or IHC-LTAg, we detected PyV in 26 out of 1,010 bladder tumors with significantly higher detection (P = 4.4 × 10-5), 25 of 554 (4.5%) in non-muscle-invasive bladder cancers (NMIBC) versus 1 of 456 (0.2%) of muscle-invasive bladder cancers (MIBC). In the NMIBC subset, the FGFR3/PIK3CA APOBEC-type driver mutations were detected in 94.7% (18/19) of PyV-positive versus 68.3% (259/379) of PyV-negative tumors (P = 0.011). BKPyV tumor positivity in the NMIBC subset with FGFR3- or PIK3CA-mutated tumors was also associated with a higher risk of progression to MIBC (P = 0.019). In conclusion, our results support smoking and BKPyV infection as risk factors contributing to bladder tumorigenesis in the general patient population through distinct molecular mechanisms.

PREVENTION RELEVANCE

Tobacco smoking likely causes one of the most common mutations in bladder tumors (FGFR3-Y375C), while viral infections might contribute to three others (FGFR3-S249C, PIK3CA-E542K, and PIK3CA-E545K). Understanding the causes of these mutations may lead to new prevention and treatment strategies, such as viral screening and vaccination.

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

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

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

SIGNIFICANCE

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

Nitrated Polycyclic Aromatic Hydrocarbon (Nitro-PAH) Signatures and Somatic Mutations in Diesel Exhaust-Exposed Bladder Tumors

BACKGROUND

Diesel exhaust is a complex mixture, including polycyclic aromatic hydrocarbons (PAH) and nitrated PAHs (nitro-PAH), many of which are potent mutagens and possible bladder carcinogens. To explore the association between diesel exposure and bladder carcinogenesis, we examined the relationship between exposure and somatic mutations and mutational signatures in bladder tumors.

METHODS

Targeted sequencing was conducted in bladder tumors from the New England Bladder Cancer Study. Using data on 797 cases and 1,418 controls, two-stage polytomous logistic regression was used to evaluate etiologic heterogeneity between bladder cancer subtypes and quantitative, lifetime estimates of respirable elemental carbon (REC), a surrogate for diesel exposure. Poisson regression was used to evaluate associations between REC and mutational signatures.

RESULTS

We observed significant heterogeneity in the diesel-bladder cancer risk relationship, with a strong positive association among cases with high-grade, nonmuscle invasive TP53-mutated tumors compared with controls [ORTop Tertile vs.Unexposed, 4.8; 95% confidence interval (CI), 2.2-10.5; Ptrend < 0.001; Pheterogeneity = 0.002]. In muscle-invasive tumors, we observed a positive association between diesel exposure and the nitro-PAH signatures of 1,6-dintropyrene (RR, 1.93; 95% CI, 1.28-2.92) and 3-nitrobenzoic acid (RR, 1.97; 95% CI, 1.33-2.92).

CONCLUSIONS

The relationship between diesel exhaust and bladder cancer was heterogeneous based on the presence of TP53 mutations in tumors, further supporting the link between PAH exposure and TP53 mutations in carcinogenesis. Future studies that can identify nitro-PAH signatures in exposed tumors are warranted to add human data supporting the link between diesel and bladder cancer.

IMPACT

This study provides additional insight into the etiology and possible mechanisms related to diesel exhaust-induced bladder cancer.

Abstract 1938: Germline mutations in breast cancer genes are associated with early age of diagnosis and triple negative disease in Guatemalan and US Hispanic women

Abstract Background: Mutations in hereditary breast cancer genes play an essential role in cancer risk. Little is known of the type and frequency of mutations in Hispanic populations in the United States and Central American countries, including Guatemala.

Methods: We used exome sequencing to identify mutations in blood DNA from unselected Hispanic breast cancer cases from community recruitment and from two hospitals each in Texas and Guatemala. Data from a structured questionnaire was used to compare mutation carriers of medium and high penetrance genes. Variants were annotated with ClinVar and VarSome.

Results: We recruited 262 Hispanic US women, 37 (14%) from community-based recruitment and 225 (86%) from hospitals in Texas. In addition, we ascertained 633 patients from two hospitals in Guatemala City. A total of 91 out of 895 subjects (10%) had a variant classified as pathogenic in a gene with known high or medium penetrance for inherited breast cancer. The most frequently mutated genes were the high penetrance BRCA1 (44/895, 4.9%) followed by BRCA2 (23/895, 2.6%), PALB2 (5/895, 0.6%), CHEK2 (5/895 0.6%), ATM (6/895, 0.7%) and TP53 (6/895, 0.7%). Pathogenic variants were also detected in the moderate penetrance genes BARD1 and MSH6, and rare pathogenic variants detected in the low penetrance genes AXIN2, FH, MLH1, MSH2, MUTYH, NF1, and SDHB. The high ratio of BRCA1/BRCA2 mutations is due to the presence of two potential founder mutations, BRCA1 c.212+1G&gt;A splice mutation (18 cases) and BRCA1 c.799delT (9 cases) in Guatemala. Compared to all others, cases with pathogenic mutations had a significantly earlier age at diagnosis (45 vs. 52 years, P&lt;0.001) and were more likely to have a diagnosis before menopause. A higher percentage of mutation carriers had a relative with any cancer (51% vs. 37%, P=0.038) or breast cancer (33% vs. 15%, P&lt;0.001). Patients with pathogenic mutations had a significantly higher percentage of triple-negative disease (60% vs. 13%, P=0.00046). Finally, mammography usage was less frequent in women with a lower socioeconomic status, indicating that this group is less likely to be screened for breast cancer (P&lt;0.001).

Conclusions: Guatemalan women have rates of hereditary breast cancer mutations like other populations, and these women are more likely to have early age at diagnosis, triple-negative disease, and

family history. Patient recruitment was higher using hospital-based versus community enrollment. This data supports genetic testing in breast cancer patients and those at high risk as part of a strategy to reduce breast cancer mortality in Hispanic women.

Citation Format: Jesica Godinez Paredes, Isabel Rodriguez, Megan Ren, Anali Orozco, Jeremy Ortiz, Anaseidy Albanez, Catherine Jones, Zeina Nahleh, Lilian Barreda, Lisa Garland, Dongjing Wu, Jiahui Wang, Victor Argueta, Roberto Orozco, Eduardo Gharzouzi, Michael Dean. Germline mutations in breast cancer genes are associated with early age of diagnosis and triple negative disease in Guatemalan and US Hispanic women [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 1938.

Abstract 4289: Long-read, assembly-based characterization of rearranged cancer karyotypes

Introduction: Recent pan-cancer whole-genome sequencing studies revealed the rich landscape of structural variants (SV), from simpler indels to complex events involving multiple breakpoints and sequence gain/loss. SVs may contribute to tumorigenesis through direct modification of coding sequence or deregulation from copy number alterations, enhancer hijacking, or topological domain modification.

A substantial part of the variation in the human genome is not accessible to short reads due to mapping ambiguities. Recent benchmarking studies reported that the best short-read methods only have 30-70% SV sensitivity. Long-read sequencing (such as PacBio or Oxford Nanopore) can overcome the limitations of short reads, however the current methods were not designed for the analysis of rearranged cancer genomes with complex copy number profiles.

Methods: We developed BGA (Breakpoint Graph Assembler), a method that combines the ideas from long-read assembly and breakpoint graph frameworks. BGA detects abnormally mapped reads and builds a breakpoint graph that characterizes the structure of derived cancer karyotypes. Complex events with multiple breakpoints form connectivity clusters and are classified based on the subgraph properties. BGA also takes advantage of phased haplotypes and can incorporate multiple related datasets (such as in tumor-normal comparison or multi-site tumor sampling). BGA is freely available at: https://github.com/KolmogorovLab/BGA.

Results: We first analyzed three cancer cell lines and corresponding matching normal DNAs (HCC1954, H2009, and COLO829). In each cell line, we identified 8-56 somatic rearrangement clusters involving more than two breakpoints and at least 1kb of sequence. In H2009, we identified a chromoplexy event involving chr13 and chr1, consistent with previous FISH experiments. COLO829 showed the lowest number of somatic rearrangement clusters (n=8), including translocation and inversion events between chr3, chr10, and chr12 within the RARB, BICC1, and TRHDE genes. Homologous recombination deficient HCC1954 has the highest number of complex events including the chr17q arm which hosts ERBB2, chromoplexy between chr8, chr5, and chromothripsis in chr21.

In addition, we analyzed three other HPV-infected cell lines (CaSki, SCC152, SNU1000). In each of them, we observed complex clusters of HPV-HPV and HPV-human breakpoints that formed cycles, suggesting extrachromosomal amplification. The HPV fragments had many interactions with chromosomal DNA in CaSki and SC152, but not in mostly episomal SNU1000 cells. We also observed karyotype-scale changes that did not involve HPV sequences, such as the simultaneous exchange of six chromosome arms of chr2, chr7, and chr17 in CaSki.

Citation Format: Ayse Keskus, Tanveer Ahmad, Ataberk Donmez, Yi Xie, Isabel Rodriguez, Rose Milano, Nicole Rossi, Hong Lou, Laksh Malik, Kimberley Billingsley, Cornelis Blauwendraat, Michael Dean, Mikhail Kolmogorov. Long-read, assembly-based characterization of rearranged cancer karyotypes. [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 4289.

Evidence for virus-mediated oncogenesis in bladder cancers arising in solid organ transplant recipients

A small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have an overall poorer outcomes, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of the tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyV-positive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyV-positive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomology-mediated end joining and coincided with focal amplifications of the tumor genome similar to other virus-mediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases, with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with an antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virus-associated.

OC-053 LONG-TERM FOLLOW-UP OF A SLOWLY RESORBABLE BIOSYNTHETIC MESH IN VHWG GRADE 3 HERNIA REPAIR

Background

Incisional hernia recurrence occurs frequently after initial repair. In potentially contaminated hernia, these numbers can rise up to 50%. Recently, resorbable meshes have been used to prevent infection in complicated incisional hernia and prevent long term mesh complications. Present study focuses on the long term outcomes of a resorbable mesh.

Methods

Patients included in the priorly conducted Phasix trial with Ventral Hernia Working Group (VHWG) grade 3 hernia were invited for an abdominal CT scan and physical examination. Primary outcome was hernia recurrence, secondary outcomes comprised of abdominal wall configuration, mesh resorption and long term mesh complications.

Results

Of the 84 patients included in the 2-year follow-up, 56 were available for long term follow-up up to five years. Mean follow-up time was 4.2 years. CT scan was made in 36 (64.3%) of patients. Kaplan-Meier analysis showed a recurrence rate of 16.7% after five years. Four (7.1%) new recurrences were found between 2-year and 4.2 years follow-up of which three by CT scan. In 50% of CT scans, no signs of postoperative distortions of the abdominal wall were observed. No long term mesh complications were reported.

Conclusion

After long term follow-up, VHWG 3 hernia repair with biosynthetic mesh demonstrated good performance with regard to hernia recurrence. Possibly, tissue reconstruction took place in the first two years after implantation by the use it or lose it principle, which led to sufficient native abdominal wall strength to prevent recurrences after resorption of the mesh.

The human ATP-binding cassette (ABC) transporter superfamily

The ATP-binding cassette (ABC) transporter superfamily comprises membrane proteins that efflux various substrates across extra- and intracellular membranes. Mutations in ABC genes cause 21 human disorders or phenotypes with Mendelian inheritance, including cystic fibrosis, adrenoleukodystrophy, retinal degeneration, cholesterol, and bile transport defects. To provide tools to study the function of human ABC transporters we compiled data from multiple genomics databases. We analyzed ABC gene conservation within human populations and across vertebrates and surveyed phenotypes of ABC gene mutations in mice. Most mouse ABC gene disruption mutations have a phenotype that mimics human disease, indicating they are applicable models. Interestingly, several ABCA family genes, whose human function is unknown, have cholesterol level phenotypes in the mouse. Genome-wide association studies confirm and extend ABC traits and suggest several new functions to investigate. Whole-exome sequencing of tumors from diverse cancer types demonstrates that mutations in ABC genes are not common in cancer, but specific genes are overexpressed in select tumor types. Finally, an analysis of the frequency of loss-of-function mutations demonstrates that many human ABC genes are essential with a low level of variants, while others have a higher level of genetic diversity.

Aflatoxin levels and prevalence of TP53 aflatoxin-mutations in hepatocellular carcinomas in Mexico

Objective. To determine the exposure to aflatoxin B1 (AFB1) in southern Mexico and the presence of the aflatoxin signature mutation in hepatocellular carcinoma (HCC) tissue from patients from a cancer referral center. Materials and methods. We estimated the prevalence and distribution of AFB1 in a representative sample of 100 women and men from Chiapas using the National Health and Nutrition Survey 2018-19. We also examined the presence of the aflatoxin signature mutation in codon 249 (R249S), and other relevant mutations of the TP53 gene in HCC tissue blocks from 24 women and 26 men treated in a national cancer referral center. Results. The prevalence of AFB1 in serum samples was 85.5% (95%CI 72.1-93.1) and the median AFB1 was 0.117 pg/μL (IQR, 0.050–0.350). We detected TP53 R249S in three of the 50 HCCs (6.0%) and observed four other G>T transversions potentially induced by AFB1. Conclusion. Our analysis provides evidence that AFB1 may have a relevant role on HCC etiology in Mexico.

Associations between aflatoxin B 1 ‐albumin adduct levels with metabolic conditions in Guatemala: A cross‐sectional study

Background and Aims

Metabolic conditions such as obesity, type 2 diabetes, metabolic syndrome, and nonalcoholic fatty liver disease (NAFLD) are highly prevalent in Guatemala and increase the risk for a number of disorders, including hepatocellular carcinoma (HCC). Aflatoxin B₁ (AFB₁) levels are also notably elevated in the population and are known to be associated with HCC risk. Whether AFB₁ also contributes to the high prevalence of the metabolic disorders has not been previously examined. Therefore, the purpose of this study was to assess the association between AFB₁ and the metabolic conditions.

Methods

Four‐hundred twenty‐three individuals were included in the study, in which AFB₁‐albumin adduct levels were measured in sera. Metabolic conditions included diabetes, obesity, central obesity, metabolic syndrome, and NAFLD. Crude and adjusted prevalence odds ratios (PORs) and 95% confidence intervals (95% CI) were estimated for the associations between the metabolic conditions and AFB₁‐albumin adduct levels categorized into quartiles.

Results

The study found a significant association between AFB₁‐albumin adduct levels and diabetes (Q4 vs Q1 POR = 3.74, 95%CI: 1.71‐8.19; P‐trend .003). No associations were observed between AFB₁‐albumin adduct levels and the other conditions.

Conclusions

As diabetes is the metabolic condition most consistently linked to HCC, the possible association between AFB₁ exposure and diabetes may be of public health importance. Further studies are warranted to replicate the findings and examine potential mechanisms.

HLA-A*03 and response to immune checkpoint blockade in cancer: an epidemiological biomarker study

BACKGROUND

Predictive biomarkers could allow more precise use of immune checkpoint inhibitors (ICIs) in treating advanced cancers. Given the central role of HLA molecules in immunity, variation at the HLA loci could differentially affect the response to ICIs. The aim of this epidemiological study was to determine the effect of HLA-A*03 as a biomarker for predicting response to immunotherapy.

METHODS

In this epidemiological study, we investigated the clinical outcomes (overall survival, progression free survival, and objective response rate) after treatment for advanced cancer in eight cohorts of patients: three observational cohorts of patients with various types of advanced tumours (the Memorial Sloan Kettering Integrated Mutation Profiling of Actionable Cancer Targets [MSK-IMPACT] cohort, the Dana-Farber Cancer Institute [DFCI] Profile cohort, and The Cancer Genome Atlas) and five clinical trials of patients with advanced bladder cancer (JAVELIN Solid Tumour) or renal cell carcinoma (CheckMate-009, CheckMate-010, CheckMate-025, and JAVELIN Renal 101). In total, these cohorts included 3335 patients treated with various ICI agents (anti-PD-1, anti-PD-L1, and anti-CTLA-4 inhibitors) and 10 917 patients treated with non-ICI cancer-directed therapeutic approaches. We initially modelled the association of HLA amino-acid variation with overall survival in the MSK-IMPACT discovery cohort, followed by a detailed analysis of the association between HLA-A*03 and clinical outcomes in MSK-IMPACT, with replication in the additional cohorts (two further observational cohorts and five clinical trials).

FINDINGS

HLA-A*03 was associated in an additive manner with reduced overall survival after ICI treatment in the MSK-IMPACT cohort (HR 1·48 per HLA-A*03 allele [95% CI 1·20-1·82], p=0·00022), the validation DFCI Profile cohort (HR 1·22 per HLA-A*03 allele, 1·05-1·42; p=0·0097), and in the JAVELIN Solid Tumour clinical trial for bladder cancer (HR 1·36 per HLA-A*03 allele, 1·01-1·85; p=0·047). The HLA-A*03 effect was observed across ICI agents and tumour types, but not in patients treated with alternative therapies. Patients with HLA-A*03 had shorter progression-free survival in the pooled patient population from the three CheckMate clinical trials of nivolumab for renal cell carcinoma (HR 1·31, 1·01-1·71; p=0·044), but not in those receiving control (everolimus) therapies. Objective responses were observed in none of eight HLA-A*03 homozygotes in the ICI group (compared with 59 [26·6%] of 222 HLA-A*03 non-carriers and 13 (17·1%) of 76 HLA-A*03 heterozygotes). HLA-A*03 was associated with shorter progression-free survival in patients receiving ICI in the JAVELIN Renal 101 randomised clinical trial for renal cell carcinoma (avelumab plus axitinib; HR 1·59 per HLA-A*03 allele, 1·16-2·16; p=0·0036), but not in those receiving control (sunitinib) therapy. Objective responses were recorded in one (12·5%) of eight HLA-A*03 homozygotes in the ICI group (compared with 162 [63·8%] of 254 HLA-A*03 non-carriers and 40 [55·6%] of 72 HLA-A*03 heterozygotes). HLA-A*03 was associated with impaired outcome in meta-analysis of all 3335 patients treated with ICI at genome-wide significance (p=2·01 × 10-8) with no evidence of heterogeneity in effect (I2 0%, 95% CI 0-0·76) INTERPRETATION: HLA-A*03 is a predictive biomarker of poor response to ICI. Further evaluation of HLA-A*03 is warranted in randomised trials. HLA-A*03 carriage could be considered in decisions to initiate ICI in patients with cancer.

FUNDING

National Institutes of Health, Merck KGaA, and Pfizer.

Sequence Variant in the TRIM39-RPP21 Gene Readthrough is Shared Across a Cohort of Arabian Foals Diagnosed with Juvenile Idiopathic Epilepsy

Juvenile idiopathic epilepsy (JIE) is a self-limiting neurological disorder with a suspected genetic predisposition affecting young Arabian foals of the Egyptian lineage. The condition is characterized by tonic-clonic seizures with intermittent post-ictal blindness, in which most incidents are sporadic and unrecognized. This study aimed to identify genetic components shared across a local cohort of Arabian foals diagnosed with JIE via a combined whole genome and targeted resequencing approach: Initial whole genome comparisons between a small cohort of nine diagnosed foals (cases) and 27 controls from other horse breeds identified variants uniquely shared amongst the case cohort. Further validation via targeted resequencing of these variants, that pertain to non-intergenic regions, on additional eleven case individuals revealed a single 19bp deletion coupled with a triple-C insertion (Δ19InsCCC) within the TRIM39-RPP21 gene readthrough that was uniquely shared across all case individuals, and absent from three additional Arabian controls. Furthermore, we have confirmed recent findings refuting potential linkage between JIE and other inherited diseases in the Arabian lineage, and refuted the potential linkage between JIE and genes predisposing a similar disorder in human newborns. This is the first study to report a genetic variant to be shared in a sub-population cohort of Arabian foals diagnosed with JIE. Further evaluation of the sensitivity and specificity of the Δ19InsCCC allele within additional cohorts of the Arabian horse is warranted in order to validate its credibility as a marker for JIE, and to ascertain whether it has been introduced into other horse breeds by Arabian ancestry.

Antimicrobial use in agriculture: critical review of the factors influencing behaviour

Antimicrobial resistance (AMR) is a global health emergency affecting humans and animals, diminishing the effectiveness of medication used to treat illness. The agri-food sector has attracted increased attention for imprudent antimicrobial use (AMU) and its contribution to AMR. Thus, ascertaining farmers' and veterinarians' behaviours surrounding AMU is essential to address imprudent AMU and generate behaviour change within the agri-food sector. Therefore, the aim of this critical review is to investigate, assess and collate the current body of evidence to identify psychosocial factors including knowledge, understanding, perceptions, attitudes and behaviours surrounding AMU. Database searches were limited to articles utilizing qualitative and quantitative methodologies, available in English with no restriction on publication year. Of the 1156 articles identified, 103 were retained for this review. Findings on the psychosocial aspects were thematically analysed. Five key themes emerged from the data: (i) knowledge and awareness of antimicrobials; (ii) attitudes towards antimicrobials; (iii) influential relationships; (iv) resources; and (v) factors influencing AMU. Results indicated that to overcome barriers experienced by key stakeholders, a carefully considered, evidence-based approach, incorporating behaviour change theory, is required when designing intricate interventions/strategies, in order to elicit successful and sustained AMU behaviour change.

Cervical Visual Inspection with Acetic Acid (VIA) and Oncogenic Human Papillomavirus Screening in Rural Indigenous Guatemalan Women: Time to Rethink VIA

Single-visit “screen-and-treat” strategies using visual inspection with acetic acid (VIA) and cryotherapy (liquid nitrous oxide ablation) in low-resource settings are commonly used to detect and treat precancerous lesions for cervical cancer prevention. This study compared VIA sensitivity and specificity in rural indigenous Guatemalan communities, to that of oncogenic human papillomavirus (HPV) testing for detection of precancerous changes, using cytology as the reference standard. Between 3–8 September 2017, trained nurses examined 222 women aged 23–58 years with VIA. Specimens for liquid-based cytology and HPV testing were obtained prior to VIA with a cytobrush and transported in PreservCyt to a US clinical laboratory. VIA and HPV test sensitivities were assessed as proportions of women with abnormal cytology that had abnormal VIA or HPV results, respectively, and specificities, as proportions with normal cytology with normal VIA or negative HPV tests. Of 222 women, 18 (8.1%) had abnormal cytology (1 carcinoma in a participant who received VIA-based cryotherapy in 2015, 4 high- and 5 low-grade squamous intraepithelial lesions, and 8 atypical squamous cells of undetermined significance (ASCUS)). Excluding ASCUS, sensitivities of VIA and HPV were 20.0% and 100%, respectively. VIA-based screening may not be acceptable for detecting precancerous lesions, and field cryotherapy for preventing malignancy. The World Health Organization recommended in 2021 “…using HPV DNA detection as the primary screening test rather than VIA or cytology”.

Joint Clustering of Single-Cell Sequencing and Fluorescence In Situ Hybridization Data for Reconstructing Clonal Heterogeneity in Cancers

Aneuploidy and whole genome duplication (WGD) events are common features of cancers associated with poor outcomes, but the ways they influence trajectories of clonal evolution are poorly understood. Phylogenetic methods for reconstructing clonal evolution from genomic data have proven a powerful tool for understanding how clonal evolution occurs in the process of cancer progression, but extant methods so far have limited the ability to resolve tumor evolution via ploidy changes. This limitation exists in part because single-cell DNA-sequencing (scSeq), which has been crucial to developing detailed profiles of clonal evolution, has difficulty in resolving ploidy changes and WGD. Multiplex interphase fluorescence in situ hybridization (miFISH) provides a more unambiguous signal of single-cell ploidy changes but it is limited to profiling small numbers of single markers. Here, we develop a joint clustering method to combine these two data sources with the goal of better resolving ploidy changes in tumor evolution. We develop a probabilistic framework to maximize the probability of latent variables given the pre-clustered datasets, which we optimize via Markov chain Monte Carlo sampling combined with linear regression. We validate the method by using simulated data derived from a glioblastoma (GBM) case profiled by both scSeq and miFISH. We further apply the method to two GBM cases with scSeq and miFISH data by reconstructing a phylogenetic tree from the joint clustering results, demonstrating their synergistic value in understanding how focal copy number changes and WGD events can collectively contribute to tumor progression.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

Outcomes following the delayed management of open tibial fractures

AIMS

National guidelines set standards for the definitive management of open fractures within 72 h. This study aims to investigate our outcomes where this timeline was unachievable for most cases due to a split-site orthoplastic service.

PATIENTS & METHODS

116 consecutive Gustilo-Anderson grade IIIB & IIIC open tibial fractures presenting to our major trauma centre (MTC) between September 2012 and April 2018 were reviewed. The mean follow up was 46 months (17 to 88). 110 (95%) were grade IIIB and 6 (5%) grade IIIC. The most common injury mechanism included road traffic accidents (59%) and falls (28%). Primary outcomes were recorded according to; timing of initial debridement and definitive cover, rates of superficial and deep infection, non-union and amputation. Subgroups were statistically analysed according to time to initial debridement, definitive soft-tissue cover and injury severity score (ISS).

RESULTS

The mean time to initial debridement was 11.3 h (2.9 to 38.9) and definitive soft-tissue cover 9.9 days (0 to 37). We recorded rates of: superficial infection; 42 cases (36%), deep infection; 14 cases (12%) and non-union requiring revision; 19 cases (16%). There were 20 amputations (17%) with 9 (8.6%) performed early and 11 (9.5%) delayed. Subgroup analysis showed higher rates of superficial infection (50%, p = 0.002) and amputation (26.6%, p = 0.01) for those debrided <12 h. A greater presenting ISS related to a delay to definitive cover >7 days (p = 0.05). Primary outcomes trended worse for those covered >7 days but did not reach significance.

CONCLUSION

Major trauma patients are particularly vulnerable to poor outcomes resulting from the delay in definitive management of open fractures. MTC's need resources and a co-located orthoplastic service to achieve national standards and better outcomes. Current guidelines do not advise for the management of patients where a delay in definitive surgery is anticipated.

Tumor heterogeneity assessed by sequencing and fluorescence in situ hybridization (FISH) data

MOTIVATION

Computational reconstruction of clonal evolution in cancers has become a crucial tool for understanding how tumors initiate and progress and how this process varies across patients. The field still struggles, however, with special challenges of applying phylogenetic methods to cancers, such as the prevalence and importance of copy number alteration (CNA) and structural variation (SV) events in tumor evolution, which are difficult to profile accurately by prevailing sequencing methods in such a way that subsequent reconstruction by phylogenetic inference algorithms is accurate.

RESULTS

In the present work, we develop computational methods to combine sequencing with multiplex interphase fluorescence in situ hybridization (miFISH) to exploit the complementary advantages of each technology in inferring accurate models of clonal CNA evolution accounting for both focal changes and aneuploidy at whole-genome scales. By integrating such information in an integer linear programming (ILP) framework, we demonstrate on simulated data that incorporation of FISH data substantially improves accurate inference of focal CNA and ploidy changes in clonal evolution from deconvolving bulk sequence data. Analysis of real glioblastoma data for which FISH, bulk sequence, and single cell sequence are all available confirms the power of FISH to enhance accurate reconstruction of clonal copy number evolution in conjunction with bulk and optionally single-cell sequence data.

AVAILABILITY

Source code is available on Github at https://github.com/CMUSchwartzLab/FISH_deconvolution.

Abstract 48: Characterization of Human Papillomavirus (HPV) Sequence and Integration and Analysis of Immune Response Genes for Application in Guatemalan Cervical Cancer

Purpose: Cervical cancer is a significant cause of mortality in Guatemala due to inefficient screening and late diagnosis. We are developing new technologies to improve the analysis of HPV-positive cervical tumors in Guatemala. The goal is to enhance understanding of HPV infection, tumor progression, and therapeutics development.

Methods: We established a cohort of 700 Guatemalan cervical cancer patients with blood, fresh-frozen tumors, and clinical data. To optimize techniques, we extracted DNA and RNA from 20 cervical cancer cell lines. Long-read whole genome and transcriptome sequencing, and long-range PCR, were performed on cervical cancer DNA. HPV and HLA gene sequences were aligned to reference genomes.

Results: Complete HPV sequences were compiled, HPV subtypes determined, and integration sites identified in cell lines with single and multiple integration events. HLA types of the class I genes were determined from the complete gene sequence, mutations detected, and HLA homozygosity. RNA was used for long-read whole transcriptome analysis of HPV and cellular genes. These results establish the use of long-read sequencing technology to characterize both HLA and host genomes. Paired Guatemalan tumor and blood DNA was successfully used for full HLA gene sequencing and mutation analysis. In addition, the complete HPV genome sequence was resolved. Binding of HPV peptides was predicted using the HPV and HLA sequences. These peptides identify the specific HPV epitopes potentially presented on the cell surface.

Conclusion: Long-read RNA and DNA sequencing technology allow the determination of the complete HPV sequence, integration, and HLA class I typing and mutation analysis. HLA typing provides essential information on an individual's immune response. These techniques and analyses are being applied to tumor samples from women in Guatemala to aid in the development of targeted treatments. These methods are portable and could be employed in low- and middle-income countries.

Citation Format: Nicole Rossi, Michael Dean. Characterization of Human Papillomavirus (HPV) Sequence and Integration and Analysis of Immune Response Genes for Application in Guatemalan Cervical Cancer [abstract]. In: Proceedings of the 9th Annual Symposium on Global Cancer Research; Global Cancer Research and Control: Looking Back and Charting a Path Forward; 2021 Mar 10-11. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2021;30(7 Suppl):Abstract nr 48.

Abstract 861: Diesel engine exhaust is associated with TP53 mutations and high-risk non-muscle invasive bladder cancer

Background: Second to cigarette smoking, occupational exposures are a leading risk factor for the development of bladder cancer. Occupational exposure to diesel exhaust in particular has been linked to increased bladder cancer risk in several epidemiologic studies. To provide insight into possible mechanisms that influence this association, we examined the relationship between diesel exhaust and somatic mutations in bladder tumors.

Methods: Targeted sequencing of frequently mutated genes in bladder cancer was conducted in formalin-fixed paraffin-embedded bladder tumors from a population-based case-control study in New England. Lifetime occupational histories combined with exposure-oriented questions were used to estimate cumulative exposure to respirable elemental carbon (REC), a primary index for diesel exhaust. REC exposure was modeled both continuously (per 100 μg/m3-yrs) and categorically (by tertiles). Two-stage polytomous logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) to evaluate etiologic heterogeneity by tumor subtypes using data on 797 cases and 1,418 controls.

Results: We observed significant heterogeneity in the relationship between diesel exhaust and risk of bladder cancer by TP53 mutation status (p-heterogeneity=0.013). When comparing patients with TP53+ tumors to controls, there was a strong association with increasing cumulative REC exposure (ORTertile 1 vs. Unexposed=1.4, 95% CI=0.7, 2.8; ORTertile 2 vs. Unexposed=1.3, 95% CI=0.6, 2.7; ORTertile 3 vs. Unexposed=2.8, 95% CI=1.4, 5.7; p-trend=0.003), while there was no apparent association between cumulative REC exposure and bladder cancer among patients with wild-type TP53 tumors (p-trend=0.603). Subtype analyses that incorporate clinical features (stage: non-muscle invasive/muscle-invasive and grade: low-grade/high-grade) with TP53 mutation, also indicate etiologic heterogeneity in the relationship between diesel exhaust and bladder cancer (p-het=0.002), with the strongest association observed between increasing cumulative REC and risk of bladder cancer among cases with non-muscle invasive, high-grade tumors (including carcinoma in situ) with TP53 mutations (OR per 100 μg/m3-yrs =1.14, 95% CI=1.00, 1.29).

Conclusions: Our results suggest that diesel engine exhaust may influence bladder cancer risk by impacting the cell cycle, especially among patients with high-risk non-muscle invasive tumors that are more likely to progress to invasive disease. More work will be needed to understand the impact of diesel-generated oxidative stress on tumor suppressor genes like TP53 and cell cycle regulation in the development of bladder cancer.

Citation Format: Nicole Gonzalez, Nina Rao, Michael Dean, Ludmila Prokunina-Olsson, Dalsu Baris, Molly Schwenn, Alison Johnson, Nathaniel Rothman, Debra Silverman, Stella Koutros. Diesel engine exhaust is associated with TP53 mutations and high-risk non-muscle invasive bladder cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 861.

Germline variants in hereditary breast cancer genes are associated with early age at diagnosis and family history in Guatemalan breast cancer

Purpose

Mutations in hereditary breast cancer genes play an important role in the risk for cancer.

Methods

Cancer susceptibility genes were sequenced in 664 unselected breast cancer cases from Guatemala. Variants were annotated with ClinVar and VarSome.

Results

A total of 73 out of 664 subjects (11%) had a pathogenic variant in a high or moderate penetrance gene. The most frequently mutated genes were BRCA1 (37/664, 5.6%) followed by BRCA2 (15/664, 2.3%), PALB2 (5/664, 0.8%), and TP53 (5/664, 0.8%). Pathogenic variants were also detected in the moderate penetrance genes ATM, BARD1, CHEK2, and MSH6. The high ratio of BRCA1/BRCA2 mutations is due to two potential founder mutations: BRCA1 c.212 + 1G > A splice mutation (15 cases) and BRCA1 c.799delT (9 cases). Cases with pathogenic mutations had a significantly earlier age at diagnosis (45 vs 51 years, P < 0.001), are more likely to have had diagnosis before menopause, and a higher percentage had a relative with any cancer (51% vs 37%, P = 0.038) or breast cancer (33% vs 15%, P < 0.001).

Conclusions

Hereditary breast cancer mutations were observed among Guatemalan women, and these women are more likely to have early age at diagnosis and family history of cancer. These data suggest the use of genetic testing in breast cancer patients and those at high risk as part of a strategy to reduce breast cancer mortality in Guatemala.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06305-5.

Aflatoxin B1 exposure and liver cirrhosis in Guatemala: a case-control study

Objective

In Guatemala, cirrhosis is among the 10 leading causes of death, and mortality rates have increased lately. The reasons for this heavy burden of disease are not clear as the prevalence of prominent risk factors, such as hepatitis B virus, hepatitis C virus and heavy alcohol consumption, appears to be low. Aflatoxin B₁ (AFB₁) exposure, however, appears to be high, and thus could be associated with the high burden of cirrhosis. Whether AFB₁ increases the risk of cirrhosis in the absence of viral infection, however, is not clear.

Design

Cirrhosis cases (n=100) from two major referral hospitals in Guatemala City were compared with controls (n=200) from a cross-sectional study. Logistic regression was used to estimate the ORs and 95% CIs of cirrhosis and quintiles of AFB₁ in crude and adjusted models. A sex-stratified analysis was also conducted.

Results

The median AFB₁ level was significantly higher among the cases (11.4 pg/mg) than controls (5.11 pg/mg). In logistic regression analyses, higher levels of AFB₁ was associated with cirrhosis (quintile 5 vs quintile 1, OR: 11.55; 95% CI 4.05 to 32.89). No attenuation was observed with adjustment by sex, ethnicity, hepatitis B virus status, and heavy alcohol consumption. A significantly increasing trend in association was observed in both models (p trend <0.01). Additionally, the cirrhosis–AFB₁ association was more prominent among men.

Conclusions

The current study found a significant positive association between AFB₁ exposure and cirrhosis. Mitigation of AFB₁ exposure and a better understanding of additional risk factors may be important to reduce the burden of cirrhosis in Guatemala.

Genome-wide homozygosity and risk of four non-Hodgkin lymphoma subtypes

AIM

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

METHODS

We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.

RESULTS

We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10-6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10-5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.

CONCLUSION

Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

Integrative molecular characterisation of gallbladder cancer reveals micro-environment-associated subtypes

BACKGROUND & AIMS

Gallbladder cancer (GBC) is the most common type of biliary tract cancer, but the molecular mechanisms involved in gallbladder carcinogenesis remain poorly understood. In this study, we applied integrative genomics approaches to characterise GBC and explore molecular subtypes associated with patient survival.

METHODS

We profiled the mutational landscape of GBC tumours (whole-exome sequencing on 92, targeted sequencing on 98, in total 190 patients). In a subset (n = 45), we interrogated the matched transcriptomes, DNA methylomes, and somatic copy number alterations. We explored molecular subtypes identified through clustering tumours by genes whose expression was associated with survival in 47 tumours and validated subtypes on 34 publicly available GBC cases.

RESULTS

Exome analysis revealed TP53 was the most mutated gene. The overall mutation rate was low (median 0.82 Mut/Mb). APOBEC-mediated mutational signatures were more common in tumours with higher mutational burden. Aflatoxin-related signatures tended to be highly clonal (present in ≥50% of cancer cells). Transcriptome-wide survival association analysis revealed a 95-gene signature that stratified all GBC patients into 3 subtypes that suggested an association with overall survival post-resection. The 2 poor-survival subtypes were associated with adverse clinicopathologic features (advanced stage, pN1, pM1), immunosuppressive micro-environments (myeloid-derived suppressor cell accumulation, extensive desmoplasia, hypoxia) and T cell dysfunction, whereas the good-survival subtype showed the opposite features.

CONCLUSION

These data suggest that the tumour micro-environment and immune profiles could play an important role in gallbladder carcinogenesis and should be evaluated in future clinical studies, along with mutational profiles.

LAY SUMMARY

Gallbladder cancer is highly fatal, and its causes are poorly understood. We evaluated gallbladder tumours to see if there were differences between tumours in genetic information such as DNA and RNA. We found evidence of aflatoxin exposure in these tumours, and immune cells surrounding the tumours were associated with survival.

Lack of transgenerational effects of ionizing radiation exposure from the Chernobyl accident

Effects of radiation exposure from the Chernobyl nuclear accident remain a topic of interest. We investigated germline de novo mutations (DNMs) in children born to parents employed as cleanup workers or exposed to occupational and environmental ionizing radiation after the accident. Whole-genome sequencing of 130 children (born 1987-2002) and their parents did not reveal an increase in the rates, distributions, or types of DNMs relative to the results of previous studies. We find no elevation in total DNMs, regardless of cumulative preconception gonadal paternal [mean = 365 milligrays (mGy), range = 0 to 4080 mGy] or maternal (mean = 19 mGy, range = 0 to 550 mGy) exposure to ionizing radiation. Thus, we conclude that, over this exposure range, evidence is lacking for a substantial effect on germline DNMs in humans, suggesting minimal impact from transgenerational genetic effects.

Targeted Deep Sequencing of Bladder Tumors Reveals Novel Associations between Cancer Gene Mutations and Mutational Signatures with Major Risk Factors

PURPOSE

Exome- and whole-genome sequencing of muscle-invasive bladder cancer has revealed important insights into the molecular landscape; however, there are few studies of non-muscle-invasive bladder cancer with detailed risk factor information.

EXPERIMENTAL DESIGN

We examined the relationship between smoking and other bladder cancer risk factors and somatic mutations and mutational signatures in bladder tumors. Targeted sequencing of frequently mutated genes in bladder cancer was conducted in 322 formalin-fixed paraffin-embedded bladder tumors from a population-based case-control study. Logistic regression was used to calculate odds ratios (OR) and 95% confidence intervals (CI), evaluating mutations and risk factors. We used SignatureEstimation to extract four known single base substitution mutational signatures and Poisson regression to calculate risk ratios (RR) and 95% CIs, evaluating signatures and risk factors.

RESULTS

Non-silent KDM6A mutations were more common in females than males (OR = 1.83; 95% CI, 1.05-3.19). There was striking heterogeneity in the relationship between smoking status and established single base substitution signatures: current smoking status was associated with greater ERCC2-Signature mutations compared with former (P = 0.024) and never smoking (RR = 1.40; 95% CI, 1.09-1.80; P = 0.008), former smoking was associated with greater APOBEC-Signature13 mutations (P = 0.05), and never smoking was associated with greater APOBEC-Signature2 mutations (RR = 1.54; 95% CI, 1.17-2.01; P = 0.002). There was evidence that smoking duration (the component most strongly associated with bladder cancer risk) was associated with ERCC2-Signature mutations and APOBEC-Signature13 mutations among current (P _trend = 0.005) and former smokers (P = 0.0004), respectively.

CONCLUSIONS

These data quantify the contribution of bladder cancer risk factors to mutational burden and suggest different signature enrichments among never, former, and current smokers.

Aflatoxin and the Etiology of Liver Cancer and Its Implications for Guatemala

During the 60 years since the first scientific reports about a relation between aflatoxin exposure and adverse health consequences, both in animals and humans, there has been a remarkable number of basic, clinical and population science studies characterizing the impact of this mycotoxin on diseases such as liver cancer. Many of these human investigations to date have focused on populations residing in Asia and Africa due to the high incidence of liver cancer and high exposures to aflatoxin. These studies formed the basis for the International Agency for Research on Cancer to classify the aflatoxins as Group 1 known human carcinogens. In addition, aflatoxin contamination levels have been used in international commodity trade to set the price of various staples such as maize and groundnuts. While there have been many case-control and prospective cohort studies of liver cancer risk over the years there have been remarkably few investigations focused on liver cancer in Latin America. Our interdisciplinary and multiple institutional collaborative has been developing a long-term strategy to characterize the role of aflatoxin and other mycotoxins as health risk factors in Guatemala and neighboring countries. This paper summarizes a number of the investigations to date and provides a roadmap of our strategies for the near term to discern the emergent etiology of liver cancer in this region. With these data in hand public health-based prevention strategies could be strategically implemented and conducted to lower the impact of these mycotoxins on human health.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in the Childhood Cancer Survivor Study

Background

Pediatric cancers are the leading cause of death by disease in children despite improved survival rates overall. The contribution of germline genetic susceptibility to pediatric cancer survivors has not been extensively characterized. We assessed the frequency of pathogenic or likely pathogenic (P/LP) variants in 5451 long-term pediatric cancer survivors from the Childhood Cancer Survivor Study.

Methods

Exome sequencing was conducted on germline DNA from 5451 pediatric cancer survivors (cases who survived ≥5 years from diagnosis; n = 5105 European) and 597 European cancer-free adults (controls). Analyses focused on comparing the frequency of rare P/LP variants in 237 cancer-susceptibility genes and a subset of 60 autosomal dominant high-to-moderate penetrance genes, for both case-case and case-control comparisons.

Results

Of European cases, 4.1% harbored a P/LP variant in high-to-moderate penetrance autosomal dominant genes compared with 1.3% in controls (2-sided P = 3 × 10-4). The highest frequency of P/LP variants was in genes typically associated with adult onset rather than pediatric cancers, including BRCA1/2, FH, PALB2, PMS2, and CDKN2A. A statistically significant excess of P/LP variants, after correction for multiple tests, was detected in patients with central nervous system cancers (NF1, SUFU, TSC1, PTCH2), Wilms tumor (WT1, REST), non-Hodgkin lymphoma (PMS2), and soft tissue sarcomas (SDHB, DICER1, TP53, ERCC4, FGFR3) compared with other pediatric cancers.

Conclusion

In long-term pediatric cancer survivors, we identified P/LP variants in cancer-susceptibility genes not previously associated with pediatric cancer as well as confirmed known associations. Further characterization of variants in these genes in pediatric cancer will be important to provide optimal genetic counseling for patients and their families.

Origins, Admixture Dynamics, and Homogenization of the African Gene Pool in the Americas

The Transatlantic Slave Trade transported more than 9 million Africans to the Americas between the early 16th and the mid-19th centuries. We performed a genome-wide analysis using 6,267 individuals from 25 populations to infer how different African groups contributed to North-, South-American, and Caribbean populations, in the context of geographic and geopolitical factors, and compared genetic data with demographic history records of the Transatlantic Slave Trade. We observed that West-Central Africa and Western Africa-associated ancestry clusters are more prevalent in northern latitudes of the Americas, whereas the South/East Africa-associated ancestry cluster is more prevalent in southern latitudes of the Americas. This pattern results from geographic and geopolitical factors leading to population differentiation. However, there is a substantial decrease in the between-population differentiation of the African gene pool within the Americas, when compared with the regions of origin from Africa, underscoring the importance of historical factors favoring admixture between individuals with different African origins in the New World. This between-population homogenization in the Americas is consistent with the excess of West-Central Africa ancestry (the most prevalent in the Americas) in the United States and Southeast-Brazil, with respect to historical-demography expectations. We also inferred that in most of the Americas, intercontinental admixture intensification occurred between 1750 and 1850, which correlates strongly with the peak of arrivals from Africa. This study contributes with a population genetics perspective to the ongoing social, cultural, and political debate regarding ancestry, admixture, and the mestizaje process in the Americas.