Multiple outgroups can cause random rooting in phylogenomics.

Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.

Morphological Characters Can Strongly Influence Early Animal Relationships Inferred from Phylogenomic Data Sets.

There are considerable phylogenetic incongruencies between morphological and phylogenomic data for the deep evolution of animals. This has contributed to a heated debate over the earliest-branching lineage of the animal kingdom: the sister to all other Metazoa (SOM). Here, we use published phylogenomic data sets ($\sim $45,000-400,000 characters in size with $\sim $15-100 taxa) that focus on early metazoan phylogeny to evaluate the impact of incorporating morphological data sets ($\sim $15-275 characters). We additionally use small exemplar data sets to quantify how increased taxon sampling can help stabilize phylogenetic inferences. We apply a plethora of common methods, that is, likelihood models and their "equivalent" under parsimony: character weighting schemes. Our results are at odds with the typical view of phylogenomics, that is, that genomic-scale data sets will swamp out inferences from morphological data. Instead, weighting morphological data 2-10$\times $ in both likelihood and parsimony can in some cases "flip" which phylum is inferred to be the SOM. This typically results in the molecular hypothesis of Ctenophora as the SOM flipping to Porifera (or occasionally Placozoa). However, greater taxon sampling improves phylogenetic stability, with some of the larger molecular data sets ($>$200,000 characters and up to $\sim $100 taxa) showing node stability even with $\geqq100\times $ upweighting of morphological data. Accordingly, our analyses have three strong messages. 1) The assumption that genomic data will automatically "swamp out" morphological data is not always true for the SOM question. Morphological data have a strong influence in our analyses of combined data sets, even when outnumbered thousands of times by molecular data. Morphology therefore should not be counted out a priori. 2) We here quantify for the first time how the stability of the SOM node improves for several genomic data sets when the taxon sampling is increased. 3) The patterns of "flipping points" (i.e., the weighting of morphological data it takes to change the inferred SOM) carry information about the phylogenetic stability of matrices. The weighting space is an innovative way to assess comparability of data sets that could be developed into a new sensitivity analysis tool. [Metazoa; Morphology; Phylogenomics; Weighting.].

Antimicrobial sensing coupled with cell membrane remodeling mediates antibiotic resistance and virulence in Enterococcus faecalis.

Bacteria have developed several evolutionary strategies to protect their cell membranes (CMs) from the attack of antibiotics and antimicrobial peptides (AMPs) produced by the innate immune system, including remodeling of phospholipid content and localization. Multidrug-resistant Enterococcus faecalis, an opportunistic human pathogen, evolves resistance to the lipopeptide daptomycin and AMPs by diverting the antibiotic away from critical septal targets using CM anionic phospholipid redistribution. The LiaFSR stress response system regulates this CM remodeling via the LiaR response regulator by a previously unknown mechanism. Here, we characterize a LiaR-regulated protein, LiaX, that senses daptomycin or AMPs and triggers protective CM remodeling. LiaX is surface exposed, and in daptomycin-resistant clinical strains, both LiaX and the N-terminal domain alone are released into the extracellular milieu. The N-terminal domain of LiaX binds daptomycin and AMPs (such as human LL-37) and functions as an extracellular sentinel that activates the cell envelope stress response. The C-terminal domain of LiaX plays a role in inhibiting the LiaFSR system, and when this domain is absent, it leads to activation of anionic phospholipid redistribution. Strains that exhibit LiaX-mediated CM remodeling and AMP resistance show enhanced virulence in the Caenorhabditis elegans model, an effect that is abolished in animals lacking an innate immune pathway crucial for producing AMPs. In conclusion, we report a mechanism of antibiotic and AMP resistance that couples bacterial stress sensing to major changes in CM architecture, ultimately also affecting host-pathogen interactions.

Pre-detection history of extensively drug-resistant tuberculosis in KwaZulu-Natal, South Africa.

