
==== Front
J Gastroenterol
J Gastroenterol
Journal of Gastroenterology
0944-1174
1435-5922
Springer Nature Singapore Singapore

39060521
2135
10.1007/s00535-024-02135-6
Original Article―Liver, Pancreas, and Biliary Tract
Profile of plasma microRNAs as a potential biomarker of Wilson’s disease
Sánchez-Monteagudo Ana 1
Ripollés Edna 12
Murillo Oihana 3
Domènech Sofia 12
Álvarez-Sauco María 4
Girona Eva 5
Sastre-Bataller Isabel 26
Bono Ariadna 27
García-Villarreal Luis 8
Tugores Antonio 8
García-García Francisco 29
González-Aseguinolaza Gloria 310
Berenguer Marina 2711
http://orcid.org/0000-0003-4435-1809
Espinós Carmen cespinos@cipf.es

121213
1 grid.418274.c 0000 0004 0399 600X Unit of Rare Neurodegenerative Diseases, Valencia Biomedical Research Foundation-Centro de Investigación Príncipe Felipe (CIPF), Calle Eduardo Primo Yúfera No. 13, 46012 Valencia, Spain
2 grid.476458.c 0000 0004 0427 8560 Rare Diseases Joint Unit, CIPF-IIS La Fe, Valencia, Spain
3 https://ror.org/02rxc7m23 grid.5924.a 0000 0004 1937 0271 DNA@RNA Medicine Division, Centro de Investigación Médica Aplicada (CIMA), University of Navarra, Pamplona, Spain
4 https://ror.org/02ybsz607 grid.411086.a 0000 0000 8875 8879 Department of Neurology, Hospital General Universitari d’Elx, Alicante, Spain
5 https://ror.org/02ybsz607 grid.411086.a 0000 0000 8875 8879 Department of Internal Medicine, Hospital General Universitari d’Elx, Alicante, Spain
6 https://ror.org/01ar2v535 grid.84393.35 0000 0001 0360 9602 Department of Neurology, Hospital Universitari i Politècnic La Fe, Valencia, Spain
7 https://ror.org/01ar2v535 grid.84393.35 0000 0001 0360 9602 Hepatology-Liver Transplantation Unit, Digestive Medicine Service, IIS La Fe and CIBER-EHD, Hospital Universitari i Politècnic La Fe, Valencia, Spain
8 grid.411322.7 0000 0004 1771 2848 Research Unit, Complejo Hospitalario Universitario Insular Materno-Infantil, Las Palmas de Gran Canaria, Spain
9 grid.418274.c 0000 0004 0399 600X Unit of Bioinformatics and Biostatistics, Valencia Biomedical Research Foundation-Centro de Investigación Príncipe Felipe (CIPF), Valencia, Spain
10 Vivet Therapeutics S.L., Pamplona, Spain
11 https://ror.org/043nxc105 grid.5338.d 0000 0001 2173 938X Department of Medicine, Universitat de València, Valencia, Spain
12 https://ror.org/01460j859 grid.157927.f 0000 0004 1770 5832 Biotechnology Department, Universitat Politècnica de València, Valencia, Spain
13 grid.413448.e 0000 0000 9314 1427 Rare Diseases Networking Biomedical Research Centre (CIBERER), Instituto de Salud Carlos III (ISCIII), Madrid, Spain
26 7 2024
26 7 2024
2024
59 10 921931
16 4 2024
17 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
Background

Wilson’s disease (WD) is a rare condition resulting from autosomal recessive mutations in ATP7B, a copper transporter, manifesting with hepatic, neurological, and psychiatric symptoms. Timely diagnosis and appropriate treatment yield a positive prognosis, while delayed identification and/or insufficient therapy lead to a poor outcome. Our aim was to establish a prognostic method for WD by characterising biomarkers based on circulating microRNAs.

Methods

We conducted investigations across three cohorts: discovery, validation (comprising unrelated patients), and follow-up (revisiting the discovery cohort 3 years later). All groups were compared to age- and gender-matched controls. Plasma microRNAs were analysed via RNA sequencing in the discovery cohort and subsequently validated using quantitative PCR in all three cohorts. To assess disease progression, we examined the microRNA profile in Atp7b−/− mice, analysing serum samples from 6 to 44 weeks of age and liver samples at three time points: 20, 30, and 40 weeks of age.

Results

In patients, elevated levels of the signature microRNAs (miR-122-5p, miR-192-5p, and miR-885-5p) correlated with serum activities of aspartate transaminase, alanine aminotransferase and gamma-glutamyl transferase. In Atp7b−/− mice, levels of miR-122-5p and miR-192-5p (miR-885-5p lacking a murine orthologue) increased from 12 weeks of age in serum, while exhibiting fluctuations in the liver, possibly attributable to hepatocyte regenerative capacity post-injury and the release of hepatic microRNAs into the bloodstream.

Conclusions

The upregulation of the signature miR-122-5p, miR-192-5p, and miR-885-5p in patients and their correlation with liver disease progression in WD mice support their potential as biomarkers of WD.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00535-024-02135-6.

Keywords

miR-122-5p
miR-192-5p
miR-885-5p
Prognosis
Drug monitoring
Fundació Per Amor a l'ArtFPAA-Wilson Espinós Carmen http://dx.doi.org/10.13039/501100004587 Instituto de Salud Carlos III PI18/00147 PI21/00103 PI21/00313 INT20/00061 CB06/04/0065 Murillo Oihana Berenguer Marina Espinós Carmen http://dx.doi.org/10.13039/501100004837 Ministerio de Ciencia e Innovación PID2021-124455OB-I00 González-Aseguinolaza Gloria http://dx.doi.org/10.13039/501100014849 Conselleria de Cultura, Educación y Ciencia, Generalitat Valenciana CIACIF/2021/252 Ripollés Edna issue-copyright-statement© Japanese Society of Gastroenterology 2024
==== Body
pmcIntroduction

