Isolation and Identification of a Large Green Alga Virus (Chlorella Virus XW01) of Mimiviridae and Its Virophage (Chlorella Virus Virophage SW01) by Using Unicellular Green Algal Cultures.

Virophages are a group of small double-stranded DNA viruses that infect protist hosts and parasitize the viral factory of host giant/large viruses to propagate. Here, we discover a novel cell-virus-virophage (CVv) tripartite interaction system by using unicellular micro-green algae (Chlorella sp.) as eukaryotic hosts for the first time. Viral particles, resembling known virophages and large alga viruses, are detected in culture supernatants and inside algal cells. Complete genomic sequences of the virophage (Chlorella virus virophage SW01 [CVv-SW01]; 24,744 bp) and large virus (Chlorella virus XW01 [CV-XW01]; 407,612 bp) are obtained from the cocultures. Both genomic and phylogenetic analyses show that CVv-SW01 is closely related to virophages previously found in Dishui Lake. CV-XW01 shares the greatest number of homologous genes (n = 82) with Cafeteria roenbergensis virus (CroV) and phylogenetically represents the closest relative to CroV. This is the first report of a large green alga virus being affiliated with a heterotrophic zooplankton-infecting Cafeteriavirus of the family Mimiviridae. Moreover, the codon usage preferences of CV-XW01 and CVv-SW01 are highly similar to those of CroV and its virophage Mavirus, respectively. The discovery of such a novel CVv system with the green alga Chlorella sp. as the single cellular eukaryotic host paves a way to further investigate the potential interaction mechanism of CVv and its significance in the ecology of green algae and the evolution of large/giant viruses and their parasitic viruses. IMPORTANCE Parasitic virophages are small unicellular eukaryotic dsDNA viruses that rely on the viral factories of coinfecting giant/large dsDNA viruses for propagation. Presently, the identified eukaryotic hosts of isolated virophages were restricted to a free-living amoeba, Acanthamoeba polyphaga, and a widespread marine heterotrophic flagellate, Cafeteria roenbergensis. In this study, we successfully discovered and identified a novel tripartite interaction system comprised of a micro-green alga (Chlorella sp.), Mimiviridae large green alga virus, and virophage at the coculture level, with Chlorella sp. as the eukaryotic host, based on combination analysis of infection, morphotype, genome, and phylogeny. The large green alga virus CV-XW01 represents the closest relative to the Mimiviridae giant virus Cafeteria roenbergensis virus, host virus of the virophage Mavirus, as well as a novel large virus of Mimiviridae that infects a non-protozoan protist host. The virophage CVv-SW01 highly resembles Mavirus in its codon usage frequency and preference, although they are phylogenetically distantly related. These findings give novel insights into the diversity of large/giant viruses and their virophages.

Functional Genomic Analyses Reveal an Open Pan-genome for the Chloroviruses and a Potential for Genetic Innovation in New Isolates.

Chloroviruses (family Phycodnaviridae) are large double-stranded DNA (dsDNA) viruses that infect unicellular green algae present in inland waters. These viruses have been isolated using three main chlorella-like green algal host cells, traditionally called NC64A, SAG, and Pbi, revealing extensive genetic diversity. In this study, we performed a functional genomic analysis on 36 chloroviruses that infected the three different hosts. Phylogenetic reconstruction based on the DNA polymerase B family gene clustered the chloroviruses into three distinct clades. The viral pan-genome consists of 1,345 clusters of orthologous groups of genes (COGs), with 126 COGs conserved in all viruses. Totals of 368, 268, and 265 COGs are found exclusively in viruses that infect NC64A, SAG, and Pbi algal hosts, respectively. Two-thirds of the COGs have no known function, constituting the "dark pan-genome" of chloroviruses, and further studies focusing on these genes may identify important novelties. The proportions of functionally characterized COGs composing the pan-genome and the core-genome are similar, but those related to transcription and RNA processing, protein metabolism, and virion morphogenesis are at least 4-fold more represented in the core genome. Bipartite network construction evidencing the COG sharing among host-specific viruses identified 270 COGs shared by at least one virus from each of the different host groups. Finally, our results reveal an open pan-genome for chloroviruses and a well-established core genome, indicating that the isolation of new chloroviruses can be a valuable source of genetic discovery. IMPORTANCE Chloroviruses are large dsDNA viruses that infect unicellular green algae distributed worldwide in freshwater environments. They comprise a genetically diverse group of viruses; however, a comprehensive investigation of the genomic evolution of these viruses is still missing. Here, we performed a functional pan-genome analysis comprising 36 chloroviruses associated with three different algal hosts in the family Chlorellaceae, referred to as zoochlorellae because of their endosymbiotic lifestyle. We identified a set of 126 highly conserved genes, most of which are related to essential functions in the viral replicative cycle. Several genes are unique to distinct isolates, resulting in an open pan-genome for chloroviruses. This profile is associated with generalist organisms, and new insights into the evolution and ecology of chloroviruses are presented. Ultimately, our results highlight the potential for genetic diversity in new isolates.