Antimicrobial-resistant (AMR) infections pose a major threat to global public health. Similar to other AMR pathogens, both historical and ongoing drug-resistant tuberculosis (TB) epidemics are characterized by transmission of a limited number of predominant Mycobacterium tuberculosis (Mtb) strains. Understanding how these predominant strains achieve sustained transmission, particularly during the critical period before they are detected via clinical or public health surveillance, can inform strategies for prevention and containment. In this study, we employ whole-genome sequence (WGS) data from TB clinical isolates collected in KwaZulu-Natal, South Africa to examine the pre-detection history of a successful strain of extensively drug-resistant (XDR) TB known as LAM4/KZN, first identified in a widely reported cluster of cases in 2005. We identify marked expansion of this strain concurrent with the onset of the generalized HIV epidemic 12 y prior to 2005, localize its geographic origin to a location in northeastern KwaZulu-Natal ∼400 km away from the site of the 2005 outbreak, and use protein structural modeling to propose a mechanism for how strain-specific rpoB mutations offset fitness costs associated with rifampin resistance in LAM4/KZN. Our findings highlight the importance of HIV coinfection, high preexisting rates of drug-resistant TB, human migration, and pathoadaptive evolution in the emergence and dispersal of this critical public health threat. We propose that integrating whole-genome sequencing into routine public health surveillance can enable the early detection and local containment of AMR pathogens before they achieve widespread dispersal.

Multidrug-Resistant Methicillin-Resistant Staphylococcus aureus Associated with Bacteremia and Monocyte Evasion, Rio de Janeiro, Brazil.

We typed 600 methicillin-resistant Staphylococcus aureus (MRSA) isolates collected in 51 hospitals in the Rio de Janeiro, Brazil, metropolitan area during 2014-2017. We found that multiple new clonal complex (CC) 5 sequence types had replaced previously dominant MRSA lineages in hospitals. Whole-genome analysis of 208 isolates revealed an emerging sublineage of multidrug-resistant MRSA, sequence type 105, staphylococcal cassette chromosome mec II, spa t002, which we designated the Rio de Janeiro (RdJ) clone. Using molecular clock analysis, we hypothesized that this lineage began to expand in the Rio de Janeiro metropolitan area in 2009. Multivariate analysis supported an association between bloodstream infections and the CC5 lineage that includes the RdJ clone. Compared with other closely related isolates, representative isolates of the RdJ clone more effectively evaded immune function related to monocytic cells, as evidenced by decreased phagocytosis rate and increased numbers of viable unphagocytosed (free) bacteria after in vitro exposure to monocytes.

Transmission of Extensively Drug-Resistant Tuberculosis in South Africa.

BACKGROUND: Drug-resistant tuberculosis threatens recent gains in the treatment of tuberculosis and human immunodeficiency virus (HIV) infection worldwide. A widespread epidemic of extensively drug-resistant (XDR) tuberculosis is occurring in South Africa, where cases have increased substantially since 2002. The factors driving this rapid increase have not been fully elucidated, but such knowledge is needed to guide public health interventions. METHODS: We conducted a prospective study involving 404 participants in KwaZulu-Natal Province, South Africa, with a diagnosis of XDR tuberculosis between 2011 and 2014. Interviews and medical-record reviews were used to elicit information on the participants' history of tuberculosis and HIV infection, hospitalizations, and social networks. Mycobacterium tuberculosis isolates underwent insertion sequence (IS)6110 restriction-fragment-length polymorphism analysis, targeted gene sequencing, and whole-genome sequencing. We used clinical and genotypic case definitions to calculate the proportion of cases of XDR tuberculosis that were due to inadequate treatment of multidrug-resistant (MDR) tuberculosis (i.e., acquired resistance) versus those that were due to transmission (i.e., transmitted resistance). We used social-network analysis to identify community and hospital locations of transmission. RESULTS: Of the 404 participants, 311 (77%) had HIV infection; the median CD4+ count was 340 cells per cubic millimeter (interquartile range, 117 to 431). A total of 280 participants (69%) had never received treatment for MDR tuberculosis. Genotypic analysis in 386 participants revealed that 323 (84%) belonged to 1 of 31 clusters. Clusters ranged from 2 to 14 participants, except for 1 large cluster of 212 participants (55%) with a LAM4/KZN strain. Person-to-person or hospital-based epidemiologic links were identified in 123 of 404 participants (30%). CONCLUSIONS: The majority of cases of XDR tuberculosis in KwaZulu-Natal, South Africa, an area with a high tuberculosis burden, were probably due to transmission rather than to inadequate treatment of MDR tuberculosis. These data suggest that control of the epidemic of drug-resistant tuberculosis requires an increased focus on interrupting transmission. (Funded by the National Institute of Allergy and Infectious Diseases and others.).

