Association Between Circular RNAs and Intracranial Aneurysm Rupture Under the Synergistic Effect of Individual Environmental Factors

Introduction: To study the association between specific circular RNAs and rupture of intracranial aneurysm. To explore its clinical diagnostic significance and synergistic effects with individual environmental influencing factors.

Methods: Three hundred and forty seven cases and controls were included in this study. Multivariate analysis was used to explore the main individual environmental factors. Intracranial aneurysm rupture related circular RNAs screened based on sequencing was verified in peripheral blood by PCR. ROC curve, logistic regression model and fork analysis were used to study the association, diagnostic values, and synergistic effects of circular RNA with intracranial aneurysms and individual environmental factors.

Results: Smoking, hair dyeing, sitting time ≥6 h/day, single animal oil intake and hypertension are the main risk factors for intracranial aneurysm rupture; People with higher education, sleeping time ≥7 h/day, tea drinking, diabetes, higher levels of (hemoglobin, low density lipoprotein, serum calcium, and apolipoprotein-A1) have a low risk of intracranial aneurysm rupture. Hsa_circ_0008433 and hsa_circ_0001946 are closely related to intracranial aneurysm rupture and have certain clinical diagnostic significance (AUC = 0.726; 95% CI: 0.668~0.784). Hsa_circ_0008433 (OR = 0.497, 95% CI: 0.338~0.731), hsa_circ_0001946 (OR = 0.682, 95% CI: 0.509~0.914) were independent epigenetic factors affecting intracranial aneurysm rupture, and have a multiplicative interaction with age (OR = 3.052, 95% CI: 1.006~9.258).

Conclusions: Low expressions of hsa_circ_0008433 and hsa_circ_0001946 are risk factors for intracranial aneurysms rupture and have good clinical diagnostic value. There was a multiplicative interaction between epigenetic score and age. The older and the higher the epigenetic score was, the more likely to have intracranial aneurysm rupture.

Introduction

Intracranial aneurysm (IA) is a localized abnormal bulging caused by high blood flow due to congenital weakness of the cerebral artery wall or acquired lesions (1). About 0.7 to 1.9% of IA can rupture under various inducements, leading to aneurysmal subarachnoid hemorrhage. The annual incidence of IA rupture worldwide is about 9.1 per 100,000, accounting for 3 to 5% of all acute strokes. The age of onset is between 40 and 60 years old, the mortality rate is close to 40%, and 46% of survivors can be disabled or long-term cognitive impairment due to multiple complications (2–5). Therefore, it is called the “time bomb” in the brain, which seriously threatens human life and health (6). Current research showed that the traditional risk factors for IA rupture are smoking, alcoholism, high blood pressure, gender, drugs, and blood lipid levels, etc (7–13). However, under the current social and economic development, the tremendous changes in people's life and habits will inevitably be accompanied by new individual environmental impact factors, and these new impact factors are not yet very clear.

Circular RNA is a type of novel nucleic acid molecule that does not have a characteristic 5'cap end and 3'poly (A) end, and exists in the form of a covalent loop (14). It is widely distributed, highly conservative in structure, and has specificity of temporal expression and spatial distribution. At present, circRNA has rapidly become a hot spot in the field of cerebrovascular disease research, and has been deeply explored as an emerging biological marker and therapeutic target (15). Recent studies have found that circRNA can competitively bind to specific miRNA through sponge adsorption to regulate the expression of downstream target genes. The related mechanism is called competitive endogenous RNA (ceRNA). A large number of subsequent studies have confirmed that the ceRNA mechanism involves many important pathological processes such as inflammation, SMC phenotype transformation and extracellular matrix, which play an important role in epigenetic regulation in vascular diseases (16, 17). However, there are still few studies on the relationship between circRNA and IA rupture recently and its synergistic effects with individual environmental factors are even rarer.

In this study, a questionnaire survey was conducted on 347 cases of IA rupture and 347 healthy examinees. Logistic regression analysis was used to study the individual environmental factors that affect the rupture of IAs. At the same time, based on the previous high-throughput sequencing results, the selected IA rupture-related circRNA were verified in the population peripheral blood, the potential correlation of these specific circRNA expression patterns in IA tissues and peripheral blood was further explored.

Materials

Study Population

A total of 347 patients with IA rupture and 347 healthy persons who were admitted to the Neurosurgery and Physical Examination Center of the Second Affiliated Hospital of Fujian Medical University from June 2017 to September 2019 were selected as case and control groups for individual environmental factors study. At the same time, 140 subjects with matching baseline data were selected from the case and the control group, their peripheral blood was collected for PCR experiments. Case group inclusion criteria: ① Age> 18 years old, and first onset; ② IA rupture confirmed by computed tomography angiography (CTA), magnetic resonance angiography (MRA), or digital subtraction angiography (DSA); ③ There was no significant systemic disease that might affect the indexes; The medical record is complete and willing to cooperate with the questionnaire. Case group exclusion criteria: ① Those with the first head CT scan or negative lumbar puncture; ② Patients with a clear history of subarachnoid hemorrhage or a family history of severe cerebrovascular disease; ③ Those with previous stroke, brain tumor, or history of craniocerebral surgery; ④ Patients with coagulopathy, fever, or pregnancy. Control group inclusion criteria: ① Immediate relatives of IA or a family history of a clear subarachnoid hemorrhage; ② Patients with previous stroke, brain tumor, or history of craniocerebral surgery; ③ Recently infected or pregnant; ④ Patients with systemic severe chronic diseases; ⑤ People with incomplete clinical data or refuse to participate in this study. All research subjects signed an informed consent form and approved by the Ethics Committee of the Second Affiliated Hospital of Fujian Medical University (2018-50).

