Next Article in Journal
Association of Macrolide Resistance Genotypes and Synergistic Antibiotic Combinations for Combating Macrolide-Resistant MRSA Recovered from Hospitalized Patients
Next Article in Special Issue
A Predictive Model for Severe COVID-19 in the Medicare Population: A Tool for Prioritizing Primary and Booster COVID-19 Vaccination
Previous Article in Journal
Oxidative Stress and Mitochondrial Damage in Dry Age-Related Macular Degeneration Like NFE2L2/PGC-1α -/- Mouse Model Evoke Complement Component C5a Independent of C3
Previous Article in Special Issue
Symptom and Age Homophilies in SARS-CoV-2 Transmission Networks during the Early Phase of the Pandemic in Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unpredictable, Counter-Intuitive Geoclimatic and Demographic Correlations of COVID-19 Spread Rates

by
Hervé Seligmann
1,2,*,
Nicolas Vuillerme
1 and
Jacques Demongeot
1
1
Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical & Labcom CNRS/UGA/OrangeLabs Telecom4Health, Faculty of Medicine, University Grenoble Alpes (UGA), 38700 La Tronche, France
2
The National Natural History Collections, The Hebrew University of Jerusalem, Jerusalem 91404, Israel
*
Author to whom correspondence should be addressed.
Biology 2021, 10(7), 623; https://doi.org/10.3390/biology10070623
Submission received: 15 January 2021 / Revised: 16 June 2021 / Accepted: 29 June 2021 / Published: 5 July 2021
(This article belongs to the Special Issue Theories and Models on COVID-19 Epidemics)

Abstract

:

Simple Summary

Rates of viral spread during first and second waves of the COVID-19 pandemic for USA states, and for consecutive nonoverlapping periods of 20 days for the USA and 51 countries across the globe associate with mean temperature, elevation, population density and age. Some associations switch directions when comparing different periods. Even population density, which presumably should always increase viral spread, at some periods seems to decrease spread rates. We also observed systematic inversions between spread rates estimated at 80–100 day intervals. These patterns remain unexplained and suggest difficulties in managing and predicting the pandemic, in particular, negative correlations between population density and spread rates, which were observed in independent samples and at different periods. Putatively, confinements could produce these patterns, by selecting viral strains with longer contagiousness and/or latent periods.

Abstract

We present spread parameters for first and second waves of the COVID-19 pandemic for USA states, and for consecutive nonoverlapping periods of 20 days for the USA and 51 countries across the globe. We studied spread rates in the USA states and 51 countries, and analyzed associations between spread rates at different periods, and with temperature, elevation, population density and age. USA first/second wave spread rates increase/decrease with population density, and are uncorrelated with temperature and median population age. Spread rates are systematically inversely proportional to those estimated 80–100 days later. Ascending/descending phases of the same wave only partially explain this. Directions of correlations with factors such as temperature and median age flip. Changes in environmental trends of the COVID-19 pandemic remain unpredictable; predictions based on classical epidemiological knowledge are highly uncertain. Negative associations between population density and spread rates, observed in independent samples and at different periods, are most surprising. We suggest that systematic negative associations between spread rates 80–100 days apart could result from confinements selecting for greater contagiousness, a potential double-edged sword effect of confinements.

1. Introduction

Spread parameters of daily new confirmed COVID-19 cases (calculated as the slope of their logarithmic regression curve) estimates viral contagiousness. Previously, it was shown that first wave spread, when comparing different countries, decreases with mean annual temperature [1], and the opposite trend with temperature occurs for second wave spread parameters [2].
This is in line with observations on variation in spread across different regions of Italy, in March 2020 (first wave period, negative correlation with temperature) and in May 2020 (second wave period, positive correlation with temperature) [3].
First and second wave spreads differ also in terms of other factors: The spread of the first wave increases with the median population age and decreases the later the date of its onset. On the contrary, the spread of the second wave decreases with the median population age and increases the later the date of its onset. First and second waves also differ because we detected no associations between second wave slopes and mean country elevation, while first wave slopes increase with elevation up to 900m and decrease beyond that approximate altitude [2].
Inversion of trends for these independent covariates is difficult to explain. One could invoke different explanations for each factor. A common, most parsimonious explanation could involve deterministic mutation dynamics, resulting in parallel evolution of distinct virus populations [4,5].
This working hypothesis expects cycles in the epidemiological behavior of the virus, where pattern inversions occur each time a sufficient number of mutations cumulated to cause a switch in secondary structure of the single-stranded RNA coronavirus. This would explain the inversion in patterns between first and second waves. Because mutations cumulate on average proportionally to time, one expects that epidemiological pattern inversions occur after a fixed time interval.
Here, we calculate spread rates for daily new confirmed COVID-19 cases for consecutive, nonoverlapping time windows of 20 days since the start of the pandemic in 51 different countries, and in each of the 50 states of the USA (and the district of Columbia) until the end of 2020, and examine these spread rates for inversions in patterns of correlations with some environmental factors (temperature, elevation) and population properties (density, median age), and across spread rates for different time windows. This approach, which calculates spread rates for each consecutive period of 20 days, avoids difficulties inherent to the objective definition of the start and end of different pandemic waves.

2. Materials and Methods

We used the methods as in previous analyses [1,2]. The coefficients (slopes) of regression analyses are considered as estimates of viral spread rates. We adjust the exponential model y = a × exp(b × x), where y is the daily number of new confirmed COVID-19 cases, x is the number of days since wave onset, a is a constant and b is the slope. The log-transformed version of this is ln y = ln a + b × x. Daily numbers of new cases per countries are from the site www.worldometers.info/coronavirus/ (accessed on 25 December 2020).
Spread rates are determined for running windows of 20 consecutive days, where consecutive time windows are nonoverlapping. Hence, if the first new case in a country occurred 280 days ago, 14 spread rates were calculated, one for each of the nonoverlapping 20-day periods. Spread rates are estimated for these running windows as described above, the slope of the exponential model of new daily cases with time.
For countries, sources for mean elevation, mean temperature, population density and median population age and numbers of genomic variants are as previously described [2]. For the 50 states of the USA, we used the following sources in Table 1:
We also analyze daily new confirmed case data for each of the 50 states of the USA and the district of Colombia, estimating first and second wave spread rates as done previously for countries across the globe [1,2], using for the USA the exact same methods as described previously [1,2].
Note that the slopes we evaluate from daily data in daily new cases estimate the acceleration in increase of new cases: in the first differential model proposed for the COVID-19 outbreak [6], the variable is the number of cumulative cases at time t, denoted Xt, the velocity is the number of daily new cases, denoted Yt = (Xt − Xt − 1)/(t − (t − 1)), and the acceleration is the slope of the curve of these daily new cases, denoted Zt = (Yt − Yt − 1)/(t − (t − 1)). This acceleration as observed in simulation of the differential model, has a first phase where daily new cases grow exponentially, until an inflection point is reached with acceleration zero before the saturation phase.
We tested for normality of the distribution of variables used in correlation tests with the Kolmogorov–Smirnov test, as available at Kolmogorov–Smirnov Calculator (Test of Normality) (socscistatistics.com, accessed 20 December 2020). For variables statistically significantly diverging from normality (p < 0.05), we calculated both parametric Pearson and nonparametric Spearman rank correlation coefficients at Spearman’s Rho Calculator (Correlation Coefficient) (socscistatistics.com, accessed on 20 December 2020)).

3. Results

3.1. Mutations Cumulate over Time during the COVID-19 Outbreak

Figure 1 plots the total number of genomic variants described for a country since the start of the pandemic, on 20 December 2020, as a function of the numbers of days since 1 January 2020 when the cumulated number of confirmed cases in a country reached 100 cases. This method for estimating the duration of the pandemic is not adequate for states/countries for which the pandemic presumably stopped before end 2020. These are excluded from analyses. Overall, the longer the pandemic is ongoing in a country (few days since 1st of January), the more mutants have been described in that country. The negative correlation confirms that numbers of mutations are proportional to time. A similar result was found on 31 May 2020 [2]. Hence, this phenomenon is stable over time.
It is probable that numbers of variants reported from a country reflect also the sequencing efforts in that country, which vary among countries. However, the association in Figure 1 is too strong to be only due to this factor, and probably reflects also, or even more, an actual natural cumulation of mutant variants. Indeed, along that rationale, population size should also increase variant numbers, but no such effect is observed, as small and large countries fit equally well the curve.

