# Multiorgan and Vascular Tropism of SARS-CoV-2

Cédric Hartard, Ahlam Chaqroun, Nicla Settembre, Guillaume Gauchotte, Benjamin Lefevre, Elodie Marchand, Charles Mazeaud, Duc Trung Nguyen, Laurent Martrille, Isabelle Koscinski, et al.

► To cite this version:

Cédric Hartard, Ahlam Chaqroun, Nicla Settembre, Guillaume Gauchotte, Benjamin Lefevre, et al. Multiorgan and Vascular Tropism of SARS-CoV-2. Viruses, 2022, 14 (3), pp.515. 10.3390/v14030515. hal-03765141

HAL Id: hal-03765141

https://hal.univ-lorraine.fr/hal-03765141

Submitted on 4 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Case Report

# Multiorgan and Vascular Tropism of SARS-CoV-2

Cédric Hartard<sup>1,2,●</sup>, Ahlam Chaqroun<sup>2</sup>, Nicla Settembre<sup>3,●</sup>, Guillaume Gauchotte<sup>4,●</sup>, Benjamin Lefevre<sup>5,6</sup>, Elodie Marchand<sup>7</sup>, Charles Mazeaud<sup>8,●</sup>, Duc Trung Nguyen<sup>9</sup>, Laurent Martrille<sup>7,●</sup>, Isabelle Koscinski<sup>10,●</sup>, Sergueï Malikov<sup>3</sup> and Evelyne Schvoerer<sup>1,2,*</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hartard, C.; Chaqroun, A.; Settembre, N.; Gauchotte, G.; Lefevre, B.; Marchand, E.; Mazeaud, C.; Nguyen, D.T.; Martrille, L.; Koscinski, I.; et al. Multiorgan and Vascular Tropism of SARS-CoV-2. *Viruses* **2022**, *14*, 515. https://doi.org/10.3390/v14030515

Academic Editors: Luis Martinez-Sobrido and Fernando Almazan Toral

Received: 21 January 2022

Accepted: 1 March 2022

Published: 3 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Although the respiratory tract is the main target of SARS-CoV-2, other tissues and organs are permissive to the infection. In this report, we investigated this wide-spectrum tropism by studying the SARS-CoV-2 genetic intra-host variability in multiple tissues. The virological and histological investigation of multiple specimens from a post-mortem COVID-19 patient was performed. SARS-CoV-2 genome was detected in several tissues, including the lower respiratory system, cardio-vascular biopsies, stomach, pancreas, adrenal gland, mediastinal ganglion and testicles. Subgenomic RNA transcripts were also detected, in favor of an active viral replication, especially in testicles. Ultra-deep sequencing allowed us to highlight several SARS-CoV-2 mutations according to tissue distribution. More specifically, mutations of the spike protein, i.e., V341A (18.3%), E654 (44%) and H655R (30.8%), were detected in the inferior vena cava. SARS-CoV-2 variability can contribute to heterogeneous distributions of viral quasispecies, which may affect the COVID-19 pathogeny.

**Keywords:** SARS-CoV-2; viral quasispecies; ultra-deep sequencing

## 1. Introduction

SARS-CoV-2, causing coronavirus disease 2019 (COVID-19), was detected for the first time in December 2019. Since then, this new pathogen has rapidly spread, with more than 300 million of infected people worldwide and more than 5 million of deaths. The clinical signs involve mainly the respiratory tract and the subsequent complication results in an acute respiratory distress syndrome (ARDS). As the COVID-19 pandemic continues to progress, knowledge about its wide-spectrum clinical manifestations and associated complications is advancing. Other organs failures are thus reported, such as cardiovascular involvement, with the description of myocarditis, vascular endothelial injury or thromboembolic events [1-3].

The ubiquitous distribution of the ACE2 receptors in several tissues may explain the multi-organ dysfunction linked to COVID-19, as well as the inflammatory and immune

<footnote n="1">1 Laboratory of Virology, Department of Microbiology, CHRU Nancy, Université de Lorraine, F-54000 Nancy, France; c.hartard@chru-nancy.fr</footnote>

<footnote n="2">2 CNRS, LCPME, Université de Lorraine, F-54000 Nancy, France; ahlam.chaqroun@univ-lorraine.fr</footnote>

<footnote n="3">3 Department of Vascular & Endovascular Surgery, CHRU Nancy, Université de Lorraine, INSERM 1116, F-54000 Nancy, France; n.settembre@chru-nancy.fr (N.S.); s.malikov@chru-nancy.fr (S.M.)</footnote>

<footnote n="4">4 Department of Biopathology, Centre de Ressources Biologiques, Université de Lorraine, F-54000 Nancy, France; g.gauchotte@chru-nancy.fr</footnote>

<footnote n="5">5 Department of Infectious and Tropical Diseases, CHRU Nancy, Université de Lorraine, F-54000 Nancy, France; b.lefevre@chru-nancy.fr</footnote>