Wilson’s disease (WD), a rare disorder (MIM 277900; ORPHA:905), with a prevalence ranging from 1/30,000 to 1/50,000 and a relatively high carrier frequency of 1/90, is characterised by autosomal recessive mutations in ATP7B, a copper transporter [1]. Improper functioning of ATP7B causes copper to accumulate in the liver, where it is released and enters the circulation. As a result, toxic copper accumulates in other organs, mainly in the brain, and also in the kidneys, eyes, heart, muscles, and bones [2]. Hence, the clinical presentation is diverse, with systemic damage stemming from copper accumulation in various tissues, particularly the liver and the brain. Liver damage can manifest as acute or chronic hepatitis, fulminant liver failure, or cirrhosis, with normal or elevated liver enzyme levels. Neurological symptoms may include parkinsonism, tremor, and dystonia, while psychiatric signs can involve behavioural changes and cognitive disorders. Fortunately, effective therapy exists, involving daily intake of copper chelators or zinc salts to promote urine copper excretion or reduce copper absorption, respectively, ultimately rebalancing copper levels. Early diagnosis and lifelong treatment adherence are crucial to prevent irreversible organ damage and disease progression.

Despite its significance, WD remains underdiagnosed. In France, the carrier frequency is 1/31, with a prevalence of around 3/200,000 [3]. Diagnosis relies on the Leipzig scale [4]. Incomplete penetrance, variable expressivity, and other disorders mimicking WD contribute to under diagnosis [5, 6]. Biochemical tests, like ceruloplasmin and copper assessments in blood, urine or liver, may be inconclusive [7].

Circulating microRNAs (miRNAs) have emerged as potential biomarkers for liver damage, serving as early indicators of disease and its progression in various pathologies, including metabolic dysfunction-associated steatotic liver disease (MASLD) and metabolic dysfunction-associated steatohepatitis (MASH) [8, 9]. MiRNAs regulate gene expression post-transcriptionally through mRNA degradation or translational repression [10, 11]. Notably, there is a strong correlation between circulating and tissue miRNAs in liver diseases [11].

This study explored plasma miRNAs as potential WD biomarkers. MiR-122-5p, miR-192-5p, and miR-885-5p exhibited altered levels in distinct patient cohorts. Findings were corroborated in an Atp7b−/− mouse model [12–15] to study disease progression. Our findings support that detected miRNA signature may serve as a valuable biomarker for monitoring WD patients.

Materials and methods

Subjects

The research performed with human samples was approved by the Ethics Committee of the Hospital Universitari i Politècnic La Fe (Valencia, Spain; protocol code: 2019/0052; 22/05/2019) and was conducted in compliance with the principles of the Declaration of Helsinki, Good Clinical Practice guidelines, and local regulatory requirements. Informed consent was obtained from all participants prior to research initiation. All data were anonymised.

All of the experimental procedures involving animals were conducted in accordance with the European Union Guidelines for the Care (European Union Directive, 2010/63/EU) and the guidelines for the use of laboratory animals. The experimental design was approved by the Ethical Committee for Animal Testing of the University of Navarra (Pamplona, Spain; protocol code: CEEA/066-22; 02/09/2022).