Questionnaire

A unified structured questionnaire is adopted, and the research objects are investigated by strictly trained staff with certain epidemiological and clinical work experience. The content includes general conditions such as gender, age, occupation, education, marital status, personal medical history, and so on. Individual behaviors include smoking, drinking alcohol and tea, chemical poison exposure, taking medicine, sitting and sleeping time per day, physical exercise and labor intensity. Eating habits mainly include eating habits and edible oil types. Define smokers who smoke ≥1 cigarette per day for >6 consecutive months or cumulatively smoke ≥100 cigarettes; Define a single drinking volume of ≥50 g as drinking; Drinking tea ≥1 cup per week for more than 6 months is defined as drinking tea; Actively participate in sports activities at least once a week and last at least 20 min each time is defined as physical exercise (18). The daily sleep time of adults over 18 years old <7 h is defined as insufficient sleep, and the daily sleep time ≥9 h is defined as excessive sleep (19). Sedentary or semi-recumbent time >6 h was definited as sedentary (20).

RNA Extraction

Two milliliters peripheral venous blood was extracted and centrifuged. Lymphocyte separation solution was added to the blood cell layer, the suspension layer of white blood cells was further extracted. Then absorb the leukocyte suspension layer and add Trizol to fully lyse, further extract the total RNA of the peripheral blood according to the kit instructions. The NanoDrop^® ND-2000 nucleic acid quantifier was used to determine the RNA content. The specimens with OD260/280 values between 1.80 and 2.00 were considered qualified samples.

cDNA Synthesis and PCR

Strictly follow the instructions of PrimeScript^TM RT reagent Kit-RR037 kit for cDNA synthesis of circular RNA. According to the instructions of the TB Green^TM Premix EX Taq^TM-RR820 kit, the Real-Time PCR System amplifier (Light Cycler 480) was used for qRT-PCR detection, establish 20 μl system. The reaction conditions are: Pre-denaturation 95°C (30 s), 1 cycle; PCR reaction 95°C (5 s), 60°C (34 s), 45 cycles; melting reaction 95°C (15 s), 60°C (1 min), 95°C (15 s), 1 cycle. Take GAPDH as internal reference.

Statistical Analysis

Two independent sample t-tests were used to compare the difference between normal distribution quantitative data. Two independent sample rank sum tests (Mann-Whitney U) were used to analyze quantitative data that did not obey the normal distribution. Chi-square test was used to infer the overall rate or composition ratio. Multivariate logistic regression was used to analyze the odds ratio (OR) of each variable and the risk of IA rupture and its 95% confidence interval (95% CI). Real-time fluorescence quantitative PCR results are expressed by -ΔΔCT method. To study the interaction between epigenetic indicators and environmental factors by using fork analysis, and further calculate the interaction index (S), attribution ratio (AP), and excess relative risk (RERI).

Result

General Demographic Characteristics

A total of 347 cases and controls were included using the propensity score method. Among them, 392 were male (56.48%) and 302 were female (43.51%), with an average age of 51.91 ± 10.15 years. The general demographic characteristics of the case group and the control group are balanced and comparable (p > 0.05) (Supplementary Table 1).

Univariate Analysis of Individual Environmental Factors in IA Rupture

Correlation Between Behavioral Factors and IA Rupture

The analysis of behavioral factors showed that there were statistical differences between the case and control group in education level, exposure to chemical poisons, daily sitting and sleeping time, physical exercise, tea drinking, and smoking (p < 0.05) (Supplementary Table 2).

Correlation Between Eating Habits and IA Rupture

The analysis of the eating habits showed that the salty and light diet, edible oil types were statistically different between the case and the control group (p < 0.05) (Supplementary Table 3).

Correlation Between Specific Physiological Indicators and IA Rupture

The analysis of the physiological indicators showed that diastolic blood pressure and the pulse pressure were statistically different between the case and the control group (p < 0.05) (Supplementary Table 4).

Correlation Between Specific Biochemical Indicators and IA Rupture

The analysis of the biochemical indicators showed that hemoglobin, low density lipoprotein, triglyceride, cholesterol, serum calcium, apolipoprotein A1, and apolipoprotein B were statistically different between the case and the control group (p < 0.05) (Supplementary Table 5).

Correlation Between Disease History and IA Rupture

The analysis of the biochemical indicators showed that hypertension and stroke disease and family history are statistically different between the case group and the control group (p < 0.05) (Supplementary Table 6).