<footnote n="6">6 APEMAC, Université de Lorraine, F-54000 Nancy, France</footnote>

<footnote n="7">7 Department of Legal Medicine, CHU Nancy, Université de Lorraine, F-54000 Nancy, France; e.marchand2@chru-nancy.fr (E.M.); laurent.martrille@chu-montpellier.fr (L.M.)</footnote>

<footnote n="8">8 Department of Urology, CHRU Nancy, Université de Lorraine, F-54000 Nancy, France; c.mazeaud@chru-nancy.fr</footnote>

<footnote n="9">9 Service ORL, CHRU Nancy, Université de Lorraine, F-54000 Nancy, France; dt.nguyen@chru-nancy.fr</footnote>

<footnote n="10">10 Laboratoire de Biologie de la Reproduction/CECOS Lorraine, CHRU Nancy, INSERM U1256 NGERE, F-54000 Nancy, France; koscinski.isa@gmail.com</footnote>

<footnote>* Correspondence: e.schvoerer@chru-nancy.fr</footnote>
components [4]. Hence, the presence of the SARS-CoV-2 genome has already been reported outside the respiratory tract, most notably in the kidney, gastrointestinal system, nervous system and blood vessels, with some arguments in favor of a local replication [5–11].

To better understand the reasons for the wide-spectrum clinical features, it is essential to investigate the viral extrapulmonary dissemination of SARS-CoV-2. Thus, autopsy studies can be particularly relevant in order to characterize possible viral- and/or host-related characteristics linked to differential tissue affinity. In the present study, we explore this issue via multiple post-mortem histological and viral analyses, including ultra-deep sequencing, from specimens collected from a COVID-19 patient during the second wave of the pandemic (October 2020, France) in order to study intra-host variability.

## 2. Materials and Methods

### 2.1. *Clinical Data and Autopsy*

A 70-year-old man was admitted to Nancy University Hospital in October 2021. He presented a medical history of chronic heart failure due to atrial fibrillation and coronary disease, peripheral vascular disease with stroke and cognitive sequelae, arterial hypertension and active smoking. SARS-CoV-2 infection, confirmed by nasopharyngeal RT-PCR, began on 12 October 2021 with confusion, cough and pharynx erythema. After 6 days, he was admitted for confusion, dehydration and ARDS, associated with a typical aspect of lung infection on a chest CT scan. Biological analysis showed thrombopenia (80 G/L), lymphopenia (0.2 G/L), hypernatremia (155 mmol/L), acute kidney injury KDIGO 3 (blood creatinine: 218 µmol/L), a high level of C-reactive protein (146 mg/L) and hyperlactatemia without acidosis (arterial lactate: 2.8 mmol/L). Despite adapted treatment in the infectious disease unit with high-flow oxygen therapy, progressive hypernatremia correction and dexamethasone, the lung function worsened. Exclusive palliative care was started on 23 October, and the patient died on 29 October of a severe ARDS.

A medical autopsy was performed two days after his death. The family consented to medical autopsy and the analysis of specimens. Lungs were congestive and globally densified. The abdominal aortic lumen was occluded by a voluminous thrombus, developed on an ulcerated atheromatous plaque. Coronary and carotid arteries showed atherosclerosis, without complication.

Multiple post-mortem specimens were collected, including olfactory clefts, trachea, lungs, pleura, mediastinal lymph nodes, stomach, liver, pancreas, spleen, duodenum, transverse colon, thyroid, brain, kidneys, adrenal glands, testicles, myocardium, inferior vena cava, mesenteric artery, carotid artery, basilar artery, aorta, blood, abdominal wall muscle, healthy and damaged skin. For each anatomical localization, tissue samples were frozen for viral analyses, and other fragments were immersed in 4% formalin solution for histopathology.

### 2.2. *Histopathology*

After the sampling of the different organs, 4 µm-thick tissue sections were obtained and stained with hematoxylin, eosin and saffron. An immunohistochemical assay was performed with an anti-SARS-CoV-2 nucleo-capsid protein antibody (SARS-CoV Nucleocapsid Protein Antibody, rabbit polyclonal, NB100-56576 Novus Biologicals, 1/50; Flex+ Envision revelation system, Dako Omnis) on lung, myocardium, kidney, pancreas, olfactory cleft, olfactory bulb and testicle tissue samples. Positive and negative controls were used for immunohistochemistry assay validation (infected and non-infected cell lines). A positive immunostaining was defined by a granular staining of the cell cytoplasm.

### 2.3. *Viral Analyses*

After viral RNA extraction using Nuclisens® easyMAG® technology, each sample was analyzed for SARS-CoV-2 detection using RT-qPCR targeting IP2 and IP4 of the RdRp gene, according to the protocol developed by the Pasteur Institut [12]. In the case of a positive result, ultra-deep sequencing (DeepCheck® Whole Genome SARS-CoV-2 Genotyping
V1, ABL SA) was performed in accordance with the manufacturer recommendation on an Illumina iSeq 100 system to study the SARS-CoV-2 genetic variability according to tissue distribution. Briefly, the RT-PCR multiplex test, which consisted of 98 amplicons of approximatively 400 bp in length, spanned the entire genome of SARS-CoV-2. Subsequently, the amplicons were used for next-generation sequencing, and the downstream analysis ABL DeepCheck® software was used to list SARS-CoV-2 mutations.

