# Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing

Nathalie Wurtz, Maelle Boussier, Louis Souville, Gwilherm Penant, Alexandre Lacoste, Philippe Colson, Bernard La Scola, Sarah Aherfi

► To cite this version:

Nathalie Wurtz, Maelle Boussier, Louis Souville, Gwilherm Penant, Alexandre Lacoste, et al., Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing. Viruses, 2023, 15 (2), pp.268. 10.3390/v15020268. hal-04082353

HAL Id: hal-04082353

https://amu.hal.science/hal-04082353

Submitted on 21 Jun 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
Article

# Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing

Nathalie Wurtz <fnref n="1" />,<fnref n="2" />, Maelle Boussier <fnref n="1" />,<fnref n="2" />, Louis Souville <fnref n="1" />,<fnref n="2" />, Gwilherm Penant <fnref n="3" />, Alexandre Lacoste <fnref n="4" />, Philippe Colson <fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID, Bernard La Scola <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,* ID and Sarah Aherfi <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*

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

check for
updates

**Citation:** Wurtz, N.; Boussier, M.; Souville, L.; Penant, G.; Lacoste, A.; Colson, P.; La Scola, B.; Aherfi, S. Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing. *Viruses* **2023**, *15*, 268. https://doi.org/10.3390/v15020268

Academic Editor: Benjamin Liu

Received: 24 November 2022

Revised: 22 December 2022

Accepted: 16 January 2023

Published: 17 January 2023

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

**Copyright:** © 2023 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:** Detecting and monitoring viruses in wastewater samples have been reported as useful ways of tracking SARS-CoV-2 epidemic trends. However, there is currently no unanimously recognised method of processing samples to identify and quantify SARS-CoV-2 variants in wastewater. We aimed to implement a method that was as simple as possible in order to be used universally. In a study performed between January 2022 and June 2022 in the city of Marseille, France, we first evaluated the impact of the sample preservation strategy. We then compared ultracentrifugation to ultrafiltration and several steps of filtration to determine the optimal approach for virus concentration. As a proof-of-concept, the definitive protocol was applied to next-generation sequencing of SARS-CoV-2 in wastewater to monitor the emergence of the Omicron variant in the city. For sewage water to be processed in the week following the sampling, storage at +4 °C is sufficient, with less than 1 Ct loss. Filtration with a 5 µm syringe filter, then with a 0.8 µm filtration unit, followed by ultrafiltration was the optimal protocol, leading to an average increase of 3.24 Ct when the starting Ct was on average 38 in the wastewater. This made it possible to observe the emergence of the Omicron 21L/BA.2 variant after Omicron 21K/BA.1 by genome sequencing over a period ranging from 20 February to 10 April 2022 in agreement with observations based on patient data. To conclude, by using a simple method requiring only basic filters and a centrifuge as equipment, it is possible to accurately track the relative incidence rates and the emergence of SARS-CoV-2 variants based on sewage samples.

**Keywords:** SARS-CoV-2; variant; Omicron; sewage; protocol; qPCR; next-generation sequencing; genomics

## 1. Introduction

