Integrated analysis of mRNA and extrachromosomal circular DNA profiles to identify the potential mRNA biomarkers in breast cancer

Extrachromosomal circular DNAs (eccDNAs) have been proved an inseparable relationship with cancer, based on the biological mechanisms of its biogenesis and impact on tumorigenesis, but still lacked of methods to analyze its function on the pathogenesis and progression of breast cancer (BC). The mRNA and eccDNA from BC cell samples (MDA-MB-453 and MCF-12A) were extracted with the removal of rRNA and linear DNA, respectively. High-throughput sequencing and bioinformatics analysis were performed to explore their expression level and molecular characterization of eccDNA. A total number of 161,062 eccDNA ranging from 33 bp to 54229 bp were detected with a median size of 1143 bp, distributed on all chromosomes and enriched on chromosome 20 the most. EccDNAs located in exons, upstream and downstream 2 kb regions were significantly increased compared with background. Analysis of eccDNA-related differentially expressed genes (eccDEGs) showed that FAT2 properly separated the two cells. CTNNB1, CACNA2D2 and CACNA1D were the hub genes with higher degrees in critical modules. All these four genes were significantly differentially expressed between breast invasive carcinoma (BRCA) tissues and normal ones. FAT2 and CTNNB1 correlated with significantly different overall survival (OS) when differentially expressed. The four genes showed a strong correlation with each other significantly and changed between tumor and normal samples. The results showed the potential of FAT2, CTNNB1, CACNA2D2 and CACNA1D as biomarkers with analysis of both DEGs and eccDEGs, which might assist in clinical medical treatment.