Mutational signatures in esophageal squamous cell carcinoma from eight countries with varying incidence

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

Dynamics of the COVID-19 Clinical Findings and the Serologic Response

The factors affecting the dynamics of lengthening of symptoms and serologic responses are not well known. In order to see how the serologic responses change in relation to the clinical features, we selected a group of 472 adults with a positive IgM/IgG antibody test result from a baseline study of the anti-SARS-CoV-2 seropositivity, assessed their COVID-19 and past medical histories, and followed them up in about 3 months. Nearly one-fourth of the subjects were asymptomatic at the baseline; 12.8% subjects became symptomatic at the follow-up (FU) when 39.8% of the subjects had some persisting symptoms. At the baseline, 6.1% showed anti-SARS-CoV-2 IgM positive, 59.3% only for IgG, and 34.5% for both. At the FU, these figures declined to 0.6, 54.0, and 4.4%, respectively, with the mean IgM and IgG levels declining about 6.3 and 2.5 folds. Blood group A was consistently linked to both sustaining and flipping of the gastrointestinal (GI) and respiratory symptoms. The baseline IgM level was associated with GI symptoms and pre-existing cirrhosis in multivariate models. Both of the baseline and FU IgG levels were strongly associated with age, male, and lung involvement seen in chest computed tomography (CT)-scan. Finally, as compared with antibody decayers, IgM sustainers were found to be more anosmic [mean difference (MD): 11.5%; P = 0.047] with lower body mass index (BMI) (MD: 1.30 kg/m²; P = 0.002), while IgG sustainers were more commonly females (MD: 19.2%; P = 0.042) with shorter diarrhea duration in the FU (MD: 2.8 days; P = 0.027). Our findings indicate how the anti-SARS-CoV-2 serologic response and COVID-19 clinical presentations change in relation to each other and basic characteristics.

SARS-CoV-2 antibody seroprevalence in the general population and high-risk occupational groups across 18 cities in Iran: a population-based cross-sectional study

Background

Rapid increases in cases of COVID-19 were observed in multiple cities in Iran towards the start of the pandemic. However, the true infection rate remains unknown. We aimed to assess the seroprevalence of antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 18 cities of Iran as an indicator of the infection rate.

Methods

In this population-based cross-sectional study, we randomly selected and invited study participants from the general population (from lists of people registered with the Iranian electronic health record system or health-care centres) and a high-risk population of individuals likely to have close social contact with SARS-CoV-2-infected individuals through their occupation (from employee lists provided by relevant agencies or companies, such as supermarket chains) across 18 cities in 17 Iranian provinces. Participants were asked questions on their demographic characteristics, medical history, recent COVID-19-related symptoms, and COVID-19-related exposures. Iran Food and Drug Administration-approved Pishtaz Teb SARS-CoV-2 ELISA kits were used to detect SARS-CoV-2-specific IgG and IgM antibodies in blood samples from participants. Seroprevalence was estimated on the basis of ELISA test results and adjusted for population weighting (by age, sex, and city population size) and test performance (according to our independent validation of sensitivity and specificity).

Findings

From 9181 individuals who were initially contacted between April 17 and June 2, 2020, 243 individuals refused to provide blood samples and 36 did not provide demographic information and were excluded from the analysis. Among the 8902 individuals included in the analysis, 5372 had occupations with a high risk of exposure to SARS-CoV-2 and 3530 were recruited from the general population. The overall population weight-adjusted and test performance-adjusted prevalence of antibody seropositivity in the general population was 17·1% (95% CI 14·6–19·5), implying that 4 265 542 (95% CI 3 659 043–4 887 078) individuals from the 18 cities included were infected by the end of April, 2020. The adjusted seroprevalence of SARS-CoV-2-specific antibodies varied greatly by city, with the highest estimates found in Rasht (72·6% [53·9–92·8]) and Qom (58·5% [37·2–83·9]). The overall population weight-adjusted and test performance-adjusted seroprevalence in the high-risk population was 20·0% (18·5–21·7) and showed little variation between the occupations included.

Interpretations

Seroprevalence is likely to be much higher than the reported prevalence of COVID-19 based on confirmed COVID-19 cases in Iran. Despite high seroprevalence in a few cities, a large proportion of the population is still uninfected. The potential shortcomings of current public health policies should therefore be identified to prevent future epidemic waves in Iran.

Funding

Iranian Ministry of Health and Medical Education.

Translation

For the Farsi translation of the abstract see Supplementary Materials section.

Abstract 1318: APOBEC/ADAR aberrations are potentially implicated in certain pancancer hypermutation patterns

Background:

Known endogenous mutators include certain APOBEC family members, whose mRNA levels have been implicated in various cancers. Other important nucleic acid deaminases include ADAR family of proteins, with some recent evidence supporting their potential mutator activity. Some external mutators have also been reported to act through endogenous ones. Although the ADAR and particularly APOBEC mutations are infrequent among cancer genomes, their potential contributions in cancer mutagenesis might not be underestimated.

