Introduction

The worldwide outbreak of the new disease called coronavirus disease 2019 (COVID-19) arises from the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has reached a pandemic condition in March 2020 [1, 2]. Existing data have indicated that SARS-CoV-2 is not usually limited to the respiratory system alone but might also attack other tissues in the body [3]. Reportedly, SARS-CoV-2 can also affect the male reproductive system, so it is observed in the semen and leads to damage in the spermatogenesis process [4, 5]. In the seminiferous tubules of the mammalian testis, Sertoli cells by providing nourishment and mechanical support, play a pivotal role in post-meiotic germ cells’ development during spermatogenesis. These cells also construct the blood–testis barrier (BTB) in the seminiferous epithelium that divides the seminiferous epithelium into the adluminal and the basal compartments [6], so there are spermatogonia and preleptotene spermatocytes in the basal compartment and other primary and secondary spermatocytes, round spermatids, and elongated spermatids in the adluminal compartment [7]. Ultrastructurally, the BTB consists of co-existing tight junctions (TJs), desmosome-like junctions, ectoplasmic specializations (ESs) and gap junctions which collectively play a crucial role in the maintenance of BTB integrity that is essential for spermatogenesis and male fertility [8]. In the testis, TJs are the most important component of the blood–testis barrier [7]. Claudins are transmembrane proteins which along with occludin, are essential for TJ formation [9]. Gap junctions and their constitutive proteins, connexins are another type of junction present in BTB. The most predominant connexins which are mainly found at the BTB are connexins-43, which probably leads to the synchronization of Sertoli cells [10].

BTB works not only as a physical barrier but rather allows the spermatogenesis process to take place in an immune-privileged environment and helps to inhibit testicular inflammation and immune response. But some viruses, such as the Mumps virus, Zika virus and HIV, can still penetrate the barriers and be seen in semen and induce inflammation, apoptosis and oxidative damage to the testis. Available data suggest that SARS-CoV-2 causes the release of several pro-inflammatory cytokines, such as interleukin 6 (IL-6), IL-1 and tumor necrosis factor-α (TNF-α) [11]. As reported in previous studies, inflammation is associated with disruption of spermatogenesis and testicular function [12].

All available studies indicate that COVID-19 could influence male reproductive health and it is suggested that the effect of SARS-CoV-2 in the male reproductive system should be considered [13]. Accordingly, this study was designed to elucidate the impacts of SARS-CoV-2 on the male reproductive system by investigating damage of BTB and histological change in human testis.

Materials and methods

Patient selection and patient criteria

Ten patient samples and ten healthy samples were included in the study (Table 1). For the control group, people in the age range (45–58 years) and reproductive health status who died due to accident, electric shock, or carbon monoxide poisoning and had no disease were used. The control group consisted of male adults without any exposure to SARS-CoV-2 infection. To prevent post-mortem delay, all control and COVID-19 groups were analyzed within a minimum period of time after death. None of the control group had a history of infertility or reproductive disorders. In addition, all of them had children.

Table 1 Clinical characteristics of COVID-19 patients

Postmortem examination and sampling of the testes

Postmortem human testes were collected from the Iranian Legal Medicine Organization. This study was conducted by the principles of the Declaration of Helsinki and the guidelines of the Chinese National Health Commission and was approved by the Ethics Committee at Shahid Beheshti University of Medical Sciences, (IR.SBMU.MSP.REC 1399.779). All cases were checked for their medical history of death, death certificate, autopsy report. Cases with a history of reproductive system disorders were excluded from the present study. According to the standard protocol, COVID-19 diagnosis was confirmed by positive nucleic acid testing of nasopharyngeal swabs, radiography (chest X-ray), computer tomography (CT) scan features of viral pneumonia, and clinical symptomatology. The deceased COVID-19 and control cases were transported by refrigerator from the hospital to forensic medicine. Postmortem examinations were carried out after informed consent from family members and were performed approximately 8 to 10 h after death. For ten deceased COVID-19 and ten control, a tissue sample of 3 × 3 cm2 was obtained via incisional autopsy for histological, cellular and molecular study.

Tissue preparation

