High Throughput Determination of Infectious Virus Titers by Kinetic Measurement of infection-Induced Changes in Cell Morphology

1. Introduction

Engineered viruses have a great potential for therapeutic applications that is reflected by the increasing number in clinical trials and the approval of virus-based therapeutics in various clinical indications including cancer, infectious diseases, and gene therapy [1]. The high structural and biological complexity of these advanced therapeutic medicinal products (ATMPs) represents a challenge for pharmaceutical development, production, and quality control. In most cases, the infectivity of virus-based therapeutics is of crucial importance, as it is the critical quality attribute that reflects viral potency and often represents the parameter defining the dose administered to the patient. A fast, precise, and accurate method to measure infectious virus titers is therefore one of the key requirements to successfully develop and manufacture such virus-based drugs in a cost-effective manner and finally to ensure patients’ safety. Physico-chemical analytical methods such as nanoparticle tracking analysis (NTA), microfluidic resistive pulse sensing (MRPS), enzyme-linked immunosorbent assays (ELISA) or PCR-based methods are usually fast and easy to apply and can be used to determine concentrations of viral particles or genomic titers, respectively [2]. These types of assays, however, cannot measure viral infectivity as they do not reflect viral functionality and fitness in the complex interplay that allows the virus to enter a target cell, replicate, and spread its progeny. Usage of cell-based assays is therefore the gold standard to determine infectious virus titers.

The plaque assay and the TCID50 assay represent two frequently used examples of such traditional cell-based infectivity assays. In the plaque assay, adherent cells are infected with different concentrations of the virus sample, before adding a viscous overlay medium, which limits virus spread to directly neighboring cells. This results in localized infection clusters, originating from a single infection event. These so-called plaques, which are identified by the presence of cytopathic effects (CPE) or by virus-specific staining, are then counted to determine the infectious sample titer. Although the TCID50 assay is performed on adherent target cells, as well, it follows a different assay principle. The TCID50 is an end-point dilution assay in which serial dilutions of the virus sample are added to the culture wells in multiple replicates. The term endpoint dilution refers to the point within the serial sample dilution at which only a single (or few) or no infectious particles are left. After sufficient time, during which viral replication and spread through the well occurs, the presence or absence of virus is evaluated, often by microscopic assessment of infection-induced CPE. This reduced binary information is then used in the statistical models of Spearman/Kärber or Reed/Münch [3,4,5] to determine the TCID50 (tissue culture infectious dose 50%), describing the sample dilution factor, at which half of the inoculated wells are infected.

The benefits of these traditional infectivity assays are that they are easy to perform, do not require specialized technological equipment and can be adapted to a large variety of virus-cell combinations. Low precision is, however, a common issue and depending on the tested virus and readout method, variability of ±0.5 log₁₀ is expected [6,7]. The low precision can often be compensated by combining results from multiple sample dilutions or independent assay repetitions, although these increase assay hands-on time. Furthermore, long incubation times required for observing the endpoint make the TCID50 assay low throughput. Despite the described limitations, the traditional infectivity assays have been used for decades. However, due to the growing complexity and sample size associated with development of virus-based advanced therapy medicinal products (ATMPs), there is an increasing need for improved concepts for infectivity assays.

During the execution of traditional infectivity assays, a large part of the hands-on time is spent on the process of plaque counting or microscopic evaluation of CPE. This is not only a time-consuming procedure, but also a source of variability due to subjective operator-dependent judgements. We have previously shown that implementation of benchtop pipetting robots in combination with an automated procedure for acquisition and analysis of label-free bright field images results in a highly standardized and ~4x faster TCID50 assay [8]. Similarly, it was reported that image-based automated counting of plaques increases the robustness and speed of plaque assays [9]. Furthermore, novel technologies offer different concepts for identification and characterization of infected cells. Laser force cytology, for example, uses the principle of optical tweezers combined with microfluidics to analyze changes of biophysical properties in infected cells. This allows one to correlate different infection metrics with the applied virus dose and to calculate sample titers relative to a calibration curve [10,11]. Dodkins and colleagues applied convolutional neural networks (CNNs), in a machine learning approach, to quantify cellular phenotypes specific to viral infection in label-free bright-field images. This AI-based image processing allowed precise determination of infectious virus titers for a variety of viruses from different families. At the same time, using subtle early cellular changes as a readout for infected cells significantly shortened the total assay duration, which represents a significant bottleneck in traditional infectivity assays [12]. The need for measurements at a single, late assay endpoint is also alleviated by the xCELLigence Real-Time Cell Analysis (RTCA) system. Kinetically recorded measurements of impedance, which corelates with the integrity of a cell monolayer in multi-well plates with embedded micro-electronic biosensors, enable tracking of virus-induced cell lysis over time. It has been shown that this time-dependent information can be used to quantify infectious Hepatitis A virus (HAV), while reducing the overall assay duration by more than 50 % compared to traditional infectivity assays [13].