Moreover, to explore the evidence of an active viral replication, an additional subgenomic RNA RT-qPCR was performed, as described previously [13].

### 3. Results

Histopathological analyses showed the presence of severe lung diffuse alveolar damages, characterized by the presence of dystrophic pneumocytes and numerous endoalveolar hyaline membranes, associated with a moderate inflammatory infiltrate (Figure 1A). The latter was mainly composed of macrophages and lymphocytes. A mild perivascular lymphocytic infiltrate was observed in the myocardium, without cardiomyocytes alteration. A recent mixt thrombus was seen in the abdominal aorta. The examination of the central nervous system revealed the presence of ancient frontal ischemic sequelae, without signs of recent necrosis or inflammation.

A strong staining of endoluminal macrophages associated with a mild straining of epithelial cells was found in the olfactory cleft (Figure 1B). In the lung, only endo-alveolar macrophages showed a significant positivity. In the testicles, a moderate cytoplasmic staining of seminiferous tubules cells was observed (Figure 1C). No significant staining was observed in the other organs (Figure 1D), but the analyses were limited by post-mortem autolysis (Figure 1E).

The SARS-CoV-2 genome was detected in several tissues, including the lower respiratory system, various cardio-vascular biopsies, the stomach, the pancreas, the adrenal gland mediastinal ganglion and the testicles (Figure 2, Table 1). Subgenomic RNA transcripts, in favor of an active SARS-CoV-2 replication site, were detected in all tissues positive for genome detection, except in the pancreas, aorta and vena cava. Concentrations were overall 1000-fold lower compared to genomic RNA, except in the nasopharynx (collected 9 days before death) and testicles, where they were only 10-fold lower compared to genomic RNA (Table 1), in favor of a strong viral replication.

All SARS-CoV-2 genomes detected showed that the virus belonged to the clade 20A.EU2, with the presence of the characteristic S477N and the D614G mutations on the spike protein. However, whole genome sequencing allowed the identification of numerous mutations, depending on the type of biopsy.

Regarding the spike protein, specific mutations were observed in vascular tissues, the H655R and the V341A in vena cava (30.8% and 18.3%, respectively), as well as the E654A (44.0%), which was also poorly detected in other tissues in which the E654G was highly represented (trachea, myocardium, stomach lining, carotid and mesenteric artery) (Table 1).

For the spike protein, the Q526* mutation (ORF1a) was only detected in the mesenteric artery, and the F3271L, absent in the initial nasopharyngeal swab, was detected in pleura (14.4%), trachea (12.5%), myocardium (13.3%), stomach lining (10.5%) and testicle (9.6%).