Methods:

We assessed 10,294 exomes from 32 TCGA cancer types for mutations among fourteen APOBEC/ADAR family members. The mean number of single nucleotide variants (SNVs) and insertions/deletions (indels) was compared among those exomes with nonsilent APOBEC/ADAR mutations and those with wild-type APOBEC/ADAR or silent mutations. Multivariate analyses identified the relative impacts of APOBEC/ADAR mutations and mRNA levels as well as DNA repair/replication defects, as determined by mutations in KEGG pathway genes, on main variant classes across pancancer exomes.

Results:

Preliminary analysis identified higher mean number of both SNVs and indels among those exomes with any of the fourteen APOBEC/ADAR genes mutated compared to those with wild-type ones. Further analyses showed higher mean number of certain indel variants among those exomes with nonsilent APOBEC1 (A1; 6.1 fold), A3B (5.1 fold), A3C (20.5 fold), or ADAR (2.2 fold) mutations, and higher mean number of T>N SNVs (2.8 fold) among those with nonsilent ADARB1 mutations compared to silent ones, suggesting their true implication in cancer hypermutation. Besides DNA repair/replication defects, multivariate analyses revealed that C>G variants were determined mostly by APOBEC (A3B, A3A, A3G) mRNA levels, while other SNV classes determined mainly by various APOBEC/ADAR mutations. Indels were found to be determined mainly by mutations in A1, A3B, ADAR, and ADARB1, with a small impact of A1 mRNA level. The impacts of DNA repair/replication defects were more prominent on indels as well as C>G SNVs.

Conclusions:

These results provide a landscape of cancer variants as correlated with endogenous mutators and DNA repair defects. They propose the implication of certain APOBEC/ADAR mutations in most SNV classes as well as indels, in contrast to C>G SNVs mostly determined by APOBEC overexpression. These potentially reflect the gain/ loss of function impacts of various mutations and/or oligomerization capacities of both families. Further analyses linking APOBEC/ADAR mutations and mRNA levels with distinct mutational signatures are shown in the accompanying abstract “Distinct pancancer mutational signatures are determined by APOBEC/ADAR aberrations”.

