The Tissue Distribution of SARS-CoV-2 in Transgenic Mice With Inducible Ubiquitous Expression of hACE2

The novel coronavirus disease COVID-19 has become one of the most socially significant infections. One of the main models for COVID-19 pathogenesis study and anti-COVID-19 drug development is laboratory animals sensitive to the virus. Herein, we report SARS-CoV-2 infection in novel transgenic mice conditionally expressing human ACE2 (hACE2), with a focus on viral distribution after intranasal inoculation. Transgenic mice carrying hACE2 under the floxed STOP cassette [(hACE2-LoxP(STOP)] were mated with two types of Cre-ERT2 strains (UBC-Cre and Rosa-Cre). The resulting offspring with temporal control of transgene expression were treated with tamoxifen to induce the removal of the floxed STOP cassette, which prevented hACE2 expression. Before and after intranasal inoculation, the mice were weighed and clinically examined. On Days 5 and 10, the mice were sacrificed for isolation of internal organs and the further assessment of SARS-CoV-2 distribution. Intranasal SARS-CoV-2 inoculation in hACE2-LoxP(STOP)×UBC-Cre offspring resulted in weight loss and death in 6 out of 8 mice. Immunostaining and focus formation assays revealed the most significant viral load in the lung, brain, heart and intestine samples. In contrast, hACE2-LoxP(STOP) × Rosa-Cre offspring easily tolerated the infection, and SARS-CoV-2 was detected only in the brain and lungs, whereas other studied tissues had null or negligible levels of the virus. Histological examination revealed severe alterations in the lungs, and mild changes were observed in the brain tissues. Notably, no changes were observed in mice without tamoxifen treatment. Thus, this novel murine model with the Cre-dependent activation of hACE2 provides a useful and safe tool for COVID-19 studies.

Introduction

A recent outbreak of the novel coronavirus disease (COVID-19), which spread rapidly worldwide in just a few months, has become one of the most socially significant infections, along with influenza and smallpox (McEntire et al., 2021). The causative agent of COVID-19 is a positive single-stranded RNA virus of Coronaviridae found in the genus Betacoronavirus, which also includes severe acute respiratory syndrome-related coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) (with 80 and 50% nucleotide homology, respectively) (Kim et al., 2020). The diameter of viral particles corresponds to a wide range of particle sizes described in the literature, 60–140 nm, and the range of values can be explained by differences in the measurement methods used (Laue et al., 2021). Mature virions include several types of protein molecules divided into structural proteins [specifically, nucleoprotein (N), envelope protein (E), transmembrane glycoprotein (M), and glycoprotein (S)] and nonstructural proteins (including RNA-dependent RNA polymerase) (Kim et al., 2020). The interaction of the SARS-CoV-2 virus with the ACE2 cell surface receptor is mediated through the S glycoprotein embedded in the envelope. In most cases, this S-protein is cleaved by host proteases into S1 and S2 subunits, which are responsible for receptor recognition and membrane fusion, respectively (Lu et al., 2015; Zhang et al., 2020). Human tissue studies have shown an abundant presence of ACE2 receptors in the epithelia of the lung, small intestine, arterial and venous endothelial cells and arterial smooth muscle cells in all organs studied, including the brain (Hamming et al., 2004).

COVID-19 is characterized by high contagiousness and the development of complications that pose a threat to life. During the first days of infection, SARS-CoV-2 in the upper respiratory tract replicates without severe clinical symptoms. A few days after infection, metabolic changes occur in the cells, leading to the apoptosis of alveocytes and the possible occurrence of acute respiratory distress syndrome (ARDS). Infiltration of the lung tissue by rod-nuclear leukocytes leads to an increase in the secretion of proteolytic enzymes, which, together with the kallikrein-kinin system, the hemocoagulation system and the activation of components of the complement system, lead to secondary alteration of the lung tissue and expand the area of damage (Fan et al., 2020; Hanff et al., 2020; Lipcsey et al., 2021). Neurological symptoms of SARS-CoV-2 infection might be associated with systemic inflammatory response syndrome, sepsis, multiorgan failure or postinfectious, immune-mediated complications when the organism responds to the infection. It is believed that this virus can also affect the central nervous system directly or indirectly. Some neurological symptoms, such as encephalitis and anosmia–hyposmia, might be caused by the direct invasion of the central nervous system by the virus (Natoli et al., 2020; Pezzini and Padovani, 2020; Liu et al., 2021).

The pathophysiological patterns of the disease can only be studied by replicating tissue-specific and systemic virus-host interactions. Among the main models for studying COVID-19 pathogenesis and anti-COVID-19 drug development are sensitive laboratory animals. Animal models susceptible to infection with the SARS-CoV-2 virus include lower primates, hamsters, ferrets and mice (Takayama, 2020). The main criterion for the selection of model animals is sensitivity to the virus, which is directly proportional to the degree of conformity of the amino acid sequence of ACE2 in humans and animals. Thus, different homologs of the ACE2 receptor can also lead to fixation of mutations in the virus genome (Damas et al., 2020). The limitations of the use of animal models for certain studies of the virus may be due to an important feature. RNA viruses such as SARS-CoV-2 tend to evolve rapidly, which is associated with a high mutation rate. This genetic variability is modulated by both internal processes and external conditions. Among the internal factors, one can single out the lack of a corrective function of viral RNA-dependent RNA-polymerase in the process of replication, together with spontaneous mutations due to physical and chemical mutagens, which can also be harbored in the genome (Benedetti et al., 2020). External conditions that lead to a change in the virus include variants of receptors with which it interacts in the host organism at the cellular level. However, upon exposure to nonhuman variants of ACE2, the receptor that allows SARS-CoV-2 to enter cells, SARS-CoV-2 may change its properties due to an adaptation process that can lead to virus mutation (Welkers et al., 2021). To reduce SARS-CoV-2 variability and gain the ability to learn details of human SARS-CoV-2 infection, a few transgenic humanized mouse strains with hACE2 (human ACE2) expression have been produced.

