viruses Article Modular Evolution of Coronavirus Genomes Yulia Vakulenko 1,2 , Andrei Deviatkin 3 , Jan Felix Drexler 1,4,5 and Alexander Lukashev 1,3,* 1 Martsinovsky Institute of Medical Parasitology, Tropical and Vector Borne Diseases, Sechenov First Moscow State Medical University, 119435 Moscow, Russia;
[email protected] (Y.V.);
[email protected] (J.F.D.) 2 Department of Virology, Faculty of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia 3 Laboratory of Molecular Biology and Biochemistry, Institute of Molecular Medicine, Sechenov First Moscow State Medical University, 119435 Moscow, Russia;
[email protected] 4 Institute of Virology, Charité-Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, 10117 Berlin, Germany 5 German Centre for Infection Research (DZIF), Associated Partner Site Charité, 10117 Berlin, Germany * Correspondence:
[email protected] Abstract: The viral family Coronaviridae comprises four genera, termed Alpha-, Beta-, Gamma-, and Deltacoronavirus. Recombination events have been described in many coronaviruses infecting humans and other animals. However, formal analysis of the recombination patterns, both in terms of the involved genome regions and the extent of genetic divergence between partners, are scarce. Common methods of recombination detection based on phylogenetic incongruences (e.g., a phylogenetic com- patibility matrix) may fail in cases where too many events diminish the phylogenetic signal. Thus, an approach comparing genetic distances in distinct genome regions (pairwise distance deviation matrix) was set up. In alpha, beta, and delta-coronaviruses, a low incidence of recombination between closely related viruses was evident in all genome regions, but it was more extensive between the spike gene and other genome regions.