Our goal was to combine the respective advantages of these previous advances in a fast infectivity assay with high precision and sample throughput, while avoiding weaknesses, such as complex data analysis or dependence on proprietary technology provided by a single manufacturer. For our assay development activities, we used VSV-GP, a virus from our own oncolytic therapy platform [14,15,16]. We show that infection of BHK-21 cells with VSV-GP induces cell rounding in a time- and dose-dependent manner. The proportion of rounded cells in kinetically acquired bright-field images can be used to determine infectious virus titers relative to a reference standard curve. The resulting assay yields high precision (CV < 20 %), shortens the time to results to only 24 h, and increases the sample throughput ~10-fold compared to traditional TCID50 assays. Importantly, the principle of the kinetic infectious virus titer (KIT) assay can be applied to different virus-cell combinations and could therefore be of interest for academic virology and pharmaceutical industry.

2. Results

2.1. Cell rounding Occurs Early after Infection and Correlates with the Applied Virus Dose

Despite significant advancements in automation of the assay-readout [8], traditional endpoint titration methods still suffer from limited precision and/or low sample throughput. To overcome this analytical bottleneck, we sought ways to replace the existing readout methods, typically based on detection of cell lysis, with an assay setup benefiting from the detection of more subtle and earlier virus-induced cellular changes. Furthermore, we implemented a kinetic readout scheme, with acquisition of label-free brightfield images over regular time intervals. The cumulative analysis of time-dependent changes increases the precision and accuracy for each sampled well and alleviates the need to incubate cells until a quantifiable infection endpoint is reached. As a model virus for assay development, we used VSV-GP. VSV-GP is an engineered oncolytic virus, based on vesicular stomatitis virus (VSV), in which the original glycoprotein was replaced by the glycoprotein (GP) of lymphocytic choriomeningitis virus (LCMV) [17]. Due to their robust growth characteristics and high susceptibility for productive VSV-GP infection, BHK-21 were chosen as model host cells. Cells were infected at different multiplicities of infection (MOI) with VSV-GP and imaged at regular time intervals using an automated imaging plate reader. Among the measured cellular parameters, the increase over time in the proportion of rounded cells demonstrated the strongest dose-dependency and greatest signal-to-noise ratio (Figure 1). In each of our analyses, the proportion of rounded cells in relation to the absolute cell number was determined to compensate for time-dependent changes in cell number arising from proliferation and/or virus-induced cytopathic effects. The baseline of rounded cells, which is shared by uninfected cells (mock) and infected cells, prior to onset of infection, depends on the definition of a rounded cell. For our assay, cells in which the ratio of the smallest cell diameter to the largest cell diameter exceeds 0.3 were classified as rounded. VSV-GP infection induced a rapid and dose-dependent increase in the proportion of rounded cells that continued until a common upper asymptote was reached for all MOIs. This progression of virus-induced cell rounding was observed independent of the software algorithm (Gen5 or Columbus) used for image analysis (Figure 1b; Figure S1a). Furthermore, it resembles the changes in cellular fractions determined by a linear classifier, which integrates multiple parameters, such as roundness, intensity, area, and surface texture (Figure S1b). This suggests that cell rounding is not only an easily measurable parameter but also representative for the multitude of cellular changes occurring upon infection. A slight increase in rounded cells was also observed in mock infected wells. This is likely a consequence of growing cell density and not critical, as its extent and time-dependent progression does never reach the levels of infection-induced cell rounding.