3.2. Spread Parameters for First and Second Waves in the USA

Table 2 presents for each of the 50 states of the USA onset dates and slopes of first and second waves, as detected by visual inspection as previously [2]. At mid-August, 42 among 50 states have a second wave. Second wave slopes are lower than first wave slopes in all but three states: Oklahoma (slopes basically identical), and Kansas and Ohio (second wave spread greater than first wave spread).
Hence, overall, spread decreased from the first to the second wave, from mean first wave slope = 0.1902 to mean second wave slope = 0.0585.
In the USA, first wave slopes decrease with time since first wave onset (r = −0.5306, two-tailed p = 0.000074, Spearman rank correlation coefficient rs = 0.415, p = 0.003). This is similar to the decrease in slopes previously reported for comparisons across countries [2]. Comparing countries, second wave slopes increase with time since second wave onset [2]. However, for the USA, no such increase could be detected, the trend might fit with that observed for the first wave (r = −0.29, p = 0.06, Spearman rank correlation coefficient rs = −0.365, p = 0.019). Also contrasting with previous analyses of variation between countries, there is no association in the USA between first and second wave slopes and mean annual temperature. A positive association exists between first wave slope and time since second wave onset (r = 0.4469, two-tailed p = 0.003), which could indicate that high initial spread rates contribute to early onset of ulterior waves. Note that this effect does not explain the absence of second waves in several states with particularly high first waves, such as Illinois, New Jersey and New York.
First wave slopes increase with population density (r = 0.3777, one-tailed p = 0.0034, Figure 2A). This association is expected, but was not observed when comparing different countries [1,2]. Six New England states have relatively low slopes considering their densities. Excluding them from calculations increases the strength of this correlation with density (r = 0.5538, one-tailed p = 0.00005).
Density in principle increases spread rate. This was not observed for first and second wave slopes when comparing countries [1,2]. However, this density principle is confirmed for first wave slopes from the USA (Figure 2A). Unexpectedly, second wave slopes decrease with population density (r = −0.459, two-tailed p = 0.0022, Figure 2B). Unexplained pattern inversions between first and second waves have been reported for covariates such as temperature, time since wave onset and median population age when comparing different countries [2]. Decrease of slopes with density in the USA is foremost astounding because it does not at all fit with any knowledge on disease spread.

3.3. Visual Inspection vs. Objective Statistical Analysis

Arguably, determining onsets of waves from visual examinations of graphs is subjective. This issue has been raised in the past [7,8], but has no obvious simple solutions and requires extensive simulation analyses curtailed to each specific dataset, meaning in this case for each state [9,10]. We therefore use a simplified method. Statistical analyses based on calculating Pearson correlation coefficients r for a running window of 20 days determine a local maximum for r within five days of the second wave onset date determined visually for 80 percent of the countries examined [2]. We applied this method for the 42 second waves detected for the USA (Dis, Table 2). Most of the onset dates determined by running windows (62%) are within five days of the date determined visually, in line with similar previous analyses for other countries [2]. The spread rates calculated for nonoverlapping, consecutive windows of 20 days presented below do not suffer difficulties inherent to defining wave onset.

3.4. Spread Rates for Consecutive Nonoverlapping 20-Day Windows

We estimated spread rates for 14 and 15 consecutive, nonoverlapping 20-day periods for the 50 states of the USA and the district of Colombia (Table 3), and 51 countries across the globe (Table 4). Note that because the start of the pandemic varies among states and countries, spread rates are not exactly for the same period when comparing states/countries. Numbers of states and countries where rates are positive are highest for the earliest period (approximately end of February to beginning of March, day 20). This number is lowest during the warmer spring and summer months. This is in line with an overall decrease of spread with high temperatures.
Table 5 presents the pairwise Pearson correlation coefficients for the spread rate data from Table 3 and Table 4 (USA, above diagonal, data from Table 3; countries across the globe, below diagonal, data from Table 4). These 91 and 105 r’s include 16 cases (eight cases for the USA) where the correlation between spread rates for two periods are positive and statistically significant (p < 0.05, 2-tailed tests), and 14 and 19 cases, respectively, where r is negative and statistically significant (p < 0.05, 2-tailed tests). Among the 16 statistically significant positive cases, a statistically significant majority, 13, are between consecutive periods (two-tailed sign test, p = 0.011). This indicates that spread rates are relatively similar across two consecutive periods of 20 days, hence for about 40 days.
The distribution of statistically significant negative correlations among spread rates from Table 3 and Table 4 in relation to time periods is greater than for positive correlations, but 13 among 14 cases are for periods separated by 40–140 days, on average 88 days with a standard deviation of 30 days (Figure 3). Two statistically significant negative correlations between spread rates are separated by 160–240 days.

3.5. Environmental Correlates of Spread Rates across Time

Table 6 presents Pearson correlation coefficients r between spread rates in Table 3 and Table 4 and environmental variables from Table 2 and Table 7. Results include inversions of directions of correlations, as also shown in Table 5 and Figure 3. For mean annual temperatures, around May (day 100), spread rates increase with temperature in the USA and in the world, as previously reported for countries across the world [2], and for the period of December (day 280, USA only). Temperature decreases spread rates at the end of summer and beginning of autumn in the USA and across the world. Results for mean altitudes show similar pattern inversions, including inversion between consecutive periods at days 180 vs. 200 (August vs. start of September, USA only). For the USA, an inversion in correlation directions also occurs with median age, with positive vs. negative statistically significant correlations at days 180 vs. 200. The negative correlation of spread rates with age corresponds to the time when schools open (end of August to beginning of September). Spread rates increase again with median age at the end of autumn (day 280, first part of December). For countries across the world, median age increases with spread rates at pandemic onset, and in summer–autumn, for the period spanning from day 160 to 240. In contrast, for the spring periods from days 60 to 100, spread rates decrease with median age when comparing countries across the world. This might result from seasonal differences between age-related behaviors, but other explanations are also possible.
The only statistically significant correlations with population density when comparing countries across the world are detected using nonparametric Spearman rank correlation analyses, with the only positive associations fitting trivial expectation on consecutive periods spanning days 180 and 200 (p < 0.05, one-tailed tests). Density is expected to increase spread rates. Such a positive association is observed at pandemic onset in the USA. This pattern is inversed for three consecutive periods in the spring in the USA and at days 80 and 280 for countries across the world.
Negative associations of viral spread rate with density are unexpected and make no sense in terms of classical epidemiological understanding, as these imply greater spread rates in populations with low density.

4. Discussion

4.1. A Potential Epidemiological Explanation for Spread Rate Inversions

