Discovery of circulating miRNAs as biomarkers of chronic Chagas heart disease via a small RNA-Seq approach

Chagas disease affects approximately 7 million people worldwide in Latin America and is a neglected tropical disease. Twenty to thirty percent of chronically infected patients develop chronic Chagas cardiomyopathy decades after acute infection. Identifying biomarkers of Chagas disease progression is necessary to develop better therapeutic and preventive strategies. Circulating microRNAs are increasingly reliable biomarkers of disease and therapeutic targets. To identify new circulating microRNAs for Chagas disease, we performed exploratory small RNA sequencing from the plasma of patients and performed de novo miRNA prediction, identifying potential new microRNAs. The levels of the new microRNAs temporarily named miR-Contig-1519 and miR-Contig-3244 and microRNAs that are biomarkers for nonchagasic cardiomyopathies, such as miR-148a-3p and miR-224-5p, were validated by quantitative reverse transcription. We found a specific circulating microRNA signature defined by low miR-Contig-3244, miR-Contig-1519, and miR-148a-3 levels but high miR-224-5p levels for patients with chronic Chagas disease. Finally, we predicted in silico that these altered circulating microRNAs could affect the expression of target genes involved in different cellular pathways and biological processes, which we will explore in the future.

manifestations 2 .Only 20-30% of patients develop symptoms after years to decades of clinical latency (indeterminate stage), showing no detectable parasitemia in microscopic evaluations at this stage.Myocardial damage in Chagas disease is a progressive process that can be classified according to the degree of myocardial involvement 3 .
MicroRNAs (miRNAs) are a class of small noncoding RNAs (20-25 nucleotides) involved in gene regulation and expression at the posttranscriptional level produced by cells that can be found in extracellular vesicles and blood 4,5 .miRNA transcription depends on RNA polimerase II, which generates a primary transcript (pri-miRNA) with hairpin structure.This pri-miRNA is processed by the Drosha complex to form the pre-miRNA that will be transported to the cytoplasm through the exportin-5.Then, it will be recognized and cut by Dicer and RNA-induced silencing complex (RISC) proteins to create the mature miRNA 6 .miRNAs are key for a set of molecular signaling pathways and cellular pathophysiological effects of different but interrelated disorders, such as cardiovascular diseases, cardiac hypertrophy, coronary heart disease, myocardial infarction [7][8][9][10][11] , infections, type 1 diabetes, inflammatory, infectious and autoimmune processes, and cardiovascular risk factors [12][13][14][15] .An advantage is that changes in miRNA levels in body fluids occur earlier than conventional biomarkers 16 .This makes them promising potential candidates as prognostic biomarkers of disease.miRNAs may play a major role in the control of gene expression in key pathological processes in chronic Chagas cardiomyopathy (CC) 5 .However, more studies are needed to understand the role of miRNAs in CC and there are only a few studies on miRNAs as biomarkers in this disease.These studies used quantitative reverse transcription PCR (RT-qPCR) for specific miRNAs and suggested a role in the pathogenesis, diagnosis, and prognosis of Chagas disease 17,18 .Circulating miR-208a, a heart-specific miRNA playing a critical role in heart failure, was increased during the chronic indeterminate phase when compared to chronic cardiac clinical forms in Chagas disease individuals 17 .Additionally, circulating miR-19a-3p, miR-21-5p, and miR-29b-3p, previously associated with heart failure, cardiac fibrosis, and hypertrophy, were increased in chronic cardiac Chagas disease and correlated with cardiac injury 18 .In addition, other authors described the dysregulation of different miRNAs in the cardiac tissue of patients with chronic Chagas cardiomyopathy 19 .
Considering the urgent need of biomarkers in CC, and the absence of studies addressing the role of miRNAs in CC, in the present study, we aimed for the first time to identify new circulating miRNAs using an exploratory small RNA sequencing (Small RNA-Seq) approach and posterior RT-qPCRs in the plasma of chagasic patients at various stages of the disease.Therefore, these miRNAs could be used as biomarkers of chronic Chagas cardiomyopathy and likely involved in its development.