Lung specimens were fixed in 10% formalin for 72 h and submitted for standard tissue processing, sectioning and hematoxylin and eosin (H&E) staining for histopathology evaluations. Evaluation of tissue sections was done by an expert histologist. The testes samples were kept in Bouin's for 48 h and then transfer to formalin 10% before placing them in paraffin blocks. We then made serial Sections 5 µm and 20 µm using a microtome (Leica RM2125 RTS, Germany) consistent with the stereological techniques. Systematic Uniform Random Sampling (SURS) was used to choose 10 sections in each sample by picking a random figure in the range of 1 to 10. They were then stained using H&E staining (Sigma, USA). It is worth noting that the testis cells were distinct in terms of morphology.

Analysis of TNF-α, IL1β, IL6, occludin, claudin-11 and connexin-43 expression using real-time PCR

After extracting the total RNA samples, they were treated with DNase I (Roche, Basel, Switzerland) to eliminate genomic DNA contamination. We used a commercial kit (Fermentas, Lithuania) to synthesize cDNA at 42 °C for 60 min in compliance with the protocols described in the manufacturer’s instructions. We used real-time PCR (TaqMan) based on QuantiTect SYBR Green RT-PCR kit (Takara Bio Inc, Japan) to quantify relative gene expression. All pairs of forward and reverse primers were designed by Primer 3 Plus software in an exon–exon junction manner to separate cDNA from genomic DNA. Before that, we tested PCR primers using the Primer-Blast tool available at the website, www.ncbi.nlm.nih.gov/tools/primer-blast (Table 2).

Table 2 Primers design

Number of testis cells

The testis cells’ number was estimated using the optical dissector method [14]. We computed the numerical density (Nv) of testis cells using the equation below:

$$N_{v} = \frac{\sum Q}{{\sum P \times h \times \frac{a}{f}}} \times \frac{t}{{{\text{BA}}}}$$

where (ΣQ) is the number of cells and (ΣP) is the number of counting frame grid in all fields; (a/f) is the area of frame; (h) is the dissector height; (t) is the real thickness of the section and (BA) is the microtome section thickness. The total number of testis cells was obtained using the equation below:

$$N_{{{\text{total}}}} = N_{v} \times V$$

Coefficient of error (CE)

CE (V) was computed according to the following equation:

$${\text{CE}}\left( V \right) = \left( {\Sigma P^{ - 1} } \right) \times \left[ {\frac{1}{240}\left( {3\Sigma P_{i}^{2} + \Sigma P_{i} P_{i + 2} - 4\Sigma P_{i} P_{i + 1} } \right)} \right]^{{{1}/{2}}}$$

where (B) and (A) are section boundary length and the section area, respectively. We calculated CE for the volume of testis and the number of testis cells using the equation below [14]:

$${\text{CE}}\left( N \right) = \left[ {{\text{CE}}^{2} \left( {N_{v} } \right) + {\text{CE}}^{2} \left( V \right)} \right]^{{{1}/{2}}}$$
$${\text{CE}}\left( {N_{v} } \right) = \left[ {\left( {\frac{n}{n - 1}} \right) \times \left[ {\left( {\frac{{\Sigma \left( Q \right)^{2} }}{{\left( {\Sigma Q} \right)^{2} }}} \right) + \left( {\frac{{\Sigma \left( P \right)^{2} }}{{\left( {\Sigma P} \right)^{2} }}} \right) - \left( {\frac{{2\Sigma \left( {QP} \right)}}{\Sigma Q\Sigma P}} \right)} \right]} \right]^{{{1}/{2}}}$$

CEs are depicted in Table 3.

Table 3 Coefficients of error (CE) for total number of sertoli cells

Covariance function

The covariance function was measured using the following equation:

$$C\left( r \right)X = \frac{{\sum {\text{DP}}\left( {{\text{XY}} r} \right)}}{{\sum {\text{DP}}\left( {{\text{ref}} r} \right)}}$$

The both end points of dipoles (DP) of class size r = 1 (equivalent to 4.3 µm). To estimate “Vv”, “C(r)”, and “g(r)”, the distance between the points (DP) ranged from r = 0 (equivalent to 0 µm) to r = 49; so, the total distance was 127.4 µm (49 × 4.3 = 210.7 µm).

Pair correlation function

The pair-correlation function is the normalized covariance function obtained by dividing the covariance by the reference value (squared volume fraction):

$$g\left( r \right) = \frac{C\left( r \right)}{{Vv2}}$$