Unmasking cellular response of a bloom-forming alga to viral infection by resolving expression profiles at a single-cell level.

Infection by large dsDNA viruses can lead to a profound alteration of host transcriptome and metabolome in order to provide essential building blocks to support the high metabolic demand for viral assembly and egress. Host response to viral infection can typically lead to diverse phenotypic outcome that include shift in host life cycle and activation of anti-viral defense response. Nevertheless, there is a major bottleneck to discern between viral hijacking strategies and host defense responses when averaging bulk population response. Here we study the interaction between Emiliania huxleyi, a bloom-forming alga, and its specific virus (EhV), an ecologically important host-virus model system in the ocean. We quantified host and virus gene expression on a single-cell resolution during the course of infection, using automatic microfluidic setup that captures individual algal cells and multiplex quantitate PCR. We revealed high heterogeneity in viral gene expression among individual cells. Simultaneous measurements of expression profiles of host and virus genes at a single-cell level allowed mapping of infected cells into newly defined infection states and allowed detection specific host response in a subpopulation of infected cell which otherwise masked by the majority of the infected population. Intriguingly, resistant cells emerged during viral infection, showed unique expression profiles of metabolic genes which can provide the basis for discerning between viral resistant and susceptible cells within heterogeneous populations in the marine environment. We propose that resolving host-virus arms race at a single-cell level will provide important mechanistic insights into viral life cycles and will uncover host defense strategies.

The mutual interplay between calcification and coccolithovirus infection.

Two prominent characteristics of marine coccolithophores are their secretion of coccoliths and their susceptibility to infection by coccolithoviruses (EhVs), both of which display variation among cells in culture and in natural populations. We examined the impact of calcification on infection by challenging a variety of Emiliania huxleyi strains at different calcification states with EhVs of different virulence. Reduced cellular calcification was associated with increased infection and EhV production, even though calcified cells and associated coccoliths had significantly higher adsorption coefficients than non-calcified (naked) cells. Sialic acid glycosphingolipids, molecules thought to mediate EhV infection, were generally more abundant in calcified cells and enriched in purified, sorted coccoliths, suggesting a biochemical link between calcification and adsorption rates. In turn, viable EhVs impacted cellular calcification absent of lysis by inducing dramatic shifts in optical side scatter signals and a massive release of detached coccoliths in a subpopulation of cells, which could be triggered by resuspension of healthy, calcified host cells in an EhV-free, 'induced media'. Our findings show that calcification is a key component of the E. huxleyi-EhV arms race and an aspect that is critical both to the modelling of these host-virus interactions in the ocean and interpreting their impact on the global carbon cycle.

Infection of phytoplankton by aerosolized marine viruses.