This model constitutively expresses the hACE2 protein during embryogenesis and ontogenesis (Jiang et al., 2020; Sun et al., 2020). However, the constitutive expression of hACE2 sensitizes the model to SARS-CoV-2, while the mice are maintained in an animal facility or during transportation: these animals are accessible to technicians who may be infected with COVID-19 during the pandemic. This results in a risk of model animals being infected and producing antibodies, which may lead to incorrect conclusions during experiments involving SARS-CoV-2.

To avoid such difficulties during the research process, a conditional hACE2 switch-on mouse model was created and characterized by conditional expression of hACE2 (Bruter et al., 2021). The basic hallmark of this strain is the ability to spatially and temporally control hACE2 expression provided by the tamoxifen-induced Cre-recombinase-LoxP (STOP) machinery (Cre-ERT2). In brief, Cre recombinase is an enzyme that catalyzes site-specific recombination between two LoxP sites, resulting in the deletion of the LoxP-limited sequence (Kim et al., 2018). Thus, if the strain is mated with Cre mice, Cre recombinase is expressed. After tamoxifen administration, Cre recombinase gains access to the nuclear compartment and removes the STOP cassette, which prevents the expression of STOP-limited hACE2. Another hallmark of this model is that hACE2 is coexpressed with green fluorescent protein (GFP), making it possible to detect the levels of transgene expression without immunostaining or other complex protocols.

In the present work, we aimed to examine the SARS-CoV-2 sensitivity of the F1 generation obtained after crossing hACE2 (LoxP-Stop) with mice expressing tamoxifen-dependent Cre recombinase (Cre-ERT2). Herein, we report severe and mild SARS-CoV-2 infection in our model depending on which Cre-ERT2 strain was utilized for transgene activation.

Materials and Methods

Animals

All animal procedures were carried out on 34 CBAxC57Bl6J mice of both sexes in accordance with the European Convention for the Protection of Vertebrate Animals Used for Experimental and Other Scientific Purposes (CETS No. 123). The experimental design was approved by the Ethical Committee of SRC VB “Vector” Rospotrebnadzor (protocol # 2 GNC VB VECTOR/02-04.2021, May 28, 2021). Groups of hACE2-harboring transgenic mice were obtained by breeding previously described hACE2(LoxP-Stop) mice (Bruter et al., 2021) and two strains of Cre-ERT2-mice: Ndor1^{Tg(UBC−cre/ERT2)1Ejb} and Gt(ROSA)26Sor^{tm1(cre/ERT2)Tyj} (Jackson Laboratory, United States). The Ndor1^{Tg(UBC−cre/ERT2)1Ejb} strain is characterized by the expression of the Cre/ERT2 gene under the control of the UBC promoter. This transgene is integrated into chromosome 2 at the Ndor1 locus. The Gt(ROSA)26Sor^{tm1(cre/ERT2)Tyj} strain is characterized by the expression of the Cre/ERT2 gene under the control of the CAG promoter. This transgene is integrated into chromosome 6 at the Rosa26 locus. Both strains are characterized by the nontissue-specific ubiquitous expression of Cre-ERT2 (Ruzankina et al., 2007; Ventura et al., 2007). Tamoxifen administration was carried out according to a previously published article (Bruter et al., 2021).

Hereinafter, hACE2(LoxP-Stop)✕Ndor1^{Tg(UBC−cre/ERT2)1Ejb} offspring are referred to as UBC-ACE2, and hACE2(LoxP-Stop)✕Ndor1^{Tg(UBC−cre/ERT2)1Ejb} offspring are referred to as Rosa-ACE2. Cre mice were genotyped according to protocols #27737 and #31150 provided by Jackson Laboratory with the forward 5′-AAA​CTG​TTC​AAT​ATG​CTG​AGG​CT-3′ and reverse 5′-AGT​AGT​TGA​GCA​GTG​GCC​TTA​C-3′ primers for hACE2 carried on the DNA cassette.

According to the experimental design, animals were divided into 7 groups: 1) infected UBC-ACE2 mice without tamoxifen administration [UBC-ACE2, (n = 2)]; 2) infected UBC-ACE2 mice treated with tamoxifen (UBC-ACE2^Tx, n = 8); 3) uninfected UBC-ACE2 mice treated with tamoxifen as a histological control [UBC-ACE2^Tx (Uninfected), n = 2], 4) infected Rosa-ACE2 mice without tamoxifen administration [Rosa-ACE2, (n = 2)]; 5) infected Rosa-ACE2 mice treated with tamoxifen (Rosa-ACE2^Tx, n = 9); 6) infected wild-type mice treated with tamoxifen (WT, n = 9); and 7) uninfected WT mice treated with tamoxifen as a histological control (n = 2—see Supplementary Table S1).

SARS-CoV-2 Infection

All mice except 2 WT (Uninfected) and 2 UBC-ACE2^Tx (Uninfected) mice (Figure 1A and Supplementary Table S1) were inoculated with SARS-CoV-2 (SARS-CoV-2/human/AUS/VIC01/2020) via the intranasal route at a dose of 4.5 lg FFU (focus-forming unit) as previously described (Agranovski et al., 2010). All experiments involving the use of live SARS-CoV-2 particles were carried out in a biosafety level-3 facility.

FIGURE 1

FIGURE 1. Evaluation of SARS-CoV-2 sensitivity in mice with Cre-dependent human (h)ACE2 expression. (A) Graphical representation of the experimental design. All the mice in the UBC-ACE2^Tx group died before the 10th day (Created with BioRender.com). (B) Dynamics of weight change after virus challenge; (C) Kaplan–Meier curve showing decreased survival in infected UBC-ACE2^Tx mice.

After viral inoculation, mice were observed daily. Weight was measured (NAVIGATOR™ NVT-2201, OHAUS, United States) before and after (5th and 10th days or just after death) inoculation. The 5th and 10th days were chosen as the endpoints for histological examination and postmortem studies of SARS-CoV-2 distribution and mutagenesis (Jiang et al., 2020) (Figure 1A).

Evaluation of Virus Production