Voronoi tessellation for distribution of the testicular cells

Area of space with closely spaced cells. Thus, the polygon region represents the space occupied by a cell. Testicular cells were mapped using the Image J Voronoi Plugin, which involves drawing a polygonal area around each cell. The area and number of closest Voronoi polygons were obtained from different parts of the testis using the microscopic images of testicular tissue, an objective lens of 40 × was obtained. The data were then analyzed using "ImageJ". The variability of polygonal regions was easily analyzed by their variance. The coefficient of variation or CV (standard deviation of polygonal regions/mean × 100) provides an indicator for the spatial distribution of Sertoli cells: CV 33–64% is associated with a random distribution of testicular cells. CVs less than 33% have a regular pattern and those with more than 64% are considered a cluster distribution. This CV is a classification and not a statistical comparison.

Immunohistochemical staining

After deparaffinization and rehydration, the sections were retrieval with 10 mM citrate buffer and then incubated with primary antibodies (Sigma-Aldrich Corporation, St. Louis, Missouri), (diluted 1 in 100). Then, the appropriate secondary IgG antibody (diluted 1 in 200), the sections were incubated with 3'-diaminobenzidine (DAB) (Dako, Glostrup, Denmark). After the immunohistochemical reaction, sections were counter-stained with hematoxylin (Merck, White house, NJ) mounted, and observed under a light microscope.

Statistical analyses

In this study, we used the Pfaffl method and t tests for statistical analysis. A significant level of p ≤ 0.05 was assumed for all data.

Results

Histopathological feature of lung tissue

Histopathological examination of the pulmonary system showed a spectrum of diffuse alveolar damage in patients, which was evidenced by the presence of intra-alveolar fibrin, hyaline membranes in the alveolar septal walls. Airways and alveolar spaces contained large, reactive mononuclear inflammatory cells. Microscopic hemorrhage was identified with diffuse alveolar damage in patients. Most patients showed variable degrees of chronic interstitial inflammation, with some having more prominent perivascular lymphocytic inflammation (Fig. 1A and B).

Fig. 1
figure 1

The effect of SARS-CoV-2 on the lung tissue in the human post-mortem. A In the control group, the lung parenchyma, including the alveolar sac and alveolar walls, has a normal appearance. B In the COVID-19 group, infiltration of lung tissue by mononuclear inflammatory cells (black arrow), along with desquamation of alveolar epithelium and formation of hyaline membrane together alveolar wall thickening was observed (red arrow)

SARS-CoV-2 increased the expression levels of genes involved in inflammation

We normalized and quantified the relative mRNA expression levels of TNF-α, IL1β and IL6, in study groups. In comparison to the control group, the transcripts for TNF-α, IL1β and IL6 in the COVID-19 group were significantly upregulated (p < 0.05, p < 0.01 and p < 0.01, respectively) (Fig. 2A).

Fig. 2
figure 2

The effect of SARS-CoV-2 on the inflammatory cytokines and junctional proteins in the human post-mortem testicular tissue. A and B Real-time PCR analyses of testes. mRNA expression levels of TNF-α, IL1β, IL6, claudin-11, occludin and connexin-43 from control and COVID-19 groups. Mean ± SD of the mRNA expression levels of TNF-α, IL1β, IL6, claudin-11, occludin and connexin-43 of testis in the study groups (*p < 0.05 and **p < 0.01)

SARS-CoV-2 decreased the expression levels of genes involved in BTB

To investigate the impact of SARS-CoV-2 on the impairment of BTB at the molecular level, the number of transcripts for the three genes contributing to the structure of BTB (occluding, claudin and connexin-43) was analyzed. Furthermore, gene expression for occluding, claudin and connexin-43 was significantly (p < 0.05) lower in the COVID-19 group compared to that of the control group (Fig. 2B).

SARS-CoV-2 decreased in the total number of Sertoli cells

At the cellular level, the total number of Sertoli cells using optical dissector was estimated. The assessment of testis tissue indicated a significant decrease in the number of Sertoli cells in the COVID-19 group compared to the normal control group (p < 0.001) (Fig. 3A–C).

Fig. 3
figure 3