To confirm specificity of the observed morphological alterations for VSV-GP infection, cell rounding was tested under conditions in which viral replication was inhibited at different stages. In the presence of a neutralizing antibody, which binds to the viral glycoprotein (GP) and inhibits viral entry, no VSV-GP-dependent rounding of BHK-21 cells could be observed (Figure 2a). Similarly, inhibition of viral RNA transcription by the polymerase inhibitor Vesiculopolin A [18] fully prevented virus-induced cell rounding (Figure 2a). The delayed rounding kinetics in cells infected with the genetically modified VSV-GP ΔM51 are consistent with previous reports [19] showing M protein functions contribute to cell rounding (Figure 2a). The ΔM51 mutation reduces the association of the viral M protein with nuclear core complexes and its ability to inhibit nucleocytoplasmic transport, which might be related to M protein mediated cytotoxicity [20,21,22]. The connection between progression of viral gene expression and cell rounding was further assessed in cells infected with a VSV-GP variant encoding the red fluorescent Katushka protein. The increase in the proportion of Katushka-positive cells was detectable ~7 hours after the onset of cell rounding. However, the dependence of cell rounding on viral gene expression shown above (Figure 2a), together with the finding that the exponential viral growth phase preceded detectability of Katushka (Figure 2b) strongly suggests that this is rather a sensitivity issue than the actual sequence of biological events. Importantly, cell rounding and Katushka expression, which reflects viral protein production, show similar curve progression (Figure 2b), independent of the MOI used for infection (Figure S2a). While at a higher MOI all kinetics shift towards earlier timepoints, the temporal relation between cell rounding, the emergence of Katushka positive cells and release of progeny virions, remains unaffected (Figure S2b). Altogether, these findings show that entry and viral gene expression are a prerequisite for virus-induced cell rounding, which reflects viral growth kinetics in infected cells.

2.2. Titer Determination with the Kinetic Infectious Virus Titer (KIT) Assay

The specificity for viral replication, its strong dose- and time-dependency, and the simple measurement based on label-free imaging qualified cell rounding kinetics as attractive readout parameter for our development of an improved setup of infectivity assays. To optimize method precision and robustness, while maintaining sufficient experimental flexibility for routine use of the assay, we employed a Design of Experiment (DoE) approach [23,24]. Factors chosen to be optimized were cell seeding density, media composition, cell pre-culture time, incubation time between cell seeding and infection and the working range (details and results will be published elsewhere). The setup of the resulting kinetic infectious virus titer (KIT) assay is outlined in Figure 3. Briefly, BHK-21 cells, which are maintained in Glasgow’s MEM (GMEM) supplemented with 6.1 % fetal calf serum (FCS) and 4.3 % tryptose phosphate broth (TPB), are seeded into 96-well plates at a density of 18,000 cells per well 22±4 hours before infection.

To achieve high intermediate precision and accuracy over time, titers are determined in a relative manner, normalized to a reference standard which is included on each plate. The reference standard consists of well-characterized and representative virus material, which is aliquoted and stored under conditions maintaining stability over a long period of time. In our setup the standard curve consists of six MOIs with ½ log₁₀ increment, resulting in a 2.5 log₁₀ assay range (MOI 31.6 – MOI 0.1). Each sample is applied at a single dose, maximizing the throughput. As, of course, the exact titer of the sample is unknown, the applied dose, i.e., the sample volume resulting in the intended target MOI, are estimated based on previous knowledge about the sample type or production process. The intended target MOI is 1.78, the middle of the working range of the assay, to allow valid measurements within the working range of the assay even if the assumed sample titer differs ±1 log10 from the actual titer. If the sample titer cannot be estimated with sufficient certainty, different dilutions of the sample can be applied to ensure that at least one of them falls within the assay working range. The challenge to dilute the samples to a fixed target MOI would be the need to pipet specific volumes for each of the samples. The material and methods section provides a detailed solution, which helps to avoid this error prone procedure and allows automated sample preparation.

Finally, 100 µl of the serially diluted standard as well as the assay dilutions of samples and assay control are transferred to triplicate wells of the 96-well plate containing the cells seeded in 100 µl per well on the previous day. Wells in which 100 µl assay medium are added are included as uninfected control (mock). Even though consistent results were generated irrespective of the positioning on the plate, applying standards and samples in alternating order (Figure 3b) further helps to reduce possible influence of plate effects, also called plate non-uniformity. Immediately after infection, incubation of the plate with imaging of all wells in regular 1 h intervals using an automated imaging reader is started, with a total duration of 24 h.

Figure 3. Workflow for performing the kinetic infectious virus titer (KIT) assay. (a) Experimental outline for determination of infectious virus titers using the KIT assay. (b) Exemplary plate layout for the KIT assay. The plate contains in alternating orientation triplicates of (1) six wells infected with a reference standard virus at ½-log₁₀ serially diluted MOIs (2) uninfected control wells (mock) (3) an assay control which is included on each assay plate and serves for trending purposes of critical assay performance parameters (4) up to 24 samples. (c) Data analysis workflow to calculate infectious virus titers based on kinetics of infection-induced cell rounding. Further details are given in the text and in the Material and Methods section.