An Analysis of the Epidemic of Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae: Convergence of Two Evolutionary Mechanisms Creates the "Perfect Storm".

Background: Carbapenem resistance is a critical healthcare challenge worldwide. Particularly concerning is the widespread dissemination of Klebsiella pneumoniae carbapenemase (KPC). Klebsiella pneumoniae harboring blaKPC (KPC-Kpn) is endemic in many areas including the United States, where the epidemic was primarily mediated by the clonal dissemination of Kpn ST258. We postulated that the spread of blaKPC in other regions occurs by different and more complex mechanisms. To test this, we investigated the evolution and dynamics of spread of KPC-Kpn in Colombia, where KPC became rapidly endemic after emerging in 2005. Methods: We sequenced the genomes of 133 clinical isolates recovered from 24 tertiary care hospitals located in 10 cities throughout Colombia, between 2002 (before the emergence of KPC-Kpn) and 2014. Phylogenetic reconstructions and evolutionary mapping were performed to determine temporal and genetic associations between the isolates. Results: Our results indicate that the start of the epidemic was driven by horizontal dissemination of mobile genetic elements carrying blaKPC-2, followed by the introduction and subsequent spread of clonal group 258 (CG258) isolates containing blaKPC-3. Conclusions: The combination of 2 evolutionary mechanisms of KPC-Kpn within a challenged health system of a developing country created the "perfect storm" for sustained endemicity of these multidrug-resistant organisms in Colombia.

A Prospective Cohort Multicenter Study of Molecular Epidemiology and Phylogenomics of Staphylococcus aureus Bacteremia in Nine Latin American Countries.

Staphylococcus aureus is an important pathogen causing a spectrum of diseases ranging from mild skin and soft tissue infections to life-threatening conditions. Bloodstream infections are particularly important, and the treatment approach is complicated by the presence of methicillin-resistant S. aureus (MRSA) isolates. The emergence of new genetic lineages of MRSA has occurred in Latin America (LA) with the rise and dissemination of the community-associated USA300 Latin American variant (USA300-LV). Here, we prospectively characterized bloodstream MRSA recovered from selected hospitals in 9 Latin American countries. All isolates were typed by pulsed-field gel electrophoresis (PFGE) and subjected to antibiotic susceptibility testing. Whole-genome sequencing was performed on 96 MRSA representatives. MRSA represented 45% of all (1,185 S. aureus) isolates. The majority of MRSA isolates belonged to clonal cluster (CC) 5. In Colombia and Ecuador, most isolates (≥72%) belonged to the USA300-LV lineage (CC8). Phylogenetic reconstructions indicated that MRSA isolates from participating hospitals belonged to three major clades. Clade A grouped isolates with sequence type 5 (ST5), ST105, and ST1011 (mostly staphylococcal chromosomal cassette mec [SCCmec] I and II). Clade B included ST8, ST88, ST97, and ST72 strains (SCCmec IV, subtypes a, b, and c/E), and clade C grouped mostly Argentinian MRSA belonging to ST30. In summary, CC5 MRSA was prevalent in bloodstream infections in LA with the exception of Colombia and Ecuador, where USA300-LV is now the dominant lineage. Clonal replacement appears to be a common phenomenon, and continuous surveillance is crucial to identify changes in the molecular epidemiology of MRSA.

Transferable vancomycin resistance in a community-associated MRSA lineage.

We report the case of a patient from Brazil with a bloodstream infection caused by a strain of methicillin-resistant Staphylococcus aureus (MRSA) that was susceptible to vancomycin (designated BR-VSSA) but that acquired the vanA gene cluster during antibiotic therapy and became resistant to vancomycin (designated BR-VRSA). Both strains belong to the sequence type (ST) 8 community-associated genetic lineage that carries the staphylococcal chromosomal cassette mec (SCCmec) type IVa and the S. aureus protein A gene (spa) type t292 and are phylogenetically related to MRSA lineage USA300. A conjugative plasmid of 55,706 bp (pBRZ01) carrying the vanA cluster was identified and readily transferred to other staphylococci. The pBRZ01 plasmid harbors DNA sequences that are typical of the plasmid-associated replication genes rep24 or rep21 described in community-associated MRSA strains from Australia (pWBG745). The presence and dissemination of community-associated MRSA containing vanA could become a serious public health concern.