The inversion of patterns with time and other variables (altitude, temperature, median population age) for slopes of pandemic waves renders predictions particularly dubious. In addition, the negative correlation between second wave slopes and USA density (Figure 2B) contradicts a proven and accepted epidemiological principle. The inversion of directions of correlations between spread rates and several environmental factors (Table 6) confirms previous observations regarding temperature [2]. Periods until inversions vary among environmental factors. For example, associations of spread rates with age flip directions from increasing to decreasing with median age at the transition between summer vacations and the beginning of the school year, at least in the USA. One could find for each of the inversions in directions in statistically significant r values in Table 6 a specific explanation, such as that stated above about school openings. However, the almost systematic inversion of correlation directions for most environmental factors suggests a common cause for most of these pattern inversions.
This point for a more general cause for pattern inversions is strengthened by correlations between spread rates at different periods, for USA states and for countries across the world. The spread rate hierarchy among US states and among countries is systematically inverted after 80–90 days (Table 4 and Figure 3). The cause for these inversions is unclear, but this inversion pattern is systematic, as spread rates at any period will be inversely proportional to those at an ulterior time, typically 80–90 days later. Note that 80–90-day intervals between slopes are too long for these slopes to be part of the ascending and descending parts of the same wave, as proven in the next section.
A first potential explanation is that spread rates vary according to cycles of maxima, intermediate and minima in daily case numbers, and that countries with the highest maxima have the lowest minima and those with the lowest maxima have the highest minima, about 80–100 days later. This hypothesis assumes that the same rules also apply to periods with intermediate spread rates, and that ups and downs in spread rates are synchronized between most countries. The latter assumption is adequate for the first wave in late winter 2020, but could not hold throughout the complete period analyzed, also because countries and states vary in confinement periods, lengths and efficiencies, which causes decreases in spread rates. Hence, additional analyses are required to understand the causes of these pattern inversions; however, the systematic aspect of these inversions in spread rates is such that these should be accounted in the predictions and policymaking. This is because results of correlation analyses hint at the possibility that steep decreases in spread rates cause steep ulterior increases. In that respect, it might be optimal to mitigate variations in spread rates by avoiding drastic policies suddenly increasing or decreasing spread rates.

4.2. Negative Correlations between Spread Rates Are Not between Ascending and Descending Parts of the Same Wave

We test here the hypothesis that negative correlations between spread rates from different periods reflect negative associations between the steepness of the rate during the ascending and the descending parts of the same wave can be tested. This is easily tested using data in Table 3 and Table 4, because spread rates in the ascending phase of a wave are positive, and those in the descending phase of a wave are by definition negative. So, if the hypothesis is correct, for statistically significant negative correlations in Table 5, we should see systematically the opposite sign for rates in the two periods considered. If this inversion of signs is not observed, the ascending/descending phase hypothesis is incorrect.
We tested this for each of the negative correlations with p < 0.05 in Table 5, as presented below for a specific example. We considered the negative association between spread rates for the USA between periods of 20 days starting at days 100 and 240 after the start of the outbreak. For these periods, among the 51 states and district of Colombia, there are 26 and 48 positive spread rates. Hence, the probability of random assortment of identical signs between these two periods expects 26 × 48/51 + 25 × 3/51 = 25.94 where spread rates in the two periods have the same sign, and 51 – 25.94 = 25.06 cases where signs are inverted between the two periods. The observed number of inverted signs between these periods are 28. This slight increase as compared to the expected has p = 0.41 according to a chi-square test.
This analysis was done for all negative correlations with p < 0.05 in Table 5. There are 10 (USA, 5; countries across the world, 5) negative associations between spread rates where sign inversion between periods is significantly greater than random predictions, and 23 cases where this effect is not statistically significant as in the above example. This means that the hypothesis of ascending/descending phases of a wave could contribute in some cases to the phenomenon of negative associations observed between spread rates from different periods, but is not the main cause of these negative associations.

4.3. Does Confinement Increase Ulterior Spread Rates?

One mechanism that could be invoked in this context is that strict confinement policies select for highly contagious viral strains and/or long contagion periods that are more likely to overcome confined conditions. As case numbers decrease during confinements, reopening societies to normal activity unleashes these more contagious viral strains while potentially competing strains with lower infection abilities disappeared, resulting in higher spread rates at the level of the whole population. This dynamic would presumably be visible at the level of analyses such as those done here at time lags of 80–90 days, and would presumably cause the systematic inversions in spread rates among countries and regions described here. Note that the natural trajectory of viral evolution typically increases contagiousness while at the same time decreases pathogenicity, also in the absence of confinement. Confinements hasten the process of evolving greater contagiousness. Steep increases in recent months (autumn of 2020) in the UK, apparently associated with a new highly contagious viral variant, seem in line with this working hypothesis (Figure 4).
A coolheaded analysis and evaluation of above and future results is crucial for better managing the pandemic.

4.4. Intrinsic vs. Extrinsic Constraints

We discuss below additional potential, more speculative and genetics-oriented explanations for inversions of directions of correlations with spread rates.
One hypothesis is that self-hybridization of the virus’ single-stranded RNA genome protects nucleotides forming stems against mutations, while favoring mutations in loops [11]. Mutation cumulation presumably causes deterministic switches between few optimal structures which differ in their properties in relation to temperature, etc. Such switches between secondary structures have been suggested for COVID-19 after small insertion/deletion mutations [12].
These patterns remind the little-known phenomenon of negative heritability. Usually, heritabilities of traits are such that offsprings resemble their parents: if a parent is for a given trait above average, his/her offspring will on average be also above average, and parents below average for that trait have on average offspring which are below average for that trait. In short, heritabilities are correlations between traits of parents and offspring. These correlations are usually positive. Surprisingly, in some cases, negative correlations were observed, meaning negative heritabilities, where above average parents produce below average offspring for given traits, and vice versa for below average parents [13,14]. These phenomena are not statistical artefacts [15] but are difficult to reproduce empirically [16]. This is expected when assuming negative heritability results from selection under changing environmental conditions [17].
Inversion of trends, such as for negative heritabilities, could presumably result from drastic environmental changes. For example, levels of channeling of Sorghum bicolor plant populations towards developmental trajectories better adapted to NaCl salinity after early low level NaCl exposure increased with population variability [18] in the first generation exposed to salinity, but decreased with populational variability in their offspring [19,20]. This is in line with the concept that COVID-19 adapts to its new human host and to various environments inhabited by that host.
It is possible that a large part of the variation in slopes during the pandemic is due to factors intrinsic to evolutionary trajectories of the virus, genetic ones included. This does not exclude environmental effects due to host population age structure, density, temperature and altitude, among others. However, inversions of trends with these cofactors at different periods of the pandemic suggest that unknown intrinsic factors have major effects on the evolution of the pandemic.
Another approach stipulates that variations across different locations reflect to large extents the dynamics at a single location, at different times. This principle where spatial variation is tantamount to temporal variation has been observed in many different contexts, such as in astrophysics where larger distances are interpreted as reflecting more ancient times [21], and in forest stage succession [22,23,24], ecological communities [25] and biomolecules such as those involved in ribosomes [26,27,28,29,30,31]. The principle is also applied in the renown Haeckelian statement “ontogeny recapitulates phylogeny” [32,33]. Variation among individuals in directional asymmetry is interpreted as reflecting different termination timings of development [34]. The principle exists within the relationship of the genetic code and the ribosome [35]. Similarly, the average order of translation of amino acids in single proteins reflects the order of integration of amino acids in the genetic code [36]. At the level of the pandemic, this would suggest that the macroscopic spread rates correspond to viral replication cycles and/or rates of production of viral particles. This hypothesis should be considered as suggestive at best, but could be proven useful as a working hypothesis for the long-term dynamics of this still unknown disease at the level of individuals.
Pattern inversions add to the problem of uncertainty in the data and in predictions [37,38]. Analyses of spread rates estimated for consecutive nonoverlapping periods of 20 days confirm results obtained for visually determined waves. They show systematic inversions in spread rates between countries/regions, which occur approximately every 80–90 days, across the whole period studied. The causes for this are not well understood but might indicate compensation mechanisms for reduced spread rates during confinement in periods following the confinement. This reduction is shown and analyzed in many countries imposing a confinement like England, France, Germany, Iran, Italy, Netherlands, Spain, United States and China compared to countries like Sweden and South Korea, which did not implement mandatory stay-at-home confinements [39,40,41]. A stochastic modeling of the reducing effect of confinement is possible and shows a control depending on the characteristics of the countries concerned and on the early or late nature of this mitigation policy [42].
However, confinements might select for more contagious viral strains, and/or viral strains that are contagious for longer periods. The increased contagiousness and possibly pathogenicity [43,44,45] of these viruses would cause higher spread rates after deconfinement and therefore reduce the positive impact of restrictive measures such as the lockdown. This working hypothesis on effects of lockdowns arises from the results presented above and could be tested by ulterior analyses specifically designed to examine this working hypothesis.

