Introduction

The identification and rapid worldwide dissemination of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), underpinning the coronavirus disease 2019 (COVID-19), led the World Health Organization (WHO) to declare the status of a pandemic in March 2020 [1]. Since the preliminary epidemiological data has been available, a clear sex disparity has been evident, with males, although not more frequently affected, often experiencing worse outcomes when compared to women [1,2,3]. The reasons underlying this association are probably multifactorial and include both social and cultural factors supporting the higher prevalence of associated morbidities observed in men when compared to women [1]. In addition, the possible contribution of hormonal factors and, in particular, testosterone (T) has also been proposed [1, 2, 4, 5]. The angiotensin converting enzyme 2 (ACE2), and the transmembrane protease, serine 2 (TMPRSS2), both crucial for viral cell entry, are highly expressed in the male genital tract and both modulated by T activity [3, 6, 7]. Limited evidence has emerged regarding the presence of SARS-CoV-2 in the male genital tract, as well as in seminal fluid [7, 8]. Similarly, the reported testis discomfort experienced, at least in a subset of patients during COVID-19 acute phase, appears to support the potential development of SARS-COV-2-related epididymal-orchitis in a number of cases [7]. Similar data has previously been reported for other coronaviruses [7]. Taken together, these observations have clearly emphasized the necessity of a more systematic evaluation of COVID-19 affected men to rapidly identify acute (semen localization and orchitis) as well as possible chronic andrological complications (i.e., infertility and hypogonadism) related to the infection. In particular, the safety and potential consequences of the utilization of reproductive cells from SARS-CoV-2-positive subjects represents a critical challenge to couples and clinicians involved in assisted reproductive care. It is noteworthy that few studies have investigated the short- and long-term andrological consequences of COVID-19. Furthermore, the lack of data derived from phase III trials related to fertility safety of COVID-19 vaccines, along with the aforementioned possible andrological consequences, still represents one the main reasons for accepting vaccination in men seeking fertility treatments [9].

The aim of the present study is to systematically review and meta-analyze all available data regarding possible short- and long-term andrological effects of COVID-19. In particular, the primary outcome is the detection rate (DR) of SARS-Cov-2 mRNA in the male genital tract and semen of infected subjects. Secondary outcomes include the comparison of semen and hormonal parameters between COVID-19 patients and controls. Finally, information regarding the safety of the COVID-19 vaccines on sperm quality was investigated.

Methods

This meta-analysis was performed in line with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) reporting guideline [see Supplementary file 1]. The protocol of this study is published on the website of the University of York (Centre for Reviews and Dissemination) https://www.crd.york.ac.uk/PROSPERO (CRD42021275185).

Search strategy

Two different extensive Medline, Embase and Cochrane searches using MeSH terms were performed.

The first search was focused on the selection of all studies evaluating the presence SARS-Cov-2 mRNA in seminal fluid of infected subjects including the following keywords: (“sarscov 2”[MeSH Terms] OR “sarscov 2”[All Fields] OR “covid”[All Fields] OR “covid 19”[MeSH Terms] OR “covid 19”[All Fields]) AND (“sperm s”[All Fields] OR “spermatozoa”[MeSH Terms] OR “spermatozoa”[All Fields] OR “sperm”[All Fields] OR “sperms”[All Fields]).

The second search was mainly focused on the selection of all studies comparing hormonal and seminal sperm parameters in infected subjects compared to controls, including the following keywords: (“covid 19”[All Fields] OR “covid 19”[MeSH Terms] OR “covid 19 vaccines”[All Fields] OR "covid 19 vaccines”[MeSH Terms] OR “covid 19 serotherapy”[All Fields] OR “covid 19 serotherapy”[Supplementary Concept] OR “covid 19 nucleic acid testing”[All Fields] OR “covid 19 nucleic acid testing”[MeSH Terms] OR “covid 19 serological testing”[All Fields] OR “covid 19 serological testing”[MeSH Terms] OR “covid 19 testing”[All Fields] OR “covid 19 testing"[MeSH Terms] OR “sarscov 2”[All Fields] OR “sarscov 2”[MeSH Terms] OR “severe acute respiratory syndrome coronavirus 2”[All Fields] OR “ncov”[All Fields] OR “2019 ncov”[All Fields] OR ((“coronavirus”[MeSH Terms] OR “coronavirus”[All Fields] OR “cov”[All Fields]) AND 2019/11/01:3000/12/31[Date—Publication])) AND (“testosterone”[MeSH Terms] OR “testosterone”[All Fields] OR “testosterone”[All Fields] OR “testosterones”[All Fields] OR “testosterone s”[All Fields]).