Genomic epidemiology of Lineage 4 Mycobacterium tuberculosis subpopulations in New York City and New Jersey, 1999-2009.

BACKGROUND: Whole genome sequencing (WGS) has rapidly become an important research tool in tuberculosis epidemiology and is likely to replace many existing methods in public health microbiology in the near future. WGS-based methods may be particularly useful in areas with less diverse Mycobacterium tuberculosis populations, such as New York City, where conventional genotyping is often uninformative and field epidemiology often difficult. This study applies four candidate strategies for WGS-based identification of emerging M. tuberculosis subpopulations, employing both phylogenomic and population genetics methods. RESULTS: M. tuberculosis subpopulations in New York City and New Jersey can be distinguished via phylogenomic reconstruction, evidence of demographic expansion and subpopulation-specific signatures of selection, and by determination of subgroup-defining nucleotide substitutions. These methods identified known historical outbreak clusters and previously unidentified subpopulations within relatively monomorphic M. tuberculosis endemic clone groups. Neutrality statistics based on the site frequency spectrum were less useful for identifying M. tuberculosis subpopulations, likely due to the low levels of informative genetic variation in recently diverged isolate groups. In addition, we observed that isolates from New York City endemic clone groups have acquired multiple non-synonymous SNPs in virulence- and growth-associated pathways, and relatively few mutations in drug resistance-associated genes, suggesting that overall pathoadaptive fitness, rather than the acquisition of drug resistance mutations, has played a central role in the evolutionary history and epidemiology of M. tuberculosis subpopulations in New York City. CONCLUSIONS: Our results demonstrate that some but not all WGS-based methods are useful for detection of emerging M. tuberculosis clone groups, and support the use of phylogenomic reconstruction in routine tuberculosis laboratory surveillance, particularly in areas with relatively less diverse M. tuberculosis populations. Our study also supports the use of wider-reaching phylogenomic and population genomic methods in tuberculosis public health practice, which can support tuberculosis control activities by identifying genetic polymorphisms contributing to epidemiological success in local M. tuberculosis populations and possibly explain why certain isolate groups are apparently more successful in specific host populations.

Induction of type I interferon signaling determines the relative pathogenicity of Staphylococcus aureus strains.

The tremendous success of S. aureus as a human pathogen has been explained primarily by its array of virulence factors that enable the organism to evade host immunity. Perhaps equally important, but less well understood, is the importance of the intensity of the host response in determining the extent of pathology induced by S. aureus infection, particularly in the pathogenesis of pneumonia. We compared the pathogenesis of infection caused by two phylogenetically and epidemiologically distinct strains of S. aureus whose behavior in humans has been well characterized. Induction of the type I IFN cascade by strain 502A, due to a NOD2-IRF5 pathway, was the major factor in causing severe pneumonia and death in a murine model of pneumonia and was associated with autolysis and release of peptidogylcan. In contrast to USA300, 502A was readily eliminated from epithelial surfaces in vitro. Nonetheless, 502A caused significantly increased tissue damage due to the organisms that were able to invade systemically and trigger type I IFN responses, and this was ameliorated in Ifnar⁻/⁻ mice. The success of USA300 to cause invasive infection appears to depend upon its resistance to eradication from epithelial surfaces, but not production of specific toxins. Our studies illustrate the important and highly variable role of type I IFN signaling within a species and suggest that targeted immunomodulation of specific innate immune signaling cascades may be useful to prevent the excessive morbidity associated with S. aureus pneumonia.

Rickettsial endosymbiont in the "early-diverging" streptophyte green alga Mesostigma viride.

