Global analysis of more than 50,000 SARS-CoV-2 genomes reveals epistasis between eight viral genes

Genome-wide epistasis analysis is a powerful tool to infer gene interactions, which can guide drug and vaccine development and lead to deeper understanding of microbial pathogenesis. We have considered all complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes deposited in the Global Initiative on Sharing All Influenza Data (GISAID) repository until four different cutoff dates, and used direct coupling analysis together with an assumption of quasi-linkage equilibrium to infer epistatic contributions to fitness from polymorphic loci. We find eight interactions, of which three are between pairs where one locus lies in gene ORF3a, both loci holding nonsynonymous mutations. We also find interactions between two loci in gene nsp13, both holding nonsynonymous mutations, and four interactions involving one locus holding a synonymous mutation. Altogether, we infer interactions between loci in viral genes ORF3a and nsp2, nsp12, and nsp6, between ORF8 and nsp4, and between loci in genes nsp2, nsp13, and nsp14. The paper opens the prospect to use prominent epistatically linked pairs as a starting point to search for combinatorial weaknesses of recombinant viral pathogens.

General Models of Multilocus Evolution

In 1991, Barton and Turelli developed recursions to describe the evolution of multilocus systems under arbitrary forms of selection. This article generalizes their approach to allow for arbitrary modes of inheritance, including diploidy, polyploidy, sex linkage, cytoplasmic inheritance, and genomic imprinting. The framework is also extended to allow for other deterministic evolutionary forces, including migration and mutation. Exact recursions that fully describe the state of the population are presented; these are implemented in a computer algebra package (available on the Web at http://helios.bto.ed.ac.uk/evolgen). Despite the generality of our framework, it can describe evolutionary dynamics exactly by just two equations. These recursions can be further simplified using a “quasi-linkage equilibrium” (QLE) approximation. We illustrate the methods by finding the effect of natural selection, sexual selection, mutation, and migration on the genetic composition of a population.

REMARKS ON THE EVOLUTIONARY EFFECT OF NATURAL SELECTION

ABSTRACT The so-called "Fundamental Theorem of Natural Selection", that the mean fitness of a population increases with time under natural selection, is known not to be true, as a mathematical theorem, when fitnesses depend on more than one locus. Although this observation may not have particular biological relevance, (so that mean fitness may well increase in the great majority of interesting situations), it does suggest that it is of interest to find an evolutionary result which is correct as a mathematical theorem, no matter how many loci are involved. The aim of the present note is to prove an evolutionary theorem relating to the variance in fitness, rather than the mean: this theorem is true for an arbitrary number of loci, as well as for arbitrary (fixed) fitness parameters and arbitrary linkage between loci. Connections are briefly discussed between this theorem and the principle of quasi-linkage equilibrium.