Differential expression and analysis of extrachromosomal circular DNAs as serum biomarkers in pulmonary arterial hypertension

Abstract

Extrachromosomal circular DNAs (eccDNAs) have been reported to play a key role in the occurrence and development of various diseases. However, the characterization and role of eccDNAs in pulmonary arterial hypertension (PAH) remain unclear. In the discovery cohort, we first explored eccDNA expression profiles by Circle-sequencing analysis. The candidate eccDNAs were validated by routine polymerase chain reaction (PCR), TOPO-TA cloning and Sanger sequencing. In the validation cohort, 30 patients with PAH and 10 healthy controls were recruited for qPCR amplification to detect the candidate eccDNAs. Datas at the baseline were collected, including clinical background, biochemical variables, echocardiography and hemodynamic factors. Receiver operating characteristic curve was used to investigate the diagnostic effect of the eccDNA. We identified a total of 21,741 eccDNAs in plasma samples of 3 IPAH patients and 3 individuals in good health, and the expression frequency, GC content, length distribution, and genome distribution of the eccDNAs were thoroughly characterized and analyzed. In the validation cohort, 687 eccDNAs were differentially expressed in patients with IPAH compared with healthy controls (screening threshold: |FC|≥2 and P < 0.05). Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the specific eccDNAs in IPAH were significantly enriched in calcium channel activity, the mitogen-activated protein kinase pathway, and the wnt signaling pathway. Verification queue found that the expression of eccDNA-chr2:131208878–131,424,362 in PAH was considerably higher than that in healthy controls and exhibited a high level of accuracy in predicting PAH with a sensitivity of 86.67% and a specificity of 90%. Furthermore, correlation analysis disclosed a significant association between serum eccDNA-chr2:131208878–131,424,362 and mean pulmonary artery pressure (mPAP) (r = 0.396, P = 0.03), 6 min walking distance (6MWD) (r = -0.399, P = 0.029), N-terminal pro-B-type natriuretic peptide (NT-proBNP) (r = 0.685, P < 0.001) and cardiac index (CI) (r = − 0.419, P = 0.021). This is the first study to identify and characterize eccDNAs in patients with PAH. We revealed that serum eccDNA-chr2:131208878–131,424,362 is significantly overexpressed and can be used in the diagnosis of PAH, indicating its potential as a novel non-invasive biomarker.