Small RNA-Seq alignment and analysis
To explore the circulating miRNA transcriptomic profile in chagasic patients, we performed exploratory small RNA-Seq of pools of miRNAs extracted from the plasma of healthy controls (HC), indeterminate Chagas disease (ICD), and mild and moderate Chagas cardiomyopathy (MCC) patients (n = 5 each) attending to the Hospital Provincial del Centernario (HPC, Rosario, Argentina).We aligned the reads to the human genome (H.sapiens, GRch38, version 25/03/2019) using Bowtie 2. The mean percent alignment of total reads aligned at least once to the human genome was 98.95% (Table S2).Then, reads aligning to known miRNAs were determined using miRBase annotations (http:// www.mirba se.org/ ftp.shtml) and normalized (Supp.Data S2).
An initial analysis was performed to check whether we could find altered levels of miRNAs in the chagasic patients compared to the HCs.We performed a preliminary PCA study with normalized reads (Fig. S1).PC components explained 100% of the variance of the data.As expected, ICD and MCC grouped separately from HC individuals, indicating that we would probably find candidate miRNAs for biomarkers.

Identification of new miRNAs
To identify new miRNAs in the small RNA-Seq, we used miRDeep2.We identified several potential new miRNAs (Table S3, column 2).miRDeep2 considers the loop sequence and structure compatible with Dicer processing as well as the conservation of the seed sequence (among other stringent parameters) to predict sequences of potentially new miRNAs and structures.
Potential new miRNAs identified in chagasic patients with at least 80% similarity were grouped into clusters, which are defined as potential precursors of miRNAs and/or mature miRNAs separated in the genome from other clusters by a minimum of 30 nucleotides 20 (Table S3, column 3).
Then, data were filtered according to star sequence, precursor, and loop structures, and only those with significant ranfold p-values were considered (Table S4).These short sequences showed a total score and secondary structure consistent with valid mature miRNA sequences, suggesting that they could be considered potential new miRNAs.Moreover, using miRDeep2, we predicted the secondary structure of these 4 new miRNAs and validated them with the precursor structure (Fig. 1).
Most canonical miRNAs are generated from long noncoding regions, and a miR gene can generate more than one mature miRNA able to target different genes (isomiRs).However, although less frequently, miRNAs can also be generated from coding regions 21 .Three of the new miRNAs mapped to the following regions of the genome: miR-Contig-5468 to Homo sapiens CACTIN antisense RNA 1 (CACTIN-AS1), long noncoding RNA (from RefSeq NR_038865); miR-Contig-2269 to Homo sapiens SREK1 interacting protein 1 (SREK1IP1), mRNA (from RefSeq NM_173829); and miR-Contig-3244 to RNA, U2 small nuclear 29, pseudogene (from HGNC RNU2-29P), while miR-Contig-1519 mapped to a noncoding region.

Validation by RT-qPCR of the new circulating miRNAs
To validate the presence of these newly predicted miRNAs in humans, we studied their expression in serum/ plasma samples from patients from Fundación Jiménez Díaz (FJD, Madrid, Spain) and Vall d'Hebrón University Hospital (HVH, Barcelona, Spain).Only 2 out of 4 were confirmed: miR-Contig-1519 and miR-Contig-3244 were detected, while miR-Contig-2269 and miR-Contig-5468 were not detected by RT-qPCR (Fig. 2A).Both miR-Contig-1519 and miR-Contig-3244 showed a marked decrease in chagasic asymptomatic and symptomatic patients of all the hospitals compared to healthy controls.Considering this result, we performed an analysis to study whether these new miRNAs could serve to discriminate patients with chronic Chagas disease from healthy individuals (Fig. 2B).The results showed that both ICD and MCC individuals grouped and separated from the HC subjects.Thus, all these results point to miR-Contig-1519 and miR-Contig-3244 as possible biomarkers of chronic Chagas disease, reducing their expression in patients.

Levels of known circulating miRNAs of interest by RT-qPCR
After the analysis of new miRNAs, we also decided to perform a study of known miRNAs, with more than 100 normalized reads, searching for previous relationships with cardiovascular diseases in databases to complete our

Target genes and pathways predicted for miR-Contig-1519, miR-Contig-3244, miR-148a-3p and miR-224-5p
From the above results, we defined a potential signature defined by low miR-Contig-3244, miR-Contig-1519, and miR-148a-3p and high miR-224-5p levels in chronic Chagas disease and cardiomyopathy (ICD and MCC).However, to explore the possible value of these miRNAs in terms of their contribution to the pathology of chronic Chagas disease, we used bioinformatics tools to predict the possible target genes of these miRNAs and the corresponding pathways in which they could be involved.Using miRDB, we predicted the target genes for these miRNAs (http:// www.mirdb.org) (Supporting Data S3).We used the Target Search tool to determine the targets of the known miR-148a-3p and miR-224-5p that are present in the database of miRDB.For the new miRNAs, miR-Contig-1519 and miR-Contig-3244, we used the Custom Prediction tool that, according to the miRNA sequence, predicts the more likely targets.Then, we considered those targets with a score equal to or higher than 90 out of 100 and performed enrichment analysis using the EnrichR website (BioPlanet and Reactome databases).
The EnrichR web tool allows enrichment analysis according to different databases to associate gene lists with specific diseases or pathways and returns an adjusted p value of the Fisher exact test after Benjamini-Hochberg correction 22,23 .We used the BioPlanet database since it integrates pathway annotations from available, manually curated sources with consistency and redundancy cross-evaluation 24 and the Reactome database, which provides a broad range of physiological and pathological biological processes in humans manually curated from the primary literature and peer-reviewed manuscripts 25 .The complete list of predicted significant pathways or biological processes likely dependent on these miRNAs of the signature is shown in Supporting Data S4.
The BioPlanet database results showed that the principal targets of these four miRNAs belong to different pathways, but most of them were related to transforming growth factor (TGF)-β signaling and the transcriptional activity of the downstream heterotrimer suppressor of mothers against decapentaplegic (SMAD)2/3-SMAD4 (Fig. 5A).Additionally, some of the targets were part of ERBB signaling and T-cell signal transduction.Considering the Reactome database results (Fig. 5B), signaling of TGF-β members and receptors were the most significant biological processes.The forkhead box O (FOXO)-mediated transcription was the third most significant, and we also detected a remarkable regulation of transcriptional pathways such as those related to RNA polymerase II.Finally, these four miRNAs also regulated clathrin-mediated endocytosis, one of the ways for the parasite to enter host cells, and signaling by the serine/threonine-protein kinase B-Raf (BRAF) and the rapidly accelerated fibrosarcoma (RAF)1.

Discussion
Previous research about the role of specific miRNAs in T. cruzi infection have mainly used RT-qPCR analysis of blood or cardiac tissue.In our study we have performed for the first time an exploratory small RNA-Seq analysis in the plasma of chagasic patients with different degrees of cardiac pathology, and we identified and validated by RT-qPCR new circulating miRNAs and a miRNA signature as a potential biomarker of chronic Chagas disease.This signature is characterized by low levels of miR-Contig-1519, miR-Contig-3244, and miR-148a-3p and high levels of miR-224-5p.
Of the four novel circulating miRNAs identified, we were able to detect miR-Contig-1519 and miR-Contig-3244, and both levels significantly decreased in all the chagasic cohorts analyzed.In the case of miR-Contig-2269 and miR-Contig-5468, either technical problems involving primer design hampered their detection or because they may be short mRNAs lacking miRNA activity.Further studies are needed to probe if they are real miRNAs and/or see their effect during the T. cruzi infection.Considering miR-Contig-1519 and miR-Contig-3244, we confirmed their use to discriminate healthy individuals from chronic chagasic subjects, independently if they are indeterminate or patients with cardiomyopathy.Therefore, this result suggests a role of these two new miRNAs in humans that is not related to the progression of the disease but to the pathological response associated with T. cruzi, that should be further elucidated in posterior research considering their targets and the affected biological pathways to provide mechanistic insights into their functions.
To complete our results and establish a more powerful signature of chronic chagasic patients, we also analyzed three known miRNAs that had been described for cardiomyopathies of other etiology but not for Chagas disease: miR-148a-3p, miR-181a-5p, and miR-224-5p.miR-148a-3p is implicated in immunopathologic conditions and could be a potential biomarker for different cardiac pathologies.According to some authors miR-148a-3p is a novel repressor of IKBKB, which encodes the inhibitor of nuclear factor kappa B (IκB) kinase β, that represses the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signalling pathway that plays a key role in the inflammatory cascade 26 .Also, overexpression of miR-148a-3p inhibited inflammatory factors and functional injury in vascular endothelial cells 27 .Therefore, according to our data, a decrease in miR-148a-3p levels in chagasic patients could be related to an increase in endothelial and cardiac inflammation.
Moreover, miR-148a-3p affects macrophage polarization and cytokine production: it regulates the activation and polarization of mouse macrophages targeting the phosphatase and tensin homolog (PTEN) through the phosphatidylinositol 3-kinase (PI3K) and protein kinase B (AKT) pathway in response to Schistosoma japonicum infection 28 .Indeed, macrophage polarization is a key in the immune response against of T. cruzi infection 29 .Also, this miRNA is a novel downstream molecule of Notch signalling, promoting the differentiation of monocytes into macrophages by the presence of granulocyte macrophage colony-stimulating factor (GM-CSF).Specifically, it enhances M1 macrophage (activated by the classic pathway) differentiation, inhibits M2 macrophage polarization (activated by the alternative pathway) and increases the ability to engulf and kill bacteria, mediated by excessive production of reactive oxygen species (ROS) through the PTEN/AKT pathway 30 .Therefore, the reduction in miR-148a-3p in some chronic chagasic patients could inhibit M1 macrophage polarization and ROS production, both related to the elimination of the parasite and its persistence in immune cells.We did not detect significant differences in the levels of miR-181a-5p between any of the groups, even though this miRNA is involved in cardiomyocyte apoptosis, a general problem of several cardiomyopathies 31 .
We found increased levels of miR-224-5p in ICD and MCC patients compared to the HC group.According to previous studies, this miRNA is a potential target molecule for treating and diagnosing atherosclerosis 32 .Additionally, it is increased in heart failure 33 , coronary artery disease 34 , and circulating extracellular vesicles of patients with reduced coronary flow reserve, which is related to endothelial dysfunction 35 .The upregulation of this miRNA in chronic chagasic patients, even in the absence of cardiac symptoms, is a new revelation that should be studied in the future to understand its specific role in the progression of Chagas disease.
Enrichment pathway analysis considering miR-Contig-3244, miR-Contig-1519, miR-148a-3p, and miR-224-5p targets was performed through the EnrichR webtool using the BioPlanet and Reactome databases.This analysis predicted that most of the targets are related to TGF-β signaling, including the transcriptional activity of the heterotrimer SMAD2/3-SMAD4.After their phosphorylation by TGF-β, these proteins form a heterotrimer that translocates to the nucleus to mediate intracellular TGF-β signal transduction 36 .In Chagas disease, TGF-β is involved the progression of the CC onset, promoting the heart cell invasion and intracellular replication by T. www.nature.com/scientificreports/cruzi, host immune and the inflammation response, and posterior cardiac fibrosis.In fact, inhibitors of TGF-β signaling revert these effects and improve cardiac electric parameters in infected mice, although more studies in infected human cells or humanized animal models are needed to clarify its potential as therapeutic approach.Additionally, Chagas disease patients with higher TGF-β in serum display a worse clinical outcome suggesting a predictive value as a surrogate biomarker [37][38][39] .Our results support these ideas, and further studies are necessary to understand the specific role and regulation that these four miRNAs may exert on TGF-β signaling and their implication as possible biomarkers and/or therapeutic approaches for Chagas disease.Another biological process regulated by these miRNAs was T-cell signal transduction.T cells are key for the immune response against the parasite.CD8+ T cells take part in the elimination of the parasite, although some evidence suggests that they are also involved in some of the clinical manifestations of Chagas disease, such as tissue damage and inflammatory processes.The balance between CD4+ Th1 and Th2 cells is also important during infection.Th1 cytokines contribute to the control and elimination of the parasite but produce an enhanced immune response, while Th2 cytokines are implicated in protective responses, although they allow host and parasite cell survival.However, some of the authors suggest that more research about the specific parasite targets of T cell memory response may contribute to the conceptualization and development of prophylactic or therapeutic vaccines [40][41][42][43] .Therefore, studies focusing on these specific types of cells could help to discover the role of these four miRNAs in the dilemma of T-cell implication in Chagas disease.
Additionally, our miRNAs regulate targets linked to clathrin-mediated endocytosis and RAF1/BRAF signaling.Clathrin-mediated endocytosis is one of the main ways of entry of T. cruzi into the host cell from the initial stages of infection until the formation of the parasitophorous vacuole, and its inhibition drastically reduces the internalization of trypomastigotes in phagocytic and epithelial cells 44 .Therefore, a critical change in this pathway may alter all the endocytic processes of the parasite and its posterior intracellular replication, although more studies with different T. cruzi strains and other cellular types that are infected by the parasite could help to confirm the relevance of the clathrin-mediated endocytosis in the infection process.
Regarding RAF1/BRAF signaling, both proteins are upstream kinases and activators of the mitogen-activated protein kinase (MAPK) cascade, which contains the mitogen-activated protein kinase kinase (MEK) and their targets, the extracellular signal-regulated kinase (ERK) family.Although there are not specific studies focused on RAF1 or BRAF proteins in T. cruzi infections, some.Some researchers studied the MAPK cascade activation and role in Chagas disease MEK/ERK pathway increase in nervous system models by the interaction of the host cells with surface molecules of T. cruzi as members of the multigene family of trans-sialidases 45,46 , while the inhibition of ERK1/2 phosphorylation produced a significant decrease in invasion and infection by this parasite in epithelial and cardiac muscle cells 47 and in cultured human umbilical vein endothelial and vascular smooth muscle cells 48 .Therefore, the regulation of targets of this pathway by these miRNAs may drastically alter the infection process of T. cruzi and should be studied in the future to discover its real potential as a treatment for Chagas disease.
Interestingly, our results complement previous observations on miRNAs in different contexts of T. cruzi infection by other authors.For instance, studies in patients pointed to miR-146a as a biomarker of the indeterminate phase of Chagas disease 49 , and miRNA-155 deficiency exacerbated T. cruzi infection in mice, suggesting that miR-155 is an important immune regulatory molecule critical for the control of T. cruzi in the acute phase of the infection 50 .In addition, it has been reported that T. cruzi induces a differential miRNA profile in human placental explants that may be potential targets for the therapeutic control of congenital Chagas 51 .On the other hand, it has been reported that modulation of miR-145-5p and miR-146b-5p levels is linked to reduced parasite load in H9C2 T. cruzi-infected cardiomyoblasts, which points to a possible direct action of miRNAs in heart tissue 52 .Thus, all these reports highlight the importance of our results and previous findings by other groups about the role of miRNAs in Chagas disease.
In summary, we have identified a signature of circulating known and new miRNAs as a potential biomarker and/or therapeutic target of the chronic phase of Chagas disease.Moreover, it is important to remember that miRNA expression could vary between individuals according to both transcriptional changes in gene expression and promoter hypermethylation, and post-transcriptional changes in miRNA processing mechanisms of regulation, as well as effects of exogenous (xenobiotics) or endogenous (hormones, cytokines) compounds 53 .We recognize the limitations of our study and results (sample size, mainly), which should be considered exploratory, and therefore, these miRNAs should be further explored in a larger cohort of individuals to validate their promising use as biomarkers, therapeutic targets, or clinical treatments.In addition, more studies including severe Chagas cardiomyopathy (SCC) patients could be useful to evaluate their prognostic value.We consider it necessary for Chagas disease since it is not possible to predict whether ICD and MCC patients will develop SCC.

Blood samples from patients
This work corresponds to an observational, cross-sectional study, where the patients and healthy volunteers were recruited from 2011 to 2017 at the Chagas disease ambulatory section from the Hospital Provincial del Centenario (HPC, Rosario, Argentina), in 2018 at the Fundación Jiménez Díaz (FJD, Madrid, Spain), and in 2020 at the Vall d'Hebrón University Hospital (HVH, Barcelona, Spain).Twenty-eight chagasic patients from Spanish hospitals were from Bolivia, 1 from Brazil, and 1 from Paraguay.Exclusion criteria comprised patients with previous treatment with Benznidazole or Nifurtimox.Healthy control subjects were seronegative for the T. cruzi-specific test.The diagnosis was based on at least two positive serological tests (either by ELISA, hemagglutination, or immunofluorescence), together with clinical symptoms, chest X-ray, 12-lead resting electrocardiogram (ECG), and/or echocardiography.Routine laboratory studies were also included.
A total of 66 subjects were included in the study: healthy controls (HC, n = 26) from endemic areas, indeterminate Chagas disease (ICD, n = 19), and mild and moderate Chagas cardiomyopathy (MCC, n = 21).Table S1 compiles all subjects included in the study classified by the hospital in which they were recruited.The criteria for the aggrupation of chagasic patients from HPC were as follows: ICD, symptomless, normal ECG and chest X-ray; MCC, presented no congestive heart failure with ECG showing incomplete or complete right bundle branch block, ventricular arrhythmia or normal heart size or only mild cardiomegaly (chest X-ray cardiothoracic ratio < 0.55).For patients from FJD: ICD, without pathologic alterations; and MCC, with pathologic alterations in ECG and/or echocardiogram and/or complementary tests compatible with CC.For patients from HVH: ICD, Kuschnir classification 0 and MCC, Kuschnir classification I. Individual information of the patients is shown in Supporting Data S1.

Circulating miRNA extraction
Blood samples from patients and healthy volunteers were collected.Plasma from blood samples of HPC and HVH was obtained from EDTA-treated blood.Serum was obtained from blood samples from FJD patients.All the samples were immediately stored at − 80 °C until miRNA extraction.To avoid impacts in the circulating miRNAs profile, as other studies have demonstrated 54,55 , hemolysis processes were checked visually according to the hemolysis reference palette 56 , and all the samples used in the study were free of hemolysis.In samples from HPC and FJD, miRNAs were extracted using a QIAamp miRNeasy Serum/plasma kit (Qiagen, Germantown, USA), and samples from HVH were extracted using a miRNeasy Serum/Plasma Advance Kit (Qiagen).A spikein mix (Qiagen, RNA Spike-In Kit containing sp2, sp4, and sp5) was added as a control for sample quality and PCR efficiency.The procedures were carried out following the directions of the manufacturer.Small RNA-Seq was performed using samples from HPC, and two or three extractions of each sample were performed to obtain at least 25 ng per sample.Then, the different extractions of each sample were gathered together.The amount of RNA was quantified using RiboGreen.Pools were made for each group (HC, ICD, and MCC), and each group was composed of 5 samples.The pools were concentrated to obtain at least 25 ng in a final volume of 6 µl, which passed the quality assessment in the bioanalyzer and were sent for sequencing to the Parque Científico de Madrid (Madrid, Spain).

Small RNA-Seq analysis
The total RNA fraction corresponding to miRNAs (15-30 bases) was purified by size exclusion in polyacrylamide gels and the miRNeasy Micro Kit (Qiagen Cat.No. 217084).Subsequently, the sequencing libraries were prepared with the TruSeq kit (Illumina) for Small RNA-Seq analysis following the specific protocols of the Illumina technology (TruSeq) without modifications by the Genomics facility of the Parque Científico de Madrid (Madrid, Spain), and readings of 14-50 bases in length were generated in a single replicate using MiSeq technology (Illumina).Differential expression analysis and miRNA de novo finding were performed using the DEseq2 57 and miRDeep2 58 packages, respectively.

miRNA analysis by RT-qPCR
RT-qPCR analysis was performed at the Parque Científico de Madrid (Madrid, Spain).cDNA synthesis was performed with the miRCURY LNA RT Kit (Qiagen Cat.No. 339340) according to the manufacturer's instructions.Then, 2 µl of purified total miRNAs (10 ng) from plasma or serum was used as input RNA and mixed with 0.5 µL synthetic RNA spike-in (sp6) before reverse transcription.The obtained cDNAs were diluted 1/40 in RNAse-free water before qPCR amplification.Real-time PCR amplification was performed using the miRCURY SYBR Green PCR Kit (Qiagen Cat.No. 339347) in the presence of PCR Locked Nucleic Acid (LNA) primers specific for spike-in sp2, sp4, sp5, and sp6, miR-103a-3p, miR-Contig-5468, miR-Contig-1519, miR-Contig-2269, miR-Contig-3244, miR-148a-3p, miR-224-5p and miR-181a-5p (Qiagen).Three replicates were run for each combination of genes and samples.To obtain Cp values, we used LightCycler® 480 Software following the 2nd derivative max method.The expression levels of target miRNAs were normalized to the expression of endogenous miR-103a-3p.The results refer to HC samples.Relative changes and statistical values of miRNA expression levels were calculated by the ΔΔC t method using LightCycler® 480 Software without correction for assay efficiency.The results were expressed as log 10 RQ (RQ = 2 -ΔΔCt ).

Statistical analysis of the small RNA-Seq and the RT-qPCR
The main limitation of the present study is the lack of technical replicates in small RNA-Seq; therefore, Basemean (normalized reads) was chosen as an additional filter to the adjusted p value as a measure of miRNA abundance.miRNAs with less than 100 normalized reads in the healthy controls and the patient samples were discarded for further analysis.All PCAs were performed using the R programming language (version 4.2.2).Bar plots and statistical significance were evaluated using GraphPad Prism (Version 8.0).After the Shapiro-Wilk normality test, two-way ANOVA tests were applied as indicated in the figure legends (*p<0.05,**p<0.01,***p<0.001,****p<0.0001).

Figure 1 .
Figure 1.Predicted secondary structure of the hairpins corresponding to the new miRNAs using miRDeep2.Red: mature miRNA, blue: star sequence, orange: loop.
This study was approved by the Bioethics Committee from the Faculty of Medical Sciences, National University of Rosario, Argentina (Resolution N° 5093/2018); the Ethical Committees of Ramón y Cajal University Hospital and Fundación Jimenez Diaz in Madrid (Spain, References FUN-BEN-2007-01 and 10/2016, respectively); and the Ethical Committee of Vall d'Hebrón University Hospital in Barcelona (Spain, Reference D-RTF080).All patients signed informed consent forms.Data on human subjects were analyzed anonymously, and clinical investigations were conducted according to the Declaration of Helsinki.