The epigenetic modification of mitochondrial DNA (mtDNA) is still in controversy. To clarify this point, we applied the gold standard method for DNA methylation, bisulfite pyrosequencing, to examine human mtDNA methylation status. Before bisulfite conversion, BamHI was used to digest DNA to open the loop of mtDNA. The results demonstrated that the linear mtDNA had significantly higher bisulfite conversion efficiency compared with circular mtDNA. Furthermore, the methylation values obtained from linear mtDNA were significantly lower than that of circular mtDNA, which was verified by SEQUENOM MassARRAY. The above impacts of circular structure were also observed in lung DNA samples but not in saliva DNA samples. Mitochondrial genome methylation of blood samples and saliva samples from 14 unrelated individuals was detected. The detected regions covered 83 CpG sites across mtDNA including D-loop, 12 S rRNA, 16 S rRNA, ND1, COXI, ND3, ND4, ND5, CYTB. We found that the average methylation levels of nine regions were all less than 2% for both sample types. In conclusion, our findings firstly show that the circular structure of mtDNA affects bisulfite conversion efficiency, which leads to overestimation of mtDNA methylation values. CpG methylation in human mtDNA is a very rare event at most DNA regions.

Aging is a complicated process characterized by progressive and extensive changes in physiological homeostasis at the organismal, tissue, and cellular levels. In modern society, age estimation is essential in a large variety of legal rights and duties. Accumulating evidence suggests roles for microRNAs (miRNAs) and circular RNAs (circRNAs) in regulating numerous processes during aging. Here, we performed circRNA sequencing in two age groups and analyzed microarray data of 171 healthy subjects (17–104 years old) downloaded from Gene Expression Omnibus (GEO) and ArrayExpress databases with integrated bioinformatics methods. A total of 1,403 circular RNAs were differentially expressed between young and old groups, and 141 circular RNAs were expressed exclusively in elderly samples while 10 circular RNAs were expressed only in young subjects. Based on their expression pattern in these two groups, the circular RNAs were categorized into three classes: age-related expression between young and old, age-limited expression-young only, and age-limited expression-old only. Top five expressed circular RNAs among three classes and a total of 18 differentially expressed microRNAs screened from online databases were selected to validate using RT-qPCR tests. An independent set of 200 blood samples (20–80 years old) was used to develop age prediction models based on 15 age-related noncoding RNAs (11 microRNAs and 4 circular RNAs). Different machine learning algorithms for age prediction were applied, including regression tree, bagging, support vector regression (SVR), random forest regression (RFR), and XGBoost. Among them, random forest regression model performed best in both training set (mean absolute error = 3.68 years, r = 0.96) and testing set (MAE = 6.840 years, r = 0.77). Models using one single type of predictors, circular RNAs-only or microRNAs-only, result in bigger errors. Smaller prediction errors were shown in males than females when constructing models according to different-sex separately. Putative microRNA targets (430 genes) were enriched in the cellular senescence pathway and cell homeostasis and cell differentiation regulation, indirectly indicating that the microRNAs screened in our study were correlated with development and aging. This study demonstrates that the noncoding RNA aging clock has potential in predicting chronological age and will be an available biological marker in routine forensic investigation to predict the age of biological samples.

In forensic science, accurate estimation of the age of a victim or suspect can facilitate the investigators to narrow a search and aid in solving a crime. Aging is a complex process associated with various molecular regulations on DNA or RNA levels. Recent studies have shown that circular RNAs (circRNAs) upregulate globally during aging in multiple organisms such as mice and C.elegans because of their ability to resist degradation by exoribonucleases. In the current study, we attempted to investigate circRNAs’ potential capability of age prediction. Here, we identified more than 40,000 circRNAs in the blood of thirteen Chinese unrelated healthy individuals with ages of 20–62 years according to their circRNA-seq profiles. Three methods were applied to select age-related circRNA candidates including the false discovery rate, lasso regression, and support vector machine. The analysis uncovered a strong bias for circRNA upregulation during aging in human blood. A total of 28 circRNAs were chosen for further validation in 30 healthy unrelated subjects by RT-qPCR, and finally, 5 age-related circRNAs were chosen for final age prediction models using 100 samples of 19–73 years old. Several different algorithms including multivariate linear regression (MLR), regression tree, bagging regression, random forest regression (RFR), and support vector regression (SVR) were compared based on root mean square error (RMSE) and mean average error (MAE) values. Among five modeling methods, regression tree and RFR performed better than the others with MAE values of 8.767 years (S.rho = 0.6983) and 9.126 years (S.rho = 0.660), respectively. Sex effect analysis showed age prediction models significantly yielded smaller prediction MAE values for males than females (MAE = 6.133 years for males, while 10.923 years for females in the regression tree model). In the current study, we first used circRNAs as additional novel age-related biomarkers for developing forensic age estimation models. We propose that the use of circRNAs to obtain additional clues for forensic investigations and serve as aging indicators for age prediction would become a promising field of interest.