A bacterial endosymbiont was unexpectedly found in the "axenic" culture strain of the streptophyte green alga Mesostigma viride (NIES-995). Phylogenetic analyses based on 16S rRNA gene sequences showed that the symbiont belongs to the order Rickettsiales, specifically to the recently designated clade "Candidatus Megaira," which is closely related to the well-known Rickettsia clade. Rickettsiales bacteria of the "Ca. Megaira" clade are found in a taxonomically diverse array of eukaryotic hosts, including chlorophycean green algae, several ciliate species, and invertebrates such as Hydra. Transmission electron microscopy, fluorescence in situ hybridi-zation, and SYBR Green I staining experiments revealed that the endosymbiont of M. viride NIES-995 is rod shaped, typically occurs in clusters, and is surrounded by a halo-like structure, presumably formed by secretory substances from the bacterium. Two additional M. viride strains (NIES-296 and NIES-475), but not SAG50-1, were found to house the rickettsial endosymbiont. Analyses of strain NIES-995 transcriptome data indicated the presence of at least 91 transcriptionally active genes of symbiont origins. These include genes for surface proteins (e.g., rOmpB) that are known to play key roles in bacterial attachment onto host eukaryotes in related Rickettsia species. The assembled M. viride transcriptome includes transcripts that code for a suite of predicted algal-derived proteins, such as Ku70, WASH, SCAR, and CDC42, which may be important in the formation of the algal-rickettsial association.

Parallel Epidemics of Community-Associated Methicillin-Resistant Staphylococcus aureus USA300 Infection in North and South America.

BACKGROUND: The community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) epidemic in the United States is attributed to the spread of the USA300 clone. An epidemic of CA-MRSA closely related to USA300 has occurred in northern South America (USA300 Latin-American variant, USA300-LV). Using phylogenomic analysis, we aimed to understand the relationships between these 2 epidemics. METHODS: We sequenced the genomes of 51 MRSA clinical isolates collected between 1999 and 2012 from the United States, Colombia, Venezuela, and Ecuador. Phylogenetic analysis was used to infer the relationships and times since the divergence of the major clades. RESULTS: Phylogenetic analyses revealed 2 dominant clades that segregated by geographical region, had a putative common ancestor in 1975, and originated in 1989, in North America, and in 1985, in South America. Emergence of these parallel epidemics coincides with the independent acquisition of the arginine catabolic mobile element (ACME) in North American isolates and a novel copper and mercury resistance (COMER) mobile element in South American isolates. CONCLUSIONS: Our results reveal the existence of 2 parallel USA300 epidemics that shared a recent common ancestor. The simultaneous rapid dissemination of these 2 epidemic clades suggests the presence of shared, potentially convergent adaptations that enhance fitness and ability to spread.

Methicillin-Susceptible, Vancomycin-Resistant Staphylococcus aureus, Brazil.

We report characterization of a methicillin-susceptible, vancomycin-resistant bloodstream isolate of Staphylococcus aureus recovered from a patient in Brazil. Emergence of vancomycin resistance in methicillin-susceptible S. aureus would indicate that this resistance trait might be poised to disseminate more rapidly among S. aureus and represents a major public health threat.

Evolution of hypervirulence by a MRSA clone through acquisition of a transposable element.

Staphylococcus aureus has evolved as a pathogen that causes a range of diseases in humans. There are two dominant modes of evolution thought to explain most of the virulence differences between strains. First, virulence genes may be acquired from other organisms. Second, mutations may cause changes in the regulation and expression of genes. Here we describe an evolutionary event in which transposition of an IS element has a direct impact on virulence gene regulation resulting in hypervirulence. Whole-genome analysis of a methicillin-resistant S. aureus (MRSA) strain USA500 revealed acquisition of a transposable element (IS256) that is absent from close relatives of this strain. Of the multiple copies of IS256 found in the USA500 genome, one was inserted in the promoter sequence of repressor of toxins (Rot), a master transcriptional regulator responsible for the expression of virulence factors in S. aureus. We show that insertion into the rot promoter by IS256 results in the derepression of cytotoxin expression and increased virulence. Taken together, this work provides new insight into evolutionary strategies by which S. aureus is able to modify its virulence properties and demonstrates a novel mechanism by which horizontal gene transfer directly impacts virulence through altering toxin regulation.

The Sorghum bicolor genome and the diversification of grasses.

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

Transcriptome deep-sequencing and clustering of expressed isoforms from Favia corals.

