<figure><img src="image_1.png" /></figure>

OPEN

# Prediction of coronary heart disease incidence in a general male population by circulating non-coding small RNA sRNY1-5p in a nested case-control study

Vera L. Costa<fnref n="1" />, Jean-Bernard Ruidavets², Vanina Bongard²,⁵, Bertrand Perret³,⁴,⁶, Emanuela Repetto<fnref n="1" />, Maria G. Stathopoulou<fnref n="1" />, Fabrizio Serra<fnref n="1" />, Mohamed Benahmed<fnref n="1" />, Claire Mauduit<fnref n="1" />, Valerie Grandjean<fnref n="1" />, Jean Ferrières²,⁵, Laurent O. Martinez²,³,⁷& Michele Trabucchi<fnref n="1" />,⁷&

During the development of atherosclerotic lesion, s-RNYs (small RNAs of about 24/34 nucleotides) are derived by the processing of long Ro-associated non-coding RNAs (RNYs) in macrophages. The levels of serum s-RNYs have been found significantly upregulated in patients with coronary heart disease (CHD) compared to age-matched CHD-free individuals. The present study aimed to examine the predictive value of serum s-RNYs for CHD events in the general male population. Within the frame of nested-case-control study, the GENES study, we measured the absolute expression of a RNY-derived small RNA, the s-RNY1-5p, in the serum of individuals (without CHD at baseline) who encountered a CHD event within 12 years of follow-up (n = 30) (Cases) and compared them to individuals who remained event-free (Controls) (n = 30). The expression of s-RNY1-5p in serum was significantly upregulated in Cases compared to Controls (p = 0.027). The proportion of CHD event-free was significantly higher among individuals with serum s-RNY1-5p below the median value (631 molecules/ mL). In a multivariable model adjusted for age, smoking, hypertension, diabetes and dyslipidemia, the risk of CHD events increased more than fourfold in individuals with serum s-RNY1-5p above the median value (HR, 4.36; 95% CI 1.22–15.60). A positive association with CHD events was also observed when considering s-RNY1-5p as a continuous variable (p = 0.022). Based on our results, we conclude that serum s-RNY1-5p is an independent predictor of CHD events in a general male population and might be a relevant biomarker for early detection of cardiovascular diseases.

Atherosclerosis is a progressive disease process that generally begins in childhood and manifests clinically in middle to late adulthood. Biological, clinical, genetic and lifestyle factors represent the major causes of lesion formation, including high level of LDL-C (Low Density Lipoprotein-Cholesterol), low level of HDL-C (High Density Lipoprotein-Cholesterol), high blood pressure, smoking, unhealthy diet, old age, and family history of heart disease. Atherosclerosis starts with the accumulation of atherogenic lipids, including cholesterol from LDL, in the arterial wall, which promotes an inflammatory response by endothelial cells and monocytes/macrophages. In particular, macrophages play a major role in atherosclerosis progression by regulating the inflammatory status of the lesion<fnref n="1" />. In advanced lesions, massive macrophage apoptosis that is not cleared by phagocytic clearance leads to the necrotic plaque formation. The build-up of necrotic plaque ultimately promotes inflammation, lesion

<footnote n="1">¹Inserm U1065, C3M, Team Control of Gene Expression (10), Université Côte D'Azur, Nice, France. ²Department of Epidemiology, Health Economics and Public Health, UMR1027 INSERM, Toulouse University, Toulouse University Hospital (CHU), Toulouse, France. ³Institut National de La Santé Et de La Recherche Médicale (INSERM), Institute of Cardiovascular and Metabolic Diseases (I2MC), UMR1297, BP 84225, 31432 Toulouse Cedex 04, France. ⁴University of Toulouse, UMR1297, Paul Sabatier University, Toulouse, France. ⁵Fédération de Cardiologie, Toulouse University Hospital (CHU), Toulouse, France. ⁶Service de Biochimie, Pôle biologie, Hôpital de Purpan, CHU de Toulouse, Toulouse, France. ⁷These authors contributed equally: Laurent O. Martinez and Michele Trabucchi. &email: Laurent.Martinez@inserm.fr; michele.trabucchi@univ-cotedazur.fr</footnote>
disruption, and eventually thrombosis and cardiovascular accident¹. When this cardiovascular accident occurs in heart arteries it is called coronary heart disease (CHD). Therefore, monitoring the degree of macrophage apoptosis in the blood vessels might be a valuable parameter to predict cardiovascular accidents.