Marine viruses constitute a major ecological and evolutionary driving force in the marine ecosystems. However, their dispersal mechanisms remain underexplored. Here we follow the dynamics of Emiliania huxleyi viruses (EhV) that infect the ubiquitous, bloom-forming phytoplankton E. huxleyi and show that EhV are emitted to the atmosphere as primary marine aerosols. Using a laboratory-based setup, we showed that the dynamic of EhV aerial emission is strongly coupled to the host-virus dynamic in the culture media. In addition, we recovered EhV DNA from atmospheric samples collected over an E. huxleyi bloom in the North Atlantic, providing evidence for aerosolization of marine viruses in their natural environment. Decay rate analysis in the laboratory revealed that aerosolized viruses can remain infective under meteorological conditions prevailing during E. huxleyi blooms in the ocean, allowing potential dispersal and infectivity over hundreds of kilometers. Based on the combined laboratory and in situ findings, we propose that atmospheric transport of EhV is an effective transmission mechanism for spreading viral infection over large areas in the ocean. This transmission mechanism may also have an important ecological impact on the large-scale host-virus "arms race" during bloom succession and consequently the turnover of carbon in the ocean.

Three-year survey of abundance, prevalence and genetic diversity of chlorovirus populations in a small urban lake.

Inland water environments cover about 2.5 percent of our planet and harbor huge numbers of known and still unknown microorganisms. In this report, we examined water samples for the abundance, prevalence, and genetic diversity of a group of infectious viruses (chloroviruses) that infect symbiotic chlorella-like green algae. Samples were collected on a weekly basis for a period of 24 to 36 months from a recreational freshwater lake in Lincoln, Nebraska, and assayed for infectious viruses by plaque assay. The numbers of infectious virus particles were both host- and site-dependent. The consistent fluctuations in numbers of viruses suggest their impact as key factors in shaping microbial community structures in the water surface. Even in low-viral-abundance months, infectious chlorovirus populations were maintained, suggesting either that the viruses are very stable or that there is ongoing viral production in natural hosts.

Unveiling of the diversity of Prasinoviruses (Phycodnaviridae) in marine samples by using high-throughput sequencing analyses of PCR-amplified DNA polymerase and major capsid protein genes.

Viruses strongly influence the ecology and evolution of their eukaryotic hosts in the marine environment, but little is known about their diversity and distribution. Prasinoviruses infect an abundant and widespread class of phytoplankton, the Mamiellophyceae, and thereby exert a specific and important role in microbial ecosystems. However, molecular tools to specifically identify this viral genus in environmental samples are still lacking. We developed two primer sets, designed for use with polymerase chain reactions and 454 pyrosequencing technologies, to target two conserved genes, encoding the DNA polymerase (PolB gene) and the major capsid protein (MCP gene). While only one copy of the PolB gene is present in Prasinovirus genomes, there are at least seven paralogs for MCP, the copy we named number 6 being shared with other eukaryotic alga-infecting viruses. Primer sets for PolB and MCP6 were thus designed and tested on 6 samples from the Tara Oceans project. The results suggest that the MCP6 amplicons show greater richness but that PolB gave a wider coverage of Prasinovirus diversity. As a consequence, we recommend use of the PolB primer set, which will certainly reveal exciting new insights about the diversity and distribution of prasinoviruses at the community scale.

Genetic data generated from virus-host complexes obtained by membrane co-immobilization are equivalent to data obtained from tangential filtrate virus concentrates and virus cultures.

The sieving and immobilization of virus-host complexes using impact filtration (aka membrane co-immobilization or MCI) is a novel approach to the study of plankton viruses. One of the most interesting characteristics of the method is the possibility of generating data on potential viral hosts without the need of culturing hosts cells. MCI has demonstrated to be useful for studying viruses of picoalgae, but studies comparing data generated by MCI to data obtained by other techniques are lacking. In this work, Ostreococcus virus (OV) and Ostreococcus sp. sequences generated from virus-host complexes obtained by MCI were compared to sequences obtained from tangential filtration (TF) concentrates and virus cultures (VC). Statistical parsimony, phylogenetic analyses, pairwise distance comparisons, and analysis of molecular variance showed that the viral and host sequences obtained by the three methods were highly related to each other, indicating that MCI, TF, and VC produce equivalent results. Minor differences were observed among viral sequences obtained from VC and TF concentrates as well as among host sequences generated from VC and MCI. As discussed in the body of the paper, the divergence observed for cultured cells could be due to selective pressures exerted by culture conditions, whereas the correlate observed for the corresponding viral sequences could obey to a hitchhiking effect.

