Seminal extracellular vesicle sncRNA sequencing reveals altered miRNA/isomiR profiles as sperm retrieval biomarkers for azoospermia.

Abstract

Non-invasive molecular biomarkers for classifying azoospermia by origin into either obstructive or non-obstructive/secretory azoospermia, as well as for inferring the spermatogenic reserve of the testis of non-obstructive/secretory azoospermia patients, are of great interest for testicular sperm retrieval outcome prediction for assisted reproduction. Prior analyses of semen small non-coding RNA expression in azoospermia have focused on microRNAs, but there has been a lack of attention on other regulatory small RNA species. In this regard, studying more in-depth expression changes of small non-coding RNA subtypes in small extracellular vesicles from semen of azoospermic individuals could be useful to select additional non-invasive biomarkers with diagnostic/prognostic purposes. A high-throughput small RNA profiling analysis to determine the expression pattern of seminal small extracellular vesicle microRNAs (analyzed at the isomiR level), PIWI-interacting RNAs, and transfer RNA-derived small RNAs in normozoospermic (n = 4) and azoospermic (obstructive azoospermia because of pathological occurring obstruction in the genital tract, n = 4; secretory azoospermic individuals with positive testicular sperm extraction value, n = 5; secretory azoospermic individuals with negative testicular sperm extraction value, n = 4) individuals was carried out. Reverse transcriptase-quantitative real-time polymerase chain reaction validation analysis of selected microRNAs was additionally performed in a larger number of individuals. Clinically relevant quantitative changes in the small non-coding RNA levels contained in semen small extracellular vesicles can be used as biomarkers for the origin of azoospermia and for predicting the presence of residual spermatogenesis. In this regard, canonical isoform microRNAs (n = 185) but also other isomiR variants (n = 238) stand out in terms of numbers and fold-change differences in expression, underlining the need to consider isomiRs when investigating microRNA-based regulation. Conversely, although transfer RNA-derived small RNAs are shown in our study to represent a high proportion of small non-coding RNA sequences in seminal small extracellular vesicle samples, they are not able to discriminate the origin of azoospermia. PIWI-interacting RNA cluster profiles and individual PIWI-interacting RNAs with significant differential expression were also not able to discriminate. Our study demonstrated that expression values of individual and/or combined canonical isoform microRNAs (miR-10a-5p, miR-146a-5p, miR-31-5p, miR-181b-5p; area under the receiver operating characteristic curve >0.8) in small extracellular vesicles provide considerable clinical value in identifying samples with a high likelihood of sperm retrieval while discriminating azoospermia by origin. Although no individual microRNA showed sufficient discriminating power on its own to identify severe spermatogenic disorders with focal spermatogenesis, multivariate microRNA models in semen small extracellular vesicles have the potential to identify those individuals with residual spermatogenesis. Availability and adoption of such non-invasive molecular biomarkers would represent a great improvement in reproductive treatment decision protocols for azoospermia in clinical practice.