Recently, we found that small RNAs of about 24/34 nucleotides (nt) in length derived by the processing of the non-coding YRNAs (RNYs) are significantly upregulated in the serum of patients with CHD compared to age-matched individuals without CHD². We have referred to these small RNAs as s-RNYs. RNY genes account for four copies in human genome (RNY1, 3, 4 and 5). Analyses in a cross-sectional study on CHD, comparing 263 CHD-patients with paired controls, positioned s-RNYs as relevant novel independent diagnostic biomarkers for CHD, which are associated with atherosclerosis burden². In the context of atherosclerotic lesion and development, serum s-RNYs would be mainly produced by macrophages². RNYs are characterized by extensive base pairing of the 5' and 3' regions and by the association with the proteins Ro60 and La/SSB, which are often targeted by the immune system in several autoimmune diseases, including the systemic lupus erythematosus and Sjogren's syndrome³. Circulating s-RNYs are associated with Ro60 protein to form a stable ribonucleoprotein complex that promotes inflammation and apoptosis of monocytes/macrophages by activating the Toll-like receptor (TLR) 7⁴, suggesting that upregulation of serum s-RNYs promotes the atherosclerosis progression.

Here, we tested the hypothesis that measurement of serum s-RNY levels could be useful in the prediction of CHD events for primary prevention of cardiovascular diseases. To explore this hypothesis, we measured the absolute expression of an sRNYs, the s-RNY1-5p, in the serum of individuals (all CHD-free at baseline) who developed CHD during the 12 years median follow-up from the blood draw and compared them to control individuals who never encountered CHD in the same period. Among the different RNY-derived small RNAs previously detected in human serum², s-RNY1-5p was selected to study the potential value of RNY-derived small RNAs to predict cardiovascular accidents. In fact, s-RNY1-5p is the most characterized RNY-derived small RNA as marker of CHD status, and like microRNAs, is very stable in the serum².

## Methods