To determine viral biodistribution, the PCR-based determination of SARS-CoV-2 genome equivalent number was performed. Samples were diluted in Opti-MEM tissue culture medium (Gibco, United Kingdom) and homogenized with a TissueLyser LT tissue homogenizer (Qiagen, United Kingdom) according to the manufacturer’s instructions (the weight of tissue samples are presented in Supplementary Table S3). The samples were centrifuged for 5 min at 10,000 g, and the resulting supernatant was collected. The purification of total RNA from 100 μl samples was performed using a “RIBO-sorb” Extraction Kit (AmpliSens, Russia) according to the manufacturer’s recommendations. An internal control exogenous RNA (InterLabService, Russia) was added to each sample before purification to control the isolation, reverse transcription and PCR steps. Reverse transcription, followed by real-time PCR, was performed with a “Reverta-L” reagent kit (InterLabService, Russia) and Vector PCRrv-2019-nCoV-RG test system (SRC VB “Vector” Rospotrebnadzor, Russia) using the following primers: 5′-GTT​GCA​ACT​GAG​GGA​GCC​TTG-3' (forward), 5′-GAG​AAG​AGG​CTT​GAC​TGC​CG-3' (reverse) and 5′-FAM-TACACCAAAAGATCACATTGGCACCCG-BHQ1-3' (probe) according to the manufacturer’s recommendations. PCR was carried out using a RotorGene 6000 thermocycler (Qiagen, United States) with the program recommended in the PCR kit manufacturer’s manual. Sample Ct values were obtained in two fluorescence channels for viral cDNA and the internal control. The data obtained were analyzed using Rotor Gene 6000 Series Software 1.8.17.5. Ct values were recalculated into genome equivalents in samples based on the obtained standard dilution curve, with points from 10⁻¹ to 10⁻⁵ of the initial viral stock used for infection. The standard curve had a high correlation coefficient (R² > 0.99).

Focus Formation Assay

To determine the concentration of live SARS-CoV-2 viral particles in the brain and lung tissues, FFU analysis of the Vero E6 cell line was performed (Supplementary Table S3). A cultured cell monolayer in 96-well plates was washed 2 times with tissue culture medium [DMEM with 300 μg/ml L-glutamine (BioloT, Russia), 100 U/ml penicillin G and 100 μg/ml streptomycin mix (Biolot, Russia), 10 μg/ml DEAE-dextran (Sigma, United States) and 5 μg/ml trypsin TPCK (Thermo FS, United States)]. In the wells of the 96-well plate, a 10-fold dilution series of the investigated viral suspension (tissue homogenate previously obtained with TissueLyser LT (Qiagen, United Kingdom)) was added. Then, 100 μl of each dilution of the virus was added to 4 wells. After that, 200 μl of cooled acetone (80%) was added for 15–20 min, then the acetone was removed and the plate was dried. Fifty microliters of a 1:1000 dilution of human monoclonal anti-SARS-CoV-2 antibody (COVID-19, 2019-nCoV) and nucleoprotein/NP monoclonal antibodies (FPZ0562, Fapon Biotech Inc., China) in PBS (Thermo FS, United States) was added to each well for 60 min. After that, the wells were washed two times and incubated with secondary rabbit antibodies against human IgG conjugated with HRP (AP309P, Sigma, United States) for 1 h. Signal detection was performed using an AEC substrate (Sigma, United States) according to the manufacturer’s protocol. The number of focus-forming units was counted in each well using an Olympus CK40 microscope (Olympus, United States). The virus titer was calculated as previously described (Sloutskin and Goldstein, 2014).

Histological Analysis

For histological analysis, tissue samples were fixed in a 4% paraformaldehyde solution (Sigma, United States) for 48 h, followed by washing in water for 10 min. The samples were sequentially dehydrated in ethanol at increasing concentrations followed by processing with a xylol-paraffin mixture. Next, the samples were fixed in paraffin blocks using the Tissue-Tek VIP 6 AI Vacuum Infiltration Processor (Sakura Finetek, Japan) according to the manufacturer’s protocols. Paraffin sections 4–5 μm thick were prepared using a Microm HM 360 Automated Microtome (Marshall Scientific, United States). Sections were stained with hematoxylin and eosin using an automated Multiple Stainer Tissue-Tek Prisma tissue stainer (Sakura Finetek, Japan) according to the manufacturer’s protocol. Light-optical examination and microphotography were performed with an Imager Z1 microscope (Zeiss, Germany). An Olympus VS200 slide scanner was used to obtain overview images and ensure the same imaging parameters. Images were captured using a UPLXAPO 20×/0.80 lens (Olympus, United States).

Immunohistochemical Analysis

For immunochemical analysis, slides of brain tissues were deparaffinized in xylene and rehydrated in a graded ethanol series. As an additional control, brain samples of SARS-CoV-2-inoculated hamsters were used. Antigen retrieval was conducted in a humidified chamber in a 37°C incubator using 0.05% trypsin (Sigma, United States) solution for 10 min. After cooling to room temperature, the slides were incubated overnight at 4°C with a nucleoprotein/NP monoclonal antibody (FPZ0562, Fapon Biotech Inc., China) diluted 1:500 in blocking buffer (1% BSA (BioloT, Russia) in TBS with Tween-20 (Abcam, United States)). After three washes in TBS (Abcam, United States), the slides were treated with a secondary Cy3-conjugated anti-human antibody (AB_2337718, Jackson ImmunoResearch, United States) diluted at 1:200 in blocking buffer. Additionally, slides were stained with DAPI and DABCO for nuclear visualization. Visualization was performed with an EVOS™ XL Core Imaging System (Thermo Fisher Scientific, United States).

High-Throughput Sequencing and Data Analysis