Furthermore, the E787V mutation (ORF1b) was detected in the totality of the analyzed tissues (8.0 to 25.6% of the viral quasispecies), but was undetected in the initial nasopharyngeal swab. Various mutations were present in the stomach lining biopsy, especially the Q57H (ORF3a), reaching 98.6%, but were not detected in other tissues.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (A) Lung tissue, standard histology: diffuse alveolar damages with hyaline membranes (arrows) and moderate inflammatory infiltrate (hematoxylin, eosin and saffron, ×200). Immunohistochemistry, SARS-CoV-2 nucleo-capsid protein: (B) olfactory cleft: strong cytoplasmic staining of endoluminal macrophages (arrows) and mild staining of the epithelium (arrowheads); (C) testicular parenchyma: moderate cytoplasmic staining of germ cells in the seminiferous tubules (arrows); (D) negative staining in renal parenchyma; (E) non-significant weak and diffuse staining in autolytic pancreatic parenchyma (×200).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Anatomical diagram of variant distribution according to autopsied patient tissues. The amino acid colored in blue matches SARS-CoV-2 Wuhan-Hu-1 sequence (NC 045512.2: 21563-25384); amino acids with other colors match the mutated sequence (designed by biorender).</figcaption></figure>
<table><caption>**Table 1.** Ultra-deep sequencing of SARS-CoV-2 genome isolated from post-mortem biopsies. The gradient color represents the notable proportion of mutations present inside the quasispecies. ND: not determined.</caption><thead><tr><th>Specimen</th><th>Nasopharynx<br>(9 Days before<br>Death)</th><th>Trachea</th><th>Lung</th><th>Pleura</th><th>Mediastenal<br>Lymph<br>Nodes</th><th>Stomach</th><th>Pancreas</th><th>Adrenal<br>Glands</th><th>Testicle</th><th>Myocardium</th><th>Inferior<br>Vena Cava</th><th>Mesenteric<br>Artery</th><th>Carotid<br>Artery</th><th>Aorta</th></tr></thead><tbody><tr><td>SARS-Co-2 genome<br>detection (IP4-Cq)</td><td>27.7</td><td>24.6</td><td>23.5</td><td>24.2</td><td>32.4</td><td>23.6</td><td>32</td><td>27.3</td><td>30</td><td>25.5</td><td>28.8</td><td>22.8</td><td>20.4</td><td>31.9</td></tr><tr><td>Subgneomic RNA<br>detection (Cq)</td><td>31.3</td><td>33.9</td><td>32.1</td><td>33.6</td><td>37.1</td><td>31.7</td><td>ND</td><td>34.7</td><td>32.1</td><td>34.7</td><td>ND</td><td>32.3</td><td>29.1</td><td>ND</td></tr><tr><td colspan="15">Mutations</td></tr><tr><td>ORF1a</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>R43H</td><td>0.00%</td><td>0.00%</td><td>ND</td><td>0.00%</td><td></td><td>14.89%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>ND</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>Q526*</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>34.00%</td><td>0.00%</td><td></td></tr><tr><td>A537V</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>18.75%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>9.15%</td><td>3.03%</td><td>0.00%</td><td>22.82%</td><td></td></tr><tr><td>V1968M</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>ND</td><td>0.00%</td><td>21.66%</td><td></td></tr><tr><td>M3087I</td><td>72.73%</td><td>96.77%</td><td>100.00%</td><td>97.12%</td><td></td><td>98.39%</td><td></td><td></td><td>98.58%</td><td>97.81%</td><td>ND</td><td>ND</td><td>ND</td><td></td></tr><tr><td>F3271L</td><td>0.00%</td><td>12.50%</td><td>0.00%</td><td>14.41%</td><td></td><td>10.48%</td><td></td><td></td><td>9.55%</td><td>13.35%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>Y3300H</td><td>97.81%</td><td>96.56%</td><td>98.89%</td><td>97.58%</td><td></td><td>98.65%</td><td></td><td></td><td>98.27%</td><td>97.68%</td><td>97.54%</td><td>99.48%</td><td>99.24%</td><td></td></tr><tr><td>D3703H</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>ND</td><td>ND</td><td>24.44%</td><td></td></tr><tr><td>G4374D</td><td>26.62%</td><td>8.41%</td><td>0.00%</td><td>0.00%</td><td></td><td>27.49%</td><td></td><td></td><td>10.25%</td><td>0.00%</td><td>ND</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>ORF1b</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>A176S</td><td>100.00%</td><td>0.00%</td><td>ND</td><td>100.00%</td><td></td><td>0.00%</td><td></td><td></td><td>100.00%</td><td>92.86%</td><td>ND</td><td>ND</td><td>ND</td><td></td></tr><tr><td>P314L</td><td>100.00%</td><td>92.86%</td><td>ND</td><td>100.00%</td><td></td><td>100.00%</td><td></td><td></td><td>76.47%</td><td>100.00%</td><td>ND</td><td>ND</td><td>ND</td><td></td></tr><tr><td>V767L</td><td>98.25%</td><td>97.01%</td><td>98.40%</td><td>98.68%</td><td></td><td>98.20%</td><td></td><td></td><td>98.93%</td><td>98.71%</td><td>99.07%</td><td>99.30%</td><td>99.22%</td><td></td></tr><tr><td>E787V</td><td>0.00%</td><td>15.20%</td><td>15.13%</td><td>12.32%</td><td></td><td>14.08%</td><td></td><td></td><td>8.31%</td><td>10.50%</td><td>8.04%</td><td>15.87%</td><td>25.62%</td><td></td></tr><tr><td>K1141R</td><td>100.00%</td><td>92.06%</td><td>100.00%</td><td>94.44%</td><td></td><td>100.00%</td><td></td><td></td><td>89.47%</td><td>96.00%</td><td>ND</td><td>100.00%</td><td>98.80%</td><td></td></tr><tr><td>E1184D</td><td>90.91%</td><td>100.00%</td><td>100.00%</td><td>84.00%</td><td>failed</td><td>97.01%</td><td>failed</td><td>failed</td><td>93.75%</td><td>92.93%</td><td>ND</td><td>100.00%</td><td>100.00%</td><td>failed</td></tr><tr><td>Spike</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>V341A</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>18.28%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>S477N</td><td>100.00%</td><td>98.76%</td><td>100.00%</td><td>97.98%</td><td></td><td>98.94%</td><td></td><td></td><td>100.00%</td><td>98.98%</td><td>100.00%</td><td>97.30%</td><td>100.00%</td><td></td></tr><tr><td>D614G</td><td>98.82%</td><td>98.29%</td><td>98.77%</td><td>98.98%</td><td></td><td>99.29%</td><td></td><td></td><td>98.50%</td><td>99.02%</td><td>98.98%</td><td>99.60%</td><td>99.51%</td><td></td></tr><tr><td>E654A</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>5.08%</td><td></td><td>2.20%</td><td></td><td></td><td>0.00%</td><td>2.79%</td><td>44.00%</td><td>0.00%</td><td>2.68%</td><td></td></tr><tr><td>E654G</td><td>9.14%</td><td>16.30%</td><td>0.00%</td><td>0.00%</td><td></td><td>18.18%</td><td></td><td></td><td>7.05%</td><td>15.08%</td><td>0.00%</td><td>14.35%</td><td>18.21%</td><td></td></tr><tr><td>H655R</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>0.00%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>30.77%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>I1172T</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>10.91%</td><td></td><td>10.82%</td><td></td><td></td><td>10.46%</td><td>8.55%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>ORF3a</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Q57H</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td><td>98.61%</td><td></td><td></td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>P207S</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td></td><td>0.00%</td><td></td><td></td><td>ND</td><td>ND</td><td>ND</td><td>25.00%</td><td>ND</td><td></td></tr><tr><td>N</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>K169N</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td>13.10%</td><td></td><td>3.32%</td><td></td><td></td><td>5.27%</td><td>2.80%</td><td>0.00%</td><td>0.00%</td><td>0.00%</td><td></td></tr><tr><td>M234I</td><td>94.12%</td><td>97.33%</td><td>81.25%</td><td>96.64%</td><td></td><td>99.17%</td><td></td><td></td><td>98.71%</td><td>98.95%</td><td>100.00%</td><td>98.61%</td><td>99.61%</td><td></td></tr><tr><td>K256*</td><td>ND</td><td>0.00%</td><td>10.27%</td><td>10.53%</td><td></td><td>10.18%</td><td></td><td></td><td>11.81%</td><td>7.76%</td><td>13.51%</td><td>0.00%</td><td>8.76%</td><td></td></tr><tr><td>A376T</td><td>100.00%</td><td>100.00%</td><td>100.00%</td><td>98.00%</td><td></td><td>98.96%</td><td></td><td></td><td>100.00%</td><td>98.98%</td><td>100.00%</td><td>98.18%</td><td>99.25%</td><td></td></tr></tbody></table>
## 4. Discussion