**Sampling frame.** The "Genetique et Environnement en Europe du Sud" (GENES) study is a case–control study designed to assess the role of gene-environment interactions in the occurrence of CHD⁵. All participants signed an informed consent form. Biological samples and clinical/physical information were collected according to the principles expressed in the Declaration of Helsinki. The study protocol was approved by the local ethics committee (CCPPRB, Toulouse/Sud-Ouest, file #1-99-48, Feb 2000). A collection of biological samples has been constituted (declared as DC-2008-463 #1 to the Ministry of Research and to the regional Health authority)⁶. Control individuals analyzed in the frame of the present work were part of the 834 male individuals from the GENES study, free of CHD and stroke, and aged 45–75 at baseline, randomly selected from the general population using electoral rolls⁵. Those control individuals were followed-up and provided the cases and controls for the nested case–control study of the present work.

Finally, the whole study was performed according to STROBE guidelines (http://www.strobe-statement.org).

**Follow-up and case ascertainment.** Control individuals from the GENES study were contacted during the year 2016 either by mail or by telephone and completed self-report questionnaires on clinical events. Among individuals with possible events, clinical information was obtained from medical records in hospital and general practitioners. All available data were gathered concerning hospital admission and type of interventions. Deceased patients' families and practitioners were contacted to clarify death circumstances, and death certificates were checked to complement clinical and postmortem information. Coronary events (unstable angina, coronary revascularization, non-fatal myocardial infarction), stroke events and the occurrence of peripheral vascular disease were adjudicated by an independent committee. The median follow-up period was 11.9 ± 4.1 years.

**Nested case–control study.** In the present work, case patients (n = 31) were individuals who were CHD-free at baseline and developed either fatal or non-fatal acute coronary events during the follow-up. The control subjects (n = 31), matched for age (± 2 years) and date of recruitment (± 3 months), were participants who remained free of cardiovascular diseases and stroke during follow-up. One Control failed for quantification of s-RNY1-5p in baseline serum and was thus excluded with its paired Case. Analyses on s-RNY1-5p were thus based on a nested case–control population of 30 case patients and 30 control participants from the GENES study. The study workflow is illustrated in Fig. 1.

**Baseline assessment.** At the baseline survey, between 2001 and 2004, participants were interviewed and examined in the same health center administered by the national health insurance system. Age, environmental characteristics and information on cardiovascular risk factors were collected through standardized face-to-face interviews, performed by a single physician. Smoking status was classified as current smokers, smokers having quitted for more than 3 years and non-smokers. Alcohol consumption was assessed using a typical week pattern. The total amount of pure alcohol consumption was calculated as the sum of different types of drinks and was expressed as grams per day. Physical activity was investigated through a standardized questionnaire⁷ and categorized into two levels, namely, “high” as physical activity for at least 20 min and at least twice a week or “low” physical activity for less than 20 min for once a week or less. Medical history was collected. Hypertension was defined as treatment with drugs, or by systolic blood pressure ≥ 160 mmHg or diastolic blood pressure ≥ 95 mmHg. Diabetes was defined as treatment with drugs, or fasting blood glucose ≥ 7.8 mmol/L. Dyslipidemia was defined as treatment with drugs or fasting serum total cholesterol ≥ 2.40 g/L⁸.

Anthropometric measurements included waist circumference, height and body weight, and body mass index was calculated (BMI, kg/m²). Blood pressure and resting heart rate were measured with an automatic
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Study flowchart. Cases (n = 31) were male subjects from the GENES study (n = 834) who experienced a Coronary Heart Disease (CHD) events during the follow-up period (11.9 ± 4.1 years). The control subjects (n = 31), matched for age and date of recruitment, were participants who remained free of cardiovascular diseases and stroke during follow-up. One Control was excluded because quantification of s-RNY1-5p in baseline serum failed. Analyses were thus based on a nested case-control population of 30 case patients and 30 control participants from the GENES study.</figcaption></figure>

sphygmomanometer (OMRON 705 CP). Measurements were performed after a minimum of 5 min rest; average values from two different measurements were recorded for further analysis.

**Laboratory assays.** Blood was collected at inclusion after an overnight fast. Serum samples aliquots were subsequently stored at – 80 °C until biological analyses. The following biomarkers were assayed with enzymatic reagents on automated analyzers (Hitachi 912 and Cobas 8000, Roche Diagnostics, Meylan, France): serum total cholesterol (TC), HDL-C, triglycerides (TG), fasting glucose, and gamma-glutamyl transferase (γGT). High-sensitive C-reactive protein (hs-CRP) was determined on the same analyzer by immunoturbidimetry assays. Apolipoprotein (Apo) A1, ApoB, and lipoprotein (a) [Lp(a)] were also assayed with an immunoturbidimetric method on an automated analyzer (Roche Diagnostics, France). LDL-C was calculated using the Friedewald formula, with VLDL-cholesterol (VLDL-C) (g/L) = TG (g/L)/5, as long as TG concentration was below 4 g/L⁹.

**Absolute quantification of s-RNY1-5p in the serum by quantitative RT-PCR.** Absolute quantification of s-RNY1-5p was performed in baseline serum samples from all participants included in the nested case-control study. Before RNA extraction from 0.2 mL human serum, hemolysis was evaluated by measuring through spectrophotometer the oxyhemoglobin at 414 nm and referring to a standard curve made by introducing haemolysis by serially diluting lysed red blood cells in non-haemolysed plasma from a healthy donor. The serum samples showing an optical density (OD) corresponding to > 0.002% of red blood cells (v/v) were discarded. 0.5 mL of Trizol LS Reagent was added to 0.2 mL of blood serum. Then (1) 5 pg of synthetic miRNA-39 from *Caenorhabditis elegans* (miScript miRNA Mimic syn-cel-miR-39-5p, Qiagen) were added as a spike-in control for purification efficiency and (2) 2 µL of glycogen (5 mg/mL) were added to enhance the efficiency of RNA column binding. Purification of extracted total RNA was performed with High Pure miRNA Isolation kit (Roche) according to the manufacturer's instructions. RNA was eluted in 50 µL of Elution Buffer. The quality of the extracted RNA was checked by ratio between the absorbance values at 260 and 280 nm and between 260 and 230 nm using NanoDrop Technologies ND-1000 spectrophotometer (Supplemental Table 1).

Five-fold dilutions of synthetic s-RNY1-5p ranging from 10⁹ to 10⁵ molecules per RT reaction were prepared freshly in Ultrapure water (Invitrogen) before each experiment to create standard curves to use in each quantitative PCR (qPCR) plate. Reverse transcriptase (RT) reaction was performed according to Repetto et al.² for the
detection of s-RNY1-5p and with TaqMan (Life Technologies) for the cel-miR-39. For s-RNY1-5p detection, RT was performed with the GoScript Reverse Transcriptase Kit (Promega), following the manufacturer's instructions. Briefly, reagents were freshly mixed before each experiment in a total volume of 20 µL, containing 0.5 mM dNTPs, 0.5 µL reverse transcriptase, 4 µL 5× buffer, 20 units RNase inhibitor, 6.6 µL nuclease-free H2O, 2 µM RT stem-loop primer, and 4 µL of serum RNA or the synthetic s-RNY1-5p. Every step was performed on ice. The reaction was performed at 25 °C for 5 min followed by 42 °C for 60 min and 70 °C for 15 min before being held at 4 °C. For the synthetic cel-miR-39 detection, RT was performed with the TaqMan microRNA Reverse Transcription Kit (Applied Biosystems), using 1 µL of serum RNA and following the manufacturer's instructions. The resulting cDNAs were stored at -80 °C, till use for the qPCR experiments.

qPCR amplifications were performed in a total volume of 10 µL, containing 2 µL of freshly prepared RT reaction, 1 µL of a mix containing the specific s-RNY1-5p oligonucleotide primer at 1 µM and the universal oligonucleotide primer at 1 µM (see below for the sequence of the oligonucleotides), 2 µL of nuclease-free H2O, and 5 µL of Fast SYBR Green qPCR Master Mix (Applied Biosystems). Cel-miR-39 detection was performed according to the manufacturer's instructions (Applied Biosystems) in a final volume of 10 µL. Expression was considered undetectable with threshold cycle (Ct) value ≥ 40. Each sample was run in technical triplicate and a RT negative control without serum RNA was added in each qPCR plate. In addition, in each RT-qPCR plate, we performed s-RNY1-5p standard curve. For each sample, both cel-miR-39 and s-RNY1-5p RT-qPCR reactions were carried out (Supplemental Tables 2 and 3).

The cel-miR-39 spike-in measured by RT-qPCR did not vary in all samples (Supplemental Table 2), suggesting an equal yield of RNA purification. These results allowed us to normalize the s-RNY1-5p qPCR data by volume of serum (Supplemental Table 3). Copy number quantification in each sample relied on the standard curve method (Supplemental Fig. 1). Standard curve was generated by the StepOne software (ThermoFisher Scientific, MA) based on the concentration of synthetic s-RNY1-5p (Supplemental Fig. 1). The automatic threshold function was applied to determine the Ct for each sample and calculate the absolute number of copies. Data were normalized to the initial serum volume (Supplemental Table 3). None value was identified as outlier with the ROUT method (Q=0.1%). Primer sequences used in this study are shown in the Supplemental Table 4.

**Statistical analysis.** Categorical (qualitative) variables are presented as proportions. Continuous variables are displayed as means and standard deviations (SD). Kurtosis and skewness tests were used to test for the assumption of normality and log-transformation was performed for some variables to improve the normality of the distribution. We first described and compared characteristics of participants according to the case-control status. Categorical variables were compared between groups using the McNemar's test. Student's t-test was used to compare the distributions of normally distributed continuous data. Wilcoxon signed rank sum test was used when distribution departed from normality, or when homoscedasticity was rejected.

Cumulative risk of patients for coronary event (fatal and non-fatal) was determined by the Kaplan-Meier method and compared using the Log-rank test for the individual endpoints of coronary events.

The relation between baseline levels of s-RNY1-5p and the case-control status was assessed using Cox proportional hazards regression analysis. Cox regression analyses were performed with adjustment on classical cardiovascular risk factors (age, smoking, dyslipidemia, hypertension and diabetes as well as for systolic and diastolic blood pressure, total cholesterol and HDL-C). We tested the proportionality assumption using cumulative sums of martingale-based residuals. We performed regression analysis with polynomial models (quadratic and cubic) to examine for possible non-linear relationship between s-RNY1-5p and occurrence of coronary event.

All statistical analyses were carried out using the SAS statistical software package 9.4 (SAS Institute, Cary, NC) or STATA statistical software, release 14.2 (STATA Corp, College Station, TX). All tests were two-sided and were considered significant at a p value < 0.05.

## Results

The Nested Case-Control study we used here, comprised 62 men from the GENES study, CHD-free at baseline. The median follow-up period was 11.9 ± 4.1 years.

Comparison of individual's data at the recruitment in the cohort is given in Table 1. Results are presented distinguishing two groups: the "Controls group" (individuals with no coronary event during the follow-up, n = 31) and the "Cases group" (individuals who developed an acute coronary event during the follow-up, n = 31). During the follow-up period, half of the Cases (n = 16) had a fatal event while a non-fatal acute coronary syndrome such as unstable angina and non-fatal myocardial infarction was reported for the other half (n = 15). There were no significant differences between biochemical and physiological variables among the Cases and Controls, except for TG levels that were significantly higher in Cases (p = 0.03). No significant difference between Cases and Controls was noted with regards to the presence of diabetes and dyslipidemia at baseline. The proportion of Cases diagnosed with hypertension tended to be higher than in Controls (35.5% versus 19.4%, p = 0.14). Similar observations were made according to the treatment status. Overall, the present cohort population was representative of the general male population of this age range in the early 2000s in Southern France in terms of biological variables and treatments¹⁰.

One Control and its paired Case were excluded from sRNY1-5p analyses (see "Methods" section). As shown in Fig. 2 and Table 2, serum levels of s-RNY1-5p at baseline, expressed in number of molecules per mL, were significantly higher in Cases (Fig. 2 and Table 2). Accordingly, the percentage of individuals with s-RNY1-5p above the median value (631 molecules/mL) was also significantly higher in Cases compared to Controls (70.0% versus 33.3%, p = 0.022). The relation between baseline levels of s-RNY1-5p and the case-control status was assessed using Cox proportional hazards regression analysis. As reported in Table 3, in a multivariable model adjusted for age, smoking, hypertension, diabetes and dyslipidemia (Model 1), the risk of CHD event
<table><caption>**Table 1.** Clinical and biological characteristics at baseline. Variables were measured at baseline (i.e. when individuals were first included in the GENES cohort). For continuous variables, values are expressed as mean ± SD. For categorical variables, values are expressed as number with frequency (%) in parentheses. BMI body mass index, hs-CRP high-sensitivity C-reactive protein, eGFR estimated glomerular filtration rate, γGT gamma glutamyl transferase, SBP systolic blood pressure, DBP diastolic blood pressure. <sup>a</sup>Paired Student's t-test, unless otherwise stated. <sup>b</sup>Wilcoxon signed rank sum test. <sup>c</sup>McNemar's test. <sup>d</sup>“High” physical activity during 20 min at least twice a week versus “low” physical activity once a week or less. <sup>e</sup>Analyses performed on log transformed data. <sup>f</sup>Hypertension: systolic blood pressure ≥ 160 mmHg or diastolic blood pressure ≥ 95 mmHg or treatment; Diabetes: fasting blood glucose ≥ 7.8 mmol/L or treatment; Dyslipidemia: total cholesterol ≥ 2.50 g/L or treatment.</caption><thead><tr><th></th><th>Controls (n = 31)</th><th>Cases (n = 31)</th><th>p<sup>a</sup></th></tr></thead><tbody><tr><td>Age, year</td><td>63.0 ± 7.8</td><td>63.8 ± 8.3</td><td>0.09<sup>b</sup></td></tr><tr><td>Alcohol, g/day</td><td>18 ± 18</td><td>18 ± 21</td><td>0.97</td></tr><tr><td>Former and current smoker, n (%)</td><td>17 (54.8)</td><td>21 (67.7)</td><td>0.29<sup>c</sup></td></tr><tr><td>Physical activity, high level<sup>d</sup>, n (%)</td><td>9 (29.0)</td><td>10 (32.3)</td><td>0.80<sup>c</sup></td></tr><tr><td>Waist circumference, cm</td><td>96.4 ± 7.7</td><td>98.3 ± 10.4</td><td>0.44</td></tr><tr><td>BMI, kg/m<sup>2</sup></td><td>26.7 ± 2.9</td><td>27.3 ± 3.6</td><td>0.48</td></tr><tr><td>Fasting glucose, mmol/L</td><td>5.25 ± 0.89</td><td>5.59 ± 0.84</td><td>0.16</td></tr><tr><td>Total cholesterol, g/L</td><td>2.18 ± 0.31</td><td>2.29 ± 0.43</td><td>0.25</td></tr><tr><td>HDL-C, g/L</td><td>0.54 ± 0.12</td><td>0.52 ± 0.10</td><td>0.47</td></tr><tr><td>LDL-C, g/L</td><td>1.43 ± 0.28</td><td>1.50 ± 0.36</td><td>0.44</td></tr><tr><td>Triglycerides, g/L</td><td>1.04 ± 0.46</td><td>1.45 ± 1.02</td><td>0.03<sup>e</sup></td></tr><tr><td>ApoA-I, g/L</td><td>1.48 ± 0.22</td><td>1.49 ± 0.24</td><td>0.92</td></tr><tr><td>ApoB, g/L</td><td>1.05 ± 0.19</td><td>1.12 ± 0.24</td><td>0.19</td></tr><tr><td>Lp(a), g/L</td><td>0.25 ± 0.28</td><td>0.35 ± 0.39</td><td>0.32<sup>e</sup></td></tr><tr><td>hs-CRP, mg/L</td><td>4.4 ± 6.5</td><td>4.6 ± 8.6</td><td>0.38<sup>e</sup></td></tr><tr><td>γGT, IU/L</td><td>45 ± 41</td><td>56 ± 63</td><td>0.34<sup>e</sup></td></tr><tr><td>SBP, mmHg</td><td>138 ± 16</td><td>138 ± 13</td><td>0.89</td></tr><tr><td>DBP, mmHg</td><td>83 ± 7</td><td>82 ± 8</td><td>0.53</td></tr><tr><td>Resting heart rate, beat/min</td><td>62 ± 9</td><td>67 ± 12</td><td>0.06</td></tr><tr><td>Hypertension<sup>f</sup>, n (%)</td><td>6 (19.4)</td><td>11 (35.5)</td><td>0.14<sup>e</sup></td></tr><tr><td>Treatment for hypertension, n (%)</td><td>6 (19.4)</td><td>10 (32.3)</td><td>0.14<sup>e</sup></td></tr><tr><td>Diabetes<sup>f</sup>, n (%)</td><td>3 (9.7)</td><td>2 (6.5)</td><td>0.65<sup>e</sup></td></tr><tr><td>Treatment for diabetes, n (%)</td><td>3 (9.7)</td><td>2 (6.5)</td><td>0.65<sup>e</sup></td></tr><tr><td>Dyslipidemia<sup>f</sup>, n (%)</td><td>12 (38.7)</td><td>17 (54.8)</td><td>0.23<sup>e</sup></td></tr><tr><td>Treatment for dyslipidemia, n (%)</td><td>7 (22.6)</td><td>7 (22.6)</td><td>1<sup>e</sup></td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Serum levels of s-RNY1-5p in individuals from the nested Case–Control study when they were first included in the GENES cohort. Box plot showing the absolute expression of circulating s-RNY1-5p in Cases (n = 30) and Controls (n = 30). Absolute expression data derived from the standard curves created by synthetic s-RNY1-5p by RT-qPCR, normalized by 1 mL of serum, and presented as median and IQR. Wilcoxon signed rank sum test **p**<0.01.</figcaption></figure>
<table><caption>**Table 2.** Serum levels of sRNY1-5p at baseline. Serum sRNY1-5p levels were measured at baseline (i.e. when individuals were first included in the GENES cohort). Values are expressed as mean ± SD or percentage. <sup>a</sup>Paired Student's t-test. <sup>b</sup>McNemar's test.</caption><thead><tr><th></th><th>Controls (n = 30)</th><th>Cases (n = 30)</th><th>P</th></tr></thead><tbody><tr><td>Mean sRNY1-5p, molecules/mL</td><td>592 ± 401</td><td>905 ± 521</td><td>0.027<sup>a</sup></td></tr><tr><td>sRNY1-5p > median value, %</td><td>33.3</td><td>70.0</td><td>0.022<sup>b</sup></td></tr></tbody></table>

<table><caption>**Table 3.** Risk of acute coronary event as a function of sRNY1-5p levels. Model 1: adjusted for age, smoking, hypertension, diabetes and dyslipidemia. Model 2: Model 1 plus systolic and diastolic blood pressure, total cholesterol and HDL-C. HR (95% CI) hazard ratio (95% confidence interval). <sup>a</sup>sRNY1-5p median value is 631 molecules/mL of serum.</caption><thead><tr><th></th><th>sRNY1-5p<br>≤ median <sup>a</sup></th><th>sRNY1-5p<br>> median <sup>a</sup></th><th>p</th><th>Per increase of 10 molecules/mL</th><th>p</th></tr></thead><tbody><tr><td>Model 1, HR (95% CI)</td><td>1</td><td>4.33 (1.23–15.30)</td><td>0.023</td><td>1.03 (1.01–1.05)</td><td>0.022</td></tr><tr><td>Model 2, HR (95% CI)</td><td>1</td><td>6.30 (1.05–37.70)</td><td>0.044</td><td>1.04 (1.01–1.08)</td><td>0.039</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Predictive value of serum s-RNY1-5p for CHD events. Kaplan–Meier probability curves to encounter CHD events according to s-RNY1-5p level, above or below the sRNY1-5p median value (631 molecules/mL). Number at risk refers to the number of individuals who are still CHD-free every 5 years with the number of CHD events in parentheses.</figcaption></figure>

increased more than fourfold in individuals with serum s-RNY1-5p above the median value (HR = 4.33 [95% CI 1.23–15.30]). Significant association with CHD events was also observed when considering s-RNY1-5p as a continuous variable. Indeed, the adjusted HR for coronary event per increase of 10 molecules/mL was 1.03 [95% CI 1.01–1.05], p = 0.022. This positive association of s-RNY1-5p with CHD events was maintained following further adjustments for systolic and diastolic blood pressure, total cholesterol and HDL-C (Model 2, HR = 1.04 [95% CI 1.01–1.08], p = 0.039). Then, cumulative risk of patients for coronary event (fatal and non-fatal) was determined by the Kaplan–Meier method according to the median s-RNY1-5p value (631 molecules/mL) and compared, using the Log-rank test for the individual endpoints of CHD events (Fig. 3). This analysis indicated that the proportion of CHD event-free individuals was significantly lower in individuals with s-RNY1-5p level above the median value (Fig. 3, Log-rank test: Chi² = 4.8, p = 0.017). The difference between the two groups was particularly evident for the follow-up period > 7.5 years. Thus, high levels of s-RNY1-5p are associated with an increased risk of developing a first coronary event.

## Discussion

s-RNYs-5p represent a large proportion of all small RNAs expressed in human serum and plasma¹¹. They circulate as part of a complex with a mass between 100 and 300 kDa¹¹, which contains the protein Ro60². Specific subtypes of s-RNYs-5p have been found dysregulated in the serum of patients with breast cancer and have been proposed as novel biomarkers for diagnosis in this pathology¹², as well in head and neck squamous cell carcinoma¹³. In 2015, we found that s-RNY1-5p is produced by macrophages and that its relative expression levels in serum, normalized by endogenous or spike-in controls with the 2⁻ΔΔCt method, was significantly upregulated in CHD patients compared to age-matched CHD-free individuals².

The goal of the present study was to evaluate whether serum s-RNY levels could be useful for primary prevention of CHD. To this aim, for the first time, we have developed and validated a method to quantify the absolute expression of s-RNY1-5p in human serum. Based on a nested case-control study of CHD-free individuals at baseline and followed for a median duration of 12 years, we observed that serum s-RNY1-5p levels are higher in
individuals who developed an acute coronary event during the follow-up. Notably, we observed that s-RNY1-5p levels were positively associated to CHD events, independently of classical risk factors. Interestingly, s-RNY1-5p seems to predict long-term and not short-term incidence of CHD. Since the expression of s-RNY1-5p in macrophages was induced by pro-apoptotic and atherogenic stimuli, we hypothesized that s-RNY1-5p levels might reflect the local inflammatory condition of the vascular cells². Thus, s-RNY1-5p might be associated to a low-grade inflammatory state of the vascular wall having deleterious effect in the long-term.

Among the small non-coding RNAs, microRNAs (miRNAs) have also become potential biomarkers for human pathologies, including cardiovascular disorders¹⁴. miRNAs are ~20 nucleotides long endogenous non-coding RNAs abundantly expressed in virtually all human cells and associated with AGO proteins to repress gene expression at the post-transcriptional level¹⁵. miRNAs are localized in mostly all cellular organelles and are also present in the extracellular environment¹⁶. More than 100 miRNAs have been detected in serum and designed as circulating miRNAs¹⁴. Because of their extracellular stability, measurement of circulating miRNA dysregulation has been proposed as novel method for the diagnosis of several human pathologies, including CHD¹⁴. However, to date, none of the dysregulated circulating miRNAs may predict cardiovascular accidents in long-range forecast. To the best of our knowledge, the present study is the first one demonstrating the potential role of a small non-coding RNA in predicting CHD events, underlining therefore the originality and importance of our findings.

The prevalence of CHD is increasing and the burden on the healthcare system stems from late diagnosis. In this context, the importance of early identification and preventative intervention is being recognized worldwide. The recommended method for risk screening according to the European Society of Cardiology (ESC) Guidelines for the diagnosis and management of chronic coronary syndromes is the Systematic COronary Risk Evaluation (SCORE)¹⁷. SCORE assesses a combination of clinical parameters (age, gender, smoking status, systolic blood pressure) and blood-based biomarkers (lipid panel of total cholesterol and HDL-C) to estimate a 10-year risk of fatal CHD. However, both the ECS and American College of Cardiology/American Heart Association (ACC/AHA) recognize that current risk evaluation methods are unreliable and have limitations. The strongest pain point experienced by the healthcare system is diagnosis of patients classified as 'moderate risk' following SCORE evaluation. These patients represent a large portion of the population, commonly those who are 35-50 years and diabetics. Primary care physicians then need to reclassify these patients as high or low risk by referring them to clinics or hospitals for expensive cardiac diagnostic testing. As a result, often the prevalence of CHD is overestimated, wasting costly and strained healthcare resources. Therefore, a more accurate biomarker that complements the lipid panel and SCORE assessment would improve the overall predictive performance of CHD risk evaluation. In this respect, our present results indicated that sRNYs, and particularly s-RNY1-5p, could address this unmet need for an early non-invasive biomarker that increase accuracy of CHD risk assessment.

One of the strengths of our study was the use of a matched control group of individuals who did not develop CHD during the follow-up period. However, several limitations of the present study must be noted. Firstly, the small sample size of the study is a limitation. Particularly, results presented here do not allow to address how much s-RNY1-5p adds to the risk prediction of CHD beyond a clinically used risk score. Thus, further measurements are foreseen within the framework of a large prospective cohort to establish more firmly the predictive value of s-RNY1-5p determinations. Yet, our study is the first to report the positive association between s-RNY1-5p upregulated levels with the prediction of CHD events. Secondly, our study was designed only with men, which has the advantage of recording a larger number of events than in a mixed all-gender cohort. Particularly, this cohort is composed of participants living in Toulouse area (South-west France), in which there is low cardiovascular mortality rate in women before 65 years of age¹⁸, thus limiting the translatability of our results to women.

To conclude, we reported here that s-RNY1-5p was positively and independently associated with CHD events in a general male population. Those results argue in favor of s-RNY1-5p being a novel predictive molecular biomarker for cardiovascular events. The present study represents a step forward towards a critical assessment of circulating small non-coding RNAs as biomarkers for prediction of CHD. In the future, measurement of s-RNY1-5p expression levels and other 5' s-RNY, such as s-RNY4-5p², could be used in clinical practice in addition to classical risk factors to identify those high-risk individuals who might benefit from prevention medicine.

Received: 19 April 2020; Accepted: 30 December 2020

Published online: 19 January 2021

## References

1. Moore, K. J., Sheedy, F. J. & Fisher, E. A. Macrophages in atherosclerosis: A dynamic balance. *Nat. Rev. Immunol.* **13**, 709–721 (2013).

2. Repetto, E. et al. RNY-derived small RNAs as a signature of coronary artery disease. *BMC Med.* **13**, 259 (2015).

3. Wolin, S. L. & Reinisch, K. M. The Ro 60 kDa autoantigen comes into focus: Interpreting epitope mapping experiments on the basis of structure. *Autoimmun. Rev.* **5**, 367–372 (2006).

4. Hizir, Z., Bottini, S., Grandjean, V., Trabucchi, M. & Repetto, E. RNY (YRNA)-derived small RNAs regulate cell death and inflammation in monocytes/macrophages. *Cell Death Dis.* **8**, e2530 (2017).

5. Bouisset, F. et al. Prognostic usefulness of clinical and subclinical peripheral arterial disease in men with stable coronary heart disease. *Am. J. Cardiol.* **110**, 197–202 (2012).

6. Genoux, A. et al. Serum IF1 concentration is independently associated to HDL levels and to coronary heart disease: The GENES study. *J. Lipid Res.* **54**, 2550–2558. https://doi.org/10.1194/jlr.P036335 (2013).

7. Genoux, A. et al. Mitochondrial inhibitory factor 1 (IF1) is present in human serum and is positively correlated with HDL-cholesterol. *PLoS ONE* **6**, e23949 (2011).

8. Duparc, T. et al. Serum level of HDL particles are independently associated with long-term prognosis in patients with coronary artery disease: The GENES study. *Sci. Rep.* **10**, 8138. https://doi.org/10.1038/s41598-020-65100-2 (2020).

9. Resnick, H. E. et al. Relationship of high and low ankle brachial index to all-cause and cardiovascular disease mortality: The Strong Heart Study. *Circulation* **109**, 733–739 (2004).
10. Wietlisbach, V., Marques-Vidal, P., Kuulasmaa, K., Karvanen, J. & Paccaud, F. The relation of body mass index and abdominal adiposity with dyslipidemia in 27 general populations of the WHO MONICA Project. *Nutr. Metab. Cardiovasc. Dis.* **23**, 432–442.
https://doi.org/10.1016/j.numecd.2011.09.002 (2013).

11. Dhahbi, J. M. et al. 5'-YRNA fragments derived by processing of transcripts from specific YRNA genes and pseudogenes are abundant in human serum and plasma. *Physiol. Genomics* **45**, 990–998 (2013).

12. Dhahbi, J. M., Spindler, S. R., Atamna, H., Boffelli, D. & Martin, D. I. Deep Sequencing of serum small RNAs identifies patterns of 5' tRNA half and YRNA fragment expression associated with breast cancer. *Biomark. Cancer* **6**, 37–47 (2014).

13. Victoria Martinez, B. et al. Circulating small non-coding RNA signature in head and neck squamous cell carcinoma. *Oncotarget* **6**, 19246–19263 (2015).

14. Zhou, S. S. et al. miRNAs in cardiovascular diseases: Potential biomarkers, therapeutic targets and challenges. *Acta Pharmacol. Sin.* **39**, 1073–1084 (2018).

15. Bartel, D. P. Metazoan MicroRNAs. *Cell* **173**, 20–51 (2018).

16. Trabucchi, M. & Mategot, R. Subcellular heterogeneity of the microRNA machinery. *Trends Genet.* **35**, 15–28 (2019).

17. Mach, F. et al. 2019 ESC/EAS Guidelines for the management of dyslipidaemias: Lipid modification to reduce cardiovascular risk. *Eur. Heart J.* **41**, 111–188. https://doi.org/10.1093/eurheartj/ehz455 (2020).

18. Tunstall-Pedoe, H. et al. Myocardial infarction and coronary deaths in the World Health Organization MONICA Project. Registration procedures, event rates, and case-fatality rates in 38 populations from 21 countries in four continents. *Circulation* **90**, 583–612 (1994).

## Acknowledgements

We are indebted to C. Lorente for technical assistance. Thanks are due to the Toulouse university hospital biobank (CRB-TBR) for providing the samples.

## Author contributions

L.O.M. and M.T. conceived and designed the project. V.L.C. and V.B., and E.R. performed the experiments. L.O.M. and M.T. wrote the manuscript. J.B.R., B.P., S.N., and M.G.S. performed the data analysis. F.S., M.B., C.M., V.G., J.F. edited the manuscript. All authors reviewed and approved the manuscript.

## Funding

V.M-C. was supported by the Satt-SudEst and FS by the UCA-IDEX. This work was supported by the ANR through the "Investments for the Future" # ANR-11-LABX-0028-01 (LABEX SIGNALIFE), FRM (grant #DEQ20140329551), Agence de la Biomédecine (AMP, Diagnostic prénatal et diagnostic génétique), Satt-SudEst and the Société Française de Nutrition (Prix de Recherche 2017) for VG and MT.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-81221-8.

**Correspondence** and requests for materials should be addressed to L.O.M. or M.T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021