Extrachromosomal circular DNAs are common and functional in esophageal squamous cell carcinoma

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) is the leading cause of cancer-related mortality. While recent studies have documented the presence of extrachromosomal circular DNAs (eccDNAs) in various tumors, to date, there have been no studies examining the distribution and function of eccDNAs in ESCC. METHODS: The eccDNAs from three surgically matched ESCC tissue samples were extracted and amplified by rolling circle amplification after removal of linear DNA and mitochondrial circular DNA. High-throughput eccDNA sequencing and bioinformatics analysis was performed to study the distribution pattern and the level of eccDNA expression. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the genes associated with the differentially expressed eccDNAs. Five up-regulated and five down-regulated candidate eccDNAs were validated by routine polymerase chain reaction (PCR), TOPO-TA cloning and Sanger sequencing. The nucleotides flanking the eccDNA junctions were analyzed to explore the mechanisms of eccDNA formation. RESULTS: A total of 184,557 eccDNAs was identified. The overall length distribution ranged from 33 to 968,842 base pairs (bp), with the peak at approximately 360 bp. These eccDNAs mainly originated from 5'- and 3'-untranslated regions (UTRs), and rarely from exons, introns, LINE, or Alu repeat regions. The chromosome distribution, length distribution, and genomic annotation of the eccDNAs were comparable between ESCC samples and matched normal epithelium. Nevertheless, 16,031 eccDNAs were found to be differentially expressed between ESCC and matched normal epithelium, including 10,126 up-regulated eccDNAs and 5,905 down-regulated eccDNAs. GO analysis and KEGG pathway analysis showed enriched in cancer pathways, mitogen-activated protein kinase (MAPK) pathway, GTPase-related activity, and cytoskeleton function. PCR, TOPO-TA cloning, and Sanger sequencing validated the junctional sites of five up-regulated candidate eccDNAs and four other unexpected eccDNAs. A repeat nucleotide pattern between the position flanking the start site and that flanking the end site was detected. CONCLUSIONS: This study demonstrated the genome-wide presence of eccDNAs, explored the differential expression of eccDNAs, and revealed the potential mechanisms of eccDNAs in ESCC. This work provides further insights into our understanding of genome plasticity, the role of eccDNAs in ESCC, and may contribute to the development of potential clinical therapies.