BACKGROUND: Genomic and transcriptomic sequence data are essential tools for tackling ecological problems. Using an approach that combines next-generation sequencing, de novo transcriptome assembly, gene annotation and synthetic gene construction, we identify and cluster the protein families from Favia corals from the northern Red Sea. RESULTS: We obtained 80 million 75 bp paired-end cDNA reads from two Favia adult samples collected at 65 m (Fav1, Fav2) on the Illumina GA platform, and generated two de novo assemblies using ABySS and CAP3. After removing redundancy and filtering out low quality reads, our transcriptome datasets contained 58,268 (Fav1) and 62,469 (Fav2) contigs longer than 100 bp, with N50 values of 1,665 bp and 1,439 bp, respectively. Using the proteome of the sea anemone Nematostella vectensis as a reference, we were able to annotate almost 20% of each dataset using reciprocal homology searches. Homologous clustering of these annotated transcripts allowed us to divide them into 7,186 (Fav1) and 6,862 (Fav2) homologous transcript clusters (E-value ≤ 2e(-30)). Functional annotation categories were assigned to homologous clusters using the functional annotation of Nematostella vectensis. General annotation of the assembled transcripts was improved 1-3% using the Acropora digitifera proteome. In addition, we screened these transcript isoform clusters for fluorescent proteins (FPs) homologs and identified seven potential FP homologs in Fav1, and four in Fav2. These transcripts were validated as bona fide FP transcripts via robust fluorescence heterologous expression. Annotation of the assembled contigs revealed that 1.34% and 1.61% (in Fav1 and Fav2, respectively) of the total assembled contigs likely originated from the corals' algal symbiont, Symbiodinium spp. CONCLUSIONS: Here we present a study to identify the homologous transcript isoform clusters from the transcriptome of Favia corals using a far-related reference proteome. Furthermore, the symbiont-derived transcripts were isolated from the datasets and their contribution quantified. This is the first annotated transcriptome of the genus Favia, a major increase in genomics resources available in this important family of corals.

Long-Read-Based Genome Assembly Reveals Numerous Endogenous Viral Elements in the Green Algal Bacterivore Cymbomonas tetramitiformis.

The marine tetraflagellate Cymbomonas tetramitiformis has drawn attention as an early diverging green alga that uses a phago-mixotrophic mode of nutrition (i.e., the ability to derive nourishment from both photosynthesis and bacterial prey). The Cymbomonas nuclear genome was sequenced previously, but due to the exclusive use of short-read (Illumina) data, the assembly suffered from missing a large proportion of the genome's repeat regions. For this study, we generated Oxford Nanopore long-read and additional short-read Illumina data and performed a hybrid assembly that significantly improved the total assembly size and contiguity. Numerous endogenous viral elements were identified in the repeat regions of the new assembly. These include the complete genome of a giant Algavirales virus along with many genomes of integrated Polinton-like viruses (PLVs) from two groups: Gezel-like PLVs and a novel group of prasinophyte-specific PLVs. The integrated ∼400 kb genome of the giant Algavirales virus is the first account of the association of the uncultured viral family AG_03 with green algae. The complete PLV genomes from C. tetramitiformis ranged between 15 and 25 kb in length and showed a diverse gene content. In addition, heliorhodopsin gene-containing repeat elements of putative mirusvirus origin were identified. These results illustrate past (and possibly ongoing) multiple alga-virus interactions that accompanied the genome evolution of C. tetramitiformis.

One high quality genome and two transcriptome datasets for new species of Mantamonas, a deep-branching eukaryote clade.

Mantamonads were long considered to represent an "orphan" lineage in the tree of eukaryotes, likely branching near the most frequently assumed position for the root of eukaryotes. Recent phylogenomic analyses have placed them as part of the "CRuMs" supergroup, along with collodictyonids and rigifilids. This supergroup appears to branch at the base of Amorphea, making it of special importance for understanding the deep evolutionary history of eukaryotes. However, the lack of representative species and complete genomic data associated with them has hampered the investigation of their biology and evolution. Here, we isolated and described two new species of mantamonads, Mantamonas vickermani sp. nov. and Mantamonas sphyraenae sp. nov., for each of which we generated transcriptomic sequence data, as well as a high-quality genome for the latter. The estimated size of the M. sphyraenae genome is 25 Mb; our de novo assembly appears to be highly contiguous and complete with 9,416 predicted protein-coding genes. This near-chromosome-scale genome assembly is the first described for the CRuMs supergroup.