Genomic and biogeographic characterisation of the novel prasinovirus Mantoniella tinhauana virus 1.

Mamiellophyceae are a ubiquitous class of unicellular green algae in the global ocean. Their ecological importance is highlighted in studies focused on the prominent genera Micromonas, Ostreococcus, and Bathycoccus. Mamiellophyceae are susceptible to prasinoviruses, double-stranded DNA viruses belonging to the nucleocytoplasmic large DNA virus group. Our study represents the first isolation of a prasinovirus in the South China Sea and the only one to infect the globally distributed genus Mantoniella. We conducted a comparative analysis with previously identified viral relatives, encompassing morphological characteristics, host specificity, marker-based phylogenetic placement, and whole-genome sequence comparisons. Although it shares morphological and genetic similarities with established prasinoviruses, this novel virus showed distinct genetic traits, confining its infection to the species Mantoniella tinhauana. We also explored the global biogeography of this prasinovirus and its host by mapping metagenomic data and analysing their relationship with various environmental parameters. Our results demonstrate a pronounced link between the virus and its host, both found predominantly in higher latitudes in the surface ocean. By gaining an increased understanding of the relationships between viruses, hosts, and environments, researchers can better make predictions and potentially implement measures to mitigate the consequences of climate change on oceanic processes.

Draft genome sequence of the coccolithovirus Emiliania huxleyi virus 202.

Emiliania huxleyi virus 202 (EhV-202) is a member of the Coccolithoviridae, a group of viruses that infect the marine coccolithophorid Emiliania huxleyi. EhV-202 has a 160- to 180-nm-diameter icosahedral structure and a genome of approximately 407 kbp, consisting of 485 coding sequences (CDSs). Here we describe the genomic features of EhV-202, together with a draft genome sequence and its annotation, highlighting the homology and heterogeneity of this genome in comparison with the EhV-86 reference genome.

Draft genome sequence of four coccolithoviruses: Emiliania huxleyi virus EhV-88, EhV-201, EhV-207, and EhV-208.

The Coccolithoviridae are a group of viruses which infect the marine coccolithophorid microalga Emiliania huxleyi. The Emiliania huxleyi viruses (known as EhVs) described herein have 160- to 180-nm diameter icosahedral structures, have genomes of approximately 400 kbp, and consist of more than 450 predicted coding sequences (CDSs). Here, we describe the genomic features of four newly sequenced coccolithoviruses (EhV-88, EhV-201, EhV-207, and EhV-208) together with their draft genome sequences and their annotations, highlighting the homology and heterogeneity of these genomes to the EhV-86 model reference genome.

Prasinoviruses of the marine green alga Ostreococcus tauri are mainly species specific.

Prasinoviruses infecting unicellular green algae in the order Mamiellales (class Mamiellophyceae) are commonly found in coastal marine waters where their host species frequently abound. We tested 40 Ostreococcus tauri viruses on 13 independently isolated wild-type O. tauri strains, 4 wild-type O. lucimarinus strains, 1 Ostreococcus sp. ("Ostreococcus mediterraneus") clade D strain, and 1 representative species of each of two other related species of Mamiellales, Bathycoccus prasinos and Micromonas pusilla. Thirty-four out of 40 viruses infected only O. tauri, 5 could infect one other species of the Ostreococcus genus, and 1 infected two other Ostreococcus spp., but none of them infected the other genera. We observed that the overall susceptibility pattern of Ostreococcus strains to viruses was related to the size of two host chromosomes known to show intraspecific size variations, that genetically related viruses tended to infect the same host strains, and that viruses carrying inteins were strictly strain specific. Comparison of two complete O. tauri virus proteomes revealed at least three predicted proteins to be candidate viral specificity determinants.

Genomic sequence and analysis of EhV-99B1, a new coccolithovirus from the Norwegian fjords.