Figure 3. Workflow for performing the kinetic infectious virus titer (KIT) assay. (a) Experimental outline for determination of infectious virus titers using the KIT assay. (b) Exemplary plate layout for the KIT assay. The plate contains in alternating orientation triplicates of (1) six wells infected with a reference standard virus at ½-log₁₀ serially diluted MOIs (2) uninfected control wells (mock) (3) an assay control which is included on each assay plate and serves for trending purposes of critical assay performance parameters (4) up to 24 samples. (c) Data analysis workflow to calculate infectious virus titers based on kinetics of infection-induced cell rounding. Further details are given in the text and in the Material and Methods section.

For every timepoint of image acquisition, the proportion of rounded cells is determined for each well. Plotting rounded cells (%) against time post infection (h.p.i.) results in sigmoidal curves, which can be fitted by a model describing an initial lag-phase, followed by an exponential increase in the proportion of rounded cells, which culminates in an upper asymptote (Figure 3c) (see Material and Method section for details about the model). The upper asymptote reflects a stable condition in which all susceptible cells are affected by virus-induced cytopathic effects and most of them meet the classification threshold qualifying them as rounded. Based on the fit, a RT50 value is determined for each curve. The RT50 value describes the time post infection at which the half-maximal extent of cell rounding is achieved (inflection point between lower and upper asymptote of the sigmoidal curves; see blue dotted lines in left panel of Figure 3c). A standard curve can be generated by plotting the six RT50 values of the reference standard dilution series against their respective MOI. On a semi-log scale, the standard curve is described by a simple linear regression model. Based on the RT50 values of the samples, their MOIs can be interpolated from the standard curve (see orange dotted lines in the middle panel of Figure 3c). The interpolated MOI can then be converted into the original infectious virus titer of the sample (see right panel of Figure 3c). As the KIT assay is a relative method, the determined infectious titers of the samples will always be expressed in the same unit as results from the assay used to characterize the reference standard.

With the described setup of the KIT assay, infectious titers of up to 24 samples can be determined on a single 96-well plate. In our routine use of the assay, turnaround time for completion of a full plate is less than 45 minutes and results are available already 24 h after infection.

2.3. Performance Characteristics of the KIT Assay

Infectious virus titers of samples obtained from a formulation development study determined by the KIT assay (3 wells per sample) correlate well (slope = 0.99; R² = 0.78) with titers determined by TCID50 assay (results from 3 independent determinations, i.e., 3 96-well plates averaged per sample) (Figure 4a). Precision and accuracy throughout the working range of the KIT assay was tested by decreasing or increasing the intended target MOI (1.78), resulting in a 10-fold higher or lower sample dilution, respectively. In each of the five experiments performed by different operators, six independent replicate measurements were performed at each target MOI. Irrespective of the target MOI, the KIT assay showed an intra-plate precision with an average CV of 11 %. With an intermediate precision of 17 % CV and an overall recovery of 91 %, the applied setup for the KIT assay yielded accurate results throughout the complete working range (Figure 4b, Table S1, Figure 4c). This shows that with only three wells per sample, variability of the KIT assay is smaller than what can be achieved with our optimized TCID50 assay [8] (44 % CV; with a single determination, i.e., one 96-well plate measured per sample) (Figure 4c).

The KIT assay was used to measure VSV-GP subjected to increased temperature or UV irradiation. The forced sample degradation manifests itself in a delayed rounding of infected BHK-21 cells (Figure 5a). This reflects a loss of infectious virus titers, the extent of which correlates with the applied stress level (Figure 5b). These findings confirm that the KIT assay is a stability-indicating method, which can resolve also small differences in virus infectivity.

2.4. The Principle of the KIT Assay Can be Applied to Different Virus-Cell-Combinations