Citation Format: Ahmadreza Niavarani, Milad Bagheri, Joseph CF Ng, Franca Fraternali, Reza Malekzadeh. APOBEC/ADAR aberrations are potentially implicated in certain pancancer hypermutation patterns [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 1318.

Pancancer analysis identifies prognostic high-APOBEC1 expression level implicated in cancer in-frame insertions and deletions

Genome insertions and deletions (indels) show tremendous functional impacts despite they are much less common than single nucleotide variants, which are at the center of studies assessing cancer mutational signatures. We studied 8891 tumor samples of 32 types from The Cancer Genome Atlas in order to explore those genes which are potentially implicated in cancer indels. Survival analysis identified in-frame indels as the most important variants predicting adverse outcome. Transcriptome-wide association study identified 16 genes overexpressed in both tumor samples and tumor types with high number of in-frame indels, of whom four (APOBEC1, BCL2L15, FOXL1 and PDX1) were identified with gene products distributed within the nucleus. APOBEC1 emerged as the mere consistently hypomethylated gene in tumor samples with high number of in-frame indels. The correlation of APOBEC1 expression levels with cancer indels was independent of age and defects in DNA homologous recombination (HR) and/or mismatch repair. Unlike frame-shift indels, triplet repeat motifs were found to occur frequently at in-frame indel sites. The splicing variant 3, making a shorter isoform b, showed essentially all the same indel correlations as of APOBEC1. Expression levels of both APOBEC1 and variant 3 were found to be predicting adverse prognosis independent of DNA HR and mismatch repair. Not less importantly, high level of variant 3 in paired normal tissues was also proved to predict cancer outcome. Our findings propose APOBEC1 and isoform b as the potential endogenous mutators implicated in cancer in-frame indels and pave the way for their use as novel prognostic tumor markers.

The Wilms Tumor-1 (WT1) rs2234593 variant is a prognostic factor in normal karyotype acute myeloid leukemia

The single-nucleotide polymorphism (SNP) within Wilms tumor-1 (WT1) exon 7, rs16754, has been arguably reported to be implicated in acute myeloid leukemia (AML) prognosis. We assessed the potential association of selected WT1 SNPs as well as WT1 mutations in normal karyotype (NK)-AML and evaluated the prognostic value of these normal gene variants. Diagnostic samples from a series of 474 young adult NK-AML patients were used to genotype five WT1 SNPs using TaqMan assays and to directly sequence WT1 exons 7 and 9. Analysis of five WT1 gene variants showed an association of rs2234593 allele C with WT1 Ex7 mutation. Prognostic study of the same variants identified rs2234593 significantly associated with relapse and overall survival (OS). Patients with rs2234593AA/AC showed significantly higher 10-year OS (50 vs 36 %, hazard ratio (HR) = 0.69 (0.52–0.90), p = 0.006) and lower cumulative incidence of relapse (CIR) (36 vs 51 %, HR = 0.62 (0.45–0.86), p = 0.004) compared to those with rs2234593CC. The effect of AA genotype on CIR remained significant after adjustment for basic covariates including FLT3 internal-tandem duplication (FLT3-ITD) and nucleophosmin 1 (NPM1) mutations (HR = 0.60 (0.41–0.89), p = 0.009), with some evidence of improved survival (HR = 0.75 (0.55–1.03), p = 0.07). A multivariate analysis showed WT1 Ex7-mutant as the major relapse predictor, with a tendency for rs2234593-A effect after allowing for Ex7 mutation (p = 0.07). No adjusted risk benefit was found for previously reported rs16754-G. In conclusion, WT1 normal gene variant rs2234593 is associated with mutational status of WT1 Ex7 and is a further prognostic marker independent from FLT3-ITD and NPM1 mutations in NK-AML.

A 4-gene expression score associated with high levels of Wilms Tumor-1 (WT1) expression is an adverse prognostic factor in acute myeloid leukaemia

Wilms Tumor-1 (WT1) expression level is implicated in the prognosis of acute myeloid leukaemia (AML). We hypothesized that a gene expression profile associated with WT1 expression levels might be a good surrogate marker. We identified high WT1 gene sets by comparing the gene expression profiles in the highest and lowest quartiles of WT1 expression in two large AML studies. Two high WT1 gene sets were found to be highly correlated in terms of the altered genes and expression profiles. We identified a 17-probe set signature of the high WT1 set as the optimal prognostic predictor in the first AML set, and showed that it was able to predict prognosis in the second AML series after adjustment for European LeukaemiaNet genetic groups. The gene signature also proved to be of prognostic value in a third AML series of 163 samples assessed by RNA sequencing, demonstrating its cross-platform consistency. This led us to derive a 4-gene expression score, which faithfully predicted adverse outcome. In conclusion, a short gene signature associated with high WT1 expression levels and the resultant 4-gene expression score were found to be predictive of adverse prognosis in AML. This study provides new clues to the molecular pathways underlying high WT1 states in leukaemia.

APOBEC3A is implicated in a novel class of G-to-A mRNA editing in WT1 transcripts

Classic deamination mRNA changes, including cytidine to uridine (C-to-U) and adenosine to inosine (A-to-I), are important exceptions to the central dogma and lead to significant alterations in gene transcripts and products. Although there are a few reports of non-classic mRNA alterations, as yet there is no molecular explanation for these alternative changes. Wilms Tumor 1 (WT1) mutations and variants are implicated in several diseases, including Wilms tumor and acute myeloid leukemia (AML). We observed two alternative G-to-A changes, namely c.1303G>A and c.1586G>A in cDNA clones and found them to be recurrent in a series of 21 umbilical cord blood mononuclear cell (CBMC) samples studied. Two less conserved U-to-C changes were also observed. These alternative changes were found to be significantly higher in non-progenitor as compared to progenitor CBMCs, while they were found to be absent in a series of AML samples studied, indicating they are targeted, cell type-specific mRNA editing modifications. Since APOBEC/ADAR family members are implicated in RNA/DNA editing, we screened them by RNA-interference (RNAi) for WT1-mRNA changes and observed near complete reversal of WT1 c.1303G>A alteration upon APOBEC3A (A3A) knockdown. The role of A3A in mediating this change was confirmed by A3A overexpression in Fujioka cells, which led to a significant increase in WT1 c.1303G>A mRNA editing. Non-progenitor CBMCs showed correspondingly higher levels of A3A-mRNA and protein as compared to the progenitor ones. To our knowledge, this is the first report of mRNA modifying activity for an APOBEC3 protein and implicates A3A in a novel G-to-A form of editing. These findings open the way to further investigations into the mechanisms of other potential mRNA changes, which will help to redefine the RNA editing paradigm in both health and disease.

Development of Transgenic Mice Expressing Calcitonin as a Beta-lactoglobulin Fusion Protein in Mammary Gland

Expression of foreign proteins in mammalian milk is becoming a widespread strategy for high-level production of recombinant pharmaceuticals, especially those with the most complex post-translational modifications. A milk-specific ovine beta-lactoglobulin (oBLG) promoter was used to drive expression of recombinant calcitonin in mouse milk. A gene construct was generated, consisting of 10.7 kbp of the oBLG gene including its promoter and 3' flanking region with the calcitonin coding sequences inserted in-frame into the oBLG fifth exon. After microinjection, six founder mice transmitted the transgene to their progeny. RT-PCR confirmed mammary-gland specific expression of recombinant mRNA in most transgenic mice and Western blot analysis confirmed expression of chimeric protein. Calcitonin can thus be expressed under the oBLG promoter and regulatory elements in a mammary-gland specific manner.