To analyze virus nucleotide sequences, cDNA libraries were amplified from RNA samples using the ARTIС and SISPA protocols (Reyes and Kim, 1991; DNA Pipelines et al., 2020). Furthermore, the libraries were prepared with the NEBNext Ultra II DNA Library Prep Kit for Illumina (New England BioLabs, United States) and sequenced (2 × 300-bp cycles) using the Illumina MiSeq platform (Illumina, United States). Raw reads were first trimmed using cutadapt 3.4 (parameters: -error-rate 0.1 –times 10 –overlap 3 –minimum-length 20 –pair-filter = any) to remove adapter and primer sequences. fastp 0.20.0 (Chen et al., 2018) (parameters: -× -3-5 -r -l 20) was used to remove low-quality sequences and reads. The data were concatenated and de novo assembled with the coronaSPAdes module (Meleshko et al., 2021) from the SPAdes 3.15.2 assembler (Bankevich et al., 2012) using the default settings. Reads were mapped with BWA 0.7.17 (Li and Durbin, 2009) to the assembled sequence, and PCR duplicates were removed using Picard MarkDuplicates 2.25.4 (Toolkit, 2019) for SISPA-amplified libraries, followed by variant calling with iVar (Grubaugh et al., 2019). Variants were called at ≥50% in reads that were supported by a minimum of 10 reads and had a MAPQ ≥ 30. Only variants present in both ARTIC- and SISPA-amplified libraries were considered mutations.

Statistical Analyses

Statistical analyses were performed in the R software environment v4.1.1. As many samples did not display a normal distribution (Shapiro–Wilk test and the Spiegelhalter test, “normtest” package) or variance equality (Levene’s test, “lawstat” package), nonparametric methods were used. For quantitative variables, the significance of the obtained results was determined using Kruskal–Wallis one-way analysis of variance and Dunn’s test of multiple comparisons using rank sums with the Benjamini–Hochberg procedure to decrease the false discovery rate (“dunn.test” package) as post-hoc analyses to identify significant differences in intergroup comparisons. GraphPad Prism 9.2.0 software was also used to create graphic material and for survival assessment via the Kaplan–Meier method together with log-rank (Mantel–Cox) and Gehan-Breslow-Wilcoxon (GBW) tests corrected by the Benjamini–Hochberg procedure. Each statistically processed sample contained at least two values. The significance threshold for all performed tests was set to p ≤ 0.05.

Results

Weight Loss After Infection

Among all SARS-CoV-2-infected groups, only UBC-ACE2^Tx mice displayed a severe clinical picture of infection. For a more demonstrative representation of the relative weight changes shown in the plot, the following formula was used (Figure 1B and Supplementary Table S2):

$Weight Change = \frac{100 +\left(Weight on the Nth day - Weight on the 0th day\right)}{Weight on the 0th day}\ast 100$