The effect of SARS-CoV-2 on the Sertoli cells number in the human post-mortem. A The total number of Sertoli cells in the COVID-19 group decreased in comparison with the control groups. (***p < 0 .001). Photomicrograph of the testis in control (B) and COVID-19 (C) groups stained with H&E, × 40. SC (Sertoli cell)

SARS-CoV-2 alters the spatial arrangement of Sertoli cells

Assessments of g(r) for Sertoli cells and the dipole distances, were plotted versus each other (Fig. 4A and B). The estimated values from the start to the end of the curve (from r = 0 to 64.5 µm), (from r = 94.6 to 103.2 µm), (from r = 124.7 to 172 µm), (from r = 184.9 to 215 µm), revealed a significant difference among the two groups. After the gap, there are the data points for both the normal control and COVID-19 groups that were arranged randomly in longer distances (p < 0.05) (Fig. 4A and B).

Fig. 4
figure 4

The effect of SARS-CoV-2 on the spatial arrangement of Sertoli cells in the human post-mortem using the second-order stereological method. A Communication between pair-correlation function and dipole distance for Sertoli cells in control and COVID-19 groups. The mean points are g(r) in the control and COVID-19 groups. The horizontal reference axis relates to values expected for a random spatial arrangement [g(r) = 1]. B A micrograph that encodes cells to produce a matrix of a transparent lattice of points to serve as a set of dipole probes. Each row consists of 50 points and 49 equal distance intervals. For each testis, 11 tests were performed (a total of p = 550 test points). For each experiment, the nature of the tissue constituent at each test point was determined and all data were recorded in 11 × 50 matrices. The point interval (r) relates to a distance of 4.3 μm, which has a chance of being in the same Sertoli profile. If the dot was placed on Sertoli cells and connective tissue, 1 and 2 were encoded, respectively. C Convert data matrix to image type and display Sertoli cells (1) yellow and connective tissue (2) dark blue. A control group; Group B, COVID-19

Convert a data matrix into the image using MATLAB software

DM changed the testis cells’ spatial distribution, leading to the dissociation of these cells in some places. To give a clear picture of this event, we converted the data matrices of both COVID-19 and normal control groups into images by MATLAB software (Fig. 4C). Photomicrograph of the testis stained with H&E and cartoon images showed the changes of tissue structure and spatial pattern of testicular cells in the seminiferous tubules in the control and COVID-19 groups (Fig. 5A and B).

Fig. 5
figure 5

The effect of SARS-CoV-2 on the testicular cells morphology in the human post-mortem. A Photomicrograph of the testis stained with H&E, ×40. SG (spermatogonia), PS (primary spermatocyte), ST (round spermatid), SC (Sertoli cell). A and B are cartoons that illustrate the differing morphologies of testicular tubules in control (A) and COVID-19 (B) groups

The spatial distribution of the testicular cells

Voronoi tessellation of the testicular cells in the control and COVID-19 groups is shown in Fig. 6A–E. The data showed that 8% of the areas of the polygons of the testicular cells in the testis in the COVID-19 samples were located in the range of ˃ 120 μm2, while in the control samples, just 17% of the areas of the polygons were located in this range. The results also showed that 50% of the areas of the polygons of the testicular cells in the testis in the COVID-19 samples were located in the range of 120–130 μm2, while in the control samples, just 60% of the areas of the polygons were located in this range.

Fig. 6
figure 6

The effect of SARS-CoV-2 on the spatial pattern of testicular cells in the human post-mortem using Voronoi tesselation. AC Mean ± SD of Voronoi polygon area (%) and Coefficient of Variation (CV) within the testis in study groups. D and E Representative photograph, polygon area, and schematic of Voronoi tessellation of the testicular cells in the control (D) and COVID-19 (E) groups

Additionally, our results also showed that 40% of the areas of the polygons of the testicular cells in the testis in the COVID-19 samples were located in the range of ˂ 130 μm2, while in the control samples, just 37% of the areas of the polygons were located in this range. Furthermore, the COVID-19 samples showed fewer testicular cells in the range of 120–130 μm2 than control samples (Fig. 6A–E). Based on the coefficient of variation (CV) classification, the mean CV of polygon areas in both groups is located in a random range (33–64%) (Fig. 6A–E).

Expression of CD68 in testis based on immunocytochemistry

