 Primary research
 Open access
 Published:
An application of slow feature analysis to the genetic sequences of coronaviruses and influenza viruses
Human Genomics volume 15, Article number: 26 (2021)
Abstract
Background
Mathematical approaches have been for decades used to probe the structure of nucleotide sequences. This has led to the development of Bioinformatics. In this exploratory work, a novel mathematical method is applied to probe the genetic structure of two related viral families: those of coronaviruses and those of influenza viruses. The coronaviruses are SARSCoV2, SARSCoV1, and MERS. The influenza viruses include H1N11918, H1N12009, H2N21957, and H3N21968.
Methods
The mathematical method used is the slow feature analysis (SFA), a rather new but promising method to delineate complex structure in nucleotide sequences.
Results
The analysis indicates that the nucleotide sequences exhibit an elaborate and convoluted structure akin to complex networks. We define a measure of complexity and show that each nucleotide sequence exhibits a certain degree of complexity within itself, while at the same time there exists complex interrelationships between the sequences within a family and between the two families. From these relationships, we find evidence, especially for the coronavirus family, that increasing complexity in a sequence is associated with higher transmission rate but with lower mortality.
Conclusions
The complexity measure defined here may hold a promise and could become a useful tool in the prediction of transmission and mortality rates in future new viral strains.
Background
Since the early 1970s, scientists have attempted to discover some kind of order or hidden structures in nucleotide sequences. With the advent of sequencing techniques in the late 1970s, scientists had the opportunity to probe nucleic acid sequences for such order [1,2,3]. Soon, mathematical approaches were employed to shed light in this endeavor, leading to the fullblown field of Bioinformatics [4,5,6,7]. We report, for the first time, the application of slow feature analysis (SFA) to genetic sequences. SFA is a procedure for extracting slowly varying, driving signals from a given nonstationary time series and is used here to delineate signals or structure in nucleotide sequences, which would not otherwise be detected. Descriptions of this procedure, which have been successfully applied in many scientific areas, have been reported previously in detail [8,9,10].
Methods
SFA is an approach that is designed to optimally identify lowfrequency behavior in a time series, thereby delineating its complex structure more effectively. This analysis is rooted, theoretically, in the timeembedding theorems. In this method, a onedimensional time series is embedded in a multidimensional space consisting of the original time series and lagged copies thereof. SFA further uses a nonlinear expansion to map this multidimensional input signal onto an even larger feature space, and then solves a linear problem to find a linear combination of featurespace variables that minimizes their time derivative (rate of change) [11]. The objective of SFA is to find the optimally filtered signals that vary as slowly as possible but still carry significant information. To ensure this, the output signals need to be uncorrelated and have unit variance [12]. This approach has been successfully applied in many areas, including climate science [13, 14].
In mathematical terms [8], the goal of SFA is, given an ndimensional input signal x(t), to find a set of realvalued inputoutput functions g_{j}(x) such that the output signals
minimize \( \varDelta \left({y}_i\right):= <{\dot{y}}_j^2{>}_t \)
under the constraints
with <∙>_{t} and \( \dot{y} \) indicating temporal averaging and the derivative of y, respectively.
The Δvalue is a measure of the temporal slowness of the signal y(t). It is given by the mean square of the signal’s time derivative. Small Δvalues correspond to slowly varying signals. The first two constraints avoid the trivial constant solution, while the last constraint guarantees that the output functions g_{j} are distinct and hence extract different information from the input signal. For a tutorial on this method, the reader could consult reference [8] or a more recent presentation in [15]. In that tutorial, a simple example of a twodimensional input signal x_{1}(t)=sin(t)+cos(11t)^{2} and x_{2}(t)=cos(11t) is considered. Both components are quickly varying, but hidden in the signal is the slowly varying “feature” y(t)=x_{1}(t)−x_{2}(t)^{2}=sin(t), which can be extracted with a polynomial of degree two, namely h(x)=x_{1}−x_{2}^{2}.
In the situation with one observable (time series of some variable) from an unknown system where the actual state space is not known (as is the case here), embedding is necessary (and essential) to delineate the underlying dynamics much like in attractor reconstructions. The SFA algorithm can be summarized as follows. Consider a time series \( {\left\{x(t)\right\}}_{t={t}_1,\dots, {t}_n} \), where t denotes time and n indicates the length of the time series. First, we embed {x(t)} into an mdimensional state space using timedelayed copies of x(t):
where x_{1}(t) = x(t); x_{2}(t) = x_{1}(t − τ); x_{3}(t) = x_{1}(t − 2τ), and so on, τ is the delay, and N = n – m + 1. Then, nonlinear expansions (usually secondorder polynomials) are used to generate a kdimensional function state space:
which can also be written as \( \mathbf{H}(t)={\left\{{h}_1(t),{h}_2(t),\dots, {h}_k(t)\right\}}_{t={t}_1,\dots, {t}_N} \), where
The expanded signal H(t) is then centered and normalized to zero mean and unit variance. This process is referred to as whitening or sphering. Thus, we have
where
\( {h}_j^{\prime }(t)=\left[{h}_j(t)\overline{h_j}\right]/S \), and \( S=\frac{1}{k}\sqrt{\sum_{j=1}^k{\left({h}_j(t)\overline{h}\right)}^2} \).
Using the Schmidt algorithm, H^{′}(t) is orthogonized into:
where the transformed signal matrix Z is column orthogonal:
The final step of SFA is to find the set of coefficients (a_{1}, a_{2}, …, a_{k}) such that the time series
varies as slowly as possible. This set is given by the eigenvector W_{1} of the timederivative covariance matrix
corresponding to the smallest eigenvalue λ_{1}. Here
and
Using W_{1}, the optimally filtered slowfeature signal (also known as a driving force factor, which can be composed of one or more components) can be written as:
where r and c are constants derived to best match y(t) and the original time series x(t).
Once the optimally filtered (lowfrequency) SFA signal has been identified, its significant periodicities can be found from the timeaveraged wavelet power spectrum. Wavelet analysis has been widely used to analyze localized structures and spectral properties of time series. For example, [16] provides a detailed description of the wavelet analysis, along with a very useful toolkit to conduct stepbystep wavelet analysis, including a statistical significance test based on the rednoise surrogate data (see http://paos.colorado.edu/research/wavelets/). We here used the Morlet wavelet with the wavenumber set to 4 to match the smoothness of the SFAderived slowfeature signal, focusing, once again, on the spectral peaks statistically significant at the 5% level. Note also that SFA is applicable to nonstationary data, so no data preprocessing is required.
The combination of the SFA and wavelet analyses we use in the present study has been shown to be more effective in diagnosing lowfrequency periodicities in data sets of a limited length than direct spectral analysis methods. Note that the driving force may not necessarily consist of just one component, but several components, which, as we will see below, correspond to forcings or signals at certain time scales. The success of SFA in delineating these slow signals lies in the fact that embedding the time series in high enough dimensions and the subsequent dynamical procedure removes the noise and smallscale features that may obscure or suppress those slow signals, thereby delineating more accurately the complex structure of a sequence.
Analysis and results
We first analyzed the nucleotide sequences from three viruses from the same family: SARSCoV2, SARSCoV1, and MERS. Those sequences are approximately 30,450 bases long and part of the now worldinfamous coronavirus family. Since a nucleotide sequence is a string of the bases A, T (U in RNA), C, and G, we first transformed it to a time series of integers in the interval [1–4] (i.e., A➔1, T/U ➔2, C ➔3, G ➔4). Here, we need to stress that a time series represents a particular type of process, where some quantity is sampled in time, t. A nucleotide sequence is a very similar object, but the “sampling” is over space. In a time series, we are interested in the dependency of observations at different time scales, whereas in nucleotide sequences, we are interested in dependencies in different space scales. As such, the mathematical tools to identify structures in time can in principle be applied to identify structure in space, as long as t is thought as a parameter identifying the scale. Transforming a nucleotide sequence into a time series has been used in the past to identify interesting properties in nucleotide sequences (such as the wellknown period 3; see [4, 5] and references therein). Note also that the above transformation of A➔1, T/U ➔2, C ➔3, and G ➔4 may, depending on frequency distribution of A, T/U, C, and G in the sequence, result in a nonstationary time series. However, unlike other spectral methods, SFA is not affected by nonstationarity in the data.
Once we have a time series, we apply SFA, and once we have the SFA signal (which as we mentioned above may be comprised of several components, see Eq. 1), we extract the SFA components by wavelet analysis. Figure 1 shows the SFA signal for SARSCoV2 virus for m=15 and τ=1. As explained above, this signal is normalized to zero mean and unit variance. Figure 2 shows the wavelet of the time series in Fig. 1. In order to extract the peak “periods” of the driving force signal, we used the Morlet wavelet to compute the timeaveraged power spectrum of the wavelet transform [16]. The black solid line in Fig. 3 is the timeaveraged power spectrum of the wavelet transform of the driving force, and the dashed line represents the 95% confidence level, estimated using AR1 surrogate data [16]. The dots show the periods of the oscillatory components of the driving force that are significant above the 95% level.
The significant peak periodicities for SARSCoV2 are as follows^{Footnote 1}:
Given the above periodicities recovered from SARSCoV2, we next construct Table 1, which shows the ratios between these peaks. We observe the following EXACT relations between peak periods:
And the following almost exact relationships based on the criterion:
Keeping only those relationships, we remain with Table 2, which could be thought as portraying the degree of structure or complexity in the SARSCoV2 sequence. We observe in the exact relationships multiples of a power of 2 and in the almost exact relationships multiples of 19 (152=2×76=8×19) and 83. Clearly, a sophisticated and rather convoluted structure, with numerous processes embedded in the sequence, is present. Keep in mind that the factors 19 and 83 (odd numbers) will appear in the rest of the sequences studied here. We define the number of entries above the diagonal in Table 2 as the degree complexity, C. In this case, C=13.
In the Supplementary material, Figures S1, S2, S3 are similar to Figs. 1, 2, and 3, and Tables ST1 and ST2 are similar to Tables 1 and 2 but for SARSCoV1. Figures S4, S5, and S6 are similar to Figs. 1, 2, and 3, and Tables ST3 and ST4 are similar to Tables 1 and 2 but for MERS.
According to Figure S3, the peak periods for SARSCoV1 are as follows:
and according to Tables ST1 and ST2, we now have five exact periodicities
and three almost exact
Again here, we observe in the exact relationships, multiples of a power of 2, and in the almost exact between P_{7} and P_{1}, P_{7} and P_{2}, and between P_{5} and P_{2}. Note again the multiples of 19 and 83. Note also that from the almost exact relationships, it follows that P_{7}/P_{5}=4, which is one of the exact relationships. Here, the degree of complexity is C=8.
According to Figure S6, the peak periods for MERS are as follows:
and according to Tables ST3 and ST4, we now have three exact periodicities
and three almost exact relationships
P_{9}, P_{8}, and P_{6} are multiples (again in a power of 2) of P_{1}, P_{2}, P_{4} (ordered in a bottomtop “symmetric” way), P_{6}, P_{7}, and P_{8} are multiples of P_{2}, P_{3}, and P_{1} but not of 2, but again of 19 and 83 (it is interesting to note that the odd multiples of 19 and 83 appear in all three sequences). Here, the degree of complexity is C=6.
By comparing Tables 2, ST2, and ST4 (and their associated C), one may argue that there is more embedded complexity and intricate patterning in SARSCoV2 than SARSCoV1 and MERS.
Other important relationships
Keeping in mind that all three sequences belong to the same coronavirus family, there are similarities and interrelationships between the sequences. For example, it is easy to observe that:
In general, SFA reveals a consistent picture between these sequences with very intricate structure with details at many scales, indicating very elaborate and sophisticated embedded processes, with complexity increasing from MERS to SARSCoV1 to SARSCoV2.
Extension of the analysis to the influenza viruses of H1N11918, H1N12009, H2N21957, and H3N21968
In an effort to provide further support for the efficiency and consistency of SFA in the analysis of nucleotide sequences, we consider four other viral sequences from a different viral family, that of the influenza viruses or the Orthomyxoviridae family [17].
In the Supplementary material, Figures S7, S8, and S9 and Tables ST5 and ST6 correspond to H1N11918 and are similar to Figs. 1, 2, and 3 and Tables 1 and 2. Figures S10, S11, and S12 and Table ST7 and ST8 correspond to H1N12009 and are again similar to Figs. 2 and 3 and Table 2. Figures S13, S14, and S15 and Tables ST9 and ST10 correspond to H2N21957 and are similar to Figs. 2 and 3 and Table 2. The same goes for Figures S16, S17, and S18 and Tables ST11 and ST12, which correspond to H3N21968. From these figures and tables, it follows that:

a)
Peak SFA periodicities for H1N11918
Exact relationships
Almost exact relationships
Complexity measure, C=6