5. Conclusions

In the USA, first and second wave slopes are not correlated with temperature, median age or time since wave onset, but with population density. The principle of inverted trends between first and second waves upholds; density effects are positive/negative for first/second wave slopes. Negative associations between population density and viral spread rates are also observed when examining countries across the world. Such negative associations of viral spread with population density are not compatible with our present understanding of the epidemiology of infectious agents. These inversions of directions of associations between viral spread rates and environmental and populational variables are confirmed by analyses of consecutive nonoverlapping periods of 20 days. We also observe that viral spread rates at any time during the study period are inversely proportional to rates 80-90 days later. These results were observed in two independent samples, US states and 51 countries across the globe. They stress that pandemic dynamics are misunderstood and probably mismanaged.

Author Contributions

Conceptualization, H.S. and J.D.; methodology, H.S.; validation, J.D., N.V. and H.S.; formal analysis, H.S.; investigation, H.S. and J.D.; resources, N.V.; writing—original draft preparation, H.S.; writing—review and editing, J.D.; visualization, H.S.; supervision, H.S. and J.D.; funding acquisition, N.V. and J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the International Excellence Fellowship of Karlsruhe Institute of Technology (KIT) funded within the Framework of the University of Excellence Concept “The Research University in the Helmholtz Association I Living the Change”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are obtained from readily available public databases that are cited in the text.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Demongeot, J.; Flet-Berliac, Y.; Seligmann, H. Temperature decreases spread parameters of the new covid-19 cases dynamics. Biology 2020, 9, 94. [Google Scholar] [CrossRef] [PubMed]
  2. Seligmann, H.; Iggui, S.; Rachdi, M.; Vuillerme, N.; Demongeot, J. Inverted covariate effects for first vs mutated second wave Covid-19: High temperature spread biased for young. Biology 2020, 9, 226. [Google Scholar] [CrossRef] [PubMed]
  3. Malki, Z.; Atlam, E.-S.; Hassanien, A.E.; Dagnew, G.; Elhosseini, M.A.; Gad, I. Association between weather data and COVID-19 pandemic predicting mortality rate: Machine learning approaches. Chaos Solution Fractals 2020, 138, 110137. [Google Scholar] [CrossRef]
  4. Forster, P.; Forster, L.; Renfrew, C.; Forster, M. Phylogenetic network analysis of SARS-CoV-2 genomes. Proc. Natl. Acad. Sci. USA 2020, 117, 9241–9243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Esteban, D. Quasispecies dynamics in disease prevention and control. Virus Popul. 2020, 263–299. [Google Scholar]
  6. Liu, Z.; Magal, P.; Seydi, O.; Webb, G. Understanding Unreported Cases in the COVID-19 Epidemic Outbreak in Wuhan, China, and the Importance of Major Public Health Interventions. Biology 2020, 9, 50. [Google Scholar] [CrossRef] [Green Version]
  7. Lund, R.; Reeves, J. Detection of Undocumented Changepoints: A Revision of the Two-Phase Regression Model. J. Climate Sci. 2002, 15, 2547–2554. [Google Scholar] [CrossRef]
  8. Percival, D.B.; Rothrock, D.A. “Eyeballing” trends in climate time series: A cautionary tale. J. Climate Sci. 2005, 18, 886–891. [Google Scholar] [CrossRef] [Green Version]
  9. Ghorbanzadeh, D.; Picard, D. Etude asymptomatique et pratique du comportement de deux tests de détection de rupture. Statistique et Analyse de Données 1992, 16, 63–84. [Google Scholar]
  10. Ondo, J.C. Étude comparative des tests de stationnarité. Ph.D. Thesis, Québec University, Quebec City, QC, Canada, 2002; p. 288. [Google Scholar]
  11. Krishnan, N.M.; Seligmann, H.; Raina, S.D.; Pollock, D.D. Detecting gradients of asymmetry in site-specific substitutions in mitochondrial genomes. DNA Cell Biol. 2004, 23, 707–714. [Google Scholar] [CrossRef]
  12. Manzourolajdad, A.; Xu, Z.; Ebrahimi, D. Novel Polybasic Cleavage Site in SARS-CoV-2 Genome Is Likely to Induce a Major Change in the RNA Secondary Structure. Preprints 2020. [Google Scholar] [CrossRef]
  13. Janssen, G.M.; de Jong, G.; Joosse, E.N.G.; Scharloo, W. A negative maternal effect in spring tails. Evolution 1988, 42, 828–834. [Google Scholar]
  14. Haldane, J.B.S. The negative heritability of neonatal jaundice. Ann. Hum. Genet. 1996, 60, 3–5. [Google Scholar] [CrossRef] [PubMed]
  15. Steinsaltz, D.; Dahl, A.; Wachter, K.W. On negative heritability and negative estimates of heritability. Genetics 2020, 215, 343–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Stam, E.; Isaaks, A.; Ernsting, G. Negative maternal effect revisited: A test on two populations of Orchesella cincta L. (Collembola: Entomobryidae). Evolution 2017, 52, 1839–1843. [Google Scholar] [CrossRef] [PubMed]
  17. Bonduriansky, R.; Crean, A.J.; Day, T. The implications of nongenetic inheritance for evolution in changing environments. Evol. App. 2011, 5, 192–201. [Google Scholar] [CrossRef]
  18. Seligmann, H.; Amzallag, G.N. Adaptive determinism during salt-adaptation in Sorghum bicolor. Biosystems 1995, 36, 71–77. [Google Scholar] [CrossRef]
  19. Seligmann, H. Transmission of acquired adjustments to salinity in Sorghum bicolor. Biosystems 1997, 40, 257–261. [Google Scholar] [CrossRef]
  20. Seligmann H 1998 Effect of environmental complexity on salt-adaptation in Sorghum bicolor. Biosystems 1998, 45, 213–220. [CrossRef]
  21. Saadeh, D.; Feeney, S.M.; Pontzen, A.; Peiris, H.V.; McEwen, J.D. How isotropic is the universe? Phys. Rev. Lett. 2016, 117. [Google Scholar] [CrossRef] [PubMed]
  22. Horn, H.S. Forest succession. Sci. Am. 1975, 232, 90–101. [Google Scholar] [CrossRef]
  23. Wardle, P. Primary succession in Westland National park and its vicinity. New Zealand J. Bot. 1980, 18, 221. [Google Scholar] [CrossRef]
  24. Reyes, H.A.; Ferreira, P.F.A.; Silva, L.C.; da Costa, M.G.; Nobre, C.P.; Gehring, C. Arbuscular mycorrhizal fungi along secondary forest succession at the eastern periphery of Amazonia: Seasonal variability and impacts of soil fertility. App. Soil Ecol. 2019, 136, 1–10. [Google Scholar] [CrossRef]
  25. Seligmann, H. Resource partition history and evolutionary specialization of subunits in complex systems. Biosystems 1999, 51, 31–39. [Google Scholar] [CrossRef]
  26. Hsiao, C.; Mohan, S.; Kalahar, B.K.; Williams, L.D. Peeling the onion: Ribosomes are ancient molecular fossils. Mol. Biol. Evol. 2009, 26, 2415–2425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Petrov, A.S.; Gulen, B.; Norris, A.M.; Kovacs, N.A.; Bernier, C.R.; Lanier, K.A.; Fox, G.E.; Harvey, S.C.; Wartell, R.M.; Hud, N.V.; et al. History of the ribosome and the origin of translation. Proc. Natl. Acad. Sci. USA 2015, 112, 15396–15401. [Google Scholar] [CrossRef] [Green Version]
  28. Gulen, B.; Petrov, A.S.; Okafor, C.D.; Vander Wood, D.; O’Neill, E.B.; Hud, N.V.; Williams, L.D. Ribosomal small subunit domains radiate from a central core. Sci. Rep. 2016, 6, 20885. [Google Scholar] [CrossRef]
  29. Lanier, K.A.; Athavale, S.S.; Petrov, A.S.; Wartell, R.; Williams, L.D. Imprint of ancient evolution on rRNA folding. Biochemistry 2016, 55, 4603–4613. [Google Scholar] [CrossRef]
  30. Demongeot, J.; Seligmann, H. Accretion history of large ribosomal subunits deduced from theoretical minimal RNA rings is congruent with histories derived from phylogenetic and structural methods. Gene 2020, 738, 144436. [Google Scholar] [CrossRef]
  31. Demongeot, J.; Seligmann, H. Comparisons between small ribosomal RNA and theoretical minimal RNA ring secondary structures confirm phylogenetic and structural accretion histories. Sci. Rep. 2020, 10, 7693. [Google Scholar] [CrossRef]
  32. Lovtrup, S. On von Baerian and Haeckelian recapitulation. Syst. Zool. 1978, 27, 348–352. [Google Scholar] [CrossRef]
  33. Kalinka, A.T.; Tomancak, P. The evolution of early animal embryos: Conservation or divergence. Trends Ecol. Evol. 2012, 27, 385–393. [Google Scholar] [CrossRef]
  34. Seligmann, H. Evolution and ecology of developmental processes and of the resulting morphology: Directional asymmetry in hindlimbs of Agamidae and Lacertidae (Reptilia: Lacertilia). Biol. J. Linn. Soc. 2000, 69, 461–481. [Google Scholar] [CrossRef]
  35. Johnson, D.B.F.; Wang, L. Imprints of the genetic code in the ribosome. Proc. Natl. Acad. Sci. USA 2010, 107, 8298–8303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Seligmann, H. Protein sequences recapitulate genetic code evolution. Comput. Struct. Biotechnol. J. 2018, 16, 177–189. [Google Scholar] [CrossRef] [PubMed]
  37. Gaudart, J.; Landier, J.; Huiart, L.; Legendre, E.; Lehot, L.; Bendiane, M.K.; Chiche, L.; Petitjean, A.; Mosnier, E.; Kirakoya-Samadoulougou, F.; et al. Factors associated with spatial heterogeneity of Covid-19 in France: A nationwide ecological study. Lancet Public Health 2021, 6, e89. [Google Scholar] [CrossRef]
  38. Lau, H.; Khosrawipour, V.; Kocbach, P.; Mikolajczyk, A.; Schubert, J.; Bania, J.; Khosrawipour, T. The positive impact of lockdown in Wuhan on containing the COVID-19 outbreak in China. J. Travel. Med. 2020, 27, taaa037. [Google Scholar] [CrossRef] [Green Version]
  39. Bendavid, E.; Oh, C.; Bhattacharya, J.; Ioannidis, J.P.A. Assessing mandatory stay-at-home and business closure effects on the spread of COVID-19. Eur. J. Clin. Investig. 2021, e13484. [Google Scholar]
  40. Bonacini, L.; Gallo, G.; Patriarca, F. Identifying policy challenges of COVID-19 in hardly reliable data and judging the success of lockdown measures. J. Popul. Econ. 2021, 34, 275–301. [Google Scholar] [CrossRef]
  41. Oraby, T.; Tyshenko, M.G.; Maldonado, J.C.; Vatcheva, K.; Elsaadany, S.; Alali, W.G.; Longenecker, J.C.; Al-Zoughool, M. Modeling the effect of lockdown timing as a COVID-19 control measure in countries with differing social contacts. Sci. Rep. 2021, 11, 3354. [Google Scholar] [CrossRef]
  42. Timelli, L.; Girardi, E. Effect of timing of implementation of containment measures on Covid-19 epidemic. The case of the first wave in Italy. PLoS ONE 2021, 16, e0245656. [Google Scholar] [CrossRef] [PubMed]
  43. Demongeot, J.; Seligmann, H. Covid-19 and miRNA-like inhibition power. Med. Hypotheses 2020, 144C, 110245. [Google Scholar] [CrossRef] [PubMed]
  44. Demongeot, J.; Griette, Q.; Magal, P. Computations of the transmission rates in SI epidemic model applied to COVID-19 data in mainland China. R. Soc. Open Sci. 2020, 7, 201878. [Google Scholar] [CrossRef] [PubMed]
  45. Soubeyrand, S.; Demongeot, J.; Roques, L. Towards unified and real-time analyses of outbreaks at country-level during pandemics. One Health 2020, 11, 100187. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Numbers of genomic variants of SARS-COVID-19 as a function of days since 1 January 2020 until 100 confirmed cases were detected in a country.