Data from January 1 2020 up to August 31 2021 were restricted to English-language articles and studies including human participants. The identification of relevant studies was performed independently by three of the authors (F.P, W.V and A.P), and conflicts were resolved by the first investigator (G.C). All the data identified during the first analysis were checked in a second wave analysis by three of the authors (G.R, D.P and E.B). Possible further conflicts were discussed and resolved by the first investigator (G.C). We did not employ search software but hand-searched the bibliographies of retrieved papers for additional references. All the authors adequately contributed to the analysis of the paper and reviewed the final version of the manuscript. The main source of information was derived from published articles.

Study selection

All prospective and retrospective observational studies reporting information on SARS-Cov-2 mRNA semen and male genital tract detection of COVID-19 subjects, as well as those reporting data on semen analysis and hormonal parameters in infected/recovered patients as compared with controls (when available), without any arbitrary restriction were included (see Supplementary Fig. 1 and Table 1) [4, 5, 8, 10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38]. In addition, studies comparing semen parameters before and after COVID-19 vaccination were also included in the analysis [39, 40].

No country restrictions were applied. Case reports were excluded from the analysis (see Supplementary Fig. 1). Studies not specifically reporting on at least hormonal or sperm parameters were also excluded from the analysis.

Outcome and quality assessment

The primary outcome was the DR of SARS-Cov-2 mRNA in the male genital tract and semen of infected subjects. Secondary outcomes included the comparison of semen and hormonal parameters between COVID-19 patients and controls. The effect of several risk factors—including age, associated morbidities, time from diagnosis as well as disease severity—on DR were investigated. Similarly, the impact of disease duration after the appearance of first symptoms and type of control populations for semen and hormonal parameters were also analyzed. Finally, the possible impact of COVID-19 vaccination in healthy controls was evaluated. The quality of trials included was assessed using the Cochrane criteria [41]. In particular, we considered the following criteria: the weaknesses of the designs that have been used (such as noting their potential to ascertain causality), the execution of the studies through a careful assessment of their risk of bias, especially the potential for selection bias and confounding to which all observational studies are susceptible, and the potential for reporting biases, including selective reporting of outcomes.

Statistical analysis

Heterogeneity in DR was assessed using I2 statistics. Even when low heterogeneity was detected, a random-effect model was applied because the validity of tests of heterogeneity can be limited with a small number of component studies. We used funnel plots and the Begg adjusted rank correlation test to estimate possible publication or disclosure bias [42]; however, undetected bias may still be present, because these tests have low statistical power when the number of trials is small. Overall DR is expressed as mean percentage (95% confidence interval).

In addition, a meta-regression analysis was performed to test the effect of different parameters on SARS-COV-2 mRNA DR. Following on that, potential predictors of DR were included as continuous variables: age, time from diagnosis of the disease, associated morbidities (including hypertension, diabetes mellitus), as well as disease severity. All data were calculated using Comprehensive Meta-analysis Version 2, Biostat (Englewood, NJ, USA).

Results

General descriptive data

Out of 93 and 111 retrieved articles, 35 were included in the study (Table 1). Among them, 31 prospectively investigated different outcomes, whereas only four retrospective data analyses had been performed in the remaining studies (Table 1). The study flow is summarized in Supplementary Fig. 1. The characteristics of the retrieved trials (including parameters on trial quality) are reported in Table 1 and Supplementary Tables I–II. Retrieved trials included 2092 patients and 1138 controls. Mean age of included population was 44.1 ± 12.6 years, and the mean follow-up 24.3 ± 18.9 days. Finally, two studies, including overall 88 subjects (mean age 38.6 ± 6.4 years), investigated the effect of COVID-19 vaccination on semen parameters (Table 1).