b)
Peak SFA periodicities for H1N12009
$$ {\displaystyle \begin{array}{c}{P}_1=55.5956928123500\\ {}{P}_2=157.248365569507\ \\ {}\begin{array}{c}{P}_3=288.394774029129\ \\ {}{P}_4=576.789548058258\ \\ {}\begin{array}{c}{P}_5=1057.83669529524\ \\ {}\begin{array}{c}{P}_6=2515.97384911212\ \\ {}{P}_7=5984.02800504983\ \end{array}\end{array}\end{array}\end{array}} $$
Exact relationships
Almost exact relationships (note 38=2×19)
Complexity measure, C=4

c)
Peak SFA periodicities for H2N21957
$$ {\displaystyle \begin{array}{c}{P}_1=72.0986935072823\\ {}{P}_2=203.925900374759\ \\ {}\begin{array}{c}{P}_3=288.394774029129\ \\ {}{P}_4=576.789548058258\ \\ {}\begin{array}{c}{P}_5=889.531084997600\ \\ {}\begin{array}{c}{P}_6=2515.97384911212\ \\ {}{P}_7=5984.02800504983\ \end{array}\end{array}\end{array}\end{array}} $$
Exact relationships
Almost exact relationships
Complexity measure, C=4

d)
Peak SFA periodicities for H3N21968
Exact relationships
Almost exact periodicities
Complexity measure, C=2
Interrelationships
As in the case of the coronaviruses, the influenza virus sequence analysis also revealed plenty of interrelationships as expected, since the four viruses belong to the same family.
Interestingly, we found that many relationships exist between the two viral families investigated here. If we compare the results in this section to the previous section, we can infer that:
More on this is discussed next.
Discussion
If we consider the peak SFA periodicities from a sequence as nodes of a community, and their relationships as links between the nodes, then, a visualization of the results for the SARSCoV2 community would look like the top left panel of Fig. 4. Since there are 10 peak periodicities, we have 10 nodes. Then, from Eqs. 3 and 4, we have 13 (recall that C=13) links between them (showing in blue). The rest of the panels correspond to the rest of the sequences in both families. The red lines give the links between the communities within a family (from Eq. 9 for the coronavirus family and from Eq. 18 for the influenza family). The black lines are the links between the two families (Eq. 19). This picture is a perfect example of complex networks, which are often characterized by a community structure, where in each community the nodes are connected in a certain way (meaning the community obeys its own dynamics), but where there exist also some connections (or interactions) between the communities (see for example [18, 19]). We note two interesting observations: (1) the influenza virus family is much more connected (more red links) than the coronavirus family, possibly indicating that the influenza strains are less mutated than the coronavirus strains, and (2) SARSCoV1 has no direct links to the influenza family.
This result supports our claims that SFA has the potential and efficiency to delineate the complex mathematical structure of genetic sequences and that it could become a useful tool in such analyses. We need to stress here that, given the mathematics behind SFA, while we can make direct comparisons of the complexity measure “C” within a certain family (where more or less the number of bases is the same), we cannot compare complexities based on “C” between different families. This is due to the differences in nucleotide length between viral families. The coronavirus family sequence length is approximately 30,450 bases, whereas the influenza family sequence length is approximately 13,500 bases. As such, SFA may “see” longer oscillations in the coronavirus family than in the influenza family. Thus, there will be more entries above the diagonal (in tables such as Table 2), and therefore, higher complexity in the coronavirus family.
Finally, it is interesting to note that the complexity measure “C” in the case of the coronavirus family relates to mortality and severity of symptoms as well as to the rate of transmission. As “C” increases, the transmission rate to humans increases, but mortality rate decreases. It is reported that symptoms of the SARSCoV2 are milder than SARSCoV1 and MERS; however, the viral transmission rate (from humantohuman) is greater than the other family members. The mortality rate of SARSCoV2 is lower (3.4%) than that of SARSCoV1 (9.6%) and MERS (35%) [20]. This relationship is not as clear, however, in the case of the influenza virus family. Unfortunately, in this case, the outbreaks span over a century, and the actual numbers are skewed by several factors such as deaths by secondary infection (due to the unavailability of antibiotics), hygiene, lack of experience and lack of proper healthcare, especially in the early outbreaks, and other problems. For example, H1N11918 (C=6) infected 30% of the planet’s population and H1N12009 (C=4) infected 10% of the population. This is consistent with “increasing C ➔ higher infection rate”, but it is not consistent with “increasing C➔ less mortality rate”. H1N11918 killed about 8% of the infected, whereas H1N12009 killed only 0.0025% of the infected [21,22,23,24,25,26,27,28]. But how can we compare the conditions in 1918 and 2009? To complicate comparisons further, there is hardly any reliable data of infection rates for H2N2 and H3N3. In any case, the complexity measure “C” may hold a promise and could become a useful tool in the prediction of transmission and mortality rates in future new viral strains.
Conclusions
In this exploratory work, a relatively recent mathematical method (SFA) is applied to probe the structure of the nucleotide sequences of two related viral families: those of coronaviruses and those of influenza viruses. The coronaviruses are SARSCoV2, SARSCoV1, and MERS. The influenza viruses include H1N11918, H1N12009, H2N21957, and H3N21968. The analysis indicates that the nucleotide sequences exhibit an elaborate and convoluted structure akin to complex networks. We define a measure of complexity and show that each nucleotide sequence exhibits a certain degree of complexity within itself, while at the same time there exists complex interrelationships between the sequences within a family and between the two families. From these relationships, we find evidence, especially for the coronavirus family, that increasing complexity in a sequence is associated with higher transmission rate but with lower mortality. As such, the complexity measure defined here may hold a promise and could become a useful tool in the prediction of transmission and mortality rates in future new viral strains.
Availability of data and materials
All nucleotide sequences used in this analysis are public domain and can be downloaded from the National Center for Biotechnology Information https://www.ncbi.nlm.nih.gov. For convenience, we have supplied all the nucleotide sequences used here.
Change history
02 June 2021
A Correction to this paper has been published: https://doi.org/10.1186/s40246021003343
Notes
Note here that before we applied SFA to actual nucleotide sequences, and in order to test the efficiency of SFA when the time series is a string of integers, we considered artificial sequences of known periodicities. SFA was able to reproduce the known periodicities.
References
Shepherd JCW. From primeval message to present day gene. CSH Symp Quant Biol. 1982;47:1099–108.
Ohno S. Codon preference is but an illusion created by the construction principle of coding sequences. Proc Natl Acad Sci USA. 1998;85:4378–82.
Yomo T, Ohno S. Concordant evolution and noncoding regions made it possible by the universal rule of TA/CG deficiencyTG/Ct excess. Proc Natl Acad Sci USA. 1989;86:8452–6.
Tsonis AA, Elsner JB, Tsonis PA. Periodicity in DNA sequences: Implications in gene evolution. J Theor Biol. 1991;151:323–31.
Tsonis AA, Kumar P, Elsner JB, Tsonis PA. Wavelet analysis of DNA sequences. Phys Rev E. 1996;53:1828–34.
Lask AM. Introduction to Bioinformatics, 3^{rd} edition, Oxford University press; 2008. p. 474.
Pevsner J. Bioinformatics and Functional Genomics, 3^{rd} edition, WilleyBlackwell; 2015. p. 1124.
Wiskott L, Sejnowski TJ. Slow Feature Analysis: Unsupervised learning of invariance. Neural Comput. 2002;14:715–70.
Wiskott L. Estimating driving forces of nonstationary time series with slow feature analysis, (2003), http://arxiv.org/abs/condmat/0312317/
Berkes P, Wiskott L. Slow feature analysis yields a rich repertoire of complex cells. J Vis. 2005;5(6):579–602.
Blaschke T, Berkes P, Wiskott L. What is the relationship between slow feature analysis and independent component analysis? Neural Comput. 2006;18(10):2495–508. https://doi.org/10.1162/neco.2006.18.10.2495.
Franzius M, Wilbert N, Wiskott L. Invariant object recognition and pose estimation with slow feature analysis. Neural Comput. 2011;23(9):2289–323. https://doi.org/10.1162/NECO_a_00171.
Yang P, Wang G, Zhang F, Zhou X. Causality of global warming seen from observations: a scale analysis of driving force of the surface air temperature time series in the Northern Hemisphere. Clim Dyn. 2015. https://doi.org/10.1007/s0038201527614.
Tsonis AA, Pan X, Wang G, Nicolis C. On the minmax estimation of mean daily temperatures. Clim Dyn. 2019;53:1981–9. https://doi.org/10.1007/s00382019047576.
Wiskott L, et al. Slow feature analysis. Scholarpedia. 2011;6(4):5282.
Torrence C, Compo GP. A practical guide to wavelet analysis. Bull Amer Meteor Soc. 1998;79:61–78. https://doi.org/10.1175/15200477.
Ng WM, Stelfox AJ, Bowden TA. Unraveling virus relationships by structurebased phylogenetic classification. Virus Evol. 2020;6(1). https://doi.org/10.1093/ve/veaa003.
Newman MEJ, Girwan M. Finding and evaluating community structure in networks. Phys Rev E. 2004;69:026113. https://doi.org/10.1103/PhysRevE.69.26113.
Newman MEJ. Modularity and community structure in networks. Proc Natl Acad Sci USA. 2006;103:8577–82.
Fani M, Teimoori A, Ghafari S. Comparison of the COVID2019 (SARSCoV2) pathogenesis with SRSCoV and MERSCoV infections. Future Virol. 2020. https://doi.org/10.2217/fvl20200050.
Andreasen V, Viboud C, Simonsen L. Epidemiologic characterization of the 1918 influenza pandemic summer wave in Copenhagen: implications for pandemic control strategies. J Infect Dis. 2008;197(2):270–728.
Cascella M, Rajnik M, Cuomo A, Dulebohn SC, Di Napoli R. Features, evaluation and treatment Coronavirus (COVID19): StatPearls Publishing; 2020.
Dawood FS, Iuliano AD, Reed C, Meltzer MI, Shay DK, Cheng PY, et al. Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study. Lancet Infect Dis. 2012;12(9):687–95.
Hassan SA, Sheikh FN, Jamal S, Ezeh JK, Akhtar A. Coronavirus (COVID19): a review of clinical features, diagnosis, and treatment. Cureus. 2020;12(3):e7355.
Liu J, Xie W, Wang Y, Xiong Y, Chen S, Han J, et al. A comparative overview of COVID19, MERS and SARS: review article. Int J Surg. 2020;81:1–8.
Taubenberger JK. The origin and virulence of the 1918 ‘Spanish’ Influenza Virus. Proc Am Philos Soc. 2006;150(1):86–112.
Viboud C, Grais RF, Lafont BAP, Miller MA, Simonsen L, Multinational Influenza Seasonal Mortality Study Group. Multinational impact of the 1968 Hong Kong influenza pandemic: evidence for a smoldering pandemic. J Infect Dis. 2005;192(2):233–48.
Viboud C, Simonsen L, Fuentes R, Flores J, Miller MA, Chowell G. Global mortality impact of the 19571959 influenza pandemic. J Infect Dis. 2016;213(5):738–45.
Acknowledgements
Not applicable
Funding
None
Author information
Authors and Affiliations
Contributions
AAT designed the research, did some of the analysis, contributed to the interpretation of the results, and wrote the first draft of the paper. GW, LZ, and WL contributed largely to SFA analysis. AK and KDRT contributed to interpretation of the results and writing of the manuscript. All authors have approved the paper for submission. There are no competing interests.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: The word DNA has been changed to nucleotide or genetic.
Supplementary Information
Additional file 1.
Supplementary figures and tables
Additional file 2.
Nucleotide sequences used in the study
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Tsonis, A.A., Wang, G., Zhang, L. et al. An application of slow feature analysis to the genetic sequences of coronaviruses and influenza viruses. Hum Genomics 15, 26 (2021). https://doi.org/10.1186/s40246021003272
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40246021003272