Mutations that appear in the SARS-CoV-2 genome are responsible for the variability of viral strains circulating in the population. This variability must also be considered within viruses multiplying in the same individual as a viral quasispecies, which could determine, at least partially, the tissue tropism. The concept of genetic variability and its involvement in tissue tropism, as well as the severity of the disease, is not new, and has been previously observed for a large number of viruses including Human Immunodeficiency Virus (HIV) and Hepatitis E Virus (HEV) [14,15]. This tropism is at least partially determined by the permissiveness of the host cells viral receptor to the virus, meaning through the spike protein in the case of coronaviruses.

The ACE2 receptors are expressed in the nasopharynx but also in several other organs, making the infection possible in multiple cells, including endothelial cells, cardiomyocytes, enterocytes, parietal cells, Leydig cells, spermatogonia or Sertoli cells [16,17]. In this report, the detection of subgenomic RNA in various tissues previously positive for SARS-CoV-2 genome is an argument in favor of a local viral replication. Even though subgenomic RNA cannot be used as a current viability marker considering the persistence of RNA [18,19], the high genome concentration in tissue associated with the presence of this replicative intermediate at a high comparative level is in favor of an active SARS-CoV-2 replication at a given time. This supports the indication of exploring SARS-CoV-2 in other tissues than the respiratory tract, as previously proposed [20]. In addition to its detection outside the respiratory tract, a deep viral characterization may be relevant to better understand COVID-19 pathogeny. At the variant scale, Wagner et al. demonstrated a higher in vitro susceptibility of coronary artery endothelial cells in the case of infection with the B.1.1.7 lineage compared to wild type [21]. At the scale of viral quasispecies, Rueca et al. described the intra-host variability of SARS-CoV-2 in the upper and lower respiratory tract, with some mutations reaching 30% of the quasispecies in favor of variant compartmentalization [22]. Again, in the respiratory tract, intra-host variability was previously reported during a longitudinal follow-up of a patient infected with SARS-CoV-2 [23]. Despite these reports, the genetic variability of SARS-CoV-2 remains poorly explored outside the epidemiology context.