Logistic Regression Analysis of Individual Factors of IA Rupture

Taking IA rupture as the dependent variable, the single factors with statistically significant difference were further included in logistic regression model for analysis (Backward: Wald). The results showed that smoking, chemical poison exposure (hair dye), sitting time >6 h/day, single animal oil intake, hypertension, larger diastolic pressure and pulse pressure differences, higher levels of plasma globulin are the risk factors for IA rupture. Tea drinking, higher education, sleep time >7 h/day, diabetes, and higher levels of (hemoglobin, low density lipoprotein, apolipoprotein A1, and serum calcium) are the protective factors for IA rupture (Table 1).

TABLE 1

Table 1. Logistic regression analysis of individual factors of IA rupture.

Expression of IA Rupture-Related circRNAs in Peripheral Blood

Based on the previous high-throughput sequencing results, seven IA rupture-related circRNAs were selected (Supplementary Table 7), and the reverse splicing technology was used in circ RNA primer design (Table 2). Their expression in the peripheral blood of the same patient (the one who provided tissue samples for RNA-Seq) was verified (n = 4). The results showed that hsa_circ_0008433 and hsa_circ_0005571 were both highly expressed in IA tissues (p < 0.05), but the expression of them in peripheral blood was low and high, respectively (Figures 1A,B). Similar to hsa_circ_0008433, hsa_circ_0033144 was highly expressed in IA tissues and low in peripheral blood of the two groups (Figure 1C, p < 0.05). Hsa_circ_0040809 was highly expressed in IA tissues and low in peripheral blood, but the differences were only statistically significant in the tissue expression (Figure 1D, p < 0.05). Hsa_circ_0056285 was highly expressed in IA tissues and low in peripheral blood, but the differences were not statistically significant (Figure 1E, p > 0.05). Hsa_circ_0007142 and hsa_circ_0072309 showed low expression both in IA tissues and peripheral blood, the differences between cases and control group of the two samples were statistically significant (Figures 1F,G, p < 0.05).

TABLE 2

Table 2. CircRNA and internal reference gene primer sequences.

FIGURE 1

Figure 1. qRT-PCR verification of specific circRNAs in IA tissue and peripheral blood of case and control group. (A–G) The PCR results of the 7 selected circRNAs hsa_circ_ (0008433, 0005571, 0033144, 0040809, 0056285, 0007142, 0072309) in IA tissue and peripheral blood, *P < 0.05.

To further study the potential association of circRNA expression in tissue and peripheral blood, the expression of the above-mentioned specific IA rupture-related circRNAs in the same patient's IA tissues and peripheral blood was described (Figures 2A–G). Pearson correlation analysis was used to study the expression patterns and correlations of these specific circRNAs in two samples of the same patient. The results showed that the expression of hsa_circ_0008433 (r = −0.778, 95% CI: −0.958 to −0.163) and hsa_circ_0033144 (r = −0.749, 95% CI: −0.951 to −0.093) in IA tissues was negatively correlated with that in peripheral blood (Figures 2A,C). The expression of hsa_circ_0007142 (r = 0.926, 95% CI: 0.635~0.987) in IA tissues was positively correlated with that in peripheral blood (Figure 2F).

FIGURE 2

Figure 2. Correlation between specific circRNAs in IA tissue and blood. (A–G) Correlations between 7 selected circRNAs hsa_circ_ (0008433, 0005571, 0033144, 0040809, 0056285, 0007142, 0072309) in each IA and STA sample from tissue and peripheral blood, *P < 0.05.

Detection of Peripheral Blood PCR in Small Sample Population

Screening circRNA Indicators for PCR Detection

Based on the RNA-Seq results and combined with the preliminary experiments (21), we further selected 4 of the 7 circRNAs verified above for PCR detection in the population peripheral blood. Among them, the expression differences of hsa_circ_0005571, hsa_circ_0033144, and hsa_circ_0007142 in case and control group of IA tissues and peripheral blood were statistically significant (p < 0.05); hsa_circ_0008433 had a large differential expression multiple in RNA-Seq (FC = 52.077). At the same time, hsa_circ_0001946 (circ-CDR1as) and hsa_circ_0000284 (circ-HIPK3) were selected by reviewing the literature, they are two circRNAs closely related to cerebrovascular diseases that are currently being noticed (Table 3) (22, 23).

TABLE 3

Table 3. CircRNAs for PCR in peripheral blood of small sample population.

qRT-PCR Detection of circRNA

Peripheral blood of 30 patients with IA rupture and 30 healthy people were collected for qRt-PCR. The results showed that the expressions of hsa_circ_0008433, hsa_circ_0005571, hsa_circ_0007142, hsa_circ_0001946 in peripheral blood of patients with IA rupture were lower than that of the control group, and the difference was statistically significant (p < 0.05) (Figure 3).

FIGURE 3

Figure 3. Expression of specific circRNAs in peripheral blood. (A–F) Population qRT-PCR verification of hsa_circ_ (0008433, 0005571, 0007142, 0033144, 0001946, 0000284) expression in peripheral blood of IA and control group, IA group (n = 30) vs. control group (n = 30), *P < 0.05.