A significant difference in mortality was observed between the UBC-ACE2^Tx group and the other groups. In this group, we observed dramatic weight loss and rapid mortality. Weighing on the 5th day revealed an ∼24% weight decrease [p = 0.0021 and 0.0004, compared to the Rosa-ACE2^Tx and WT groups, respectively (Figure 1B: # and §)] and an ∼35% decrease up to the moment the last animal died on the 8^th day.

Intergroup comparison of survival probability revealed that UBC-ACE2^Tx animals had a significantly higher probability of death than mice from the Rosa-ACE2^Tx and WT groups (in both comparisons p = 0.002 and 0.006, for the Mantel–Cox and GBW tests, respectively) (Figure 1C). Notably, although a portion of the animals were sacrificed and censored on the 5th day, based on the remaining animals, the Kaplan–Meier curve predicted 100% mortality by the 10th day. Neither the UBC-ACE2 without tamoxifen nor the Rosa-ACE2^Tx group showed any clinical sign of infection, mortality or statistically significant weight change, suggesting no dramatic infection process in these groups.

SARS-CoV-2 Distribution in Organs

Next, PCR and FFU analyses confirmed SARS-CoV-2 invasion and clarified its tissue distribution. On the 5th day, tissue samples were collected from all sacrificed and deceased animals. Intergroup comparisons of the obtained data were carried out only between the Rosa-ACE2^Tx, UBC-ACE2^Tx and WT mice due to the small sample size in the remaining groups (Figure 2 and Supplementary Table S3).

FIGURE 2

FIGURE 2. SARS-CoV-2 genome equivalent (A) and focus-forming unit (B) analyses in different samples of mouse organs on the 5th day of infection. (A) Significant differences between WT and UBC-ACE2^Tx and between Rosa-ACE2^Tx in heart and brain; significant differences between overall UBC-ACE2^Tx, Rosa-ACE2^Tx and WT in lungs. (B) Significant differences between WT and UBC-ACE2Tx and between Rosa-ACE2^Tx in lungs and brain; Note: *p < 0.05, **p < 0.001, #p < 0.1.

Analysis of genome equivalents in lung samples revealed significant differences between all compared groups on the 5th day (Figure 2A). While the median level of viral genome equivalents in UBC-ACE2^Tx mice was ∼8 times higher than that in Rosa-ACE2^Tx mice (8465.05 vs. 1048.5, respectively, p = 0.0460), no virus was detected in the lungs of WT animals (p = 0.0004 and 0.0437 compared to UBC-ACE2^Tx and Rosa-ACE2^Tx mice, respectively). The same results were found during FFU analysis (Figure 2B): the UBC-ACE2^Tx animals again had the highest amount of viable virus particles, which was ∼91 times more than the median value in the Rosa-ACE2^Tx mice (8206.69 vs. 90.25, p = 0.0285). In contrast to the UBC-ACE2^Tx and Rosa-ACE2^Tx groups, there were no detected viable virus particles in lungs collected from WT animals. Similarly, the brain samples of UBC-ACE2^Tx and Rosa-ACE2^Tx animals were characterized by the presence of viral genome equivalents and viable virus particles. Only UBC-ACE2^Tx mice showed the presence of viral genome equivalents in the heart, liver and intestine (Figure 2).

In the brain samples of #1 Rosa-ACE2 and #2 UBC-ACE2 mice, the presence of live viral particles was detected in the FFU assay, with values of 609.14 and 547.45 FFU normalized to the weight of the sample (Supplementary Table S3). During the PCR study of the genome equivalent, we did not detect the RNA of the virus in the tissue samples of these animals. This is because infected samples are inactivated at high temperature for PCR analysis, according to the standard protocol for working with pathogenic agents, which leads to the degradation of a part of the genomic RNA and to a decrease in the sensitivity of the method. Thus, the FFU assay is the more sensitive method of analysis in this case and allows the detection of even insignificant levels of infection in laboratory mice even in the absence of transgenesis, which was previously shown (Leist et al., 2020). However, in Rosa-ACE2 and UBC-ACE2 without tamoxifen administration, the FFU values were significantly lower than those in the experiment, and such an insignificant level of viral load was not pathogenic to laboratory mice.

The same PCR and FFU analyses were also performed on the 10th day (Supplementary Table S3). However, these data are not used for statistical estimation of virus distribution in organs due to the death of UBC-ACE2^Tx mice and the absence of the viral genome in Rosa-ACE2^Tx and WT samples.

High-Throughput Sequencing Data Analysis

To analyze changes in the virus genome during mouse infection, we sequenced a virus from the original stock as well as postinfection virus samples from the brain and lung from #6 UBC-ACE2^Tx as a representative of the UBC-ACE2^Tx group. Here, using metatranscriptomic data, we obtained an almost complete genome (99% of reference genome SARS-CoV-2 coverage) with an average read depth of 277×. The sequencing reads were deposited in the NCBI SRA database under the accession number PRJNA780672. The assembled genomes were deposited in the GenBank database under accession number OL677207.1. None of the tissue-derived SARS-CoV-2 sequences contained any mutations relative to the original viral stock used for mouse infection.

SARS-CoV-2-Induced Tissue Lesions

Histological examination of animal tissue samples was conducted on the 5th day of infection (Figure 3). Lung samples from both UBC-ACE2 and Rosa-ACE2 mice without tamoxifen administration showed a normal tissue structure, as did the WT mouse group. Histological sections (Figures 3A3–4,B) show vessels, terminal respiratory bronchioles lined with ciliated columnar epithelium and alveolar ducts, sacs and alveoli lined with simple squamous epithelium that is interrupted with club cells on top of the septs. Unlike previous samples, UBC-ACE2^Tx and Rosa-ACE2^Tx mice showed pathological changes in lung tissues (Figures 3A1,2). These changes include the septal thickening of alveolar ducts and sacs caused by the diffuse hyperplasia of type II alveolar epithelial cells with mild lymphocytic infiltration and edema. In the bronchioles, the diffuse hyperplasia of the epithelium was observed. Moreover, lung tissue samples from UBC-ACE2^Tx and Rosa-ACE2^Tx mice were distinguished by multiple sludges of erythrocytes in vessels. This pattern of pathological changes corresponds to the proliferative phase of diffuse alveolar damage.

FIGURE 3

FIGURE 3. Histological alterations in organs on the 5th day after SARS-CoV-2 inoculation. (A) Hematoxylin/eosin histostained lung tissue samples from UBC-ACE2^Tx(A1), UBC-ACE2 (A3), Rosa-ACE2^Tx (A2) and Rosa-ACE2 (A4) groups. (C) Hematoxylin/eosin histostained brain samples from the UBC-ACE2^Tx (C1), UBC-ACE2 (C3), Rosa-ACE2^Tx (C2) and Rosa-ACE2 (C4) groups. The pattern of pathological changes in UBC-ACE2^Tx(A1) and Rosa-ACE2^Tx (A2) lung samples included the septal thickening of alveolar ducts and sacs caused by the diffuse hyperplasia of type II alveolar epithelial cells with mild lymphocytic infiltration and edema (green arrow), hyperplasia of the epithelium of bronchioles (blue arrow) and multiple sludges of erythrocytes in vessels (black arrow). Examination of brain samples revealed multiple sludges of erythrocytes (black arrow) only in brain vessels of UBC-ACE2^Tx (C1) mice. In both the UBC-ACE2^Tx− and Rosa-ACE2^Tx-infected groups, tissue lesions were detected only after tamoxifen administration. (B), (D) WT mice displayed no damage in the studied organs.

In contrast to lung samples, histological analysis of the brain samples did not reveal any significant changes among all groups (Figures 3C,D). Only multiple sludges of erythrocytes were observed in the vessels of UBC-ACE2^Tx mice (Figure 3C1). UBC-ACE2^Tx (Uninfected) and WT (Uninfected) mice treated with tamoxifen showed no pathological changes in tissue samples, indicating the absence of an effect of tamoxifen on the histological assay.

Immunohistochemical Analysis

Since there is no consensus on SARS-CoV-2 neurotropism and neuroinvasion ability (Natoli et al., 2020; Pezzini and Padovani, 2020; Liu et al., 2021), we further visualized the presence of viral proteins in UBC-ACE2^Tx brain samples. The most consideration was given to the dentate gyrus of the hippocampus due to its connection with the lateral entorhinal cortex that receives olfactory inputs via direct projections from the olfactory bulb (Woods et al., 2020). The virus has been shown to be located in neurons, which were distinguishable in slides as cells with soma and dendrites (Figure 4, Cy3). In mice without tamoxifen treatment (UBC-ACE2), infected WT mice and infected hamsters, specific signals indicating the presence of the SARS-CoV-2 N-protein were not observed (Supplementary Figures S1A–C).

FIGURE 4

FIGURE 4. Immunohistochemical detection of SARS-CoV-2 N-protein in UBC-ACE2^Tx brain samples. Viral N-protein was labeled with Cy3-conjugated secondary antibody (red field). Cell nuclei and DNA were stained with DAPI (blue field).

Discussion

In this article, we present the results of a study investigating the susceptibility of switch-on mice with inducible hACE2 expression to SARS-CoV-2 infection (Bruter et al., 2021). As the LoxP-STOP system provides temporal and spatial control of transgene expression, SARS-CoV sensitivity in ACE2-LoxP(STOP) mice might be precisely regulated by breeding with various Cre strains. Indeed, we show here that the UBC-ACE2 and Rosa-ACE2 strains are characterized by different severities of infection course. Although SARS-CoV-2 invasion was confirmed in both strains, only UBC-ACE2 displayed acute weight loss and mortality. However, weight loss can be associated with either a viral infection or renin-angiotensin system since ACE2 needs to balance AngII levels. Its changes can lead to mouse feeding behavior alterations and consequently weight loss (Shorning et al., 2012). Unequal SARS-CoV-2 sensitivity might be related to different levels of transgene expression, consistent with previous data on UBC-driven and ROSA-26 locus-integrated promoter activity (Chen et al., 2011). Thus, two strains express the same Cre-ERT2 enzyme; however, the difference in its expression enables the resulting offspring to serve various research purposes.

When modeling a human disease using animal systems, we must estimate and consider the similarities and differences in the pathogenetic change patterns, for example, when developing antiviral therapy methods. To do this, we studied pathological changes in tissues during infection with the virus. Although SARS-CoV-2 causes severe pathology of the respiratory system, there is evidence that it affects the central nervous system in humans (Natoli et al., 2020; Pezzini and Padovani, 2020; Liu et al., 2021). For example, SARS-CoV-2 mRNA has been recovered from cerebrospinal fluid, suggesting that it can cross the blood–brain barrier (Moriguchi et al., 2020). Moreover, the virus has been found in brain tissues (Ren et al., 2021). It has also been shown that the penetration of the virus into the brain causes cerebral circulation (Helms et al., 2020). Some animal models do not show the presence of the virus in the brain. For example, the virus was not found in brain samples from the Syrian hamster. This has been described previously but was also confirmed in this study by immunohistochemistry. Our research showed that in the UBC-ACE2 line, the virus penetrated directly into neurons in the gyrus dentate area. Considering that this part of the hippocampus receives olfactory inputs (Woods et al., 2020) and SARS-CoV-2 enters the central nervous system through the olfactory tract (Meinhardt et al., 2021), model mice can be used to study the interplay between COVID-19 and the nervous system. Although Cre-ERT2 is ubiquitously expressed, the analysis of SARS-CoV-2 distribution revealed the most significant viral load in brain and lung tissues. Apparently, these findings are related to the intranasal route of viral inoculation and may be similar to SARS-CoV-2 distribution in humans (Chen et al., 2021). We also provided evidence that mice without tamoxifen administration are resistant to SARS-CoV-2, suggesting that there is no undesired leakage in hACE2 expression; hence, the model is biosafe until an experiment is started.

It is also important to study whether the virus undergoes mutations during the mouse infection process. This information is necessary to understand in which studies these mice can be used, for example, antiviral drug research in which two factors are important: maintaining viral replication in mice and an adequate infection pattern. Both study types involve a search for differences in different viral strain activities, and thus, the use of certain cellular receptors is necessary for the correct modeling of the infection process. RNA-containing viruses are capable of rapid changes in their primary sequence, which is associated with their increased mutagenicity, leading to worldwide spread. For SARS-CoV-2, the use of a human ACE2 receptor can likely reduce its variability due to the lack of targeted selection for other types of this virus in other animal models. Indeed, no change in the viral sequence was detected during mouse infection.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (and its supporting information files). The sequencing reads were deposited in the NCBI SRA database under the accession number PRJNA780672. The assembled genomes were deposited in the GenBank database under accession number OL677207.1.

Ethics Statement

The animal study was reviewed and approved by Experiments on laboratory animals were carried out in accordance with protocol #GNC VB VECTOR/02-04.2021 supported by the Ethical Committee of SRC VB “Vector” Rospotrebnadzor by protocol #2 May 28, 2021.

Author Contributions

Conceptualization, AAD and DY; Data curation, AAD, AG, and DY; Formal analysis, AAD, AG, OT, IG, EMS, DP, EVS, SB, ND, AM, PT, YB, GO, MK, VS, EI, TT, AR, and DY; Investigation, AAD and AG; Methodology, AAD, AG, OT, IG, DP, VS, SB, TT, and AVD; Project administration, EG, RM, and DY; Resources, MK, VS, AB, AR, RM, OT, AVD, and DY; Visualization, AAD, IG, and VS; Writing, original draft, AAD, EMS, DP, VS, IG, and AG; Writing, review and editing, DY and AVD.

Funding

This work was supported by the Ministry of Science and Higher Education of the Russian Federation (agreement # 075-15-2019-1665).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.821506/full#supplementary-material

Supplementary Figure S1 | Immunohistochemistry control. Immunohistochemical detection of SARS-CoV-2 particles in UBC-ACE2 brain samples. The N-protein of SARS-CoV-2 (red field) was not observed in either WT mice (A) or UBC-ACE2 without tamoxifen treatment (B). Additionally, we did not find N-protein-positive signals in hamster brain samples, which is consistent with previous literature data related to SARS-CoV-2 infection patterns in these animals (C). Cell nuclei were stained with DAPI (blue field).

Supplementary Table S1 | The experimental design. In the table, the number of animals, their genotypes, experimental conditions and outcomes are summarized.

Supplementary Table S2 | Changes in the weight of animals. The table discloses all primary and processed data concerning weight dynamics.

Supplementary Table S3 | SARS-CoV-2 genome equivalent number and focus-forming unit analysis. Virus production is presented as genome equivalents measured for three groups of mice. Here, we describe additional data about virus detection in different organs (brain, lungs, heart, etc.). The results of FFU analysis are also shown.

References

Agranovski, I. E., Pyankov, O. V., Pyankova, O. G., Sergeev, A. A., Sergeev, A. N., Smetannikova, M. A., et al. (2010). Development of a New Procedure for Precise Determination of Viral Aerosol Lethal Dose (ALD50) for Birds. J. Aerosol Sci. 41, 161–169. doi:10.1016/j.jaerosci.2009.10.003

CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: A New Genome Assembly Algorithm and its Applications to Single-Cell Sequencing. J. Comput. Biol. 19, 455–477. doi:10.1089/cmb.2012.0021

CrossRef Full Text | Google Scholar

Benedetti, F., Pachetti, M., Marini, B., Ippodrino, R., Ciccozzi, M., and Zella, D. (2020). SARS‐CoV‐2: March toward Adaptation. J. Med. Virol. 92, 2274–2276. doi:10.1002/jmv.26233

CrossRef Full Text | Google Scholar

Bruter, A. V., Korshunova, D. S., Kubekina, M. V., Sergiev, P. V., Kalinina, A. A., Ilchuk, L. A., et al. (2021). Novel Transgenic Mice with Cre-dependent Co-expression of GFP and Human ACE2: a Safe Tool for Study of COVID-19 Pathogenesis. Transgenic Res. 30, 289–301. doi:10.1007/s11248-021-00249-8

CrossRef Full Text | Google Scholar

Chen, C.-m., Krohn, J., Bhattacharya, S., and Davies, B. (2011). A Comparison of Exogenous Promoter Activity at the ROSA26 Locus Using a PhiC31 Integrase Mediated Cassette Exchange Approach in Mouse ES Cells. PLoS One 6, e23376. doi:10.1371/journal.pone.0023376

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R., Wang, K., Yu, J., Howard, D., French, L., Chen, Z., et al. (2021). The Spatial and Cell-type Distribution of SARS-CoV-2 Receptor ACE2 in the Human and Mouse Brains. Front. Neurol. 11, 573095. doi:10.3389/fneur.2020.573095

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an Ultra-fast All-In-One FASTQ Preprocessor. Bioinformatics 34, i884–i890. doi:10.1093/bioinformatics/bty560

PubMed Abstract | CrossRef Full Text | Google Scholar

Damas, J., Hughes, G. M., Keough, K. C., Painter, C. A., Persky, N. S., Corbo, M., et al. (2020). Broad Host Range of SARS-CoV-2 Predicted by Comparative and Structural Analysis of ACE2 in Vertebrates. Proc. Natl. Acad. Sci. USA 117, 22311–22322. doi:10.1073/pnas.2010146117

PubMed Abstract | CrossRef Full Text | Google Scholar

Dna Pipelines R&D. Farr, B., Rajan, D., Betteridge, E., Shirley, L., Quail, M., et al. (2020). COVID-19 ARTIC V3 Illumina Library Construction and Sequencing Protocol (Hinxton, United Kingdom: Wellcome Sanger Institute). Available at: dx.doi.org/10.17504/protocols.io.bibtkann.

Google Scholar

Fan, E., Beitler, J. R., Brochard, L., Calfee, C. S., Ferguson, N. D., Slutsky, A. S., et al. (2020). COVID-19-associated Acute Respiratory Distress Syndrome: Is a Different Approach to Management Warranted? Lancet Respir. Med. 8, 816–821. doi:10.1016/S2213-2600(20)30304-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Grubaugh, N. D., Gangavarapu, K., Quick, J., Matteson, N. L., De Jesus, J. G., Main, B. J., et al. (2019). An Amplicon-Based Sequencing Framework for Accurately Measuring Intrahost Virus Diversity Using PrimalSeq and iVar. Genome Biol. 20, 8. doi:10.1186/s13059-018-1618-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamming, I., Timens, W., Bulthuis, M., Lely, A., Navis, G., and van Goor, H. (2004). Tissue Distribution of ACE2 Protein, the Functional Receptor for SARS Coronavirus. A First Step in Understanding SARS Pathogenesis. J. Pathol. 203, 631–637. doi:10.1002/path.1570

CrossRef Full Text | Google Scholar

Hanff, T. C., Harhay, M. O., Brown, T. S., Cohen, J. B., and Mohareb, A. M. (2020). Is There an Association between COVID-19 Mortality and the Renin-Angiotensin System? A Call for Epidemiologic Investigations. Clin. Infect. Dis. 71, 870–874. doi:10.1093/cid/ciaa329

PubMed Abstract | CrossRef Full Text | Google Scholar

Helms, J., Kremer, S., Merdji, H., Clere-Jehl, R., Schenck, M., Kummerlen, C., et al. (2020). Neurologic Features in Severe SARS-CoV-2 Infection. N. Engl. J. Med. 382, 2268–2270. doi:10.1056/NEJMc2008597

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, R.-D., Liu, M.-Q., Chen, Y., Shan, C., Zhou, Y.-W., Shen, X.-R., et al. (2020). Pathogenesis of SARS-CoV-2 in Transgenic Mice Expressing Human Angiotensin-Converting Enzyme 2. Cell 182, 50–58. doi:10.1016/j.cell.2020.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Lee, J.-Y., Yang, J.-S., Kim, J. W., Kim, V. N., and Chang, H. (2020). The Architecture of SARS-CoV-2 Transcriptome. Cell 181, 914–921. doi:10.1016/j.cell.2020.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H., Kim, M., Im, S.-K., and Fang, S. (2018). Mouse Cre-LoxP System: General Principles to Determine Tissue-specific Roles of Target Genes. Lab. Anim. Res. 34, 147. doi:10.5625/lar.2018.34.4.147

PubMed Abstract | CrossRef Full Text | Google Scholar

Laue, M., Kauter, A., Hoffmann, T., Möller, L., Michel, J., and Nitsche, A. (2021). Morphometry of SARS-CoV and SARS-CoV-2 Particles in Ultrathin Plastic Sections of Infected Vero Cell Cultures. Sci. Rep. 11, 3515. doi:10.1038/s41598-021-82852-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Leist, S. R., Dinnon, K. H., Schäfer, A., Tse, L. V., Okuda, K., Hou, Y. J., et al. (2020). A Mouse-Adapted SARS-CoV-2 Induces Acute Lung Injury and Mortality in Standard Laboratory Mice. Cell 183, 1070–1085. doi:10.1016/j.cell.2020.09.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipcsey, M., Persson, B., Eriksson, O., Blom, A. M., Fromell, K., Hultström, M., et al. (2021). The Outcome of Critically Ill COVID-19 Patients Is Linked to Thromboinflammation Dominated by the Kallikrein/Kinin System. Front. Immunol. 12, 211. doi:10.3389/fimmu.2021.627579

CrossRef Full Text | Google Scholar

Liu, J. M., Tan, B. H., Wu, S., Gui, Y., Suo, J. L., and Li, Y. C. (2021). Evidence of central Nervous System Infection and Neuroinvasive Routes, as Well as Neurological Involvement, in the Lethality of SARS‐CoV‐2 Infection. J. Med. Virol. 93, 1304–1313. doi:10.1002/jmv.26570

CrossRef Full Text | Google Scholar

Lu, G., Wang, Q., and Gao, G. F. (2015). Bat-to-human: Spike Features Determining 'host Jump' of Coronaviruses SARS-CoV, MERS-CoV, and beyond. Trends Microbiol. 23, 468–478. doi:10.1016/j.tim.2015.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

McEntire, C. R. S., Song, K.-W., McInnis, R. P., Rhee, J. Y., Young, M., Williams, E., et al. (2021). Neurologic Manifestations of the World Health Organization's List of Pandemic and Epidemic Diseases. Front. Neurol. 12, 161. doi:10.3389/fneur.2021.634827

CrossRef Full Text | Google Scholar

Meinhardt, J., Radke, J., Dittmayer, C., Franz, J., Thomas, C., Mothes, R., et al. (2021). Olfactory Transmucosal SARS-CoV-2 Invasion as a Port of central Nervous System Entry in Individuals with COVID-19. Nat. Neurosci. 24, 168–175. doi:10.1038/s41593-020-00758-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Meleshko, D., Hajirasouliha, I., and Korobeynikov, A. (2021). coronaSPAdes: from Biosynthetic Gene Clusters to RNA Viral Assemblies. Bioinformatics 0728, 224584. doi:10.1093/bioinformatics/btab597

CrossRef Full Text | Google Scholar

Moriguchi, T., Harii, N., Goto, J., Harada, D., Sugawara, H., Takamino, J., et al. (2020). A First Case of Meningitis/encephalitis Associated with SARS-Coronavirus-2. Int. J. Infect. Dis. 94, 55–58. doi:10.1016/j.ijid.2020.03.062

CrossRef Full Text | Google Scholar

Natoli, S., Oliveira, V., Calabresi, P., Maia, L. F., and Pisani, A. (2020). Does SARS‐Cov‐2 Invade the Brain? Translational Lessons from Animal Models. Eur. J. Neurol. 27, 1764–1773. doi:10.1111/ene.14277

PubMed Abstract | CrossRef Full Text | Google Scholar

Pezzini, A., and Padovani, A. (2020). Lifting the Mask on Neurological Manifestations of COVID-19. Nat. Rev. Neurol. 16, 636–644. doi:10.1038/s41582-020-0398-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, A. L., Digby, R. J., and Needham, E. J. (2021). Neurological Update: COVID-19. J. Neurol. 268, 4379–4387. doi:10.1007/s00415-021-10581-y

CrossRef Full Text | Google Scholar

Reyes, G. R., and Kim, J. P. (1991). Sequence-independent, Single-Primer Amplification (SISPA) of Complex DNA Populations. Mol. Cell Probes 5, 473–481. doi:10.1016/S0890-8508(05)80020-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruzankina, Y., Pinzon-Guzman, C., Asare, A., Ong, T., Pontano, L., Cotsarelis, G., et al. (2007). Deletion of the Developmentally Essential Gene ATR in Adult Mice Leads to Age-Related Phenotypes and Stem Cell Loss. Cell Stem Cell 1, 113–126. doi:10.1016/j.stem.2007.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Shorning, B. Y., Jardé, T., McCarthy, A., Ashworth, A., de Leng, W. W. J., Offerhaus, G. J. A., et al. (2012). Intestinal Renin-Angiotensin System Is Stimulated after Deletion of Lkb1. Gut 61, 202–213. doi:10.1136/gutjnl-2011-300046

PubMed Abstract | CrossRef Full Text | Google Scholar

Sloutskin, A., and Goldstein, R. (2014). Infectious Focus Assays and Multiplicity of Infection (MOI) Calculations for Alpha-Herpesviruses. BIO-PROTOCOL 4, e1295. doi:10.21769/BioProtoc.1295

CrossRef Full Text | Google Scholar

Sun, S.-H., Chen, Q., Gu, H.-J., Yang, G., Wang, Y.-X., Huang, X.-Y., et al. (2020). A Mouse Model of SARS-CoV-2 Infection and Pathogenesis. Cell Host & Microbe 28, 124–133. doi:10.1016/j.chom.2020.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Takayama, K. (2020). In Vitro and Animal Models for SARS-CoV-2 Research. Trends Pharmacol. Sci. 41, 513–517. doi:10.1016/j.tips.2020.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Toolkit, P. (2019). Broad institute, GitHub Repository. Available at: http//broadinstitute. github. io/picard.

Google Scholar

Ventura, A., Kirsch, D. G., McLaughlin, M. E., Tuveson, D. A., Grimm, J., Lintault, L., et al. (2007). Restoration of P53 Function Leads to Tumour Regression In Vivo. Nature 445, 661–665. doi:10.1038/nature05541

PubMed Abstract | CrossRef Full Text | Google Scholar

Welkers, M. R. A., Han, A. X., Reusken, C. B. E. M., and Eggink, D. (2021). Possible Host-Adaptation of SARS-CoV-2 Due to Improved ACE2 Receptor Binding in Mink. Virus. Evol. 7, veaa094. doi:10.1093/ve/veaa094

PubMed Abstract | CrossRef Full Text | Google Scholar

Woods, N. I., Stefanini, F., Apodaca-Montano, D. L., Tan, I. M. C., Biane, J. S., and Kheirbek, M. A. (2020). The Dentate Gyrus Classifies Cortical Representations of Learned Stimuli. Neuron 107, 173–184. doi:10.1016/j.neuron.2020.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Penninger, J. M., Li, Y., Zhong, N., and Slutsky, A. S. (2020). Angiotensin-converting Enzyme 2 (ACE2) as a SARS-CoV-2 Receptor: Molecular Mechanisms and Potential Therapeutic Target. Intensive Care Med. 46, 586–590. doi:10.1007/s00134-020-05985-9

PubMed Abstract | CrossRef Full Text | Google Scholar