Table 1 Characteristics of trials included in the meta-analysis

Male genital tract SARS-Cov-2 mRNA detection rate

Among the retrieved trials, 19 studies analyzed the possible presence of SARS-Cov-2 mRNA in male genital tract (DR) samples. Among the available studies, 17 investigated the possible presence of SARS-Cov-2 mRNA directly in the semen, whereas two evaluated viral presence from testis autopsies of deceased subjects due to COVID-19 (see Supplementary Table 1). The I2 was 0; p = 0.665. Mean crude DR was 8 [5; 12]% (Fig. 1). A funnel plot and Begg adjusted rank correlation test (Kendall’s τ: 0.09; p = 0.576) suggested no major publication bias (Supplementary Fig. 2). Similar data were found when retrospective studies were compared to prospective ones (DR = 9 [6; 14]% vs. 2 [1; 7]%, respectively; Q = 0.05, p = 0.82) or after the exclusion of those studies evaluating the presence of SARS-Cov-2 in testis autopsy (7 [5; 11]%). Conversely, the DR was significantly lower when low quality studies were excluded from the analysis (DR = 8 [5; 12]% vs. 9 [3; 26]%, respectively; Q = 5,27, p = 0.02).

Fig. 1
figure 1

Male genitalia tract mRNA SARS-Cov-2 detection rate (%)

Meta-regression analysis showed that DR was not influenced by patient age or by disease severity or associated morbidities (Fig. 2A–D). Conversely, DR was significantly and inversely related to COVID-19 diagnosis timing (Fig. 2E). Accordingly, the DR was significantly higher in those studies assessing the viral mRNA presence in the semen less than 11 days after the diagnosis (Q = 5.611; p = 0.018; see also Supplementary Fig. 3, panel A). The latter was confirmed even when those studies evaluating the presence of SARS-Cov-2 in testis autopsy were excluded from the analysis (Q = 5.951; p = 0.015, Panel B).

Fig. 2
figure 2

Influence of age (A), disease severity (B), arterial hypertension (C), diabetes mellitus (D), and time from diagnosis (E) on male genitalia tract mRNA COVID-19 detection rate. The size of the circles indicates sample dimension

Semen parameters

Five trials evaluated the effect of COVID-19 on semen parameters, compared to healthy controls. In addition, one study compared the same data between hospitalized and not hospitalized infected patients (Supplementary Table 1). When all studies were considered, COVID-19 was associated with a significant reduction of total sperm count, sperm concentration and total sperm motility (Fig. 3A–C), whereas no difference in sperm morphology or progressive motility was observed (not shown). In addition, a lower seminal volume in men with COVID-19 was also detected, when compared to controls (Fig. 3D). Similar results were observed when the only study not considering healthy controls [25] was excluded from the analysis (Supplementary Fig. 4 Panels A–C). Finally, the exclusion of those studies considering subjects with positive oropharyngeal swab [11, 13, 36] from the analysis only partially attenuated the observed negative results (Supplementary Fig. 5, Panels A–C).

Fig. 3
figure 3

Semen parameters in COVID-19 subjects as compared to controls: total sperm count (A), sperm concentration (B), sperm total motility (C), semen volume (D)

Hormonal parameters

Among the included studies, 11 trials evaluated the effect of COVID-19 on different hormonal parameters, including FSH (n = 8), LH (n = 9) and total testosterone (n = 11) (See also Table 1 and Supplementary Tables I–II). In addition, four and five studies investigated, at baseline, the relationship between low total T level and worse COVID-19 outcomes, including probability of being admitted to ICU and mortality risk, respectively (See also Supplementary Table 1).

SARS-Cov-2-infected patients were characterized by reduced total T levels, whereas no difference in either LH or FSH levels was observed (Fig. 4A–C). Similar results were obtained when patients with severe diseases were compared to those with milder forms (Supplementary Fig. 6 Panels A–C).