Detection of Peripheral Blood PCR in Expanded Sample Population

Study Population and General Demographic Characteristics

The peripheral blood of 140 patients with IA rupture and 140 healthy people were further collected for PCR verification in expanded sample population. The general demographic characteristics of the case and the control group were balanced and comparable (p > 0.05) (Supplementary Table 8).

qRT-PCR Detection of circRNA

Select 4 IA rupture-related circRNA indicators (hsa_circ_0008433, hsa_circ_0005571, hsa_circ_0007142, hsa_circ_0007142, hsa_circ_0001946) with significant expression differences in small sample population peripheral blood for the next verification in an expanded sample population. The results showed that compared with the control group, 75.0% (105/140), 65.71% (92/140), 55.71% (78/140), and 67.14% (94/140) of patients showed down-regulation of the four circRNAs in peripheral blood (Figure 4).

FIGURE 4

Figure 4. General expression of specific circRNAs in blood of population. (A–D) General expression of hsa_circ_ (0008433, 0005571, 0007142, 0001946) in blood of population.

The results of qRT-PCR showed that compared with the control group, the expression of hsa_circ_0008433, hsa_circ_0005571, and hsa_circ_0001946 in the peripheral blood of patients with IA were down-regulated, and the difference was statistically significant (Figures 5A–D, p < 0.05).

FIGURE 5

Figure 5. Expression of specific circRNAs in peripheral blood. (A–D) Population qRT-PCR verification of hsa_circ_ (0008433, 0005571, 0001946, 0007142) expression in peripheral blood of IA and control group, IA group (n = 140) vs. control group (n = 140), *P < 0.05.

Association of circRNA and IA Rupture

A logistic regression model was established to analyze whether these specific circRNAs were independent influencing factors for IA rupture. The results showed that after adjustment of the main individual environmental factors, has_circ_0008433 (OR = 0.497, 95% CI: 0.338~0.731) and has_circ_0001946 (OR = 0.682, 95% CI: 0.509~0.914) were independent influencing factors of IA rupture. High-level expression of has_circ_0008433 and has_circ_0001946 in peripheral blood may have a protective effect on IA rupture (Table 4).

TABLE 4

Table 4. Multivariate analysis of circRNA and IA rupture.

ROC Curve

ROC curve is used to explore the significance of specific circRNA in the diagnosis of IA rupture. The results show that the area under the ROC curve (AUC) of has_circ_0008433 is 0.7028, the maximum Youden index is 0.3000, and the corresponding cutoff is −1.9051 (Figure 6A); The AUC of has_circ_0005571 is 0.617, the maximum Youden index is 0.1928, and the corresponding cutoff is −0.5889 (Figure 6B); The AUC of has_circ_0001946 is 0.664, the maximum Youden index is 0.2571, and the corresponding cutoff is −1.6451 (Figure 6C, Supplementary Table 9).

FIGURE 6

Figure 6. ROC curve of IA related circRNAs and combined diagnosis. (A–D) ROC curve of hsa_circ_ (0008433, 0005571, 0001946) and combined factor. AUC, Area under ROC curve.

The three circRNAs were further adopted into the logistic regression equation to construct new combined diagnostic factors and draw ROC curves. The results showed that the AUC of the new combined diagnostic factor was 0.726 (95% CI: 0.668~0.784), the sensitivity and specificity were 0.793 and 0.564, respectively. Further comparison test results of the ROC curves of these circRNA and new combined factors showed that, hsa_circ_0008433 has higher diagnostic recognition than has_circ_0005571, and when the three circRNAs were jointly predicted, the diagnostic value is better than that of circ_0005571 or circ_0001946 single index (p < 0.05) (Figure 6, Supplementary Table 10).

The Combined Effects of Epigenetics and Other Influencing Factors

has_circ_0008433, has_circ_0005571, and has_circ_0001946 were further adopted in logistic regression model, the case and the control group were divided into high and low risk score according to the median epigenetic risk score of the control group, and the epigenetic IA risk score was constructed (the higher the score, the greater the risk of IA rupture), so as to explore the joint effects of circRNA epigenetic risk and age. The results show that patients with high epigenetic risk and older (≥55 years old) have an increased risk of IA rupture (p < 0.05), which indicated that there is a certain multiplicative interaction between epigenetic factors and age (Table 5).

TABLE 5

Table 5. Combined effects of epigenetic risk and age on IA.

Discussion

Current studies suggest that IA is a complex disease with multiple environmental factors and genetic regulation (24). This study indicated that people with smoking, chemical poison exposure (hair dye), sitting time >6 h/day, single animal oil intake, larger diastolic pressure and pulse pressure differences, higher levels of plasma globulin have higher risk of IA rupture. On the contrary, people with tea drinking habits, higher education, sleep time >7 h/day, diabetes and higher levels of (hemoglobin, low density lipoprotein, apolipoprotein A1, and serum calcium) have a lower risk of IA rupture. After controlling these major individual environmental factors, the low expression of circular RNA has_circ_0008433 and has_circ_0001946 in peripheral blood is the independent epigenetic risk factors for IA rupture. At the same time, the epigenetic risk score constructed based on these two circular RNAs interacts with age and has a certain clinical diagnostic values for IA repture.