Morphological changes of host cells are a common phenomenon upon infection with various types of viruses. This suggests that the principle of measuring the kinetics of cell rounding to determine infectious titers is not limited to VSV-GP. To confirm this, we performed a feasibility study applying the principle of the KIT assay to determine titers of Newcastle disease virus (NDV) based on rounding of DF-1 chicken fibroblasts [25]. As observed with VSV-GP in BHK-21 cells, NDV infection induced dose- and time-dependent rounding of DF-1 cells (Figure 6a). While all curves start at a common basal value of rounded cells, the height of the upper asymptote depends on the applied MOI. It is not surprising that the relation between virus dose and cell rounding differs from the observations made with VSV-GP. Factors contributing to these characteristics include the apparently less effective spread and induction of cell death by NDV, which allows a certain level of proliferation of regular shaped DF-1 cells and limits the proportion of rounded cells, especially at lower MOIs. These virus- and host cell-dependencies are reflected by the RT50 values, which were determined using our kinetic curve fit model and exhibit a sigmoidal curve shape, when plotted against the respective MOIs (Figure 6b, Table S2). Precision and accuracy of the measurements for NDV were evaluated by interpolating MOIs of three NDV samples (M 1 – M 3) at four different concentrations. The mean CV of 13 % and 107 % recovery of the actual sample titer highlight that the KIT assay does not only allow precise and accurate determination of VSV-GP titers, but this assay principle can also be applied to completely unrelated virus-cell-combinations, for which titers can be determined with the same high throughput.

3. Discussion

Infectivity is the most critical quality attribute for many virus-based therapeutics because it determines their ability to enter target cells, deliver their therapeutic payload, and/or replicate. Therefore, the ability to monitor this fundamental mode-of action is key for evaluating the potency, safety, and efficacy of virus-based therapeutics. Pharmaceutical industry has thus a strong interest in infectivity assays that allow precise determination of infectious virus titers in a fast and cost-efficient manner. However, traditional assay principles, such as the TCID50 or plaque assay, suffer from high variability and limited sample throughput. This creates an analytical bottleneck, especially for complex process and formulation development studies. Here precise measurements of large sample numbers are required to make informed decisions, which for example allow to maximize production yields and therefore decrease the production costs. To overcome these limitations, we developed a high-throughput method to determine infectious virus titers based on automatic quantification of infection-induced cell rounding in kinetically acquired bright-field images. The kinetic infectious virus titer (KIT) assay (Figure 3) is characterized by its straightforward principle, low running costs, short time to results and superior precision despite a ~10-fold increased throughput compared to the traditional TCID50 assay (Figure 7).

The TCID50 assay has been widely utilized in various virological disciplines for many years due to its simplicity and the independence of special instrumentation. However, as an endpoint dilution assay it requires every sample to be diluted to a point at which infection is absent in the inoculated culture wells or caused by only a single (or few) infectious virus particle. Furthermore, the assay design includes infection in multiple replicates for each dilution step to calculate the TCID50, describing the sample dilution at which 50 % of the inoculated wells shows signs of infection. This high number of pipetting steps and the laborious readout procedure, which often relies on manual microscopic evaluation of infection-induced cytopathic effects (CPE), limit the sample throughput. Furthermore, the need to wait until also a single initial infection event results in an infection outcome, which by manual microscopic evaluation can be unambiguously identified as such, keeps the time to result long. We have previously shown that automation of the sample dilution procedure combined with automated image-based CPE evaluation reduces the hands-on time ~4-fold [8]. Automated image analysis can for both, plaque and TCID50 assays, be combined with virus-specific immune staining. The resulting increased sensitivity allows an earlier measurement with operator-independent results [9,26]. However, this comes at the cost of additional working steps and makes the assay dependent on a constant supply with a critical, possibly expensive, reagent. In our KIT assay, we enabled an early readout based on bright-field images without additional staining steps. To do so, a software algorithm quantifies infection-induced changes in the percentage of rounded cells. Cell rounding occurs early after infection with a variety of viruses and is likely the result of virus-induced subversion of the actin cytoskeleton [27] and/or represents the onset of CPE. Our data show that infection-induced cell rounding occurs already within few hours after infection and depends on viral entry as well as viral gene expression (Figure 1 and Figure 2). The delayed cell rounding kinetics observed for the VSV-GP ΔM51 variant emphasize the involvement of M protein functions in this process and are in alignment with previous reports showing, that mutations disrupting the expression of the alternative translation products M2 and M3 (M33,51A) reduce rounding of infected BHK-21 cells [19]. The specificity of cell rounding for viral infection is further confirmed by the strong time-dependent progression, which correlates with the applied virus dose.

The algorithm-based image analysis allows reliable assessment of subtle infection-induced cellular changes, much earlier than possible by eye. Implementation of this readout procedure could therefore also be an option for traditional infectivity assays, where presence/absence of CPE or number of plaques is determined as at single timepoint. Several examples of using artificial intelligence (AI) to detect CPEs in bright field images have been presented previously [12,28,29]. However, it is often not fully understood, on which properties the AI bases its decision. Using these sophisticated software tools in infectivity assay therefore entails the challenge of robust validation, quality control, and regulatory compliance, especially when applying the assay under GMP-conditions. This is more straightforward if simple and transparent readout parameters, such as cell rounding, are used.