In this study, ultra-deep sequencing allowed the observation of several SARS-CoV-2 mutations in various proportions according to the tissue distribution. Ultra-deep sequencing has to be applied outside the respiratory tract to obtain data concerning intra-host SARS-CoV-2 quasispecies variability, differential tropism or pathogeny. When focusing on the spike protein, signature mutations associated with vascular tropism may be suggested in vena cava (i.e., V341A, E654A and H655R) (Figure 1, Table 1). These mutations have not yet been characterized in the literature; however, position 655 is also mutated in the 20J/501Y.V variant (H655Y) and has been correlated to antibody neutralization escape and to the modulation of the interaction between the spike glycoprotein and ACE2 receptors [24]. Moreover, the H655 residue is closed to the E654A mutation, and both are near the furin cleavage site of the spike protein (S1/S2 site) involved in SARS-CoV-2 cell entry. The E654A mutation suggests a significant change in the physicochemical properties of the amino acid sequence, since a glutamic acid (E) is polar and negatively charged with a higher molecular weight (MW = 128 g/mol) compared to an alanine (A), which is nonpolar and uncharged (71 g/mol). The major impact of such changes has already been described for the D614G mutation in which a polar, negatively charged aspartic acid (D) with a higher MW (114 g/mol) is substituted with a nonpolar uncharged glycine (G) with a lower MW (57 g/mol). The difference in D614G physicochemical properties has a significant impact on the three-dimensional structure of the spike protein, giving the RBD (Receptor Binding Domain) a more open conformation than the wild type strain (Wuhan-Hu-1) and, therefore, an improvement in virus interaction with the host's cell receptors [25]. Henceforth, a potential impact on the protein functions could be suggested concerning the implication of the E654A mutation regarding its interaction with endothelial cell receptors.

Variability was not limited to the vascular system since the significant compartmentalization of viral quasi-species was also observed in other organs (i.e., the stomach lining and
pleura), in association with spike mutations but also regarding other regions of the genome (ORF1a, ORF1b and ORF3a).

Additionally, viral analyses allowed the detection of SARS-CoV-2 in the testicular parenchyma. Early in the pandemic, the expression of the ACE2 receptors by germ cells, Sertoli cells and Leydig cells suggested a potential testicular reservoir for SARS-CoV-2 with possible sexual transmission [26], as well as a potential impact on male fertility [27] or even on offspring [28]. Consequently, the presence of the virus has been studied in semen from patients with infections of varying severity and expression, with conflicting results [29,30]. Further investigation of the testicular consequences of COVID-19 has highlighted that SARS-CoV-2 can directly impair male fertility through binding to the ACE2 receptor and the invasion of all testicular cell types [11]. Despite the documented infections, little is known about mutations occurring in the testes, their selective advantage and their consequences. This case report demonstrated a high level of SARS-CoV-2 replication in the testes, but the viral genome pattern was not testis specific, and to our knowledge, there is no evidence in the literature to support a link between the clinic and a testis-specific variant.

Finally, from a technical aspect, RT-qPCR appeared more sensitive than immunohistochemical assay to highlight the presence of SARS-CoV-2 in tissues. Ultra-deep sequencing, essential to study quasispecies distribution, is however, more uncertain since a low coverage may be observed for some regions of the genome (i.e., vascular samples in this study), or may even fail, despite a significant viral load in tissue (e.g., adrenal glands). The degradation of post-mortem samples may be the cause of the low sensitivity. In addition, the technique has to be improved to explore the genetic variability of SARS-CoV-2 in samples that are more complex than the classical nasopharyngeal swabs, especially in the case of low viral load, as observed here in the aorta, pancreas or mediastinal lymph nodes (i.e., Cq values ≥ 32).

**5. Conclusions**

To summarize, SARS-CoV-2 was detected in multiple tissues collected from a death related to COVID-19. A local replication was observed in various tissues, including vessels and testicular parenchyma, leading to the relatively heterogeneous distribution of viral quasispecies associated with compartmentalization. Larger studies are now needed to explore the SARS-CoV-2 variability during the course of the infection, and in different clinical conditions, to better understand the dynamics and impact of the minority SARS-CoV-2 variants on pathogeny.

**Author Contributions:** Conceptualization, E.S. and S.M.; methodology, C.H., N.S., G.G., L.M., S.M. and E.S.; formal analysis, C.H. and A.C.; data curation, C.H., A.C., N.S., G.G., B.L., E.M., C.M., D.T.N., L.M., I.K. and S.M.; writing—original draft preparation, C.H., A.C., G.G., B.L. and I.K.; writing—review and editing, C.H., N.S., G.G., B.L., E.M., C.M., D.T.N., I.K. and E.S.; supervision, C.H., E.S. and S.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Informed consent was obtained from the family of the subject involved in the study.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Ackermann, M.; Verleden, S.E.; Kuehnel, M.; Haverich, A.; Welte, T.; Laenger, F.; Vanstapel, A.; Werlein, C.; Stark, H.; Tzankov, A.; et al. Pulmonary Vascular Endothelialitis, Thrombosis, and Angiogenesis in Covid-19. *N. Engl. J. Med.* **2020**, *383*, 120–128. [CrossRef] [PubMed]

2. Gauchotte, G.; Venard, V.; Segondy, M.; Cadoz, C.; Esposito-Fava, A.; Barraud, D.; Louis, G. SARS-CoV-2 fulminant myocarditis: An autopsy and histopathological case study. *Int. J. Legal Med.* **2021**, *135*, 577–581. [CrossRef] [PubMed]
3. Basso, C.; Leone, O.; Rizzo, S.; De Gaspari, M.; Van Der Wal, A.C.; Aubry, M.C.; Bois, M.C.; Lin, P.T.; Maleszewski, J.J.; Stone, J.R. Pathological features of COVID-19-associated myocardial injury: A multicentre cardiovascular pathology study. *Eur. Heart J.* **2020**, *41*, 3827–3835. [CrossRef] [PubMed]