Fig. 4
figure 4

Hormonal parameters in COVID-19 subjects as compared to controls: total testosterone (A), follicular stimulating hormone (FSH; B), luteinizing hormone (LH; C)

By performing sensitivity analyses, the negative effects of SARS-Cov-2 infection on T levels were confirmed when only those studies that included patients in the acute phase were investigated (mean difference in total T levels − 2.19 [− 7.08; − 1.20] nmol/l; p = 0.01) but not in the only study [36] including subjects in the recovering phase (mean difference in total T levels − 2.19 [− 7.80; 3.41] nmol/l; p = 0.44) (See also supplementary Fig. 7 Panel A). Similar to what was observed for seminal parameters, the exclusion of those studies not considering healthy controls only partially attenuated the results (see also Supplementary Fig. 7, Panel B). Similar data were observed considering LH and FSH levels (not shown).

Finally, when the effects of reduced T levels at baseline on COVID-19 outcomes were investigated, low T resulted in up to four- and fivefold increased risk to be admitted to the Intensive Care Unit (ICU) or to die, after the adjustment for confounders (Fig. 5A and B).

Fig. 5
figure 5

Fully adjusted risk to be admitted to Intensive Care Unit (A) or to die (B) due to COVID-19 according to low baseline T levels

Vaccination and semen parameters

Two studies evaluated the possible impact of COVID-19 vaccination on semen parameters in healthy volunteers. Among the subjects included, the vast majority received BNT162b2 mRNA Covid-19 vaccine whereas a small fraction received mRNA-1273 COVID-19 vaccine (see Supplementary Table 1). No significant effect of vaccination in sperm motility and semen volume was observed (Supplementary Fig. 8, Panels B and C), whereas a positive effect on sperm concentration was detected (Supplementary Fig. 8, Panel A). No further data for analysis of other semen parameters were available.

Discussion

Our data suggest that SARS-CoV-2 mRNA can be detected in the semen only in the acute phase of COVID-19 infection but data are not conclusive. The possible testicular localization of the virus can result in impaired sperm and T production. Low T observed in the acute phase of the disease is intimately associated with an increased risk of worsening outcomes. Finally, the use of mRNA COVID-19 vaccines does not affect semen parameters.

Paoli et al. [43] did not find SARS-CoV-2 mRNA in seminal fluid and urine in 31-year-old men in the recovered phase of the infection. The same authors confirmed this result in two mild COVID-19 patients in the acute phase with a positive nasopharyngeal swab and in recovered patients with a negative swab [44]. Similar results were reported by Pan et al. [14]. The Italian Society of Andrology and Sexual Medicine (SIAMS), based on the available evidence, produced a specific Position Statement not supporting the presence of SARS-COV-2 in semen and the risk of sexual transmission [6]. Following these observations, several authors have further investigated the presence of SARS-CoV-2 in semen but to date only two studies documented semen viral mRNA detection in a relevant number of patients [8, 34]. Li et al. [8] in May 2020 reported for the first time the detection of SARS-CoV-2 in the semen in a limited number of patients after a relative short interval from the onset of COVID-19 symptoms. The latter study was presumably conducted in patients with a severe form and during the acute phase of the disease. This condition may correspond to a higher blood viral load and, thus, a higher chance of reaching other organs and body fluids including semen; at the same time, a confined environment with severe cases of viral disease is more likely to be susceptible to contamination. Likewise, the study of Saylam et al. [34] found four patients with SARS-CoV-2-positive semen during the viremia period of the disease. However, they were recovered in an intensive care unit and, despite a positive swab, little evidence was provided for the presence of viremia. The presence of fever and the severity of the infection have been proposed as possible facilitatory mechanisms for detection of the virus in the semen [7]. Our data indicates that the stage of the disease is the only positive predicting factor. It can be speculated that the seminal identification of the virus in the early phase of the infection can be the consequence of the alteration of Blood–Testis Barrier or secondary to its excretion in the seminal fluid. The majority of available studies provide limited information regarding the method of semen collection and preparation. Hence, the possibility that SARS-CoV-2 semen localization could reflect possible contamination from the feces, urine, hands or respiratory droplets cannot be excluded. Accordingly, DR was higher in low quality studies when compared to higher quality reports. In line with the latter hypothesis, a recent study, applying a new RT-PCR, specifically validated to detect SARS-CoV-2 RNA in sperm, did not document any positive test in patients evaluated from 6 to 181 days after a positive SARS-CoV-2 nasopharyngeal PCR test [45]. Similarly, SARS-CoV-2 mRNA was not found in 111 semen samples cryopreserved in the Florentine bank [46].