In addition to its user-friendliness, the major benefit of the KIT assay lies within the principle of kinetic measurements. The curve fit of the time-dependent progression of cell rounding takes into consideration the values obtained over the complete assay duration. Thus, a curve parameter such as the RT50 has a much higher accuracy and precision in comparison to a measurement at a single timepoint. This makes the measurement robust, less sensitive to outliers and yields precise results with a CV < 20 % using only 3 wells per sample (Figure 4). Our TCID50 assay, which depends on CPE evaluation at a single timepoint, exhibits a CV of about 45 % with a single determination per sample. This is already quite good for such a traditional infectivity assay, where a variability of ±0.5 log₁₀ is not uncommon [6,7]. However, to reach the precision of the KIT assay, data from six measurements, equaling six 96-well plates, would have to be averaged (Figure 7) [8]. Of note, we have successfully up-scaled the assay to 384-well plates without major adaptations (data not shown). If higher assay precision is required in the KIT assay, this could therefore be achieved by increasing the number of replicates, with only a minor reduction of sample throughput. Due to the kinetic measurement, it is also not necessary to define a fixed timepoint for the readout as images can simply be acquired until all curves reach the upper asymptote, represented by stable percentage of rounded cells. This results in a broad assay range and in most cases abolishes the need to test samples at different concentrations, which has a positive effect on number of samples which can be measured in parallel. As prolonged measurements negatively affect the time to results, a reasonable and practical assay range and thus duration for image acquisition should be chosen. Altogether, the presented design of the KIT assay allows to measure up to 24 samples on a single 96-well plate within 24 hours, while keeping the variability below a CV of 20 %.

As the kinetics of cell rounding reflect productive viral replication, the KIT assay is not only able to measure differences in concentration but also in viral functionality. The KIT assay is thus a stability-indicating method (Figure 5). Due to the determination of sample titers relative to a reference standard, the influence of day-to-day variability, which is often caused by cultured cells, is minimized, what makes the KIT assay a powerful tool for long-term stability studies. Moreover, it is expected that the collection of kinetic data enables the KIT assay to resolve differences in viral fitness, which are reflected by altered infection and replication kinetics. These differences can for example be caused by viral mutations and might be neglected by traditional infectivity assays, which depend on a single endpoint readout.

The power of kinetically acquired data is also appreciated in studies using the xCELLigence Real-Time Cell Analysis (RTCA) system for determination of infectious virus titers [13]. However, due to the need of using specialized plates, measuring cellular impedance causes relatively high running costs and the instrument is provided by a single manufacturer. In contrast, automated imaging readers are commonly available and can be combined either with commercially available or custom-built automated plate incubators. Importantly, independent algorithms with different ways of evaluating cell rounding yielded highly comparable quantifications of the proportion of rounded cells (Figure S1), showing that this readout method does not dependent on a special feature of a specific software. The flexibility regarding the used instrumentation and software is especially of advantage if the assay needs to be aligned between multiple sites with different pre-existing lab configurations as it is often the case in the pharmaceutical industry.

We used our engineered oncolytic virus VSV-GP as a model for the initial assay development. However, it is likely that the KIT assay can also be applied and might thus be of interest for the growing number of other VSV-based therapeutics or vaccines [30,31,32,33,34]. Importantly, we could also show that the KIT assay allows determining infectious titers of unrelated viruses, such as Newcastle disease virus (NDV) (Figure 6). While the process of image acquisition and analysis was identical to VSV-GP, the replication characteristics of NDV in DF-1 cells made it necessary to prolong the overall assay duration and to modify the fit model for generation of the standard curve. Such biological differences regarding the growth rate of host cells, effectiveness of virus spread, or the strength of infection induced cytopathicity, will certainly require slight adaptations in the analysis procedure also for other virus-cell combinations. With the wide range of conventional or AI-based image analysis tools, it should, however, be possible to find a suitable infection-specific morphological measurement parameter with strong time- and dose-dependency for most types of viruses. This makes the kinetic infectivity assay a versatile tool and therefore also attractive for academic research, where in most cases a high level of flexibility is required. Also, the additional level of information regarding virus kinetics might be a bigger advantage for academic research in comparison to pharmaceutical industry, where functional differences or genetic variants will routinely be identified by orthogonal analytics to determine, for example, the proportion of infectious particles or genomic integrity.