It is worth mentioning that, in addition to the traditionally recognized influencing factors (smoking, drinking, biochemical indicators changes, etc.), this study also has some new findings. For example, people with higher diastolic pressure and pulse pressure differences are susceptible to IA rupture. Our data shows that for every 1 mmHg increase in diastolic blood pressure, the risk of IA rupture increases by about 4.6% (OR = 1.046, 95% CI: 1.007~1.033). For every 1 mmHg increase in pulse pressure difference, the risk of IA rupture increases by about 2.1% (OR = 1.021, 5% CI: 1.007~1.033). Studies have shown that patients with simple diastolic hypertension have a higher incidence of cerebrovascular disease. In Asian and European populations, for every 5 mmHg decrease in diastolic blood pressure, the risk of stroke is reduced by 44 and 27%, respectively (25). At the same time, the increase in pulse pressure difference means that the stress load on the blood vessel wall increases, which will lead to changes in arterial structure and decreased compliance, that constitutes a local weakening factor of the blood vessel wall (26, 27). Therefore, more attention should be paid to the monitoring and management of diastolic blood pressure and pulse pressure difference in the control of IA rupture, which may have more positive significance for its prevention.

Sedentary activity has been listed by the American Heart Association (AHA) as one of the independent risk factors for cerebrovascular diseases (28). According to research data, the all-cause mortality rate of people who sit for more than 11 h a day will increase by 40% (29). Sedentary reduces muscle contraction and blood circulation, promotes an increase in arterial pressure. In addition, pathological changes such as oxidative stress and endothelial dysfunction induced by continuous high perfusion become an important pathological basis for the occurrence of vascular diseases (30). In our data, the OR of people with sitting time> 6 h/day suffering from IA rupture was 3.382 (95% CI: 1.262~9.062). Meanwhile, hair dye exposure is probably another bad lifestyle that may induce IA rupture. P-phenylenediamine (PPD) in hair dyes has a variety of target organ toxicity, which can induce TP53 gene mutation and cause DNA damage (31, 32). Although there is no direct report that the hair dye is related to the occurrence of IA, PPD as an environmental chemical poison and a new potential influencing factor of IA still deserve our attention.

Recent studies have shown that inflammation is the most important pathological mechanism in IA rupture, and circRNA may play an epigenetic regulatory role in the pathological mechanism of IA through various inflammatory pathways (33–35). including regulation of inflammatory gene expression, recruitment, chemotactic aggregation of macrophages and regulation of key signaling pathways, etc (36). In this study, based on high-throughput sequencing to analyze the circRNA expression profiles in IA tissues and their biological functions. The results showed that these abnormally expressed IA-related circRNAs are closely related to inflammation, and some important biological processes and cell signal pathways that associate with inflammation may be regulated by circular RNA epigenetic regulation.

Among them, the expressions of hsa_circ_0008433 and hsa_circ_0005571 in IA tissues and normal blood vessels are significantly different, and they still show such expression differences in the subsequent peripheral blood verification to expand the sample size. After adjusting for individual environmental factors that are important for the rupture of IA, the circular RNAs hsa_circ_0008433, hsa_circ_0005571, and hsa_circ_0001946 are independent epigenetic influencing factors for IA rupture, and have a certain synergistic interaction with age (>55 years). The results indicated that circRNA may participate in the pathological process of IA through various of mechanisms, which provide certain references for the molecular etiology of IA rupture (37–40).

At the same time, a “star molecule” of circular RNA circ-CDR1as (hsa_circ_0001946) which is widely studied currently was selected for peripheral blood PCR verification of IA rupture. The results showed that the expression of hsa_circ_0001946 in IA patients and healthy people was significantly different (p < 0.05). Studies have shown that circ-CDR1as can specifically bind miR-7-5p through ceRNA mechanism, regulate downstream target genes expression, and participate in the pathological mechanisms of various cerebrovascular diseases (22, 41). Genes such as PARP1 and RAF1, which are important downstream targets of miR-7-5p, are closely related to various pathophysiological mechanisms. Dysregulation expression of these genes can affect the proliferation and migration of important functional cells such as microvascular endothelial and vascular smooth muscle cells, causing damage and weakness of the vessel wall, and promoting the development of atherosclerosis (42).