Coccolithoviruses are giant dsDNA viruses that infect Emiliania huxleyi, the most ubiquitous marine microalga. Here, we present the genome of the latest coccolithovirus strain to be sequenced, EhV-99B1, and compare it with two other coccolithovirus genomes (EhV-86 and EhV-163). EhV-99B1 shares a pairwise nucleotide identity of 98% with EhV-163 (the two strains were isolated from the same Norwegian fjord but in different years), and just 96.5% with EhV-86 (isolated in the same spring as EhV-99B1 but in the English Channel). We confirmed and extended the list of relevant genomic differences between these EhVs from the Norwegian fjord and EhVs from the English Channel, namely the removal/insertions of: a phosphate permease, an endonuclease, a transposase, and two specific tRNAs. As a whole, this study provided new clues and insights into the diversity and mechanisms driving the evolution of these large oceanic viruses, in particular those processes involving selfish genetic elements.

[Computational genome analysis of three marine algoviruses].

Computational analysis of genomic sequences of three new marine algoviruses: Tetraselmis viridis virus (TvV-S20 and TvV-SI1 strains) and Dunaliella viridis virus (DvV-SI2 strain) was conducted. Both considerable similarity and essential distinctions between studied strains and the most studied marine algoviruses of Phycodnaviridae family were revealed. Our data show that the tested strains are new viruses with the following features: only they were isolated from marine eukaryotic microalgae T. viridis and D. viridis, coding sequences (CDSs) of their genomes are localized mainly on one of the DNA strands and form several clusters with short intergenic spaces; there are considerable variations in genome structure within viruses and their strains; viral genomic DNA has a high GC-content (55.5 - 67.4%); their genes contain no well-known optimal contexts of translation start codones, and the contexts of terminal codons read-through; the vast majority of viral genes and proteins do not have any matches in gene banks.

Genome sequence of Ostreococcus tauri virus OtV-2 throws light on the role of picoeukaryote niche separation in the ocean.

Ostreococcus tauri, a unicellular marine green alga, is the smallest known free-living eukaryote and is ubiquitous in the surface oceans. The ecological success of this organism has been attributed to distinct low- and high-light-adapted ecotypes existing in different niches at a range of depths in the ocean. Viruses have already been characterized that infect the high-light-adapted strains. Ostreococcus tauri virus (OtV) isolate OtV-2 is a large double-stranded DNA algal virus that infects a low-light-adapted strain of O. tauri and was assigned to the algal virus family Phycodnaviridae, genus Prasinovirus. Our working hypothesis for this study was that different viruses infecting high- versus low-light-adapted O. tauri strains would provide clues to propagation strategies that would give them selective advantages within their particular light niche. Sequence analysis of the 184,409-bp linear OtV-2 genome revealed a range of core functional genes exclusive to this low-light genotype and included a variety of unexpected genes, such as those encoding an RNA polymerase sigma factor, at least four DNA methyltransferases, a cytochrome b(5), and a high-affinity phosphate transporter. It is clear that OtV-2 has acquired a range of potentially functional genes from its host, other eukaryotes, and even bacteria over evolutionary time. Such piecemeal accretion of genes is a trademark of large double-stranded DNA viruses that has allowed them to adapt their propagation strategies to keep up with host niche separation in the sunlit layers of the oceanic environment.

Revealing the Viral Community in the Hadal Sediment of the New Britain Trench.

Marine viruses are widely distributed and influence matter and energy transformation in ecosystems by modulating hosts' metabolism. The hadal trenches represent the deepest marine habitat on Earth, for which the viral communities and related biogeochemical functions are least explored and poorly understood. Here, using the sediment samples (8720 m below sea level) collected from the New Britain Trench (NBT), we investigated the viral community, diversity, and genetic potentials in the hadal sediment habitat for the first time by deep shotgun metagenomic sequencing. We found the NBT sediment viral community was dominated by Siphoviridae, Myoviridae, Podoviridae, Mimiviridae, and Phycodnaviridae, which belong to the dsDNA viruses. However, the large majority of them remained uncharacterized. We found the hadal sediment virome had some common components by comparing the hadal sediment viruses with those of hadal aquatic habitats and those of bathypelagic and terrestrial habitats. It was also distinctive in community structure and had many novel viral clusters not associated with the other habitual virome included in our analyses. Further phylogenetic analysis on its Caudovirales showed novel diversities, including new clades specially evolved in the hadal sediment habitat. Annotation of the NBT sediment viruses indicated the viruses might influence microbial hydrocarbon biodegradation and carbon and sulfur cycling via metabolic augmentation through auxiliary metabolic genes (AMGs). Our study filled in the knowledge gaps on the virome of the hadal sediment habitats and provided insight into the evolution and the potential metabolic functions of the hadal sediment virome.

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples.

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.