Immunostaining showed that CD68 was and strongly expressed in the interstitial tissue of testis in the COVID-19 group when compared with the control group (p < 0.001) (Fig. 7C). In the interstitium of testicular tissue, there was mild inflammatory infiltrate composed predominantly of macrophage (CD68-positive cells) (Fig. 7A and B) as confirmed immunohistochemistry.

Fig. 7
figure 7

The effect of SARS-CoV-2 on the macrophage biomarker in the human post-mortem testcular tissue. Photomicrograph of the testis immunohistochemistry staining for CD68, ×40. A Control group and B COVID-19 group. C Total number of CD68-positive cells in study groups. (***p < 0.001). Arrows are CD68-positive cells

Immunohistochemistry of occluding, claudin-11 and connexin-43

Figure 8A shows the expression of claudin-11, occludin and connexin-43 in testicular tissues in COVID-19 groups. A significant difference was observed in the expression of target proteins that are involved in the establishment of BTB. The expression of claudin-11, occludin and connexin-43 in testis showed a significant decrease in the COVID-19 group when compared with the normal control group (p < 0.001, p < 0.001, p < 0.001, respectively) (Fig. 8B). Immunohistochemical staining showed that SARS-CoV-2 induced disrupts BTB by decreasing BTB junction proteins, such as claudin, occludin and connexin-43. The brownish color was considered to be evidence of a positive expression of claudin-11, occludin and connexin-43 in the testis. To convert the intensity and quantification of immunoreactivity of protein expression, we used the optical density (OD), following this formula: OD = log (max intensity/Mean intensity), where max intensity = 255 for 8-bit images.

Fig. 8
figure 8

The effect of SARS-CoV-2 on the junctional proteins of BTB in the human post-mortem testicular tissue. A The patterns of immunostaining of claudin-11, occluding, and connexin-43 in the control and COVID-19 groups. Immunoreactivity pattern (OD = log(max intensity/Mean intensity) of claudin-11, occluding, and connexin-43 in testis in the study groups. B The significant difference between COVID-19 with the control groups is indicated. (***p < 0.001). Arrows indicate junctional proteins

Discussion

Some researchers have confirmed that COVID-19 could induce male reproductive dysfunction, but the exact mechanisms of COVID-19 on the male reproductive system are still not clear. The analysis within this report provides evidence that the COVID-19 could induce structural changes of the seminiferous tubules and alterations of the blood–testis barrier which gives insight into the damage to the seminiferous tubules associated with COVID-19.

Then, we first proved that COVID-19 infection would disrupt BTB integrity through producing higher levels of pro-inflammatory cytokines, including TNF-α, IL-1b and IL-6 as well as decreasing expressions of genes and proteins related to BTB, including occludin, claudin and connexin-43.

It is right to point out that based on previous studies, the spikes of the SARS-CoV-2 virus interact with angiotensin-converting enzyme 2 (ACE2) to infect the cell. Except for the lung and some other organs, ACE2 expresses extensively on the surface of most testicular cells, interstitial cells, Sertoli cells, Leydig cells and the early stage of germ cells [15,16,17]. As reported in other studies, angiotensinogen is converted via renin to angiotensin I (Ang I). ACE2 regulates negatively the function of ACE by converting Ang I to Ang 1–9 and Ang II to Ang 1–7 [5, 18, 19]. In general, the activated ACE-Ang-II-AT1R axis leads to pro-inflammatory cytokines. On the other hand, the ACE-2-Ang- (1–7)-Mas axis has an anti-inflammatory effect [20,21,22]. It is also suggested that the attachment of SARS-CoV-2 to the ACE2 receptor may trigger an inflammatory response that could lead to Sertoli and Leydig cells dysfunction [23]. Our previous results unpublished showed that the expression of the ACE2 gene and protein was a reduction in testicular cells in the COVID-19 samples compared with the healthy control samples. Therefore, the possibilities of SARS-CoV-2 virus entry to the testicular cells are unlikely. Therefore, it can be said, following the down-regulation of ACE2 and consequently reduction of Ang 1–7, SARS-CoV-2 has been able to increase the level of TNF-α, IL-1b and IL-6 gene expressions in the COVID-19 group in comparison with the control group [24]. The TNF-α, mediated by activation of the ERK1/2 signaling pathway, leads to increase tight junction permeability [25]. Our present study showed COVID-19 increased the mRNA levels of TNF-α, IL-1β and IL-6, suggesting that COVID-19 induced inflammation in testis. And we speculated that COVID-19 reduced protein levels of occludin, claudin-11 and connexin-43 associated with inflammation in testis. Several inflammatory cytokines, such as TGF-b3, TNF-α, and IL-1, play an effective role in the regulation of blood–testis barrier (BTB) dynamics [26]. Furthermore, IL6 might play a role in the downregulation of occludin expression and the modulation of BTB permeability [27]. Additionally, TNF-α temporarily facilitates permeability of the BTB during spermatogenesis to germ cell migration [28] and it has been reported that TNF-α inhibits the expression of claudin-11 in Sertoli cells [29]. Therefore, based on these studies, it can be explained, how an increase in expression of inflammatory cytokines, including TNF-α, IL-1β and IL-6, has caused a decrease in BTB proteins expression including occludin, claudin-11 and connexin-43 in the COVID-19 group. Another study has reported that the levels of pro-inflammatory chemokines and cytokines, including IL-6, TNF-α, and MCP-1, play a noteworthy role in the immunopathology of SARS-CoV-2, in semen samples of COVID-19 patients increased [30]. However, one of the mechanisms underlying COVID-19-associated barrier disruption is increasing the susceptibility of testicular cells to interactions with other immune cells, leading to increased infiltration across the barrier 31). The rapid proliferation of the virus upon arrival may lead to the secretion of immune cells and apoptosis of endothelial cells in the testes with subsequent release of inflammatory biomarkers including cytokines, chemokines, and adhesion molecules [32]. Also, inflammatory cytokines including TNF-α, IL-1β may cause oxidative stress in Sertoli cells and alter BTB integrity [33].