circRNA has good structural stability and can be released into extracellular space through exosomes, so it can be detected in body fluids such as saliva, milk, and plasma (43), these characteristics make the exploration of circRNA as a biological marker of disease widely concerned and deeply studied. In this study, we verified the results of high-throughput sequencing of circRNA simultaneously in two human specimens of cerebral artery tissue and peripheral blood. The analysis showed that certain circRNAs are differentially expressed in ruptured IA tissues and peripheral blood, there is a certain correlation between their expression patterns. For example, compared with the control group, hsa_circ_0007142 was down-regulated in IA tissues and peripheral blood, hsa_circ_0005571 and hsa_circ_0033144 were up-regulated in IA tissues, and down-regulated in IA peripheral blood. In the same patient, the expression pattern of hsa_circ_0007142 in tissue and peripheral blood were positively correlated, and the expression patterns of hsa_circ_0008433 and hsa_circ_0033144 in tissue and peripheral blood were negatively correlated. However, there was no significant difference in the expression of hsa_circ_0033144 and hsa_circ_0007142 in the follow-up population validation, which may be related to the influence factors such as the small number of tissue samples included in RNA-Seq, so the verification of further expansion of the sample size is needed.

However, these circRNAs with certain expression correlation in small samples of IA tissue and peripheral blood are still worthy of attention. This suggests that circRNA may have some correlation expression characteristics in central tissue and peripheral blood, which provides a theoretical possibility for the use of peripheral blood to diagnose IAs. In addition, we screened the three circular RNAs (has_circ_0008433, has_circ_0005571, has_circ_0001946) indicators closely related to IA rupture, and confirmed that they have certain reference significance for IA discrimination (AUC = 0.726, 95% CI: 0.668~0.784). This also indicated that circRNA may have potential values as a biomarker, and provides a new perspective for the future research of IA non-invasive diagnosis and treatment strategies. In the study, we also found that specific circRNAs also interact with individual factors such as age, thereby further confirming that IA is a complex disease caused by the combination of environment and genetic factors.

The study also has certain limitations. First of all, we used a structured questionnaire to collect the data of the research subjects in the case-control study, and tried to reflect the clinical status of the patient as objectively as possible, but there were still some deviations inevitably. For example, the information distortion caused by the subjective will or memory bias of the respondents may have a certain impact on the research results. In the follow-up research, the sample size can be further expanded, and cross-regional cooperation can be strengthened to precisely control confounding factors and make the results more scientific. Secondly, this study does not yet have the conditions for multi-center research, so there are certain deficiencies in sample representativeness. Multi-center research needs to be carried out in follow-up studies to diversify the target population, thereby improving the accuracy and accuracy of the research. Finally, due to the limitations of case-control studies, it is not yet possible to accurately determine the causal relationship between circRNA and IA rupture, so further functional tests are needed to verify.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the Second Affiliated Hospital of Fujian Medical University (2018-50). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

SW, HL, and DK contributed to the study design. QinH, YS, QiuH, and YZ performed statistical analysis, interpretation, and drafted the manuscript. QiuH, YS, SL, YC, and XX contributed to data collection and laboratory test. QinH, YS, SW, and HL revised the manuscript. All authors contributed to critical revision of the final manuscript and approved the final version of the manuscript.

Funding

This research was funded by Joint Funds for the Innovation of Science and Technology, Fujian Province (2018Y9089), the Natural Science Foundation of Fujian Province (2019J01315), and Professor Development Fund Project of Fujian Medical University (JS15002).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank the participants and participating physicians from Second Affiliated Hospital of Fujian Medical University, China, as well as investigators and staff for making this research possible.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.594835/full#supplementary-material

Abbreviations

IA, intracranial aneurysm; circRNA, circular RNA; ceRNA, competitive endogenous RNA; CTA, computed tomography angiography; MRA, magnetic resonance angiography; DSA, digital subtraction angiography; OR, odds ratio; 95%CI, 95% confidence interval; AUC, area under the ROC curve.

References

1. Marbacher S, Wanderer S, Strange F, Grüter BE, Fandino J. Saccular aneurysm models featuring growth and rupture: a systematic review. Brain Sci. (2020) 10:101. doi: 10.3390/brainsci10020101

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Vlak MH, Algra A, Brandenburg R, Rinkel GJ. Prevalence of unruptured intracranial aneurysms, with emphasis on sex, age, comorbidity, country, and time period: a systematic review and meta-analysis. Lancet Neurol. (2011) 10:626–36. doi: 10.1016/S1474-4422(11)70109-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chen J, Li M, Zhu X, Chen L, Yang S, Zhang C, et al. Atorvastatin reduces cerebral vasospasm and infarction after aneurysmal subarachnoid hemorrhage in elderly Chinese adults. Aging. (2020) 12:2939–51. doi: 10.18632/aging.102788

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Krishnamurthi RV, Ikeda T, Feigin VL. Global, regional and country-specific burden of ischaemic stroke, intracerebral haemorrhage and subarachnoid haemorrhage: a systematic analysis of the global burden of disease study 2017. Neuroepidemiology. (2020) 54:171–9. doi: 10.1159/000506396

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Florez WA, García-Ballestas E, Deora H, Agrawal A, Martinez-Perez R, Galwankar S, et al. Intracranial hypertension in patients with aneurysmal subarachnoid hemorrhage: A systematic review and meta-analysis. Neurosurg Rev. (2020) 44:203–11. doi: 10.1007/s10143-020-01248-9

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Jiang Z, Chen Y, Zeng C, Feng J, Wan Y, Zhang X. Neurosurgical clipping versus endovascular coiling for patients with intracranial aneurysm: a systematic review and meta-analysis. World Neurosurg. (2020) 138:e191–22. doi: 10.1016/j.wneu.2020.02.091