Sperm quality was significantly impaired in men with COVID-19 compared to controls. The latter results were essentially confirmed when only recovered subjects were considered. Several direct and indirect mechanisms have been proposed to explain this negative impact on sperm production. The ACE2 receptor has been documented on seminiferous duct cells including the spermatogonia, Leydig cells (LC) and Sertoli cells (SC). Although the specific role of the renin–angiotensin system in testicular function regulation has not been completely elucidated, ACE2 has been proposed as playing a putative role in steroidogenesis regulation. In addition, angiotensin II can impair sperm motility and fertility through type 1 and type 2 receptors (AT1 and AT2) leading to cell apoptosis and senescence.

Conversely, data on the expression of TMPRSS2 in the male genitalia tract are more conflicting [47]. The lack of co-expression of ACE2 and TMPRSS2 on sperm cells has been suggested to reduce the risk of viral direct damage [48]. However, it is important to recognize that, besides ACE2 and TMPRSS2, several other proteins have been shown the capacity to interact with SARS-CoV-2 during virus cell internalization. Among them, the human tyrosine-protein kinase receptor (AXL), a member of the TAM receptor family, highly expressed in Sertoli and Leydig cells can facility SARS-CoV-2 entrance into human pulmonary epithelial cells in an ACE2-independent manner [49]. Hence, after testis localization, SARS-CoV-2 can potentially enter into LC and SC, contributing to impaired T and sperm production. In addition, the inflammatory response due to viral localization within the testis can induce the development of an intense local immune reaction, supporting the development of a viral orchitis, potentially evolving into a vasculitis or to an autoimmune response [3, 7]. In addition, although never directly demonstrated, available evidence indicates that spermatozoa are able to bind to SARS-CoV-2 supporting a potential direct viral insult to spermatozoa [7]. Furthermore, local inflammation and an increased concentration of seminal leucocytes might increase oxidative stress, leading to sperm DNA breakage [7]. Finally, drugs used for COVID-19 treatment including antivirals, antibiotics and steroids can all negatively affect sperm as well as T production [7]. Despite this evidence, a recent study, including 120 SARS-CoV-2-infected subjects tested after a mean of 52.7 days from COVID-19 infection, showed that sperm parameters progressively improved according to the time lapse since virus infection [45]. Data derived from the same study also documented a strong inverse correlation between semen parameters and COVID-19 immune reaction whereas no correlation with fever was observed [45]. Taken together the available data indicate an immunologic rather than a fever-derived sperm damage. In addition, long-term sperm quality recovery was also reported [45]. Present data cannot adequately clarify the latter point due to limited follow-up of the included studies.