4. Guéant, J.; Guéant-Rodriguez, R.; Fromonot, J.; Oussalah, A.; Louis, H.; Chery, C.; Gette, M.; Gleye, S.; Callet, J.; Raso, J.; et al. Elastase and exacerbation of neutrophil innate immunity are involved in multi-visceral manifestations of COVID-19. *Allergy* **2021**, *76*, 1846–1858. [CrossRef]

5. Delgado-Gonzalez, P.; Gonzalez-Villarreal, C.A.; Roacho-Perez, J.A.; Quiroz-Reyes, A.G.; Islas, J.F.; Delgado-Gallegos, J.L.; Arellanos-Soto, D.; Galan-Huerta, K.A.; Garza-Treviño, E.N. Inflammatory effect on the gastrointestinal system associated with COVID-19. *World J. Gastroenterol.* **2021**, *27*, 4160–4171. [CrossRef]

6. McGrowder, D.A.; Miller, F.; Anderson Cross, M.; Anderson-Jackson, L.; Bryan, S.; Dilworth, L. Abnormal Liver Biochemistry Tests and Acute Liver Injury in COVID-19 Patients: Current Evidence and Potential Pathogenesis. *Diseases* **2021**, *9*, 50. [CrossRef]

7. Liu, J.; Li, Y.; Liu, Q.; Yao, Q.; Wang, X.; Zhang, H.; Chen, R.; Ren, L.; Min, J.; Deng, F.; et al. SARS-CoV-2 cell tropism and multiorgan infection. *Cell Discov.* **2021**, *7*, 2–5. [CrossRef]

8. Bhatnagar, J.; Gary, J.; Reagan-Steiner, S.; Estetter, L.B.; Tong, S.; Tao, Y.; Denison, A.M.; Lee, E.; Deleon-Carnes, M.; Li, Y.; et al. Evidence of Severe Acute Respiratory Syndrome Coronavirus 2 Replication and Tropism in the Lungs, Airways, and Vascular Endothelium of Patients with Fatal Coronavirus Disease 2019: An Autopsy Case Series. *J. Infect. Dis.* **2021**, *223*, 752–764. [CrossRef]

9. Casagrande, M.; Fitzek, A.; Spitzer, M.; Püschel, K.; Glatzel, M.; Krasemann, S.; Aepfelbacher, M.; Nörz, D.; Lütgehetmann, M.; Pfefferle, S.; et al. Detection of SARS-CoV-2 genomic and subgenomic RNA in retina and optic nerve of patients with COVID-19. *Br. J. Ophthalmol.* **2021**, 1–5. [CrossRef]

10. Campos, R.K.; Camargos, V.N.; Azar, S.R.; Haines, C.A.; Eyzaguirre, E.J.; Rossi, S.L. Sars-cov-2 infects hamster testes. *Microorganisms* **2021**, *9*, 1318. [CrossRef]

11. Ardestani Zadeh, A.; Arab, D. COVID-19 and male reproductive system: Pathogenic features and possible mechanisms. *J. Mol. Histol.* **2021**, *2*, 869–878. [CrossRef] [PubMed]

12. Etievant, S.; Bal, A.; Escuret, V.; Brengel-Pesce, K.; Bouscambert, M.; Cheynet, V.; Generenaz, L.; Oriol, G.; Destras, G.; Billaud, G.; et al. Performance Assessment of SARS-CoV-2 PCR Assays Developed by WHO Referral Laboratories. *J. Clin. Med.* **2020**, *9*, 1871. [CrossRef] [PubMed]

13. Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Müller, M.A.; Niemeyer, D.; Jones, T.C.; Vollmar, P.; Rothe, C.; et al. Virological assessment of hospitalized patients with COVID-2019. *Nature* **2020**, *581*, 465–469. [CrossRef]

14. Kamar, N.; Izopet, J.; Cintas, P.; Garrouste, C.; Uro-Coste, E.; Cointault, O.; Rostaing, L. Hepatitis E Virus-Induced Neurological Symptoms in a Kidney-Transplant Patient with Chronic Hepatitis. *Am. J. Transplant.* **2010**, *10*, 1321–1324. [CrossRef] [PubMed]

15. Dunfee, R.L.; Thomas, E.R.; Gorry, P.R.; Wang, J.; Taylor, J.; Kunstman, K.; Wolinsky, S.M.; Gabuzda, D. The HIV Env variant N283 enhances macrophage tropism and is associated with brain infection and dementia. *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 15160–15165. [CrossRef] [PubMed]

16. Verma, S.; Saksena, S.; Sadri-Ardekani, H. ACE2 receptor expression in testes: Implications in coronavirus disease 2019 pathogenesis. *Biol. Reprod.* **2020**, *103*, 449–451. [CrossRef]