CrossRef Full Text | Google Scholar

7. Pera J, Ruigrok YM. More evidence against alcohol or smoking in patients with unruptured intracranial aneurysm. Neurology. (2015) 84:442–3. doi: 10.1212/WNL.0000000000001222

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kernan WN, Viscoli CM, Brass LM, Broderick JP, Brott T, Feldmann E, et al. Phenylpropanolamine and the risk of hemorrhagic stroke. N Engl J Med. (2000) 343:1826–32. doi: 10.1056/NEJM200012213432501

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wermer MJ, van der Schaaf IC, Algra A, Rinkel GJ. Risk of rupture of unruptured intracranial aneurysms in relation to patient and aneurysm characteristics: an updated meta-analysis. Stroke. (2007) 38:1404–10. doi: 10.1161/01.STR.0000260955.51401.cd

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gross BA, Rosalind Lai PM, Frerichs KU, Du R. Aspirin and aneurysmal subarachnoid hemorrhage. World Neurosurg. (2014) 82:1127–30. doi: 10.1016/j.wneu.2013.03.072

CrossRef Full Text | Google Scholar

11. Tirschwell DL, Smith NL, Heckbert SR, Lemaitre RN, Longstreth WT Jr, Psaty BM. Association of cholesterol with stroke risk varies in stroke subtypes and patient subgroups. Neurology. (2004) 63:1868–75. doi: 10.1212/01.WNL.0000144282.42222.DA

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Tokuda Y, Stein GH. Serum lipids as protective factors for subarachnoid hemorrhage. J Clin Neurosci. (2005) 12:538–41. doi: 10.1016/j.jocn.2004.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang H, Yang J, Yang J, Fan Z, Yang C. Circular RNAs: novel rising stars in cardiovascular disease research. Int J Cardiol. (2016) 202:726–7. doi: 10.1016/j.ijcard.2015.10.051

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Rybak-Wolf A, Stottmeister C, GlaŽar P, Jens M, Pino N, Giusti S, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. (2015) 58:870–85. doi: 10.1016/j.molcel.2015.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhang Y, Liang W, Zhang P, Chen J, Qian H, Zhang X, et al. Circular RNAs: emerging cancer biomarkers and targets. J Exp Clin Cancer Res. (2017) 36:152. doi: 10.1186/s13046-017-0624-z

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang Y, Liang W, Zhang P, Chen J, Qian H, Zhang X, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. (2013) 495:333–8. doi: 10.1038/nature11928

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhang Y, Xue W, Li X, Zhang J, Chen S, Zhang JL, et al. The biogenesis of nascent circular RNAs. Cell Rep. (2016) 15:611–24. doi: 10.1016/j.celrep.2016.03.058

PubMed Abstract | CrossRef Full Text | Google Scholar

18. World Health Organization. International Guide for Monitoring Alcohol Consumption and Related Harm[M]. World Health Organization (2000).

Google Scholar

19. Gottlieb DJ, Punjabi NM, Newman AB, Resnick HE, Redline S, Baldwin CM, et al. Association of sleep time with diabetes mellitus and impaired glucose tolerance. Arch Intern Med. (2005) 165:863–7. doi: 10.1001/archinte.165.8.863

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bankoski A, Harris TB, McClain JJ, Brychta RJ, Caserotti P, Chen KY, et al. Sedentary activity associated with metabolic syndrome independent of physical activity. Diabetes Care. (2011) 34:497–503. doi: 10.2337/dc10-0987

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Huang Q, Huang QY, Sun Y, Wu S. High-throughput data reveals novel circular RNAs via competitive endogenous RNA networks associated with human intracranial aneurysms. Med Sci Monit. (2019) 25:4819–30. doi: 10.12659/MSM.917081

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Xu B, Yang T, Wang Z, Zhang Y, Liu S, Shen M. CircRNA CDR1as/miR-7 signals promote tumor growth of osteosarcoma with a potential therapeutic and diagnostic value. Cancer Manag Res. (2018) 10:4871–80. doi: 10.2147/CMAR.S178213

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Shan K, Liu C, Liu BH, Chen X, Dong R, Liu X, et al. Circular noncoding RNA HIPK3 mediates retinal vascular dysfunction in diabetes mellitus. Circulation. (2017) 136:1629–42. doi: 10.1161/CIRCULATIONAHA.117.029004

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Tromp G, Weinsheimer S, Ronkainen A, Kuivaniemi H. Molecular basis and genetic predisposition to intracranial aneurysm. Ann Med. (2014) 46:597–606. doi: 10.3109/07853890.2014.949299

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li Y, Wei FF, Wang S, Cheng YB, Wang JG. Cardiovascular risks associated with diastolic blood pressure and isolated diastolic hypertension. Curr Hypertens Rep. (2014) 16:489. doi: 10.1007/s11906-014-0489-x

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zakopoulos NA, Lekakis JP, Papamichael CM, Toumanidis ST, Kanakakis JE, Kostandonis D, et al. Pulse pressure in normotensives: a marker of cardiovascular disease. Am J Hypertens. (2001) 14:195–9. doi: 10.1016/S0895-7061(00)01268-1

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kim KM, Gwak MS, Choi SJ, Kim MH, Park MH, Heo BY. Pulse pressure variation and stroke volume variation to predict fluid responsiveness in patients undergoing carotid endarterectomy. Korean J Anesthesiol. (2013) 65:237–43. doi: 10.4097/kjae.2013.65.3.237