LC damage can contribute to impairment in T production, as documented in the current study. Ma et al. [29] demonstrated significantly higher serum LH levels and a lower T to LH ratio (T/LH) in COVID-19 patients compared to controls. Similarly, Rastrelli et al. [4] in a small sample of subjects with COVID-19 reported that a worsening clinical status was associated not only with reduced T levels but also with an increased serum LH concentration. Conversely, Salonia et al. [5] showed that secondary hypogonadism was more frequently observed in SARS-CoV-2-infected patients compared to controls. The development of secondary hypogonadism is a frequent complication of several acute and chronic illnesses. This might represent a protective mechanism, by switching off T-dependent functions (such as reproduction and/or physical and sexual activity) that are not desirable when the physical condition is ailing [50]. Present data seem to support the latter hypothesis. A secondary or at least a mixed hypogonadism is associated with COVID-19 during the acute phase, which does not persist in the recovery phase. However, Salonia et al. [51] in a 7-month follow-up of their data, reported that although total T levels increased over time, more than 50% of men who recovered from the disease still had circulating T levels suggestive of a condition of hypogonadism. Despite these considerations, we here report that low T observed at hospital admission was associated with worse disease outcomes. Similarly, low endogenous T has been associated with an increased risk of overall mortality as well as cardiovascular (CV) mortality and morbidity in the general population [52]. The role of low T in the stratification of CV risk is still conflicting. Several associated morbidities including metabolic diseases [53, 54], drugs [55], HIV [56] as well as heart failure, obstructive pulmonary disease, chronic kidney diseases, and bowel inflammatory diseases have been associated with reduced T levels [57]. Salonia et al. [51] showed that the higher the burden of comorbid conditions at presentation, the lower the probability of T level recovery over time after COVID-19. The real significance of T replacement therapy (TRT) on long-term outcomes in the latter conditions is still conflicting and the optimization of the disease status, as well as weight loss, in the case of obesity, or drug withdrawal when possible, can ameliorate T levels [57, 58]. Functional hypogonadism is the emerging term frequently used to describe the latter cases [57]. Preset data, however, indicated that low T is associated with an increased risk of worsening COVID-19 outcomes, even after the adjustment for possible confounders, supporting a primary pathogenetic role of hypogonadism. Accordingly, T has been demonstrated to have an anti-inflammatory effect in either preclinical [59, 60] or clinical studies [61]. Unfortunately, no data on the effect of TRT in COVID-19 are available [4]. In addition, genetic polymorphism of androgen receptors and ACE2 may modulate susceptibility to SARS-COV-2 infection and COVID-19 outcomes in males [62].

Phase III clinical trials showed that all the available COVID-19 vaccines had an acceptable efficacy/safety profile [63]. Up to January 28rd, 2022, a total of 9.854.237.363 vaccine doses had been administered worldwide (https://covid19.who.int/) with and excellent safety/efficacy profile [64]. A recent study performed in the US indicates that COVID-19 vaccines prevented more than 139,000 deaths during the first five months that they were available with an estimated economic value ranging from $625 billion and $1.4 trillion [65]. Since reproductive toxicity was not evaluated in the clinical trials, one of the reasons for vaccine hesitancy is represented by the potential negative effects on fertility. Gonzalez et al. [39] reported the first evidence of the safety of mRNA-derived COVID-19 vaccines on sperm quality in a small sample of volunteers. Similar results were reported in the same period by Safrai et al. [40]. By meta-analyzing this data, we can confirm the safety of COVID-19 vaccination on sperm parameters. It should be recognized that only mRNA-derived vaccines were used in the available studies, supporting the necessity of further evidence with the use of other types of vaccines to confirm these results.

Several limitations should be recognized. The vast majority of the meta-analyzed data are derived from observational studies, which present an important risk of bias due to the lack of completeness of follow-up and the accrual of missing data. A further limitation is that estimating reproduction numbers for SARS-CoV-2 presents challenges due to the high proportion of infections not correctly detected by health systems, due to paucity or even the lack of symptoms and to changes in testing policies, thus resulting in different proportions of infections being detected over time and between countries. Hence, the reproducibility of our data warrants caution. The detection of RT-PCR reflects only the presence of the viral RNA and does not necessarily indicate the presence of the virus. The follow-up period is limited. Hence, long-term andrological effects of COVID-19 infection should be investigated though larger and longer studies. Similarly, the impact of associated morbidities must be better evaluated. Finally, the possible use of anti-androgen drugs could have represented a further source of bias. However, no sufficient information was available in the considered studies.

In conclusion, present data show that COVID-19 is associated with short-term sperm and T impairment, whereas long-term consequences have still not been sufficiently clarified. Limited evidence indicates no effects of COVID-19 vaccination on sperm quality. Hence, each patient recovered from COVID-19 should be accurately monitored to rule out sperm and T abnormalities. Specific contribution of reduced T levels during the acute phase of the infection still needs to be better clarified.