On the other hand, based on previous studies, connexin-43, the predominant gap junction protein in the testis, is necessary for tight junction reassembly at the BTB [34] and is essential for the initiation and maintenance of spermatogenesis [35]. Furthermore, reduced expression of the gap junction protein connexin-43 can lead to damage to BTB. Tight junction proteins, such as JAM-A, claudin-11, occludin, ZO-1, and N-cadherin, play a key role in keeping the integrity of the structure and the stability of the cell's internal environment. The present study showed that COVID-19 infection disrupts BTB permeability and spatial arrangement by reducing the expression of junctional molecules, which is in agreement with the fact that tight junctions play a critical role in the junctions between Sertoli cells. Additionally, the current study is the first to do an investigation to evaluate the spatial arrangement of testicular cells using second-order stereology and Voronoi tessellation test by diminished populations of testicular cells in the testes after infection. [31,32,33], suggests that there is a correlation between the disruption of the spatial arrangement of the testicular cells and SARS-CoV-2. The results of this study indicate that Sertoli cells are damaged, in patients infected by SARS-CoV-2. The plots of the pair-correlation function of Sertoli cells indicated a wide range of gaps in the spatial distribution of Sertoli cells in the COVID-19 group. Our Voronoi tessellation findings also showed that changes in the spatial arrangement of testicular cells as well as we observed testicular cells were randomly distributed in the control and COVID-19 groups. The spatial arrangement of the cell is highly associated with cell-to-cell interactions. Thus, these events lead to the loss of a spermatogenic cell in the testes and will undoubtedly change the spatial pattern of the Sertoli and spermatogenic cells. The number of Sertoli cells determines efficient spermatogenesis and directs correlation of the number of germ cells with the number of functional Sertoli cells. Decreasing the number of Sertoli cells and subsequent changes in the spatial arrangement of the cells disrupt the structure and function of BTB and lead to impaired spermatogenesis. In this study, we observed that the coronavirus can disrupt the BTB. Given that coronavirus is an emerging virus that has recently negatively affected the lives of many people around the world, we still have to wait for more extensive research to better identify the chronic complications of this disease. In conclusion, the present study demonstrated that SARS-CoV-2 could induce the up-regulation of the pro-inflammatory cytokine and down-regulation of junctional proteins that cause a disruption of BTB and impairment of spermatogenesis. This implies the importance of applying further care in reproductive health in men infected by SARS-CoV-2. Screening of these patients is necessary to look over the reversibility of these impacts and improve standard therapeutic strategies to aid these patients.