17. Poma, A.M.; Bonuccelli, D.; Giannini, R.; Macerola, E.; Vignali, P.; Ugolini, C.; Torregrossa, L.; Proietti, A.; Pistello, M.; Basolo, A.; et al. COVID-19 autopsy cases: Detection of virus in endocrine tissues. *J. Endocrinol. Investig.* **2021**, *45*, 209–214. [CrossRef]

18. Verma, R.; Kim, E.; Martínez-Colón, G.J.; Jagannathan, P.; Rustagi, A.; Parsonnet, J.; Bonilla, H.; Khosla, C.; Holubar, M.; Subramanian, A.; et al. SARS-CoV-2 Subgenomic RNA Kinetics in Longitudinal Clinical Samples. *Open Forum Infect. Dis.* **2021**, *8*, ofab310. [CrossRef]

19. Alexandersen, S.; Chamings, A.; Bhatta, T.R. SARS-CoV-2 genomic and subgenomic RNAs in diagnostic samples are not an indicator of active replication. *Nat. Commun.* **2020**, *11*, 6059. [CrossRef]

20. Veronese, G.; Ammirati, E.; Chen, C.; Klingel, K.; Suzuki, M.; Okumura, T.; Maisch, B.; Zuo, H.; Ni, L.; Jiang, J.; et al. Management perspectives from the 2019 Wuhan international workshop on fulminant myocarditis. *Int. J. Cardiol.* **2021**, *324*, 131–138. [CrossRef]

21. Wagner, J.U.G.; Bojkova, D.; Shumliakivska, M.; Luxán, G.; Nicin, L.; Aslan, G.S.; Milting, H.; Kandler, J.D.; Dendorfer, A.; Heumueller, A.W.; et al. Increased susceptibility of human endothelial cells to infections by SARS-CoV-2 variants. *Basic Res. Cardiol.* **2021**, *116*, 42. [CrossRef] [PubMed]

22. Rueca, M.; Bartolini, B.; Gruber, C.E.M.; Piralla, A.; Baldanti, F.; Giombini, E.; Messina, F.; Marchioni, L.; Ippolito, G.; Di Caro, A.; et al. Compartmentalized Replication of SARS-CoV-2 in Upper vs. Lower Respiratory Tract Assessed by Whole Genome Quasispecies Analysis. *Microorganisms* **2020**, *8*, 1302. [CrossRef]

23. Jary, A.; Leducq, V.; Malet, I.; Marot, S.; Klement-Frutos, E.; Teyssou, E.; Soulié, C.; Abdi, B.; Wirden, M.; Pourcher, V.; et al. Evolution of viral quasispecies during SARS-CoV-2 infection. *Clin. Microbiol. Infect.* **2020**, *26*, 1560.e1–1560.e4. [CrossRef] [PubMed]

24. Yuan, M.; Liu, H.; Wu, N.C.; Wilson, I.A. Recognition of the SARS-CoV-2 receptor binding domain by neutralizing antibodies. *Biochem. Biophys. Res. Commun.* **2021**, *538*, 192–203. [CrossRef]

25. Mahmoudi Gomari, M.; Rostami, N.; Omidi-Ardali, H.; Arab, S.S. Insight into molecular characteristics of SARS-CoV-2 spike protein following D614G point mutation, a molecular dynamics study. *J. Biomol. Struct. Dyn.* **2021**, *1–9*. [CrossRef] [PubMed]
26. Sharun, K.; Tiwari, R.; Dhama, K. SARS-CoV-2 in semen: Potential for sexual transmission in COVID-19. *Int. J. Surg.* **2020**, *84*, 156–158. [CrossRef]

27. Omolaoye, T.S.; Adeniji, A.A.; Cardona Maya, W.D.; du Plessis, S. SARS-COV-2 (Covid-19) and male fertility: Where are we? *Reprod. Toxicol.* **2020**, *99*, 65–70. [CrossRef] [PubMed]

28. He, Y.; Wang, J.; Ren, J.; Zhao, Y.; Chen, J.; Chen, X. Effect of COVID-19 on Male Reproductive System—A Systematic Review. *Front. Endocrinol. (Lausanne)* **2021**, *12*, 677701. [CrossRef]

29. Holtmann, N.; Edimiris, P.; Andree, M.; Doehmen, C. Assessment of SARS-CoV-2 in human semen—A cohort study. *Fertil. Steril.* **2020**, *114*, 233–238. [CrossRef]

30. Vahedian-Azimi, A.; Karimi, L.; Makvandi, S.; Jamialahmadi, T.; Sahebkar, A. Does SARS-CoV-2 Threaten Male Fertility? In *Clinical, Biological and Molecular Aspects of COVID-19*; Guest, P.C., Ed.; Advances in Experimental Medicine and Biology Series; Springer: Cham, Switzerland, 2021; pp. 139–146.