A functional calcium-transporting ATPase encoded by chlorella viruses.

Calcium-transporting ATPases (Ca(2+) pumps) are major players in maintaining calcium homeostasis in the cell and have been detected in all cellular organisms. Here, we report the identification of two putative Ca(2+) pumps, M535L and C785L, encoded by chlorella viruses MT325 and AR158, respectively, and the functional characterization of M535L. Phylogenetic and sequence analyses place the viral proteins in group IIB of P-type ATPases even though they lack a typical feature of this class, a calmodulin-binding domain. A Ca(2+) pump gene is present in 45 of 47 viruses tested and is transcribed during virus infection. Complementation analysis of the triple yeast mutant K616 confirmed that M535L transports calcium ions and, unusually for group IIB pumps, also manganese ions. In vitro assays show basal ATPase activity. This activity is inhibited by vanadate, but, unlike that of other Ca(2+) pumps, is not significantly stimulated by either calcium or manganese. The enzyme forms a (32)P-phosphorylated intermediate, which is inhibited by vanadate and not stimulated by the transported substrate Ca(2+), thus confirming the peculiar properties of this viral pump. To our knowledge this is the first report of a functional P-type Ca(2+)-transporting ATPase encoded by a virus.

Isolation of prasinoviruses of the green unicellular algae Ostreococcus spp. on a worldwide geographical scale.

Ostreococcus spp. are extremely small unicellular eukaryotic green algae found worldwide in marine environments, and they are susceptible to attacks by a diverse group of large DNA viruses. Several biologically distinct species of Ostreococcus are known and differ in the ecological niches that they occupy: while O. tauri (representing clade C strains) is found in marine lagoons and coastal seas, strains belonging to clade A, exemplified by O. lucimarinus, are present in different oceans. We used laboratory cultures of clonal isolates of these two species to assay for the presence of viruses in seawater samples from diverse locations. In keeping with the distributions of their host strains, we found a decline in the abundance of O. tauri viruses from a lagoon in southwest France relative to the Mediterranean Sea, whereas in the ocean, no O. tauri viruses were detected. In contrast, viruses infecting O. lucimarinus were detected from distantly separated oceans. DNA sequencing, phylogenetic analyses using a conserved viral marker gene, and a Mantel test revealed no relationship between geographic and phylogenetic distances in viruses infecting O. lucimarinus.

From small hosts come big viruses: the complete genome of a second Ostreococcus tauri virus, OtV-1.

Ostreococcus tauri virus (OtV-1) is a large double-stranded DNA virus and a prospective member of the family Phycodnaviridae, genus Prasinovirus. OtV-1 infects the unicellular marine green alga O. tauri, the smallest known free-living eukaryote. Here we present the 191 761 base pair genome sequence of OtV-1, which has 232 putative protein-encoding and 4 tRNA-encoding genes. Approximately 31% of the viral gene products exhibit a similarity to proteins of known functions in public databases. These include a variety of unexpected genes, for example, a PhoH-like protein, a N-myristoyltransferase, a 3-dehydroquinate synthase, a number of glycosyltransferases and methyltransferases, a prolyl 4-hydroxylase, 6-phosphofructokinase and a total of 8 capsid proteins. A total of 11 predicted genes share homology with genes found in the Ostreococcus host genome. In addition, an intein was identified in the DNA polymerase gene of OtV-1. This is the first report of an intein in the genome of a virus that infects O. tauri. Fifteen core genes common to nuclear-cytoplasmic large dsDNA virus (NCLDV) genomes were identified in the OtV-1 genome. This new sequence data may help to redefine the classification of the core genes of these viruses and shed new light on their evolutionary history.