4. Materials and Methods

4.1. Cell Lines

BHK-21 cells (CLS #603126) were cultured in Glasgow’s MEM (GMEM) (Gibco #21710) supplemented with 6.1 % fetal calf serum and 4.3 % tryptose phosphate broth (Thermo Fisher #18050039). DF-1 cells (Elabscience #CL-0279) were cultured in DMEM (Gibco #11960) supplemented with 2 % Ultroser G (Sartorius), 4 mM L-Glutamin and 1 mM sodium pyruvate. Cells were passaged 2-3 times per week before reaching full confluence. For passaging, cells were washed with PBS and detached by incubation with TrypLE^TM Select Enzyme (Thermo Fisher #12563011) for 6-8 minutes. Cells were taken up in medium, counted using a NucleoCounter NC-200 (ChemoMetec) and then further cultivated or seeded for an assay. Cells are incubated at 37 °C and 5 % CO₂ in humidified atmosphere.

4.2. Kinetic Infectious Virus Titer (KIT) Assay

Sample dilution and infection procedure

BHK-21 cells were seeded into 96-well flat bottom plates (Thermo Fisher #161093) at a density of 18,000 cells in 100 µl 22±4 hours before infection. All pipetting steps for cell seeding and virus dilutions were either performed manually or in a semi-automated fashion using Integra pipetting robots (Mini96, VIAFLO96, Assist Plus). All dilutions for the reference standard, samples and assay control were prepared using supplemented cell culture medium as diluent.

In our typical assay setup for VSV-GP, the reference standard curve consists of six MOIs with ½ log₁₀ increment, resulting in a 2.5 log₁₀ assay range (e.g., MOI 31.6 – MOI 0.1). To broaden the assay range, lower MOIs can be included in the reference standard curve, but the duration of image acquisition might have to be prolonged to ensure that all curves reach the constant upper asymptote of rounded cells.

In addition to the reference standard, an assay control (optional) can be included. The assay control is a representative virus sample, which is included on every assay plate. It is treated just like all other samples, but the infectious virus titers determined for the assay control should be documented together with other critical assay details (e.g., operator, cell passage, cell pre-culture time, time between seeding and infection, batch ID of critical reagents) to monitor assay performance over time. Just like for all other samples, measurement of the assay control is performed at a single concentration, i.e., MOI. In our assays, the target MOI used for infection was 1.78, which lies in the middle of the standard curve and allows valid measurements within the working range of the assay even if the assumed sample titer differs ±1 log₁₀ from the actual titer. MOI calculations are made based on expected titers, which are informed by previous knowledge about the sample type or production process. One challenge that comes with the constant target MOI is that different sample volumes would have to be used for infection for each sample, which is error prone and represents a considerable limitation for assay automation. To bypass this issue, expected sample titers are rounded to the next integer log₁₀ value and only then used to calculate the sample volumes required to reach the target MOI (refer to Table M1 for an example of possible dilution schemes). Serial 10-fold sample dilutions are used to obtain the required sample concentration in volumes that are easy to pipet. Due to the previous rounding step, the volumes of pure sample required to infect at MOI 1.78 always differ by a factor of 10 between samples with different expected titers. Thus, only the number of 10-fold pre-dilutions must be adapted, but all pipetted volumes remain constant. This is also true for the final dilution step (in our example: dilution of 192 µl pre-diluted virus in 408 µl medium), in which the highest 10-fold pre-dilution is used to prepare the assay dilution, which is used to infect the BHK-21 cells.

Finally, 100 µl of the serially diluted standard as well as the assay dilutions of samples and assay control were transferred to replicate wells (triplicates in the standard assay setup) of the 96-well plate containing the cells seeded on the previous day. Wells in which 100 µl assay medium (no virus) were added served as uninfected control (mock). 100 µl medium was added to unused wells to maintain consistent culturing conditions throughout the whole plate.

The overall infection procedure was the same when titers of Newcastle disease virus (NDV) were determined in DF-1 cells. The only differences were a seeding density of 15,000 cells/well and different MOIs used for infection.

Automated image acquisition

Immediately after infection, incubation of the plate with imaging of all wells in regular intervals (1 h for VSV-GP, 2 h for NDV) using an the automated Cytation5 (Agilent/BioTek) multimode imaging reader in combination with the automated BioSpa (Agilent/BioTek) incubator was started. The total duration of image acquisition depended on the virus, the MOIs used for infection and the specific purpose of the experiment. Imaging was performed using a wide field of view objective with 4x magnification and laser autofocus. Expression of the red fluorescent Katushka protein was monitored at an excitation/emission wavelength of 584±20 nm/625±20 nm.

Image analysis

Analysis of cell rounding other infection-induced morphological changes was performed on bright-field images. If not stated otherwise, image analysis was performed using Gen5 (Agilent/BioTek). Cellular analysis in Gen5 was restricted to an area with 3000 µm diameter in the center of the image, avoiding the outmost areas of the well with possible lower optical quality. Within the total number of identified cells (object size 5-75 µm) the software determined the rounded cell proportion. Cells are analyzed by approximating their shape with an ellipse and defined as rounded if the ratio of its smallest diameter to the largest diameter exceeds 0.3. A different threshold value would shift the absolute value of % rounded cells of the upper and lower asymptote. Importantly, however, the time-dependent progression of cell rounding would be maintained and due to calculation of the titers relative to a reference standard, this would not significantly affect the results.

Columbus^TM (PerkinElmer) was used as a complementary software tool to determine cell rounding and other cellular changes occurring upon infection. Two approaches with different levels of complexity were applied. First, as done by the Gen5 software, cells were defined as rounded if the ratio of their smallest diameter to the largest diameter exceeded 0.3. Second, a linear classifier model based on the Columbus^TM PhenoLOGIC^TM technology was applied to discriminate elongated versus rounded cells. The linear classifier was trained with both cell populations (around 200 cells each across many different images) and then all images were analyzed. For classification of cells, Columbus integrated different parameters including roundness, width-length-ratio, cell area, intensity and cell surface texture.

Statistical analysis

In the KIT assay the outcome of the image analysis is the percentage of rounded cells in each well at a specific point of time after infection. If these two factors are plotted against each other with time post infection on the x-axis and % rounded cells on the y-axis, this results in a sigmoidal curve (Figure 3). The following kinetic curve fit model is used to determine the RT50 value, i.e., the time after infection at which the half-maximal cell-rounding is reached, from each of the sigmoidal curves.

$$y=LA\frac{UA-LA}{1+{10}^{-k(x-RT50)}}$$

In which: LA = lower asymptote. Value is dependent on the baseline of rounded cells per well. UA = upper asymptote. Value is the maximum % rounded cells, that is reached based on the defined image analysis parameters. k = slope at the inflection point of the curve. RT50 = time after infection at which half-maximal extent of cell rounding is reached, which is identical to the inflection point of the curve.

The RT50 values of the different doses of the reference standard are plotted against the corresponding known MOIs, resulting in a standard curve (Figure 3). For VSV-GP in BHK-21 cells this standard curve is reflected by a semi-log line, whereas a sigmoidal four parameter logistic fit model was applied for NDV in DF-1 cells. Based on the determined RT50 values of the assay control and the samples with unknown titer, their MOIs can be interpolated from this standard curve. The interpolated MOIs, together with the known number of cells seeded per well and the volume of the samples that was added per well, is then used to calculate the infectious titer. The calculated infectious virus titer has the same unit in which infectivity of the standard sample is expressed.

4.3. TCID50 Assay

The semi-automated version of the TCID50 assay with algorithm-based evaluation of cytopathic effects (CPE) was performed as described previously [8]. BHK-21 cells were seeded in 96-well plates. 24 h later, adherent cells were infected with 8 replicates of a ½-log10 serial sample dilution. At ~72 hours post infection (h.p.i.) bright-field images were acquired using a Cytation5 Multimode Imaging Reader (Agilent) with 4X objective. An algorithm determined the cell area in each well. Wells with a cell area below a certain threshold were classified as CPE-positive. The model of Spearman and Kärber of Spearman/Kärber or Reed/Münch [3,4] to was used to calculate the infectious virus titer.

4.4. qRT-PCR

TaqMan^TM qRT-PCR was used to monitor VSV-GP replication in infected BHK-21 cells. Genomic viral RNA was isolated from cell culture supernatants using the KingFisher Flex System (Thermo Fisher) in combination with the MagMAX^TM Viral/Pathogen Nucleic Acid Isolation Kit (Applied Biosystems #A42352). Reverse transcription of RNA and cDNA amplification were conducted in single-tube format using the TaqMan^TM Fast Virus 1-Step Master Mix (Applied Biosystems #4444432). VSV-GP-specific TaqMan probe and primers binding the m and gp gene were used. Signals were quantified relative to a synthetic RNA standard are results are expressed as genomic copies (GC) per ml.