The emergence of new SARS-CoV-2 variants, associated with either increased transmissibility, escaping antibody responses, or both, have recently attracted attention [1,2]. As of September 2022, several variants have now spread around the world [3–5]. To date, one of the most recent of these is the Omicron variant (B.1.1.529), which arrived from South Africa [6] including sub-lineage variants such as 21L (Nextstrain classification (https://clades.nextstrain.org/; access date: 16 January 2023)) or BA.2 (Pangolin classification (https://cov-lineages.org/resources/pangolin.html; access date: 16 January 2023)), which was first highly prevalent in Denmark, then spread to across the whole of Europe [7,8]. Monitoring these variants within a community is necessary to understand their transmission dynamics and take appropriate public health measures. This surveillance is currently mainly performed using real-time reverse transcription (RT)-PCR (qPCR) or whole genome sequencing from nasopharyngeal samples from infected patients. This approach is relatively expensive, labour intensive, and time consuming, and cannot be

<footnote n="1">¹ MEPHI, Institut de Recherche pour le Développement (IRD), Aix-Marseille Université, 13005 Marseille, France</footnote>

<footnote n="2">² Institut Hospitalo-Universitaire Méditerranée Infection, 19-21 Boulevard Jean Moulin, 13005 Marseille, France</footnote>

<footnote n="3">³ Assistance Publique—Hôpitaux de Marseille (AP-HM), Aix-Marseille Université, 13005 Marseille, France</footnote>

<footnote n="4">⁴ Bataillon des Marins Pompiers de la ville de Marseille, 13005 Marseille, France</footnote>

<footnote>* Correspondence: bernard.la-scola@univ-amu.fr (B.L.S.); aherfi.s@ap-hm.fr (S.A.);
Tel.: +33-4-1373-2401 (B.L.S. & S.A.)</footnote>
exhaustively implemented in all settings such as low resource environments. In addition, this approach presents a major bias, since only the genomes of symptomatic people who have attended hospital or presented at a laboratory to be tested for SARS-CoV-2 by RT-PCR and who presented well-defined criteria were sequenced and analysed. This approach thus overlooks asymptomatic infections that may represent more than one-third of SARS-CoV-2 infections, as suggested by one recent study [9].

Several studies have reported the presence of SARS-CoV-2 RNA in stool and anal/rectal swabs, not only from symptomatic but also from asymptomatic patients [10–12]. It should be noted that even if the faecal–oral transmission of SARS-CoV2 has not been clearly demonstrated, viruses in stools can rarely still remain infectious [13,14]. Therefore, the analysis of SARS-CoV-2 in wastewater would appear to be an efficient approach to monitoring the burden of disease and infection in the community. Monitoring the SARS-CoV-2 epidemic through wastewater may also make it possible to more easily monitor particular establishments such as nursing homes and hospitals, particular districts in a city, or may even make it possible to assess the effects of sanitary measures and lockdowns that may be imposed in the context of the COVID-19 pandemic [15,16].

The detection and monitoring of SARS-CoV-2 in wastewater samples has been reported in many countries [15,17–21]. Nevertheless, few studies have performed SARS-CoV-2 genome sequencing on sewage samples in order to determine the strains circulating within a community and to study their genetic diversity [16,22–29]. Unexpectedly, some sequences were not identified in clinical samples [20,22–26,30,31]. This method could, therefore, also be a more efficient, faster, and more realistic method of monitoring the emergence and spread of variants within a community.

Wastewater surveillance of SARS-CoV-2 comprises several phases, starting with the sampling of wastewater and possibly storing wastewater before processing, followed by virus concentration, RNA extraction, SARS-CoV-2 RNA detection, and quantification (qPCR) as well as the potential sequencing of any variants present. Although most studies agree that collected wastewater must be stored at low temperatures (usually at 4 °C) to preserve the viral load [27,28], a large panel of procedures is available that has been proposed to concentrate viruses. Numerous studies have evaluated procedures concentrating on enveloped viruses from wastewater including ultracentrifugation [24,32], ultrafiltration [32–37], polyethylene glycol (PEG) and salt-based precipitation [32,35,36,38], skimmed milk flocculation [36,37], sludge extraction [36], bag-mediated filtration [36], and magnesium- and aluminium-based adsorption [34,38]. In a previous study, we used the ultracentrifugation method to concentrate SARS-CoV-2 in wastewater [24]. However, this method requires expensive equipment (an ultracentrifuge), which may not be available in a standard laboratory. Moreover, as the laboratory in charge may be far from the sampling place, it is important to ensure good preservation of samples. In this study, we first evaluated the best method for preserving wastewater samples to avoid the progressive loss of viral RNA, then tried to determine the best method of concentrating the viral RNA present in sewers using a technique that could be implemented in any laboratory.

## 2. Materials and Methods

### 2.1. Wastewater Collection

The study was performed between January 2022 and June 2022 in the city of Marseille in south-eastern France. The SERAMM (Marseille Metropole Sanitation Department) collected a daily sample of 250 mL of wastewater from the separate sewer networks (indicated “RS”) using an automatic “ASP-Station 2000 RPS20B” sampler (Endress Hauser, Huningue, France). This network, referred to as the “RS network”, drains wastewater from Marseille and 16 other municipalities including nearly all hospitals in the city as well as COVID-19 dedicated units. Marseille city has more than 870,000 inhabitants. The wastewater treatment plant capacity of Marseille is one of the most important in France and Europe, servicing 1.82 million equivalent inhabitants. This type of sampling made it possible to fill a refrigerated flask with 20 litres of wastewater per 24 h, collected between
8 am and 8 am the following day. Every day, samples were transferred on ice to the nuclear, radiological, biological, chemical laboratory (the "NRBC unit") of the Bataillon de Marins-Pompiers de Marseille (BMPM), treated within one hour of collection for multiplex qPCR, then transferred directly to the University Hospital Institute Méditerranée Infection for immediate use.

## 2.2. Preservation of Wastewater Samples

A quantity of 10 mL of wastewater was immediately filtered using 5 µm sterile syringe filters (Clearline, Merck, Darmstadt, Germany). Half of the samples were kept at room temperature while the other half were placed immediately at 4 °C. The study was performed on seven wastewater samples corresponding to seven different dates (seven time points). The quantity of SARS-CoV-2 present in the wastewater was measured on the day of sampling then every week for three weeks. The loss of virus quantity compared with the initial quantity was assessed by calculating the differences between the qPCR cycle threshold values (Delta Ct).

## 2.3. RNA Extraction and RT-qPCR

Nucleic acid extraction and qPCR were performed as previously described [24]. Briefly, the extraction of viral nucleic acids was performed using the EZ1 Virus Mini Kit (Qiagen, Hilden, Germany), following the manufacturer's recommendations using 200 µL of wastewater and eluted in 60 µL of elution buffer. Real-time RT-PCRs were carried out specifically targeting the N-gene, using the primers previously described [24]. A dilution series of the calibration range of the PCR, tested in triplicate, provided an efficiency of 1.99, a R² of 0.9987, and a SD between the replicates <0.167 for each point. The RT-PCR was carried out using the Superscript III Platinum One-step Quantitative RT-PCR system with the ROX Kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer's recommendations, with a final concentration of 400 nM of primers, of 200 nM of probe, in a final volume of 20 µL with 5 µL of RNA. The RT-PCR program was performed according to the manufacturer's recommendations, and carried out on a LightCycler 480i (Roche Diagnostics, Meylan, France).

An external control based on a SARS-CoV2 strain and calibrated at 30 Ct, with a tolerance of +/- 0.5Ct for this N-gene targeting system, allowed us to control the standardization of all the PCR runs. A negative template control was also put in each PCR run.

## 2.4. Wastewater Filtration and Concentration Method

### 2.4.1. Comparison of Ultracentrifugation and Ultrafiltration

The first steps were comparable to that set out in our previous study [24]. Briefly, if large particles were present in wastewater, a low-speed centrifugation step was conducted at 4000× g for 10 min at 4 °C. Then, 150 mL of wastewater was serially filtered through a 5 µm polycarbonate membrane filter (Merck Millipore, Burlington, MA, USA) using a vacuum filtration flask, a 0.8 µm disposable filter unit (Thermo Scientific, Waltham, MA, USA), a 0.45 µM bottle-top vacuum filtration system (VWR, Radnor, PA, USA), and finally, a 0.2 µM sterile syringe filter (ClearLine, Dutscher, Bernolsheim, France). For the next steps in comparing ultracentrifugation to ultrafiltration, 50 mL of the filtered samples was centrifuged at 100,000× g for two hours using a Sorvall Discovery 90SE ultra-centrifuge. The pellet was resuspended in 200 µL of PBS. At the same time, 50 mL of filtered samples was introduced into Centricon tubes (Centricon Plus centrifugal filter, 70 Merck Millipore, Burlington, MA, USA), then centrifuged for 30 min at 3000× g. The Centricons were turned over and centrifuged for two minutes at 1000× g. About 200 µL of the eluate was recovered. These tests were carried out in quadruplicate.

### 2.4.2. Comparison of Filtration Steps

To optimise the viral concentration of SARS-CoV-2 from wastewater, tests on the usefulness of the different filtration steps, based on a previous study [24], were carried out in order to simplify the filtration procedure. Filtration steps at 5 µm, 0.8 µm, and 0.2 µm were
performed on a 33 mm filter syringe, while 0.45 µm was performed on the 500 mL filtration units. Concentrations of SARS-CoV-2 were compared using three different protocols. We first tested the filtration at 5 µm, 0.8 µm, 0.45 µm, and 0.2 µm, plus ultrafiltration versus filtration on 5 µm, 0.8 µm and 0.2 µm, plus ultrafiltration (quadruplicate), then filtration at 5 µm, 0.8 µm, and 0.2 µm, plus ultrafiltration versus filtration at 5 µm and 0.8 µm plus ultrafiltration (15 replicates). In a final assay, we compared the filtration at 5 µm, 0.8 µm, plus ultrafiltration, but replaced the 0.8 µm syringe filter with a 0.8 µm filter unit (12 replicates).

## 2.5. Sequencing of SARS-CoV-2 in Wastewater and Sequence Analysis

### 2.5.1. RT-PCR

Reverse transcription was performed using the SuperScript VILO cDNA Synthesis Kit (11754-250, Thermo Fisher, Waltham, MA, USA), with 10 µL of RNA from each wastewater concentrated sample according to the manufacturer's protocols. The programme was 25 °C for 10 min, 42 °C for 120 min, and 85 °C for five minutes. The provided cDNA was amplified with the two primer pools (ARTIC nCoV-2019 V3 Panel and 500rxn of IDT 10006788, Integrated DNA Technologies, Inc., Coralville, IA, USA) in 16 wells (8 wells for each pool). The PCR reaction mixture was prepared with 2.5 µL of cDNA template, primer pool (1.6-µM final concentration), dNTP (0.2-mM final concentration), 10 µL of 10× PCR Buffer containing 15-mM MgCl2 (1X final concentration), HotStarTaq DNA Polymerase (Qiagen 203205, Hilden, Germany) (2.5 U per reaction), and DEPC-treated water (qsp 25 µL). The amplification programme was 95 °C for 15 min, followed by 50 cycles, each cycle consisting of denaturation at 95 °C for 30 seconds, annealing at 65 °C for 30 seconds, and extension at 72 °C for 5 min. The program included a final extension step at 72 °C for 10 min. The PCR products (total volume = 400 µL) were concentrated by Nucleofast 96P (Macherey Nagel ref 743100.50, Hoerdt, France) following the manufacturer's recommendations for a total elution volume of 30 µL for each pool and were separated using electrophoresis in 2% agar gel and visualised with SYBR safe (Invitrogen S3102, Waltham, MA, USA). The Monarch® DNA Gel Extraction Kit (New England BioLabs, Ipswich, MA, USA; ref T1020L, Evry-Courcouronnes, France) was used for the purification of an interest band at 400 pb with an elution at 40 µL for each pool.

### 2.5.2. Sequencing

The libraries were prepared using the Illumina COVIDSeq protocol (Illumina Inc., San Diego, CA, USA). Each previously purified PCR product was processed for tagmentation and adapter ligation using IDT for the Illumina Nextera UD Indices. Enrichment and clean-up were performed as per the protocols provided by the manufacturer (Illumina Inc., San Diego, CA, USA). All libraries were pooled together, and the pooled samples were quantified using a Qubit 2.0 fluorometer (Invitrogen Inc., MA, USA). Fragment sizes were analysed in an Agilent Fragment analyser 5200 (Agilent Inc., Santa Clara, CA, USA). Sequencing was performed on the Illumina MiSeq instrument (Illumina Inc., San Diego, CA, USA). Negative template controls were performed for each sequencing run.

### 2.5.3. Bioinformatics

The reads generated from the wastewater samples were analysed as previously described [24]. Briefly, reads from pool1 and pool2 provided by the ARTIC procedure were mapped together against the Wuhan-Hu-1 SARS-CoV-2 isolate genome (GenBank accession number NC_045512.2) using the CLC genomics software v7.5 (Qiagen Digital Insights, Hilden, Germany) with the default parameters. To consider a position, a minimum coverage of 10 reads was considered. The non-synonymous mutations present in more than 10% of the reads were considered. A custom Python script was used to associate mutations present in the output file generated by the CLC genomics workbench with a database of mutations of SARS-CoV-2 variants built by our institute based on the Nextstrain tool [39]. The script generated the mutational events detected, their frequency, the numbers of reads covering a given nucleotide position, the frequency of mutational events for each variant,
and the list of variants associated with the given mutation. Then, these results were used to generate graphics with the plotly packages [40,41].

Proportions of reads harbouring mutations specific to the Omicron 21K/BA.1 or BA.1.1 lineage on one hand, and mutations specific to 21L/BA.2 on the other hand, were more deeply observed.

### 2.5.4. Clinical Samples

Data obtained here from the sewage samples were compared to those obtained from the patients' respiratory samples and produced in the framework of the EMERGEN consortium led by the French Ministry of Solidarity and Health and Santé Publique France (https://www.santepubliquefrance.fr/dossiers/coronavirus-covid-19/consortium-emergen; access date 16 January 2023). Briefly, these data were generated from SARS-CoV-2-diagnosed patients through screening using the TaqPath COVID-19 Kit (Thermo Fisher Scientific, Waltham, MA, USA), which targets viral genes ORF1, N (nucleocapsid), and S (spike), and using qPCR assays specific to variants that target spike mutations L452R and K417N. Negativity for L452R, and positivity for K417N with positive signals by the TaqPath COVID-19 Kit for all three genes targeted indicated a 21L/BA.2 variant; conversely, a 21K/BA.1 was indicated by the negativity of S gene detection by the TaqPath COVID-19 Kit. When performed, definitive identification of the variant was reached through analysis of the viral genome obtained by next-generation sequencing with the Illumina COVID-Seq protocol on the NovaSeq 6000 instrument (Illumina Inc.), followed by sequence read processing and genome analysis, as previously described [42].

### 2.6. Statistical Tests

In order to compare the two groups of samples with regard to the conservation of wastewater samples (4 °C and room temperature), a paired non-parametric Student T test was carried out (Wilcoxon matched pairs signed rank test). In order to compare the effects of the filtration steps, paired parametric Student T tests were carried out. A verification of the normal distribution was carried out before further analyses in order to identify the Gaussian distribution of the samples (using a Shapiro-Wilk normality test or an Agostino and Pearson normality test, depending on the number of samples). Sample groups were considered statistically different if the p-value was less than 0.05. Statistical tests were performed using GraphPad Prism 7.00 software (https://www.graphpad.com/scientific-software/prism/; access date: 16 January 2023).

## 3. Results

### 3.1. Conservation of Wastewater

The results of the quantity of virus present in each of the two conditions (4 °C and room temperature) are presented in Figure 1 to compare the effect of storage temperature.

On days 7 and 14, viral RNA was significantly better preserved at 4 °C than at room temperature (p < 0.05 at 7 and 14 days of conservation). Moreover, the viral load decreased drastically after 21 days of storage, even when the samples were stored at 4 °C.

### 3.2. Comparison of Concentration Methods

Two different concentration methods of SARS-CoV-2 were tested from wastewater: ultracentrifugation and ultrafiltration with a Centricon filter and centrifuge. The comparison of these two concentration methods is shown in Figure 2, which represents the variation in the amounts of SARS-CoV-2 in wastewater after the two concentration methods had been tested.

The method using ultrafiltration with Centricons was more efficient, with an average Delta Ct gain from the native samples of 3.24 compared to the ultracentrifugation method which provided an average Delta Ct of only 0.34. The difference between the two methods was statistically significant (p = 0.0150). Thus, concentration by ultrafiltration with Centricons was used for the subsequent tests.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Effect of temperature storage on SARS-CoV-2. NS, not significant; RT, room temperature. \* indicate significant differences (p < 0.05).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Variation in the amount of SARS-CoV-2 virus in wastewater following two concentration methods; mean Ct ranged from 38 (initial concentration) to 34.8 (ultrafiltration). \* indicate significant differences (p < 0.05).</figcaption></figure>

### 3.3. *Comparison of Filtration Steps*

We first compared the viral load obtained after filtrations using 5 µm then 0.8 µm filters, followed or not followed by a filtration step on a 0.45 µm filter, and then on a 0.2 µm filter. The results are shown in Figure 3a.

No statistically significant difference between the 0.45 µm filtered and unfiltered wastewater samples was observed (p = 0.5559). The filtration step at 0.45 µm was therefore not useful.

We then compared the viral load obtained after filtrations using 5 µm and 0.8 µm filters followed or not followed by a final step with a 0.2 µm filter. The results are shown in Figure 3b. No statistically significant difference was observed between the samples of wastewater filtered on 0.8 µm and 0.2 µm, and those filtered on 0.8 µm only (p = 0.3695). The filtration step on 0.2 µm was therefore unnecessary.

Finally, we compared the viral load obtained after filtration on a 5 µm filter, then either on a 0.8 µm syringe filter or on a sterile 0.8 µm filtration unit. The results are shown in Figure 3c. The viral load obtained using filtration units was significantly higher than those obtained with a syringe filter (difference of 3 Ct; p = 0.0147).

Hence, the only useful filtration that increased the virus concentration was the filtration step on a sterile 0.8 µm filtration unit.
<figure><img src="image_6.png" /></figure>

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

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

**Figure 3.** Variation in the SARS-CoV-2 RNA load in wastewater after different filtration steps followed by ultrafiltration: (a) comparison of the viral load obtained after filtration on a 0.45 µm filter of wastewater or for unfiltered water followed by ultrafiltration; the mean Ct ranged from 38.1 (initial concentration) to 31.9 (ultrafiltration only); (b) comparison of viral load obtained after the filtration of wastewater using a 0.8 µm filter and a 0.2 µm filter, and a 0.8 µm filter only, followed by ultrafiltration, where the mean Ct ranged from 36.2 (initial concentration) to 31.3 (0.8 µm + 0.2 µm filtration + ultrafiltration); (c) comparison of the viral load obtained after filtration using a 0.8 µm syringe filter and a 0.8 µm wastewater filtration unit, followed by ultrafiltration, where the mean Ct ranged from 37.6 (initial concentration) to 33.3 (0.8 µm filtration unit + ultrafiltration). \* indicate significant differences (p < 0.05). NS indicates no significant difference (p > 0.05).
### 3.4. Sequencing SARS-CoV-2 from Sewers

For each concentrated wastewater sample, a total of 80 µL (40 µL for each pool) of purified DNA extract was sequenced over a three-month period from February 2022 to April 2022. A total of 13 samples was collected on 20 and 22 February, 5, 6, 12, 13, 19, 26 and 27 March, and 2, 3, 9, and 10 April. Mapping on the reference genome of NGS reads provided good SARS-CoV-2 genome coverage, ranging from 87.13% to 99.44% of the genome length. Mean depth of coverage ranged from 28 (for sample collected on 6 March) to 275 (for sample collected on 9 April). The Omicron 21K/BA.1 and 21L/BA.2 variants were clearly predominant throughout the study period. When considering the mutation patterns of 21L/BA.2 that were absent in 21K/BA.1 and BA.1.1, they represented 36–50% of the total number of reads on 20 February, reaching 95–100% on 10 April (Figure 4).

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Mean proportions of specific mutations harboured by the Omicron 21K/BA.1 or BA.1.1 lineages and by the Omicron 21L/BA.2 lineage. Mutations A2832G; A11537G; G22898A; G23048A; C23202A; C24130A; and C24503T were found in 21K/BA.1 and in BA.1.1. but absent in 21L/BA.2. Mutations T670G; C2790T; G4184A; C4321T; C9344T; A9424G; C9534T; C9866T; C17410T; C19955T; and C26060T were found in 21L/BA.2 and absent in 21K/BA.1 and BA.1.1.</figcaption></figure>

### 3.5. Sequencing SARS-CoV2 in Clinical Samples

We retrieved the SARS-CoV-2 sequencing data from 9,158 clinical samples collected between 14 February and 10 April 2022 and produced in the framework of the EMER-GEN consortium (https://www.santepubliquefrance.fr/dossiers/coronavirus-covid-19/consortium-emergen, access date: 16 January 2023). On week 8, 21K/BA.1 and 21L/BA.2 were respectively present in 77.0% and 21.7% of the patients, reaching 3.5% and 96.2% on week 15, respectively. Other genotypes identified included the 21I, 21J, and 21M lineages, but in negligible proportions. These results were correlated with the variant identifications based on SARS-CoV-2 sequencing from the wastewater samples (Figure 5).
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Proportions of reads identified from the 21K/BA.1 (red) and 21L/BA.2 (blue) variants in the patients (histograms) and wastewater (curves) samples. The curves represent the proportion of reads harbouring the C23202A (red) or the G4184A (blue) mutations specific to the 21K/BA.1 and 21L/BA.2 variants, respectively.</figcaption></figure>

**4. Discussion**

Wastewater surveillance of SARS-CoV-2 involves a succession of steps, starting with wastewater collection followed by transport, storage, +/- virus concentration methods, RNA extraction, and SARS-CoV-2 detection/quantification/sequencing. All of these steps can influence the reliability of the results, as previously described [32,43,44]. A major part of the improvements into research on wastewater has focused on extraction and quantification methods [38,45]. However, a better understanding of sample storage and concentration methods will help to improve the detection and sequencing of SARS-CoV-2.

Previous studies have shown that the storage of wastewater samples at 4 °C had no significant effect, while storage at -20 °C or -80 °C led to a significant reduction in the copy numbers of SARS-CoV-2 RNA [46,47]. First, in the present study, we evaluated the effect of the storage of wastewater samples at 4 °C and at room temperature. We showed that wastewater samples should be stored at 4 °C for short-term preservation (<14 days). A refrigerator is therefore the only equipment needed to store wastewater and can be installed in any standard laboratory.

Second, wastewater samples must be sufficiently concentrated to allow for the effective detection of very few viruses in a large volume of water. Several studies have evaluated the best method of concentrating enveloped viruses in wastewater [26,35,36,38]. For SARS-CoV-2, a preliminary study carried out in our laboratory showed that ultracentrifugation was the best way to concentrate the virus [24]. However, this method requires expensive equipment (an ultracentrifuge) that is available in research laboratories but almost never in laboratories performing routine clinical biology or wastewater analyses. In this study, the ultrafiltration concentration method was able to concentrate SARS-CoV-2 RNA more efficiently than ultracentrifugation. The individual tests for each filtration step showed that filtration at 0.8 µm using a sterile filtration unit after filtration at 5 µm was the only procedure that had a real impact on the viral concentration. This simple procedure for concentrating and filtering wastewater samples can be carried out in any laboratory, since a simple centrifuge is all that is needed. A photographic summary of the different stages is presented in Figure 6.
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Photographic summary of the definitive protocol for the preparation of wastewater samples.</figcaption></figure>

Results from the whole genome sequencing of the wastewater samples enabled the detection and discrimination of the specific mutations of the two predominant variants, 21K/BA.1 and 21L/BA.2, which were circulating during this period. Analysis of the proportions of the reads harbouring the specific mutations of these two variants in the sequential wastewater samples made it possible to track the quasi extinction of the 21K/BA.1 variant and the rapid emergence of 21L/BA.2 over time. These simultaneous evolutions, with the curves crossing at week 8 of 21K/BA.1 and 21L/BA.2, were superimposable and correlated with those found based on the analysis of respiratory samples from patients. These results should probably be weighted according to the transmissibility of the viral variants. Indeed, the SARS-CoV-2 load in the stool collected from a patient infected by a highly transmissible
variant will be higher than that for a patient infected by a less transmissible variant. For the 21K/BA.1 and 21L/BA.2 variants, it has been shown that the mean Ct obtained from 21L/BA.2-infected patients was lower than those obtained from the 21K/BA.2-infected patients [48,49]. It has been previously demonstrated that rainfall may impact the SARS-CoV2 concentration in wastewater [50]. It is noteworthy that the RS network of Marseille collects only domestic wastewater; storm water being collected in a different and adapted network. Moreover, average rainfall amounts in Marseille range from 500 to 700 mm per year. Between 20 February and 10 April 2022, the rainfall amount was low, reaching only 6.4 mm distributed on 3 days (13, 20, and 30 March).

In addition to SARS-CoV-2, the proposed protocol could make it possible to simply process wastewater samples in developing countries or during occasional studies in various areas of the world without the need for very complex equipment. Similarly, this protocol should probably be tested for monitoring other respiratory viruses, gastroenteritis, or even arboviruses when they are excreted in urine and/or stool.

**Author Contributions:** Conceptualisation, B.L.S.; Methodology, B.L.S., A.L., N.W., and S.A.; Software, L.S. and S.A.; Validation, B.L.S.; Formal analysis, M.B., L.S., N.W., G.P., and S.A.; Investigation, A.L.; Resources, A.L.; Data curation, N.W., M.B., L.S., and S.A.; Writing—original draft preparation, N.W., M.B., L.S., and S.A.; Writing—review and editing, P.C., B.L.S. and S.A.; Visualisation, M.B., L.S., and S.A.; Supervision, B.L.S. and S.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the French Government under the "Investments for the Future" programme managed by the National Agency for Research (ANR), Méditerranée-Infection 10-IAHU-03.

**Institutional Review Board Statement:** The present study was approved by the ethics committee of the University Hospital Institute (IHU) Méditerranée Infection (No. 2022-001).

**Informed Consent Statement:** Not applicable. This study was conducted on sequences obtained from nasopharyngeal samples sent to our clinical microbiology/virology laboratory for the purpose of diagnosis of SARS-CoV-2 infection and from which the SARS-CoV-2 genomes had been obtained as recommended in France (https://www.santepubliquefrance.fr/dossiers/coronavirus-covid-19/consortium-emergen; access date: 16 January 2023). Data were anonymised. No additional clinical samples were specifically collected for this study.

**Data Availability Statement:** Not applicable.

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

## References

1. Altmann, D.M.; Boyton, R.J.; Beale, R. Immunity to SARS-CoV-2 Variants of Concern. *Science* **2021**, *371*, 1103–1104. [CrossRef] [PubMed]

2. Challen, R.; Brooks-Pollock, E.; Read, J.M.; Dyson, L.; Tsaneva-Atanasova, K.; Danon, L. Risk of Mortality in Patients Infected with SARS-CoV-2 Variant of Concern 202012/1: Matched Cohort Study. *BMJ* **2021**, *372*, n579. [CrossRef]

3. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D.; et al. Estimated Transmissibility and Impact of SARS-CoV-2 Lineage B.1.1.7 in England. *Science* **2021**, *372*, eabg3055. [CrossRef] [PubMed]

4. Faria, N.R.; Mellan, T.A.; Whittaker, C.; Claro, I.M.; Candido, D.d.S.; Mishra, S.; Crispim, M.A.E.; Sales, F.C.S.; Hawryluk, I.; McCrone, J.T.; et al. Genomics and Epidemiology of the P.1 SARS-CoV-2 Lineage in Manaus, Brazil. *Science* **2021**, *372*, 815–821. [CrossRef] [PubMed]

5. Volz, E.; Mishra, S.; Chand, M.; Barrett, J.C.; Johnson, R.; Geidelberg, L.; Hinsley, W.R.; Laydon, D.J.; Dabrera, G.; O'Toole, Á.; et al. Assessing Transmissibility of SARS-CoV-2 Lineage B.1.1.7 in England. *Nature* **2021**, *593*, 266–269. [CrossRef] [PubMed]

6. Kannan, S.; Shaik Syed Ali, P.; Sheeza, A. Omicron (B.1.1.529)-Variant of Concern-Molecular Profile and Epidemiology: A Mini Review. *Eur. Rev. Med. Pharmacol. Sci.* **2021**, *25*, 8019–8022. [CrossRef]

7. Mahase, E. Omicron Sub-Lineage BA.2 May Have "Substantial Growth Advantage," UKHSA Reports. *BMJ* **2022**, *376*, o263. [CrossRef]

8. Desingu, P.A.; Nagarajan, K. Omicron BA.2 Lineage Spreads in Clusters and Is Concentrated in Denmark. *J. Med. Virol.* **2022**, *94*, 2360–2634. [CrossRef]
9. Oran, D.P.; Eric, J.T. The Proportion of SARS-CoV-2 Infections That Are Asymptomatic: A Systematic Review. *Ann. Intern. Med.* **2021**, *174*, 655–662. [CrossRef]

10. Gu, J.; Han, B.; Wang, J. COVID-19: Gastrointestinal Manifestations and Potential Fecal-Oral Transmission. *Gastroenterology* **2020**, *158*, 1518–1519. [CrossRef]

11. Xiao, F.; Tang, M.; Zheng, X.; Liu, Y.; Li, X.; Shan, H. Evidence for Gastrointestinal Infection of SARS-CoV-2. *Gastroenterology* **2020**, *158*, 1831–1833.e3. [CrossRef]

12. An, P.; Song, P.; Wang, Y.; Liu, B. Asymptomatic Patients with Novel Coronavirus Disease (COVID-19). *Balk. Med. J.* **2020**, *37*, 229–230. [CrossRef]

13. Dergham, J.; Delerce, J.; Bedotto, M.; La Scola, B.; Moal, V. Isolation of Viable SARS-CoV-2 Virus from Feces of an Immunocompromised Patient Suggesting a Possible Fecal Mode of Transmission. *J. Clin. Med.* **2021**, *10*, 2696. [CrossRef]

14. Guo, M.; Tao, W.; Flavell, R.A.; Zhu, S. Potential Intestinal Infection and Faecal–Oral Transmission of SARS-CoV-2. *Nat. Rev. Gastroenterol. Hepatol.* **2021**, *18*, 269–283. [CrossRef]

15. Wurtz, N.; Lacoste, A.; Jardot, P.; Delache, A.; Fontaine, X.; Verlande, M.; Annessi, A.; Giraud-Gatineau, A.; Chaudet, H.; Fournier, P.-E.; et al. Viral RNA in City Wastewater as a Key Indicator of COVID-19 Recrudescence and Containment Measures Effectiveness. *Front. Microbiol.* **2021**, *12*, 664477. [CrossRef] [PubMed]

16. Nghiem, L.D.; Morgan, B.; Donner, E.; Short, M.D. The COVID-19 Pandemic: Considerations for the Waste and Wastewater Services Sector. *Case Stud. Chem. Environ. Eng.* **2020**, *1*, 100006. [CrossRef]

17. Wurtzer, S.; Marechal, V.; Mouchel, J.M.; Maday, Y.; Teyssou, R.; Richard, E.; Almayrac, J.L.; Moulin, L. Evaluation of Lockdown Effect on SARS-CoV-2 Dynamics through Viral Genome Quantification in Waste Water, Greater Paris, France, 5 March to 23 April 2020. *Eurosurveillance* **2020**, *25*, 2000776. [CrossRef] [PubMed]

18. Haramoto, E.; Malla, B.; Thakali, O.; Kitajima, M. First Environmental Surveillance for the Presence of SARS-CoV-2 RNA in Wastewater and River Water in Japan. *Sci. Total Environ.* **2020**, *737*, 140405. [CrossRef] [PubMed]

19. Randazzo, W.; Truchado, P.; Cuevas-Ferrando, E.; Simón, P.; Allende, A.; Sánchez, G. SARS-CoV-2 RNA in Wastewater Anticipated COVID-19 Occurrence in a Low Prevalence Area. *Water Res.* **2020**, *181*, 115942. [CrossRef]

20. Izquierdo-Lara, R.; Elsinga, G.; Heijnen, L.; Munnink, B.B.O.; Schapendonk, C.M.E.; Nieuwenhuijse, D.; Kon, M.; Lu, L.; Aarestrup, F.M.; Lycett, S.; et al. Monitoring SARS-CoV-2 Circulation and Diversity through Community Wastewater Sequencing, the Netherlands and Belgium. *Emerg. Infect. Dis.* **2021**, *27*, 1405–1415. [CrossRef]

21. La Rosa, G.; Iaconelli, M.; Mancini, P.; Bonanno Ferraro, G.; Veneri, C.; Bonadonna, L.; Lucentini, L.; Suffredini, E. First Detection of SARS-CoV-2 in Untreated Wastewaters in Italy. *Sci. Total Environ.* **2020**, *736*, 139652. [CrossRef]

22. Crits-Christoph, A.; Kantor, R.S.; Olm, M.R.; Whitney, O.N.; Al-Shayeb, B.; Lou, Y.C.; Flamholz, A.; Kennedy, L.C.; Greenwald, H.; Hinkle, A.; et al. Genome Sequencing of Sewage Detects Regionally Prevalent SARS-CoV-2 Variants. *mBio* **2021**, *12*, e02720–e02720. [CrossRef] [PubMed]

23. Nemudryi, A.; Nemudraia, A.; Wiegand, T.; Surya, K.; Buyukyoruk, M.; Cicha, C.; Vanderwood, K.K.; Wilkinson, R.; Wiedenheft, B. Temporal Detection and Phylogenetic Assessment of SARS-CoV-2 in Municipal Wastewater. *Cell Rep. Med.* **2020**, *1*, 100098. [CrossRef]

24. Wurtz, N.; Revol, O.; Jardot, P.; Giraud-Gatineau, A.; Houhamdi, L.; Soumagnac, C.; Annessi, A.; Lacoste, A.; Colson, P.; Aherfi, S.; et al. Monitoring the Circulation of SARS-CoV-2 Variants by Genomic Analysis of Wastewater in Marseille, South-East France. *Pathogens* **2021**, *10*, 1042. [CrossRef] [PubMed]

25. Rios, G.; Lacoux, C.; Leclercq, V.; Diamant, A.; Lebrigand, K.; Lazuka, A.; Soyeux, E.; Lacroix, S.; Fassy, J.; Couesnon, A.; et al. Monitoring SARS-CoV-2 Variants Alterations in Nice Neighborhoods by Wastewater Nanopore Sequencing. *Lancet Reg. Health Eur.* **2021**, *10*, 100202. [CrossRef]

26. Wurtzer, S.; Waldman, P.; Levert, M.; Cluzel, N.; Almayrac, J.L.; Charpentier, C.; Masnada, S.; Gillon-Ritz, M.; Mouchel, J.M.; Maday, Y.; et al. SARS-CoV-2 Genome Quantification in Wastewaters at Regional and City Scale Allows Precise Monitoring of the Whole Outbreaks Dynamics and Variants Spreading in the Population. *Sci. Total Environ.* **2022**, *810*, 152213. [CrossRef] [PubMed]

27. Cutrupi, F.; Cadonna, M.; Manara, S.; Foladori, P. Surveillance of SARS-CoV-2 in Extensive Monitoring of Municipal Wastewater: Key Issues to Yield Reliable Results. *Water Sci. Technol.* **2021**, *84*, 3508–3514. [CrossRef] [PubMed]

28. Medema, G.; Heijnen, L.; Elsinga, G.; Italiaander, R.; Brouwer, A. Presence of SARS-Coronavirus-2 RNA in Sewage and Correlation with Reported COVID-19 Prevalence in the Early Stage of the Epidemic in The Netherlands. *Environ. Sci. Technol. Lett.* **2020**, *7*, 511–516. [CrossRef]

29. McGowan, J.; Borucki, M.; Omairi, H.; Varghese, M.; Vellani, S.; Chakravarty, S.; Fan, S.; Chattopadhyay, S.; Siddiquee, M.; Thissen, J.B.; et al. SARS-CoV-2 Monitoring in Wastewater Reveals Novel Variants and Biomarkers of Infection. *Viruses* **2022**, *14*, 2032. [CrossRef]

30. Agrawal, S.; Orschler, L.; Lackner, S. Metatranscriptomic Analysis Reveals SARS-CoV-2 Mutations in Wastewater of the Frankfurt Metropolitan Area in Southern Germany. *Microbiol. Resour. Announc.* **2021**, *10*, e00280-21. [CrossRef]

31. La Rosa, G.; Mancini, P.; Bonanno Ferraro, G.; Veneri, C.; Iaconelli, M.; Lucentini, L.; Bonadonna, L.; Brusaferro, S.; Brandtner, D.; Fasanella, A.; et al. Rapid Screening for SARS-CoV-2 Variants of Concern in Clinical and Environmental Samples Using Nested RT-PCR Assays Targeting Key Mutations of the Spike Protein. *Water Res.* **2021**, *197*, 117104. [CrossRef]

32. Ahmed, W.; Bertsch, P.M.; Bivins, A.; Bibby, K.; Farkas, K.; Gathercole, A.; Haramoto, E.; Gyawali, P.; Korajkic, A.; McMinn, B.R.; et al. Comparison of Virus Concentration Methods for the RT-QPCR-Based Recovery of Murine Hepatitis Virus, a Surrogate for SARS-CoV-2 from Untreated Wastewater. *Sci. Total Environ.* **2020**, *739*, 139960. [CrossRef] [PubMed]
33. Forés, E.; Bofill-Mas, S.; Itarte, M.; Martínez-Puchol, S.; Hundesa, A.; Calvo, M.; Borrego, C.M.; Corominas, L.L.; Girones, R.; Rusiñol, M. Evaluation of Two Rapid Ultrafiltration-Based Methods for SARS-CoV-2 Concentration from Wastewater. *Sci. Total Environ.* **2021**, *768*, 144786. [CrossRef] [PubMed]

34. Jafferali, M.H.; Khatami, K.; Atasoy, M.; Birgersson, M.; Williams, C.; Cetecioglu, Z. Benchmarking Virus Concentration Methods for Quantification of SARS-CoV-2 in Raw Wastewater. *Sci. Total Environ.* **2021**, *755*, 142939. [CrossRef]

35. LaTurner, Z.W.; Zong, D.M.; Kalvapalle, P.; Gamas, K.R.; Terwilliger, A.; Crosby, T.; Ali, P.; Avadhanula, V.; Santos, H.H.; Weesner, K.; et al. Evaluating Recovery, Cost, and Throughput of Different Concentration Methods for SARS-CoV-2 Wastewater-Based Epidemiology. *Water Res.* **2021**, *197*, 117043. [CrossRef]

36. Philo, S.E.; Keim, E.K.; Swanstrom, R.; Ong, A.Q.W.; Burnor, E.A.; Kossik, A.L.; Harrison, J.C.; Demeke, B.A.; Zhou, N.A.; Beck, N.K.; et al. A Comparison of SARS-CoV-2 Wastewater Concentration Methods for Environmental Surveillance. *Sci. Total Environ.* **2021**, *760*, 144215. [CrossRef]

37. Rusiñol, M.; Martínez-Puchol, S.; Forés, E.; Itarte, M.; Girones, R.; Bofill-Mas, S. Concentration Methods for the Quantification of Coronavirus and Other Potentially Pandemic Enveloped Virus from Wastewater. *Curr. Opin. Environ. Sci. Health* **2020**, *17*, 21–28. [CrossRef]

38. Pérez-Cataluña, A.; Cuevas-Ferrando, E.; Randazzo, W.; Falcó, I.; Allende, A.; Sánchez, G. Comparing Analytical Methods to Detect SARS-CoV-2 in Wastewater. *Sci. Total Environ.* **2021**, *758*, 143870. [CrossRef]

39. Aksamentov, I.; Roemer, C.; Hodcroft, E.B.; Neher, R.A. Nextclade: Clade Assignment, Mutation Calling and Quality Control for Viral Genomes. *J. Open Source Softw.* **2021**, *6*, 3773. [CrossRef]

40. Sievert, C. *Interactive Web-Based Data Visualization with R, Plotly, and Shiny*; CRC Press: Boca Raton, FL, USA, 2020; ISBN 978-0-429-44727-3.

41. Hossain, S. Visualization of Bioinformatics Data with Dash Bio. In Proceedings of the 18th Python in Science Conference, SCIPY, Austin, TX, USA, 8–14 July 2019.

42. Colson, P.; Fournier, P.-E.; Chaudet, H.; Delerce, J.; Giraud-Gatineau, A.; Houhamdi, L.; Andrieu, C.; Brechard, L.; Bedotto, M.; Prudent, E.; et al. Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics. *Front. Microbiol.* **2022**, *12*, 4202. [CrossRef]

43. Pecson, B.M.; Darby, E.; Haas, C.N.; Amha, Y.M.; Bartolo, M.; Danielson, R.; Dearborn, Y.; Di Giovanni, G.; Ferguson, C.; Fevig, S.; et al. Reproducibility and Sensitivity of 36 Methods to Quantify the SARS-CoV-2 Genetic Signal in Raw Wastewater: Findings from an Interlaboratory Methods Evaluation in the U.S. *Environ. Sci.* **2021**, *7*, 504–520. [CrossRef]

44. Medema, G.; Been, F.; Heijnen, L.; Petterson, S. Implementation of Environmental Surveillance for SARS-CoV-2 Virus to Support Public Health Decisions: Opportunities and Challenges. *Curr. Opin. Environ. Sci. Health* **2020**, *17*, 49–71. [CrossRef] [PubMed]

45. Ciesielski, M.; Blackwood, D.; Clerkin, T.; Gonzalez, R.; Thompson, H.; Larson, A.; Noble, R. Assessing Sensitivity and Reproducibility of RT-DdPCR and RT-QPCR for the Quantification of SARS-CoV-2 in Wastewater. *J. Virol. Methods* **2021**, *297*, 114230. [CrossRef]

46. Islam, G.; Gedge, A.; Lara-Jacobo, L.; Kirkwood, A.; Simmons, D.; Desaulniers, J.-P. Pasteurization, Storage Conditions and Viral Concentration Methods Influence RT-QPCR Detection of SARS-CoV-2 RNA in Wastewater. *Sci. Total Environ.* **2022**, *821*, 153228. [CrossRef] [PubMed]

47. Markt, R.; Mayr, M.; Peer, E.; Wagner, A.O.; Lackner, N.; Insam, H. Detection and Stability of SARS-CoV-2 Fragments in Wastewater: Impact of Storage Temperature. *Pathogens* **2021**, *10*, 1215. [CrossRef] [PubMed]

48. Mastrorosa, I.; Cozzi-Lepri, A.; Colavita, F.; Lalle, E.; Mazzotta, V.; Cimaglia, C.; Paulicelli, J.; Matusali, G.; Fabeni, L.; Carletti, F.; et al. SARS-CoV-2 Nasopharyngeal Viral Load in Individuals Infected with BA.2, Compared to Alpha, Gamma, Delta and BA.1 Variants: A Single-Center Comparative Analysis. *J. Clin. Virol.* **2022**, *157*, 105299. [CrossRef] [PubMed]

49. Migueres, M.; Dimeglio, C.; Mansuy, J.M.; Abravanel, F.; Raymond, S.; Latour, J.; Jeanne, N.; Ranger, N.; Lhomme, S.; Saune, K.; et al. Influence of Nasopharyngeal Viral Load on the Spread of the Omicron BA.2 Variant. *Clin. Infect. Dis.* **2022**, ciac563. [CrossRef] [PubMed]

50. Janssens, R.; Hanoteaux, S.; Maloux, H.; Klamer, S.; Laisnez, V.; Verhaegen, B.; Linard, C.; Lahousse, L.; Delputte, P.; Terwagne, M.; et al. SARS-CoV-2 Surveillance in Belgian Wastewaters. *Viruses* **2022**, *14*, 1950. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.