PubMed Abstract | CrossRef Full Text | Google Scholar

28. de Rezende LF, Rodrigues Lopes M, Rey-López JP, Matsudo VK, Luiz Odo C. Sedentary behavior and health outcomes: an overview of systematic reviews. PLoS ONE. (2014) 9:e105620. doi: 10.1371/journal.pone.0105620

PubMed Abstract | CrossRef Full Text | Google Scholar

29. van der Ploeg HP, Chey T, Korda RJ, Banks E, Bauman A. Sitting time and all-cause mortality risk in 222497 Australian adults. Arch Intern Med. (2012) 172:494–500. doi: 10.1001/archinternmed.2011.2174

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Thosar SS, Johnson BD, Johnston JD, Wallace JP. Sitting and endothelial dysfunction: the role of shear stress. Med Sci Monit. (2012) 18:173–80. doi: 10.12659/MSM.883589

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Elevli M, Civilibal M, Ersoy O, Demirkol D, Gedik AH. Paraphenylene diamine hair dye poisoning: an uncommon cause of rhabdomyolysis. Indian J Pediatr. (2014) 81:709–11. doi: 10.1007/s12098-013-1074-z

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Huang YC, Hung WC, Kang WY, Chen WT, Chai CY. p-Phenylenediamine induced DNA damage in SV-40 immortalized human uroepithelial cells and expression of mutant p53 and COX-2 proteins. Toxicol Lett. (2007) 170:116–23. doi: 10.1016/j.toxlet.2007.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

33. van Rossum D, Verheijen BM, Pasterkamp RJ. Circular RNAs: novel regulators of neuronal development. Front Mol Neurosci. (2016) 9:74. doi: 10.3389/fnmol.2016.00074

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Pawlowska E, Szczepanska J, Wisniewski K, Tokarz P, Jaskólski DJ, Blasiak J. NF-kappaB-Mediated inflammation in the pathogenesis of intracranial aneurysm and subarachnoid hemorrhage. does autophagy play a role? Int J Mol Sci. (2018) 19:1245. doi: 10.3390/ijms19041245

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Dong Y, Fan C, Hu W, Jiang S, Ma Z, Yan X, et al. Melatonin attenuated early brain injury induced by subarachnoid hemorrhage via regulating NLRP3 inflammasome and apoptosis signaling. J Pineal Res. (2016) 60:253–62. doi: 10.1111/jpi.12300

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Yu L, Wang J, Wang S, Zhang D, Zhao Y, Wang R, et al. DNA methylation regulates gene expression in intracranial aneurysms. World Neurosurg. (2017) 105:28–36. doi: 10.1016/j.wneu.2017.04.064

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Giang S, La Cava A. IRF1 and BATF: key drivers of type 1 regulatory T-cell differentiation. Cell Mol Immunol. (2017) 14:652–4. doi: 10.1038/cmi.2017.38

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wen C, Xu M, Mo C, Cheng Z, Guo Q, Zhu X. JMJD6 exerts function in neuropathic pain by regulating NFkappaB following peripheral nerve injury in rats. Int J Mol Med. (2018) 42:633–42. doi: 10.3892/ijmm.2018.3613

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liu Y, Yan X, Zhou T. TBCK influences cell proliferation, cell size and mTOR signaling pathway. PLoS ONE. (2013) 8:e71349. doi: 10.1371/journal.pone.0071349

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yan L, Zhu YQ, Li MH, Tan HQ, Cheng YS. Geometric, hemodynamic, and pathological study of a distal internal carotid artery aneurysm model in dogs. Stroke. (2013) 44:2926–9. doi: 10.1161/STROKEAHA.113.002290

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Piwecka M, GlaŽar P, Hernandez-Miranda LR, Memczak S, Wolf SA, Rybak-Wolf A, et al. Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science. (2017) 357:8526. doi: 10.1126/science.aam8526

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Faltz M, Bergin H, Pilavachi E, Grimwade G, Mabley JG. Effect of the anti-retroviral drugs efavirenz, tenofovir and emtricitabine on endothelial cell function: role of PARP. Cardiovasc Toxicol. (2017) 17:393–404. doi: 10.1007/s12012-016-9397-4

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Li Y, Zheng Q, Bao C, Li S, Guo W, Zhao J, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. (2015) 25:981–4. doi: 10.1038/cr.2015.82

PubMed Abstract | CrossRef Full Text | Google Scholar