Figure 1. Numbers of genomic variants of SARS-COVID-19 as a function of days since 1 January 2020 until 100 confirmed cases were detected in a country.
Biology 10 00623 g001
Figure 2. Spread parameters (slope of exponential regression) of first (A) and second (B) waves in USA states as a function of population density. Trend inverted between first and second wave. A: the New England states with high densities (filled circles) follow a different pattern than the rest of states, with a slower increase in spread with increasing density. B: the second wave slopes (as determined in August 2020) are overall lower than first wave slopes, and counter-intuitively decrease with population density. Patterns remain statistically significant with nonparametric Spearman rank correlation tests, see text.
Figure 2. Spread parameters (slope of exponential regression) of first (A) and second (B) waves in USA states as a function of population density. Trend inverted between first and second wave. A: the New England states with high densities (filled circles) follow a different pattern than the rest of states, with a slower increase in spread with increasing density. B: the second wave slopes (as determined in August 2020) are overall lower than first wave slopes, and counter-intuitively decrease with population density. Patterns remain statistically significant with nonparametric Spearman rank correlation tests, see text.
Biology 10 00623 g002
Figure 3. Pearson correlation coefficients r between spread rates from Table 3 and Table 4 (open cirles, USA; and +, other countries, r from Table 5) as a function of numbers of days between the time frames when spread rates were estimated. Blue lines indicate crirical values of r for p < 0.05 according to two-tailed tests. Most positive r’s with p < 0.05 are for pairs of consecutive periods (20 days on x axis, green box), most negative r’s with p < 0.05 are for periods separated by approximately 80 days, orange box.3.5.
Figure 3. Pearson correlation coefficients r between spread rates from Table 3 and Table 4 (open cirles, USA; and +, other countries, r from Table 5) as a function of numbers of days between the time frames when spread rates were estimated. Blue lines indicate crirical values of r for p < 0.05 according to two-tailed tests. Most positive r’s with p < 0.05 are for pairs of consecutive periods (20 days on x axis, green box), most negative r’s with p < 0.05 are for periods separated by approximately 80 days, orange box.3.5.
Biology 10 00623 g003
Figure 4. Daily new confirmed COVID-19 cases in the UK since 15 February 2020 until the end of 2020. The dates, with months counted in roman numbers are indicated for specific datapoints.
Figure 4. Daily new confirmed COVID-19 cases in the UK since 15 February 2020 until the end of 2020. The dates, with months counted in roman numbers are indicated for specific datapoints.
Biology 10 00623 g004
Table 1. Web sites for data sources, accessed on 20 August 2020.
Table 2. Spread parameters of first and second waves in the 50 states of the USA and the district of Columbia.
Table 2. Spread parameters of first and second waves in the 50 states of the USA and the district of Columbia.
StateTEDADate 1Slope 1Date 2DisSlope 2
Alabama17.11503739.221.030.3338
Alaska−3580134.616.030.17523.0570.051
Arizona15.712502337.922.030.184125.0560.0733
Arkansas15.82002238.319.030.096710.0580.0502
California15.28809536.814.030.166210.05−40.0223
Colorado7.320702036.920.030.198911.0620.0408
Connecticut9.41502864118.030.2759
DC 50 3421.030.101130.0660.0719
Delaware12.92018740.722.030.14636.06−60.0412
Florida21.510514542.215.030.258426.05−110.0785
Georgia17.51806836.915.030.253625.0500.0293
Hawaii21.19208639.223.030.14332.05−210.0906
Idaho6.91520736.624.030.27935.06−30.0991
Illinois111808938.316.030.430815.0610.0289
Indiana10.92867137.924.030.109322.0610.0348
Iowa8.83402138.224.030.101716.0420.0893
Kansas12.46101436.923.030.087612.0410.096
Kentucky13.12304338.923.030.174310.6−70.0344
Louisiana19.1304137.715.030.202223.05−60.0382
Maine51801644.921.030.109218.0400.0373
Maryland12.311023838.825.030.121614.0400.0127
Massachusetts8.815033639.413.030.301515.06−330.0129
Michigan6.92706739.818.030.417315.0620.0612
Minnesota5.13702638.120.030.119414.04−60.1172
Mississippi17.4902437.720.030.0863
Missouri12.52403438.721.030.176229.05140.0307
Montana5.91040239.914.030.24855.0640.0965
North Carolina152107938.919.030.14584.0500.0272
North Dakota4.7580435.227.030.05458.0570.0914
Nebraska9.3790936.629.030.077519.06−70.0248
Nevada9.916801038.118.030.199223.06−20.054
New Hampshire6.6300574321.030.1426
New Jersey11.5804704016.030.3797
New Mexico11.917401738.124.030.11269.06−50.0347
New York7.43001623910.030.5503
Ohio10.426010939.419.030.236915.0410.2485
Oklahoma15.34002236.724.030.12678.0620.1285
Oregon9.110101639.419.030.10617.05−90.0564
Pennsylvania9.334011040.817.030.293514.0610.0472
Rhode Island10.16039440.121.030.14624.0510.0435
South Carolina16.91106239.618.030.153624.0500.0609
South Dakota7.3670437.129.030.1716
Tennessee14.22706138.817.030.2258
Texas18.25204034.817.030.24120.410.0281
Utah9.21860143119.030.114611.05−60.0329
Vermont6.13002642.822.030.0546
Virginia12.82908138.419.030.178915.0600.0173
West Virginia114602942.723.030.09855.0610.0565
Washington9.15204137.214.030.065714.0400.0494
Wisconsin6.23204139.618.030.091012.0400.0200
Wyoming5.6204023828.030.234629.0500.0646
P K-S0.9070.0080.0070.530 0.274 0.228
Columns are: 1. State, 2. Mean annual temperature, 3. Mean elevation, 4. Population density, 5. Median age, 6. Onset date of first wave, 7. Exponential slope of first wave, 8. Onset date of second wave, 9. Days separating onset date in column 8 vs. onset date determined by statistical analysis. P K-S is p-values along the Kolmogorov–Smirnov test for normality.
Table 3. Spread rates (×100) for COVID-19 in the 50 states of the USA and the district of Columbia for 14 consecutive nonoverlapping 20-day time periods. N > 0 are numbers of positive rates, Mean indicates mean rates. Days are since pandemic onset in that state, which differs among states. P K-S are p-values along the Kolmogorov–Smirnov test for normality.
Table 3. Spread rates (×100) for COVID-19 in the 50 states of the USA and the district of Columbia for 14 consecutive nonoverlapping 20-day time periods. N > 0 are numbers of positive rates, Mean indicates mean rates. Days are since pandemic onset in that state, which differs among states. P K-S are p-values along the Kolmogorov–Smirnov test for normality.
State\Day20406080100120140160180200220240260280
Alabama24.820.512.674.616.495.17−0.59−4.30−0.262.370.881.711.072.13
Alaska14.70−8.70−2.7614.640.526.415.180.600.351.642.813.030.29−4.71
Arizona28.61−0.072.091.798.311.97−0.58−6.02−3.773.924.054.036.062.25
Arkansas23.241.66−4.1712.993.53−0.570.39−1.081.241.19−0.71−0.150.721.30
California26.210.060.552.831.543.880.300.94−1.79−0.50−0.063.147.856.02
Colorado23.83−0.170.42−2.01−4.762.251.56−2.610.595.252.595.331.36−2.47
Connecticut32.483.46−3.92−3.20−4.993.623.101.255.22−1.332.864.440.24−0.64
DC18.431.182.00−2.08−4.002.501.52−1.83−1.700.984.374.062.961.71
Delaware28.25−4.309.34−4.42−2.685.02−0.380.44−0.161.332.042.785.103.54
Florida22.08−2.26−0.88−0.087.875.21−1.21−3.72−1.810.111.263.722.181.86
Georgia25.76−0.580.030.602.285.680.41−1.63−1.06−1.14−0.251.872.484.07
Hawaii22.85−7.89−9.60−3.8416.155.103.196.530.84−3.67−2.120.70−0.431.59
Idaho32.63−9.28−5.8810.374.656.68−0.40−1.76−1.814.391.931.501.40−0.50
Illinois23.161.722.51−1.82−5.171.471.680.610.530.472.264.61−0.67−0.52
Indiana34.250.060.50−0.26−0.383.051.49−0.260.080.833.394.901.090.23
Iowa26.936.861.39−2.82−0.964.401.350.95−1.622.782.556.94−4.71−1.99
Kansas24.211.493.91−4.313.205.06−3.57−2.44−2.711.112.230.112.03−1.89
Kentucky23.894.070.83−1.77−3.163.472.67−0.021.911.600.862.502.570.39
Louisiana29.89−7.88−2.11−0.572.969.530.65−2.440.31−1.71−4.58−0.964.921.24
Maine22.35−0.725.222.02−3.99−2.011.06−1.420.172.290.417.931.975.99
Maryland30.025.992.67−1.11−5.580.263.50−1.580.94−1.871.252.113.171.53
Massachusetts33.102.80−2.26−3.74−8.440.532.26−1.500.581.501.083.731.752.75
Michigan37.94−4.00−2.67−2.12−4.923.781.31−0.700.280.302.465.67−0.76−5.05
Minnesota19.014.6010.340.17−2.172.230.68−0.942.252.731.345.571.74−4.48
Mississippi32.254.000.992.180.034.083.05−2.570.33−1.730.05−2.050.651.93
Missouri27.29−0.48−2.610.990.452.813.594.79−0.071.111.812.43−2.17−2.36
Montana23.86−6.12−12.880.371.806.772.181.473.112.095.241.532.020.28
Nebraska19.146.554.390.15−3.171.242.96−1.57−0.873.003.525.69−3.28−4.55
Nevada26.55−0.48−1.131.493.875.541.94−2.36−3.26−0.860.931.233.481.25
New Hampshire21.280.151.25−1.06−5.27−0.77−0.26−2.452.87−1.763.974.103.682.58
New Jersey37.460.15−3.57−2.64−4.82−0.340.721.421.392.893.094.412.640.98
New Mexico14.282.461.26−1.04−1.643.721.49−3.76−1.334.515.042.873.75−1.11
New York24.05−2.39−3.75−2.75−3.60−1.07−1.390.342.43−0.172.132.712.961.72
North Carolina29.431.702.242.552.251.67−0.78−4.92−0.541.871.861.483.163.14
North Dakota16.306.56−1.08−3.77−1.376.401.570.541.574.273.214.11−1.41−4.33
Ohio28.975.10−1.13−0.190.402.70−0.40−1.752.29−1.132.333.813.651.40
Oklahoma28.27−3.560.02−1.777.423.312.65−2.220.922.810.651.260.550.97
Oregon20.27−1.791.86−4.647.452.490.80−0.93−0.161.030.643.442.56−0.62
Pennsylvania36.040.850.71−2.20−2.913.201.79−0.871.900.382.643.763.932.27
Rhode Island22.393.75−2.60−4.36−4.45−0.314.360.15−2.402.904.503.311.90−3.86
South Carolina25.64−1.690.971.625.172.57−1.47−3.292.79−0.94−0.610.542.354.14
South Dakota28.015.633.39−1.53−1.630.913.101.720.785.425.222.97−1.94−2.23
Tennessee26.846.020.080.871.115.392.35−1.620.780.352.850.38−0.823.97
Texas27.730.052.81−0.584.554.74−1.22−2.03−2.721.572.344.111.293.50
Utah27.540.260.870.111.530.46−1.55−1.201.684.670.971.450.94−0.63
Vermont22.07−8.53−0.138.59−2.92−2.88−3.218.341.07−1.869.304.999.213.59
Virginia23.593.992.081.06−4.810.551.80−0.061.31−1.401.701.573.723.33
W Virginia26.40−2.56−1.811.894.968.561.410.193.711.373.184.711.54−0.02
Washington8.22−6.721.030.590.893.150.85−1.31−2.752.272.454.662.430.46
Wisconsin23.990.014.044.10−0.823.900.40−0.821.674.200.652.05−3.06−1.76
Wyoming17.98−3.488.950.179.483.003.04−4.67−2.255.933.273.861.90−3.35
N > 05129322426443716313745484131
Mean25.38−0.0380.3230.3160.3973.0691.086−0.8310.2531.3212.0753.0331.7650.491
P K-S0.7760.3990.5290.1130.7540.9890.8830.5450.8670.9200.7940.7900.6540.670
Table 4. Spread rates (×100) for COVID-19 in the 51 countries for 15 consecutive nonoverlapping 20-day time periods. N > 0, numbers of positive rates, Mean is mean rates. Days are since pandemic onset in that state, which differs among states. P K-S are p values along the Kolmogorov–Smirnov test for normality.
Table 4. Spread rates (×100) for COVID-19 in the 51 countries for 15 consecutive nonoverlapping 20-day time periods. N > 0, numbers of positive rates, Mean is mean rates. Days are since pandemic onset in that state, which differs among states. P K-S are p values along the Kolmogorov–Smirnov test for normality.
Country\Day20406080100120140160180200220240260280300
Afghanistan15.285.007.265.32−1.29−2.26−0.48−6.00−6.08−1.2211.191.604.800.36
Albania14.49−0.84−5.79−2.205.020.971.301.88−1.001.650.892.664.131.11
Armenia16.65−0.795.727.072.12−0.68−4.33−0.870.481.453.921.69−4.83−3.31−5.01
Australia25.07−4.71−11.82−0.040.8311.069.072.56−7.88−4.05−7.59−0.08−3.59−1.247.08
Austria30.234.16−10.29−4.76−6.941.505.692.135.720.122.434.227.34−2.42−2.19
Bangladesh5.2624.624.003.712.10−0.57−0.491.25−2.11−1.020.531.102.44−1.72
Belgium29.945.88−2.79−3.27−5.26−1.731.746.97−0.492.053.529.38−8.02−4.66−0.16
Bolivia14.657.188.705.721.632.121.100.23−2.76−0.96−2.08−1.721.906.46
Bosnia18.434.10−0.38−5.195.366.150.900.321.22−0.42−2.277.65−2.09−2.28−3.57
Brazil33.738.575.186.822.172.05−0.931.46−0.73−2.52−0.860.531.652.280.04
Canada25.0414.581.87−1.82−2.58−4.040.101.970.882.094.291.542.681.281.28
Chile31.066.441.735.411.47−1.95−2.31−0.97−1.280.701.060.002.050.881.54
Czechia30.072.34−6.20−1.42−2.702.91−4.901.870.348.702.598.27−2.65−1.155.81
Georgia7.565.78−0.160.082.65−6.343.13−6.75−2.3115.536.308.634.030.63−6.27
Germany35.167.43−5.98−5.09−5.063.36−2.524.944.69−2.042.076.734.300.303.30
Greece18.28−1.52−3.81−2.056.380.931.2512.281.550.620.305.472.79−3.74−4.50
Hungary19.486.03−0.99−3.38−7.37−1.736.181.612.7310.631.614.883.251.92−4.47
India19.7814.274.255.483.253.713.242.230.401.75−0.67−1.730.01−0.20−2.43
Indonesia31.306.620.213.573.001.290.540.852.181.26−0.09−0.512.311.571.37
Iran40.406.66−4.41−0.191.85−0.870.100.02−0.48−1.722.881.144.162.47−2.79
Iraq14.828.94−5.013.219.963.941.202.043.050.49−0.95−1.32−0.96−2.68−2.93
Ireland30.396.56−5.86−7.61−10.15−4.524.669.261.824.313.793.43−3.89−1.988.65
Israel28.6112.23−2.47−11.5312.255.227.81−0.061.183.821.80−10.380.551.703.79
Italy35.436.45−1.99−5.16−5.06−2.65−2.391.624.416.080.357.276.76−0.68−2.62
Kazakhstan21.105.980.556.71−0.228.67−0.94−6.96−7.08−1.472.079.071.27−0.50
Kuwait4.176.193.047.45−0.56−2.010.330.151.490.99−1.782.36−0.29−4.74−1.18
Lithuania21.70−3.11−9.79−2.26−0.202.466.224.34−0.175.181.448.250.671.40
Luxemburg29.22−8.42−6.55−6.810.4614.025.36−5.470.302.604.257.57−2.280.16
Malaysia31.131.47−6.38−2.084.30−3.722.533.783.426.774.925.480.821.232.00
Mexico26.388.114.823.361.520.790.89−0.240.39−1.991.880.09−0.831.380.45
Moldova19.261.800.031.234.22−4.230.681.600.83−0.640.09−0.301.41−0.78
Morocco25.026.01−1.31−5.954.986.504.173.80−0.232.652.040.16−1.43−1.66
N Macedonia19.605.10−8.861.8010.530.380.31−0.69−4.683.925.905.551.80−3.23
Netherlands33.966.11−1.67−4.50−2.12−6.07−1.058.251.242.966.154.91−3.62−0.795.22
Norway32.243.04−4.34−4.61−5.32−4.10−8.3410.912.805.89−0.321.045.25−2.751.26
Pakistan32.055.036.722.797.68−2.83−5.12−5.37−3.392.12−0.643.034.910.68−3.49
Panama25.652.630.321.014.091.48−1.03−0.600.20−0.351.031.064.053.421.44
Philippines23.021.90−0.50−0.893.503.592.595.20−1.100.90−1.60−2.30−1.77−0.620.97
Poland32.906.28−0.371.451.00−2.380.934.221.28−3.036.777.001.16−4.49−0.04
Portugal32.573.86−4.33−1.822.481.430.19−2.99−0.434.721.367.843.59−4.45−1.16
Romania25.793.61−0.81−4.610.441.885.160.03−0.501.013.582.861.94−1.54−3.23
Russia24.4915.795.52−0.90−0.10−1.62−0.46−0.76−0.460.813.511.891.351.180.40
S Korea39.90−5.49−8.29−5.987.203.220.91−2.40−3.5613.72−5.90−0.991.885.773.62
Saudi Arabia26.616.476.993.681.52−0.50−2.70−3.71−1.80−2.88−2.03−0.83−1.19−1.990.01
Serbia24.038.49−5.76−1.455.409.671.79−3.79−6.400.434.4914.247.720.85−3.55
Turkey42.423.06−4.63−2.493.35−1.02−0.541.671.820.180.311.416.301.44
UAE5.1216.421.892.99−2.70−0.51−2.32−1.133.674.102.321.31−0.340.030.45
UK35.6614.620.24−1.45−4.09−2.13−3.280.732.284.964.963.900.68−3.834.09
Ukraine27.0410.09−0.42−0.493.69−0.962.532.751.241.992.883.102.34−1.29
USA35.9324.230.89−0.29−0.670.173.350.60−0.80−1.640.540.793.632.381.99
Uzbekistan14.331.901.023.814.714.621.32−0.89−3.712.43−2.68−2.47−1.171.02
N > 051.0044.0021.0021.0033.0027.0033.0033.0027.0036.0037.0039.0035.0025.0021.00
Mean25.156.02−1.20−0.351.350.880.851.06−0.152.031.582.871.31−0.330.14
P K-S0.7300.1090.9270.8710.8800.6570.5320.3240.4300.2240.9940.5990.9330.9620.868
Table 5. Pearson correlation coefficients r (×100) between spread rates. Above the diagonal, consecutive nonoverlapping 20-day periods from the USA (Table 3), and, below the diagonal, from various countries (Table 4). Highlights are for p < 0.05, two-tailed test.
Table 5. Pearson correlation coefficients r (×100) between spread rates. Above the diagonal, consecutive nonoverlapping 20-day periods from the USA (Table 3), and, below the diagonal, from various countries (Table 4). Highlights are for p < 0.05, two-tailed test.
Period20406080100120140160180200220240260280
20 8.44−20.28−14.06−12.912.06−5.351.0819.97−20.66−16.90−13.991.2522.38
40−12.19 34.97−35.06−41.88−27.1619.09−20.758.9511.979.1612.29−32.21−10.46
60−30.7046.88 −13.30−14.46−21.45−17.35−37.87−25.2427.824.1526.255.78−1.48
80−43.4819.6060.91 18.563.12−7.800.98−3.242.781.03−19.449.147.52
100−16.05−12.482.7414.58 45.82−16.25−15.34−29.57−0.86−34.56−38.750.6710.38
120−3.47−28.73−29.64−9.7130.60 10.69−14.23−11.02−2.38−31.61−31.53−16.50−14.27
140−10.69−13.90−38.06−35.5210.3339.24 16.8610.116.74−2.213.73−42.42−38.85
16018.81−3.84−31.98−36.01−28.21−25.494.07 33.33−29.2722.7113.97−7.00−7.52
18019.5313.42−11.88−35.10−34.54−36.00−9.6449.67 −24.47−2.14−2.01−10.990.32
200−1.15−18.88−28.08−37.32−7.39−25.391.72−4.1115.19 30.5428.37−32.57−49.68
2200.6517.936.56−8.7−23.01−35.36−7.58−8.6710.0011.19 46.6410.23−23.15
2405.33−19.27−32.94−13.92−32.562.86−9.33−4.04−3.3319.6740.71 −8.76−22.56
26015.6714.91−4.48−3.114.91−15.03−14.76−22.094.903.989.6610.70 56.77
28014.850.8613.41.4816.4711.0412.40−29.26−19.9811.38−20.61−31.8131.23
30043.392.23−32.48−29.66−26.214.744.0533.80−4.15−4.35−10.73−26.07−35.8412.88
Table 6. Pearson correlation coefficients r (×100) between spread rates for consecutive nonoverlapping 20-day periods from Table 3 and Table 4 with environmental variables from Table 2 and Table 6. Rows starting with USA give r between spread rates and environmental factors for the US states, rows starting with world are for countries. Rows starting with rs present nonparametric Spearman rank correlation coefficients rs for the two variables, altitude and density, that differ at p < 0.05 from normality according to Kolmogorov–Sminov tests. Bold, underlined indicate that correlations are statistically significant at p < 0.05 according to a two- and one-tailed test, respectively.
Table 6. Pearson correlation coefficients r (×100) between spread rates for consecutive nonoverlapping 20-day periods from Table 3 and Table 4 with environmental variables from Table 2 and Table 6. Rows starting with USA give r between spread rates and environmental factors for the US states, rows starting with world are for countries. Rows starting with rs present nonparametric Spearman rank correlation coefficients rs for the two variables, altitude and density, that differ at p < 0.05 from normality according to Kolmogorov–Sminov tests. Bold, underlined indicate that correlations are statistically significant at p < 0.05 according to a two- and one-tailed test, respectively.
Variable\Day20406080100120140160180200220240260280
Temp. USA25.124.70−7.16−16.2147.0621.18−21.98−19.91−21.60−40.10−50.69−49.4616.2951.98
World−19.044.5123.0334.3131.6615.947.09−6.62−5.28−23.11−38.76−38.72−10.6813.34
Altitude USA−23.09−18.78−1.666.5332.7111.811.28−19.04−29.2049.9814.481.343.90−28.26
rs−4.8−19.91.512.832.313.85−8.560.52−6.6942.0421.236.97−1.97−23.1
World−11.64−10.6129.9233.4618.73−3.51−13.12−35.17−30.11−5.927.80−24.3412.0822.81
rs−4.9−8.922.8819.7927.616.34−5.22−36.68−32.67−13.16−1.60−21.7526.3625.92
Density USA38.8517.62−17.62−36.69−41.00−31.0314.3213.5614.52−18.433.4110.1614.5812.67
rs11.5426.24−10.30−36.64−34.65−5.497.6116.8612.70−38.32−9.8610.748.7720.96
World−19.9024.746.45−4.704.56−4.781.6210.297.0013.47−13.05−10.58−11.70−14.60
rs8.63−5.25−20.88−25.967.270.914.3319.3532.8336.690.2212.44−10.41−23.4
Age USA15.17−9.38−5.15−7.11−13.10−21.44−7.1317.2532.06−33.3111.6825.9121.8430.23
World33.71−16.10−50.50−54.53−39.52−2.629.0830.1425.7429.1012.2757.837.62−22.37
Table 7. Environmental variables for 51 countries, mean annual temperature, mean elevation, population density per square kilometer, and median age. P K-S is the p-value along the Kolmogorov–Smirnov test for normality.
Table 7. Environmental variables for 51 countries, mean annual temperature, mean elevation, population density per square kilometer, and median age. P K-S is the p-value along the Kolmogorov–Smirnov test for normality.
CountryTemperatureElevationDensityAge
Afghanistan12.618854918.9
Albania11.470899.7332.9
Armenia7.1517929935.1
Australia21.65330338.7
Austria6.3591010644
Bangladesh2585117526.7
Belgium9.5518137841.4
Bolivia21.5511921024.3
Bosnia9.855006942.1
Brazil24.953202532.6
Canada−5.35487442.2
Chile8.4518712334.4
Czech7.5543313542.1
Georgia5.8143253.5138.1
Germany8.526323347.1
Greece15.449821044.5
Hungary9.7514327242.3
India23.6562141228.1
Indonesia25.8536714130.2
Iran17.2513055130.3
Iraq14.033129020
Ireland9.31187036.8
Israel19.250841729.9
Italy12.4553820045.5
Kazakhstan6.4387730.6
Kuwait25.3510864229.3
Lithuania6.21107343.7
Luxemburg8.65325237.3939.3
Malaysia25.45389928.5
Mexico2111116428.3
Moldova9.4513979.2736.7
Morocco17.19098029.3
N Macedonia9.87418137.9
Netherlands9.253042142.6
Norway1.54601739.2
Pakistan20.290027423.8
Panama25.43605629.2
Philippines25.8544236223.5
Poland7.8517312340.7
Portugal15.1537211242.2
Romania8.841481.441.1
Russia−5.1600939.6
S Korea11.528251741.8
Saudi Arabia24.656651627.5
Serbia10.5547389.0842.6
Turkey11.11132106.1230.9
UAE27149116.8730.3
UK8.4516228040.5
Ukraine8.317569.4940.6
USA8.557603438.1
Uzbekistan12.0527507328.6
P K-S0.1990.0240.0020.259
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Seligmann, H.; Vuillerme, N.; Demongeot, J. Unpredictable, Counter-Intuitive Geoclimatic and Demographic Correlations of COVID-19 Spread Rates. Biology 2021, 10, 623. https://doi.org/10.3390/biology10070623

AMA Style

Seligmann H, Vuillerme N, Demongeot J. Unpredictable, Counter-Intuitive Geoclimatic and Demographic Correlations of COVID-19 Spread Rates. Biology. 2021; 10(7):623. https://doi.org/10.3390/biology10070623

Chicago/Turabian Style

Seligmann, Hervé, Nicolas Vuillerme, and Jacques Demongeot. 2021. "Unpredictable, Counter-Intuitive Geoclimatic and Demographic Correlations of COVID-19 Spread Rates" Biology 10, no. 7: 623. https://doi.org/10.3390/biology10070623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop