While our knowledge about the roles of microbes and viruses in the ocean has increased tremendously due to recent advances in genomics and metagenomics, research on marine microbial eukaryotes and zooplankton has benefited much less from these new technologies because of their larger genomes, their enormous diversity, and largely unexplored physiologies. Here, we use a metatranscriptomics approach to capture expressed genes in open ocean Tara Oceans stations across four organismal size fractions. The individual sequence reads cluster into 116 million unigenes representing the largest reference collection of eukaryotic transcripts from any single biome. The catalog is used to unveil functions expressed by eukaryotic marine plankton, and to assess their functional biogeography. Almost half of the sequences have no similarity with known proteins, and a great number belong to new gene families with a restricted distribution in the ocean. Overall, the resource provides the foundations for exploring the roles of marine eukaryotes in ocean ecology and biogeochemistry.

Environmental genomics can describe all forms of organisms—cellular and viral—present in a community. The analysis of such eco-systems biology data relies heavily on reference databases, e.g., taxonomy or gene function databases. Reference databases of symbiosis sensu lato, although essential for the analysis of organism interaction networks, are lacking. By mining existing databases and literature, we here provide a comprehensive and manually curated database of taxonomic links between viruses and their cellular hosts.

Recent discoveries of new large DNA viruses reveal high diversity in their morphologies, genetic repertoires, and replication strategies. Here, we report the novel features of medusavirus, a large DNA virus newly isolated from hot spring water in Japan. Medusavirus, with a diameter of 260 nm, shows a T=277 icosahedral capsid with unique spherical-headed spikes on its surface. It has a 381-kb genome encoding 461 putative proteins, 86 of which have their closest homologs in Acanthamoeba, whereas 279 (61%) are orphan genes. The virus lacks the genes encoding DNA topoisomerase II and RNA polymerase, showing that DNA replication takes place in the host nucleus, whereas the progeny virions are assembled in the cytoplasm. Furthermore, the medusavirus genome harbored genes for all five types of histones (H1, H2A, H2B, H3, and H4) and one DNA polymerase, which are phylogenetically placed at the root of the eukaryotic clades. In contrast, the host amoeba encoded many medusavirus homologs, including the major capsid protein. These facts strongly suggested that amoebae are indeed the most promising natural hosts of medusavirus, and that lateral gene transfers have taken place repeatedly and bidirectionally between the virus and its host since the early stage of their coevolution. Medusavirus reflects the traces of direct evolutionary interactions between the virus and eukaryotic hosts, which may be caused by sharing the DNA replication compartment and by evolutionarily long lasting virus-host relationships. Based on its unique morphological characteristics and phylogenomic relationships with other known large DNA viruses, we propose that medusavirus represents a new family, Medusaviridae. IMPORTANCE We have isolated a new nucleocytoplasmic large DNA virus (NCLDV) from hot spring water in Japan, named medusavirus. This new NCLDV is phylogenetically placed at the root of the eukaryotic clades based on the phylogenies of several key genes, including that encoding DNA polymerase, and its genome surprisingly encodes the full set of histone homologs. Furthermore, its laboratory host, Acanthamoeba castellanii, encodes many medusavirus homologs in its genome, including the major capsid protein, suggesting that the amoeba is the genuine natural host from ancient times of this newly described virus and that lateral gene transfers have repeatedly occurred between the virus and amoeba. These results suggest that medusavirus is a unique NCLDV preserving ancient footprints of evolutionary interactions with its hosts, thus providing clues to elucidate the evolution of NCLDVs, eukaryotes, and virus-host interaction. Based on the dissimilarities with other known NCLDVs, we propose that medusavirus represents a new viral family, Medusaviridae.

Cell type–specific gene expression is driven by the interplay between lineage-specific transcription factors and cis-regulatory elements to which they bind. Adaptive immunity relies on RAG-mediated assembly of T cell receptor (TCR) and immunoglobulin (Ig) genes. Although Rag1 and Rag2 expression is largely restricted to adaptive lymphoid lineage cells, it remains unclear how Rag gene expression is regulated in a cell lineage–specific manner. Here, we identified three distinct cis-regulatory elements, a T cell lineage–specific enhancer (R-TEn) and the two B cell–specific elements, R1B and R2B. By generating mice lacking either R-TEn or R1B and R2B, we demonstrate that these distinct sets of regulatory elements drive the expression of Rag genes in developing T and B cells. What these elements have in common is their ability to bind the transcription factor E2A. By generating a mouse strain that carries a mutation within the E2A binding site of R-TEn, we demonstrate that recruitment of E2A to this site is essential for orchestrating changes in chromatin conformation that drive expression of Rag genes in T cells. By mapping cis-regulatory elements and generating multiple mouse strains lacking distinct enhancer elements, we demonstrate expression of Rag genes in developing T and B cells to be driven by distinct sets of E2A-dependent cis-regulatory modules.

Xanthomonas virus (phage) XacN1 is a novel jumbo myovirus infecting Xanthomonas citri, the causative agent of Asian citrus canker. Its linear 384,670 bp double-stranded DNA genome encodes 592 proteins and presents the longest (66 kbp) direct terminal repeats (DTRs) among sequenced viral genomes. The DTRs harbor 56 tRNA genes, which correspond to all 20 amino acids and represent the largest number of tRNA genes reported in a viral genome. Codon usage analysis revealed a propensity for the phage encoded tRNAs to target codons that are highly used by the phage but less frequently by its host. The existence of these tRNA genes and seven additional translation-related genes as well as a chaperonin gene found in the XacN1 genome suggests a relative independence of phage replication on host molecular machinery, leading to a prediction of a wide host range for this jumbo phage. We confirmed the prediction by showing a wider host range of XacN1 than other X. citri phages in an infection test against a panel of host strains. Phylogenetic analyses revealed a clade of phages composed of XacN1 and ten other jumbo phages, indicating an evolutionary stable large genome size for this group of phages.

Ralstonia solanacearum phages ΦRP12 and ΦRP31 are jumbo phages isolated in Thailand. Here we show that they exhibit similar virion morphology, genome organization and host range. Genome comparisons as well as phylogenetic and proteomic tree analyses support that they belong to the group of ΦKZ-related phages, with their closest relatives being R. solanacearum phages ΦRSL2 and ΦRSF1. Compared with ΦRSL2 and ΦRSF1, ΦRP12 and ΦRP31 possess larger genomes (ca. 280 kbp, 25% larger). The replication of ΦRP12 and ΦRP31 was not affected by rifampicin treatment (20 μg/ml), suggesting that phage-encoded RNAPs function to start and complete the infection cycle of these phages without the need of host-encoded RNAPs. In contrast, ΦRSL2 and ΦRSF1, encoding the same set of RNAPs, did not produce progeny phages in the presence of rifampicin (5 μg/ml). This observation opens the possibility that some ΦRP12/ΦRP31 factors that are absent in ΦRSL2 and ΦRSF1 are involved in their host-independent transcription.

Adaptive immunity relies on the V(D)J DNA recombination of immunoglobulin (Ig) and T cell receptor (TCR) genes, which enables the recognition of highly diverse antigens and the elicitation of antigen-specific immune responses. This process is mediated by recombination-activating gene (Rag) 1 and Rag2 (Rag1/2), whose expression is strictly controlled in a cell type-specific manner; the expression of Rag1/2 genes represents a hallmark of lymphoid lineage commitment. Although Rag genes are known to be evolutionally conserved among jawed vertebrates, how Rag genes are regulated by lineage-specific transcription factors (TFs) and how their regulatory system evolved among vertebrates have not been fully elucidated. Here, we reviewed the current body of knowledge concerning the cis-regulatory elements (CREs) of Rag genes and the evolution of the basic helix-loop-helix TF E protein regulating Rag gene CREs, as well as the evolution of the antagonist of this protein, the Id protein. This may help to understand how the adaptive immune system develops along with the evolution of responsible TFs and enhancers.