Statistical power of the data was evaluated from the normalised counts of the experimental subjects, using Bioconductor’s RnaSeqSampleSize R package (https://www.bioconductor.org/packages/release/bioc/html/RnaSeqSampleSize.html). For a cohort of 20 individuals, a probability (power) of 0.70 was obtained, considering that 10% of miRNAs with a differential expression profile may have a relevant role.

Patients

A first cohort (discovery) including 20 in-depth phenotyped WD patients (Supplementary Table S1) was screened by high-throughput small RNA-seq (RNA-sequencing). For validation purposes, a second cohort (validation) and a third cohort (follow-up) comprising 22 and 25 WD patients, respectively (Supplementary Tables S1 and S2), were screened by quantitative PCR (qPCR). The discovery and follow-up cohorts comprised the same patients studied 3 years apart. Two patients included in the discovery cohort could not be further studied because they did not attend the medical check-up, while seven new patients diagnosed during these years and not evaluated in the discovery cohort were included in the follow up cohort. Individuals from the discovery and follow-up cohorts were supervised at Hospital Universitari i Politècnic La Fe (Valencia, Spain) and at Hospital General Universitari d’Elx (Alicante, Spain). Individuals from the validation cohort were recruited at Complejo Hospitalario Universitario Insular Materno Infantil de Gran Canaria (Las Palmas, Spain). These clinical series have been partially reported [6, 16]. Differential diagnosis is based on clinical features, biochemical criteria, presence of corneal Kaiser-Fleischer (KF) ring and other hepatic disturbances. All patients included met the following inclusion criteria: (1) Leipzig score ≥ 3 at diagnosis (without considering ATP7B mutations); and (2) conclusive genetic analysis. Other chronic liver diseases (hepatitis C or B, immune mediated liver diseases or alcohol related liver disease) were ruled out. We excluded patients who had undergone liver transplantation. Supplementary Tables S1 and S2 show the clinical and biochemical features of the patients.

Peripheral blood samples for plasma isolation and biochemical analysis were obtained during clinical follow-up visits. Biochemical determinations analysed using standard protocols (upon request): AST, aspartate transaminase; ALT, alanine aminotransferase; GGT, gamma-glutamyl transferase; bilirubin; AP, alkaline phosphatase; cholesterol; and triglycerides. For miRNA analyses, EDTA-blood samples were processed for plasma separation within 2–3 h following collection. After centrifugation at 1900×g for 10 min (4 °C), plasma fraction was aliquoted and stored at – 80 °C.

Healthy subjects

Plasma from controls were provided by the Biobanco para la Investigación Biomédica y en Salud Pública de la Comunitat Valenciana (PT13/0010/0064), integrated in the Spanish National Biobanks Network and in the Valencian Biobanking Network. Samples were processed following standard operating procedures. Demographic data from controls recorded in general health and lifestyle questionnaires were provided (Supplementary Table S3). Healthy subjects were matched by sex and age with patients, and those who registered neoplasms, cardiovascular, respiratory, bone, mental or endocrine disease risk factors were discarded.

Mice

Serum and liver samples of wild-type (WT) and Atp7b−/− (WD) mice were received from Dra. González-Aseguinolaza laboratory (Centro de Investigación Médica Aplicada, Pamplona, Spain). Atp7b−/− on a C57BL/6 J background were bred and maintained under pathogen-free conditions and genotyped at 3 weeks of age as described [14, 15]. These Atp7b−/− mice show no ATP7B expression in the liver and exhibit the typical biochemical and physiopathological alterations observed in WD patients, except for the neurological signs [12, 13]. Alterations comprise high copper excretion in urine, low holoceruloplasminemia, high serum transaminase levels and increased liver copper concentration with associated hepatocellular damage [12–15].

The cohort comprised 20 or 10 WD mice for the studies using serum or liver, respectively, with the same number of WT mice. Serum samples were collected at 6–7, 10, 12, 16, 20, 30 and 40 weeks old, while liver samples were obtained at 20, 30 and 40 weeks old. Samples were immediately frozen at -80ºC and no tissue preservatives were used.

Circulating miRNA-enriched total RNA isolation

Total RNA enriched for circulating miRNA fraction was extracted from 400 µL human plasma, 10 µL mouse serum and 20 mg mouse liver per subject and column using miRNeasy Mini kit (Qiagen, Hilden, Germany), following manufacturer’s protocol (miRNeasy Mini Handbook 11/2020). Aliquots from plasma/serum samples were thawed on ice and centrifuged (16,000×g, 5 min, 4 °C) to remove cell debris prior to RNA isolation. During isolation from mouse serum, 5.1 × 108 copies of synthetic cel-miR-39 (IDT, Leuven, Belgium) were added into each sample to be used as exogenous miRNA reference.

Small RNA library preparation and high-throughput sequencing

Individual small RNA libraries were prepared from 6 µL of total RNA enriched for circulating miRNA fraction isolated from human plasma samples with NEBNext Multiplex Small RNA Library Prep kit (New England Biolabs, Ipswich, MA, USA). To improve final miRNA sequencing yield, library prep protocol was optimised according to dilution 1:5 of 3′ and 5′ SR RT primer and 16× PCR cycles for indexing. Pre- and post-gel size selection PCR-indexed libraries were purified with Nucleospin Gel and PCR Clean-Up (Macherey–Nagel, Düren, Germany) following manufacturer’s protocol (Macherey–Nagel—03/2023, Rev. 08) for PCR clean-up and DNA purification from polyacrylamide gel, respectively.

Pre- and post-gel size selection library QC (Quality Control) was performed using High Sensitivity D1000 ScreenTape in a TapeStation 4200 (Agilent Technologies, Santa Clara, CA, USA). Post-gel purified libraries were quantified by qPCR with KAPA Library Quantification Kit (Kapa Biosystems, Wilmington, MA, USA) and pooled for sequencing in Illumina HiSeq 2500 (1 × 50 bp, v4).

Bioinformatics analyses

Primary analysis

FASTQ raw data files obtained from small RNA-seq were processed for computational analysis. Before and after removing adapter sequences and selecting reads by size with Cutadapt (v2.6), QC with FastQC (v0.11.8) was performed to check if reads size distribution and quantity was compatible with mature miRNA (16–28 bp) and sufficient for differential representation analyses. Selected reads were aligned against pre-miRNA sequences (hairpin) from the latest miRBase v22 using Bowtie (v1.1.2). Counts of reads matching hairpin arms were obtained with Subread (v1.6.0) and a custom GFF (General Feature Format) file with 5p and 3p coordinates corresponding to annotated human mature miRNAs in miRBase v22.

Differential expression analysis and functional enrichment

To identify differential expressed mature miRNAs, edgeR package was used [17]. Starting from extracted mature miRNA counts, those low-represented were filtered before setting per-sample library size for normalisation by the trimmed mean of M values (TMM) method. QLF (Quasi-Likelihood F) and LRT (Likelihood Ratio Test), both considered generalised linear models (GLM), were applied to determine differentially expressed mature miRNAs in patients compared to controls (design 1), and in a second approach, including the covariables sex and age, to adjust comparisons with controls (design 2). miRNAs with a false discovery rate (FDR) < 0.05 were considered as significantly deregulated.

miRNA detection by quantitative PCR

To obtain cDNA from mature miRNAs, a 2–2.5 µL total RNA enriched for circulating miRNA fraction isolated from human plasma or murine serum or 10 ng if isolated from tissue sample was used with the TaqMan Advanced miRNA cDNA Synthesis Kit (Applied Biosystems, Foster City, CA, USA).

Levels of mature miRNAs of interest were detected by qPCR with TaqMan Advanced (Applied Biosystems, Foster City, CA, USA) probes (Supplementary Table S4) and master mix in a LightCycler 480 II (Roche, Mannheim, Germany) programmed. Pre-amplified cDNA was diluted 1/10 and 2.5 µL per reaction were used as template for qPCR, scaling reaction components to a 10 µL final volume.

Relative miRNA expression levels in case sample (human or murine) compared to control sample was calculated as 2-ΔCtCase-AverageΔCtControl, being ΔCt=CttargetmiRNA-AverageCtreferencemiRNA. As reference miRNAs for qPCR analysis in human plasma samples, endogenous levels of miR-16-5p and miR-484 were used, while for mouse serum, endogenous miR-484 and exogenous cel-miR-39 spike-in levels were considered to normalise data. Finally, for mouse liver samples, endogenous miR-16-5p was used as reference miRNA.

Statistical analyses

Shapiro–Wilk test was applied to biochemical parameters to check for normality distribution, and as all resulted negative, non-parametric Wilcoxon signed-rank test was used. Paired t test was performed to assess differences of relative miRNAs expression levels, represented as log2FC, between cases and controls. For in-group comparisons, two-sample t test was used instead.

For murine samples, Shapiro–Wilk and Levene tests were applied to check for normality distribution and homoscedasticity of miRNA expression levels represented as log2FC. Two-way mixed ANOVA (analysis of variance between-sex factor and within-age factor) and pairwise post hoc t test were used to determine differences between groups.

The Spearman’s Rho (Rs) coefficient of correlation was used to establish the association of relative miRNAs expression levels with biochemical parameters in human plasma and in murine serum, and additionally, with Leipzig score (including ATP7B mutations) and age in human plasma.

Results

Cohorts

Studies were first performed in the discovery cohort by small RNA-seq, and they were further validated in both cohorts, discovery and validation, by qPCR. Supplementary Table S5 shows the demographic and clinical characteristics of the patient cohorts. Both cohorts are quite similar regarding the number of women versus men, with mean ages between 36 and 43 years old. The hepatic presentation was more frequent than the neurologic one, whereas the mixed presentation was rare. Asymptomatic individuals were detected in mutational screenings of relatives. Biochemical determinations resulted to be quite homogeneous between both cohorts, and although ALT levels were more disparate, differences among groups were not significant (Supplementary Table S5).

Determination of miRNAs in plasma of WD patients and controls

Initially, we assessed the quantity and quality of reads obtained in small miRNA-seq, focusing on their processing. Across the 40 sequenced samples, approximately 4–8 million raw reads were obtained, reduced to 4–6 million after removing library adapters and discarding reads smaller than 15 bp (Supplementary Figure S1). Subsequently, 2–4 million reads per sample, corresponding to mature miRNAs (size 16–28 bp), were selected for alignment against precursor miRNA sequences from miRBase v.22. For differential representation analysis, counts of sequences mapped to the determined coordinates in the 5p and 3p arms of each precursor, representing mature miRNAs, were obtained. Proportionally, the number of mapped readings belonging to mature miRNAs sequences was higher in the patient group, reflecting the physiological differences between the two groups (Supplementary Figure S1). After filtering and normalisation, 2158 mature miRNAs out of 2650 were retained, excluding low-represented miRNAs.

Two different experimental designs were applied for differential representation analysis: design 2 considers as covariables sex and age (Supplementary Figure S2). The LRT method identified 42 and 37 differentially expressed miRNAs between designs 1 and 2, respectively (Supplementary Figure S2). The QLF method revealed 29 deregulated miRNAs in design 1 and 30 in design 2. Importantly, 18 miRNAs were common for the four strategies, and therefore, do not depend on the covariates age and sex. The 18 miRNAs presented a positive log2FC value, and hence, they would be overrepresented in patient plasmas. Supplementary Table S6 lists the results obtained using both methods, QLF and LRT, with the two designs.

Profile of miRNAs and correlation with biochemical parameters

From the list of 18 miRNAs, eight miRNAs were selected for further validation (Supplementary Table S6). Seven miRNAs (miR-122-5p, miR-193b-5p, miR-885-5p, miR-885-3p, miR-455-5p, miR-485-3p and miR-340-3p) were selected based on the following criteria: (1) value of log2FC > 1.4 in the differential representation analysis in both experimental designs; and (2) association with liver disease and/or metal metabolism. Of note, because miR-122b-5p, miR-122b-3p and miR-122-5p belong to the same family, we only included miR-122-5p, which is the most studied isoform related to hepatic disease. In addition, the miR-192-5p was also investigated because this miRNA is the second most abundant miRNA in the liver (after miR-122) and is widely related to liver damage [18].

Three (miR-193b-5p, miR-455-5p and miR-885-3p) out of the eight selected miRNAs were ruled out because they had an indistinguishable fluorescence signal. Therefore, five miRNAs were validated by qPCR in the discovery cohort (Supplementary Table S6, shaded in blue). Levels for miR-122-5p, miR-192-5p, miR-885-5p, miR-485-3p, and miR-340-3p were significantly higher in the patients (Fig. 1a), consistent with RNA-seq findings. Using the same approach in the validation (Fig. 1b) and follow-up cohorts (Fig. 1c), miR-122-5p, miR-192-5p, and miR-885-5p exhibited consistently increased levels in patient plasmas. In contrast, levels of miR-485-3p and miR-340-3p were significantly decreased (Fig. 1b, c), leading to inconclusive findings. Finally, no significant differences were observed when cohorts were split by sex (Fig. 2a, c, e), and age did not influence miRNA levels in plasma across the studied cohorts (Fig. 2b, d, f).Fig. 1 Analysis of the miRNA profile. a Discovery, b validation, and c follow-up cohorts. Comparative analysis of miRNA levels between patients and controls, representing IQR (interquartile range) of log2FC relative expression (2−ΔΔCt). ****P value < 0.0001; ***P value < 0.001; **P value < 0.01; *P value < 0.05

Fig. 2 Comparative analysis of the relative expression of miRNA levels by gender from case groups from the a discovery, c validation, and e follow-up cohorts. Correlation matrix including biochemical parameters, age, and levels of miRNAs profile in patients belonging to the b discovery, d validation, and f follow-up cohorts. In each cell, Spearman’s coefficient (Rs) of correlation by pair of variables is indicated. Coloured cells represent significant correlation (*P value < 0.05) in orange when positive, and in blue when negative. ns not significant, AP alkaline phosphatase, AST aspartate aminotransferase, ALT alanine aminotransferase, GGT gamma-glutamyl transferase, L-score Leipzig score

Positive correlations were detected between AST, ALT, and GGT and miR-122-5p, miR-192-5p, and miR-885-5p (Fig. 2b, d, f), supporting their association with liver disease. Triglycerides and/or cholesterol showed a positive correlation with miR-122-5p and miR-885-5p. Additionally, miR-192-5p showed a positive and significant correlation with total bilirubin, although only in the validation cohort. Notably, Rs coefficients between miR-122-5p, miR-192-5p, and miR-885-5p were consistently significant across all cohorts (Rs: 0.48–0.97) (Fig. 2b, d, f), emphasizing their participation in common pathways. Based on all the above findings, the miRNA profile -miR-122-5p, miR-192-5p, and miR-885-5p- emerged as a potential biomarker for WD.

MiR-485-3p did not show significant correlations in any cohort (Fig. 2b, d, f) while miR-340-3p displayed variable correlations with biochemical parameters (Fig. 2b, d). Additionally, neither miR-485-3p nor miR-340-3p exhibited notable associations with the other three miRNAs across the cohorts. Therefore, miR-485-3p and miR-340-3p presented conflicting findings, hindering a definitive conclusion regarding their role as WD biomarkers.

Significant correlation with the Leipzig score was only detected for miR-485-3p in the discovery cohort, and was not replicated in the remaining cohorts of patients (Fig. 2b, d, f).

Profile of circulating miRNAs in serum and liver of a mouse model of WD

The miRNA profile in the Atp7b−/− mice [12–15] was investigated to establish the presentation pattern as the disease progresses. Unfortunately, miR-885-5p lacks its murine orthologue. Since we used human plasma, we first tested miR-122-5p and miR-192-5p in plasma from five WT and WD mice at the age of 21 weeks. Both miRNAs were significantly overrepresented in the WD mice (data not shown). Next, we investigated both miRNAs in serum from WT and WD mice. Significantly increased levels of miR-122-5p and miR-192-5p were detected from the age of 12 weeks to 40 (Fig. 3a, b).Fig. 3 Validation of the miRNAs profile in serum of wild-type (WT) and Atp7b−/− mice (WD, Wilson’s disease) samples. a Comparative analysis of miR-122 levels and, b miR-192 levels between groups at different ages, representing IQR (interquartile range) of log2FC relative expression (2−ΔΔCt). c Comparative analysis of relative expression miRNA levels from Atp7b−/− group from a and b by sex. ****P value 0.0001; **P value < 0.01; *P value < 0.05

We investigated the influence of sex in the presentation of the miRNAs in serum in the Atp7b−/− mice. Significant differences were observed at the ages of 30 and 40 weeks (Fig. 3c), when disease is in advanced state: miR-122-5p and miR-192-5p had an increased presentation in males compared to females.

In WD mice, we found a positive Rs coefficient of the miRNA presentation with ALT and AST (GGT was not available in mice) starting at 16 weeks, with correlations being significant in most estimates (Supplementary Figure S3). In addition, the Rs coefficient between miR-122-5p and miR-192-5p showed a strong association (Rs > 0.58) in all ages.

Lastly, we investigated miR-122-5p and miR-192-5p in liver from mice at 20, 30 and 40 weeks. Both miRNAs presented significantly decreased levels compared to WT, mainly at 40 weeks of age (Fig. 4).Fig. 4 Validation of the liver tissue miRNAs profile in wild-type (WT) and Atp7b−/− (WD) mice samples. a Comparative analysis of miR-122 levels and, b miR-192 levels between groups at different ages, representing IQR (interquartile range) of log2FC relative expression (2−ΔΔCt)

Discussion

Our objective was to assess the utility of plasma miRNAs as biomarkers of WD patients. In human plasma across three patient cohorts, miR-122-5p, miR-192-5p, and miR-885-5p displayed significantly elevated levels compared to healthy controls. In a Atp7b−/− mouse model [12–15], miR-122-5p and miR-192-5p in serum mirrored the pattern observed in patients, while liver levels of both miRNAs decreased, especially at 40 weeks, coinciding with the progression from chronic hepatitis (30 weeks) to cirrhosis (36–44 weeks) [14]. The fluctuation in liver miRNA levels could be attributed partly to hepatocyte regenerative capacity and the release of miRNAs into circulation during disease progression [19].

MiRNAs from the liver may enter the bloodstream passively through apoptosis and necrosis or actively through exosome and viral particle secretion [20]. Therefore, hepatic miRNA levels correlate with serum/plasma levels, offering a non-intrusive method to study liver injury. MiR-122 is associated with hepatocyte differentiation, development, and homeostasis [11]. Inhibition of miR-122 leads to inflammation, necrosis, steatosis and fibrosis [21]. Circulating miRNA-122 is used as biomarker of liver injury in chronic hepatitis B or C, hepatocellular carcinoma, MASH / MASLD, drug-induced liver disease (DILI) and alcohol-related liver disease (ALD) [8, 22–24]. In Long-Evans Cinnamon (LEC) rats used as WD models after induction of fulminant hepatitis with a high copper diet, levels of miR-122 in serum resulted to be significantly increased before the elevation of hepatitis-associated serum markers ALT, AST and bilirubin [19]. Similarly, in a rat model of MASLD, miR-122 was also superior to clinical markers traditionally used to monitor hepatic disease [25]. Moreover, in a study encompassing liver diseases of different aetiology, circulating miR-122 was proposed as a disease severity-dependent biomarker at an early stage because its levels became elevated earlier than that of serum ALT [26]. In the WD model here used [12–15], there was a substantial increase in levels of miRNA-122-5p as well as miRNA-192-5p at week 12 (P value < 0.0001).

The Rs coefficients of miR-122-5p, miR-192-5p, and miR-885-5p demonstrated a significant and strong association, indicating their involvement in common pathways. These miRNAs are frequently found in profiles associated with liver injury [27–31]. Higher serum levels of a miRNA signature (which included miR-885-5p, miR-122-5p and miR-192-5p) have been found to be upregulated in chronic hepatitis B and MASLD [28, 29]. Correlation with ALT, AST and GGT was found for miR-885-5p and miR-122, and even a role regulating ALT activity was proposed for miR-122 [30, 31]. In line with this, our study identified a positive association between the miRNA signature and liver enzymes (ALT, AST, and GGT).

Significant correlation between miRNAs and the Leipzig score was only detected for miR-485-3p in the discovery cohort, discarded as a potential biomarker for WD. The absence of correlation between miRNAs and the Leipzig score could, therefore, be due to different type of data to be considered clinical signs, histopathological studies, biochemical tests and genetic analysis [4].

A role of sex was found in a study with a stronger association in men between the miR-122-5p and miR-885-5p and steatotic liver [31]. In addition, sex differences observed in the Atp7b−/− mice at 30–40 weeks of age, when the liver is cirrhotic, were described to be associated with higher survival rates in females [15]. In our study, no differences were appreciated between women and men, suggesting the need for further exploration in larger patient series.

In conclusion, the increased levels of the miRNA profile miR-122-5p, miR-192-5p and miR-885-5p in plasma of WD patients compared to healthy controls and its association with liver enzymes suggest it might be a potential non-invasive biomarker in this disease. Some of these miRNAs (isolated or with others) have also been reported to be upregulated in other hepatic conditions. A further step would be to assess the role of this miRNA profile as a biomarker determining treatment response. While we did not study this, sensitive biomarkers, such as circulating miRNAs, hold promise for both prognosis of liver disease and assessment of treatment response, vital for identifying optimal therapeutic strategies.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (PDF 721 KB)

Supplementary file2 (DOCX 66 KB)

Acknowledgements

We are in debt with the patients and their families for agreeing to participate in this study. The authors also thank Aída Regi-Cosculluela and Faustino Giménez-Felices, Spanish association of patients with Wilson disease (AEFE Wilson), always willing to bring science closer to patients and their families. This study has been funded by the Fundació Per Amor a l’Art (FPAA); by the Instituto de Salud Carlos III (ISCIII) through the projects PI18/00147 and PI21/00103 to CE, PI21/00313 to OM, and INT20/00061 to MB, and co-funded by the European Union; by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of, making Europe” through the project PID2021-124455OB-I00 to GGA; and by the Centro de Investigación Biomédica en Red (CIBER) through the project CB06/04/0065 to MB. Part of the equipment employed in this work was funded by Generalitat Valenciana and co-financed with ERDF (OP ERDF of Comunitat Valenciana 2014-2020). ER has a PhD fellowship from the Generalitat Valenciana [CIACIF/2021/252].

Author contributions

AS-M and CE designed the study. AS-M, ER, OM and SD carried out the experimental tasks. MÁ-S, EG, IS-B, AB, LG-V, AT and MB recruited and supervised the patients and acquired the clinical data. AS-M, ER and FG-G performed the statistical analyses. AS-M, ER and CE analysed the data. AS-M and CE drafted the manuscript. All authors read and approved the final manuscript.

Declarations

Conflict of interest

G. González-Aseguinolaza is Vivet therapeutics employee and shareholder. No sponsor had a role in the study design, the data collection, the analysis and interpretation of data, the writing of the paper or the decision to submit the article for publication.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ana Sánchez-Monteagudo, Edna Ripollés, Marina Berenguer, Carmen Espinós: Equal contribution.
==== Refs
References

1. Sánchez-Monteagudo A Ripolles E Berenguer M Espinós C Wilson’s disease: facing the challenge of diagnosing a rare disease Biomedicines 2021 9 1100 10.3390/biomedicines9091100 34572285
Sánchez-Monteagudo A, Ripolles E, Berenguer M, Espinós C. Wilson’s disease: facing the challenge of diagnosing a rare disease. Biomedicines. 2021;9:1100.34572285
2. Teschke R Eickhoff A Wilson disease: copper-mediated cuproptosis, iron-related ferroptosis, and clinical highlights, with comprehensive and critical analysis update Int J Mol Sci 2024 25 4753 10.3390/ijms25094753 38731973
Teschke R, Eickhoff A. Wilson disease: copper-mediated cuproptosis, iron-related ferroptosis, and clinical highlights, with comprehensive and critical analysis update. Int J Mol Sci. 2024;25:4753.38731973
3. Collet C Laplanche JL Page J High genetic carrier frequency of Wilson’s disease in France: discrepancies with clinical prevalence BMC Med Genet 2018 19 143 10.1186/s12881-018-0660-3 30097039
Collet C, Laplanche JL, Page J, et al. High genetic carrier frequency of Wilson’s disease in France: discrepancies with clinical prevalence. BMC Med Genet. 2018;19:143.30097039
4. Liver EAftSot. EASL Clinical practice guidelines: Wilson’s disease. J Hepatol. 2012;56:671–685.
5. Espinós C Ferenci P Are the new genetic tools for diagnosis of Wilson disease helpful in clinical practice? JHEP Rep 2020 2 100114 10.1016/j.jhepr.2020.100114 32613181
Espinós C, Ferenci P. Are the new genetic tools for diagnosis of Wilson disease helpful in clinical practice? JHEP Rep. 2020;2: 100114.32613181
6. Sánchez-Monteagudo A Alvarez M Sastre I Genetics of Wilson disease and Wilson-like phenotype in a clinical series from eastern Spain Clin Genet 2020 97 758 763 10.1111/cge.13719 32043565
Sánchez-Monteagudo A, Alvarez M, Sastre I, et al. Genetics of Wilson disease and Wilson-like phenotype in a clinical series from eastern Spain. Clin Genet. 2020;97:758–63.32043565
7. Czlonkowska A Rodo M Wierzchowska-Ciok A Accuracy of the radioactive copper incorporation test in the diagnosis of Wilson disease Liver Int 2018 38 1860 1866 10.1111/liv.13715 29418065
Czlonkowska A, Rodo M, Wierzchowska-Ciok A, et al. Accuracy of the radioactive copper incorporation test in the diagnosis of Wilson disease. Liver Int. 2018;38:1860–6.29418065
8. Szabo G Bala S MicroRNAs in liver disease Nat Rev Gastroenterol Hepatol 2013 10 542 552 10.1038/nrgastro.2013.87 23689081
Szabo G, Bala S. MicroRNAs in liver disease. Nat Rev Gastroenterol Hepatol. 2013;10:542–52.23689081
9. Agbu P Carthew RW MicroRNA-mediated regulation of glucose and lipid metabolism Nat Rev Mol Cell Biol 2021 22 425 438 10.1038/s41580-021-00354-w 33772227
Agbu P, Carthew RW. MicroRNA-mediated regulation of glucose and lipid metabolism. Nat Rev Mol Cell Biol. 2021;22:425–38.33772227
10. Doghish AS Elballal MS Elazazy O The role of miRNAs in liver diseases: Potential therapeutic and clinical applications Pathol Res Pract 2023 243 154375 10.1016/j.prp.2023.154375 36801506
Doghish AS, Elballal MS, Elazazy O, et al. The role of miRNAs in liver diseases: Potential therapeutic and clinical applications. Pathol Res Pract. 2023;243: 154375.36801506
11. Szelenberger R Kacprzak M Saluk-Bijak J Plasma MicroRNA as a novel diagnostic Clin Chim Acta 2019 499 98 107 10.1016/j.cca.2019.09.005 31499022
Szelenberger R, Kacprzak M, Saluk-Bijak J, et al. Plasma MicroRNA as a novel diagnostic. Clin Chim Acta. 2019;499:98–107.31499022
12. Buiakova OI Xu J Lutsenko S Null mutation of the murine ATP7B (Wilson disease) gene results in intracellular copper accumulation and late-onset hepatic nodular transformation Hum Mol Genet 1999 8 1665 1671 10.1093/hmg/8.9.1665 10441329
Buiakova OI, Xu J, Lutsenko S, et al. Null mutation of the murine ATP7B (Wilson disease) gene results in intracellular copper accumulation and late-onset hepatic nodular transformation. Hum Mol Genet. 1999;8:1665–71.10441329
13. Lutsenko S Atp7b-/- mice as a model for studies of Wilson’s disease Biochem Soc Trans 2008 36 1233 1238 10.1042/BST0361233 19021531
Lutsenko S. Atp7b-/- mice as a model for studies of Wilson’s disease. Biochem Soc Trans. 2008;36:1233–8.19021531
14. Murillo O Luqui DM Gazquez C Long-term metabolic correction of Wilson’s disease in a murine model by gene therapy J Hepatol 2016 64 419 426 10.1016/j.jhep.2015.09.014 26409215
Murillo O, Luqui DM, Gazquez C, et al. Long-term metabolic correction of Wilson’s disease in a murine model by gene therapy. J Hepatol. 2016;64:419–26.26409215
15. Murillo O Moreno D Gazquez C Liver expression of a MiniATP7B gene results in long-term restoration of copper homeostasis in a Wilson disease model in mice Hepatology 2019 70 108 126 10.1002/hep.30535 30706949
Murillo O, Moreno D, Gazquez C, et al. Liver expression of a MiniATP7B gene results in long-term restoration of copper homeostasis in a Wilson disease model in mice. Hepatology. 2019;70:108–26.30706949
16. Garcia-Villarreal L Hernandez-Ortega A Sánchez-Monteagudo A Wilson disease: revision of diagnostic criteria in a clinical series with great genetic homogeneity J Gastroenterol 2021 56 78 89 10.1007/s00535-020-01745-0 33159804
Garcia-Villarreal L, Hernandez-Ortega A, Sánchez-Monteagudo A, et al. Wilson disease: revision of diagnostic criteria in a clinical series with great genetic homogeneity. J Gastroenterol. 2021;56:78–89.33159804
17. Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data Bioinformatics 2010 26 139 140 10.1093/bioinformatics/btp616 19910308
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.19910308
18. Gu Y Wei X Sun Y miR-192-5p silencing by genetic aberrations is a key event in hepatocellular carcinomas with cancer stem cell features Cancer Res 2019 79 941 953 10.1158/0008-5472.CAN-18-1675 30530815
Gu Y, Wei X, Sun Y, et al. miR-192-5p silencing by genetic aberrations is a key event in hepatocellular carcinomas with cancer stem cell features. Cancer Res. 2019;79:941–53.30530815
19. Siaj R Sauer V Stoppeler S Longitudinal analysis of serum miR-122 in a rat model of Wilson’s disease Hepatol Int 2012 6 770 777 10.1007/s12072-012-9348-5 23125884
Siaj R, Sauer V, Stoppeler S, et al. Longitudinal analysis of serum miR-122 in a rat model of Wilson’s disease. Hepatol Int. 2012;6:770–7.23125884
20. Novellino L Rossi RL Bonino F Circulating hepatitis B surface antigen particles carry hepatocellular microRNAs PLoS ONE 2012 7 e31952 10.1371/journal.pone.0031952 22470417
Novellino L, Rossi RL, Bonino F, et al. Circulating hepatitis B surface antigen particles carry hepatocellular microRNAs. PLoS ONE. 2012;7: e31952.22470417
21. Sendi H Mead I Wan M miR-122 inhibition in a human liver organoid model leads to liver inflammation, necrosis, steatofibrosis and dysregulated insulin signaling PLoS ONE 2018 13 e0200847 10.1371/journal.pone.0200847 30024933
Sendi H, Mead I, Wan M, et al. miR-122 inhibition in a human liver organoid model leads to liver inflammation, necrosis, steatofibrosis and dysregulated insulin signaling. PLoS ONE. 2018;13: e0200847.30024933
22. Cermelli S Ruggieri A Marrero JA Circulating microRNAs in patients with chronic hepatitis C and non-alcoholic fatty liver disease PLoS ONE 2011 6 e23937 10.1371/journal.pone.0023937 21886843
Cermelli S, Ruggieri A, Marrero JA, et al. Circulating microRNAs in patients with chronic hepatitis C and non-alcoholic fatty liver disease. PLoS ONE. 2011;6: e23937.21886843
23. Waidmann O Bihrer V Pleli T Serum microRNA-122 levels in different groups of patients with chronic hepatitis B virus infection J Viral Hepat 2012 19 e58 65 10.1111/j.1365-2893.2011.01536.x 22239527
Waidmann O, Bihrer V, Pleli T, et al. Serum microRNA-122 levels in different groups of patients with chronic hepatitis B virus infection. J Viral Hepat. 2012;19:e58-65.22239527
24. Starkey Lewis PJ Dear J Platt V Circulating microRNAs as potential markers of human drug-induced liver injury Hepatology 2011 54 1767 1776 10.1002/hep.24538 22045675
Starkey Lewis PJ, Dear J, Platt V, et al. Circulating microRNAs as potential markers of human drug-induced liver injury. Hepatology. 2011;54:1767–76.22045675
25. Yamada H Ohashi K Suzuki K Longitudinal study of circulating miR-122 in a rat model of non-alcoholic fatty liver disease Clin Chim Acta 2015 446 267 271 10.1016/j.cca.2015.05.002 25958847
Yamada H, Ohashi K, Suzuki K, et al. Longitudinal study of circulating miR-122 in a rat model of non-alcoholic fatty liver disease. Clin Chim Acta. 2015;446:267–71.25958847
26. Zhang Y Jia Y Zheng R Plasma microRNA-122 as a biomarker for viral-, alcohol-, and chemical-related hepatic diseases Clin Chem 2010 56 1830 1838 10.1373/clinchem.2010.147850 20930130
Zhang Y, Jia Y, Zheng R, et al. Plasma microRNA-122 as a biomarker for viral-, alcohol-, and chemical-related hepatic diseases. Clin Chem. 2010;56:1830–8.20930130
27. de Paredes AG Manicardi N Tellez L Molecular profiling of decompensated cirrhosis by a novel MicroRNA signature Hepatol Commun 2021 5 309 322 10.1002/hep4.1642 33553977
de Paredes AG, Manicardi N, Tellez L, et al. Molecular profiling of decompensated cirrhosis by a novel MicroRNA signature. Hepatol Commun. 2021;5:309–22.33553977
28. Tan Y Ge G Pan T Serum MiRNA panel as potential biomarkers for chronic hepatitis B with persistently normal alanine aminotransferase Clin Chim Acta 2015 451 232 239 10.1016/j.cca.2015.10.002 26483130
Tan Y, Ge G, Pan T, et al. Serum MiRNA panel as potential biomarkers for chronic hepatitis B with persistently normal alanine aminotransferase. Clin Chim Acta. 2015;451:232–9.26483130
29. Tan Y Ge G Pan T A pilot study of serum microRNAs panel as potential biomarkers for diagnosis of nonalcoholic fatty liver disease PLoS ONE 2014 9 e105192 10.1371/journal.pone.0105192 25141008
Tan Y, Ge G, Pan T, et al. A pilot study of serum microRNAs panel as potential biomarkers for diagnosis of nonalcoholic fatty liver disease. PLoS ONE. 2014;9: e105192.25141008
30. Pirola CJ Fernandez T Castano GO Circulating microRNA signature in non-alcoholic fatty liver disease: from serum non-coding RNAs to liver histology and disease pathogenesis Gut 2015 64 800 812 10.1136/gutjnl-2014-306996 24973316
Pirola CJ, Fernandez T, Castano GO, et al. Circulating microRNA signature in non-alcoholic fatty liver disease: from serum non-coding RNAs to liver histology and disease pathogenesis. Gut. 2015;64:800–12.24973316
31. Raitoharju E Seppala I Lyytikainen LP Blood hsa-miR-122-5p and hsa-miR-885-5p levels associate with fatty liver and related lipoprotein metabolism-The Young Finns Study Sci Rep 2016 6 38262 10.1038/srep38262 27917915
Raitoharju E, Seppala I, Lyytikainen LP, et al. Blood hsa-miR-122-5p and hsa-miR-885-5p levels associate with fatty liver and related lipoprotein metabolism-The Young Finns Study. Sci Rep. 2016;6:38262.27917915
