The giant Cafeteria roenbergensis virus that infects a widespread marine phagocytic protist is a new member of the fourth domain of Life

Background

A recent work has provided strong arguments in favor of a fourth domain of Life composed of nucleo-cytoplasmic large DNA viruses (NCLDVs). This hypothesis was supported by phylogenetic and phyletic analyses based on a common set of proteins conserved in Eukarya, Archaea, Bacteria, and viruses, and implicated in the functions of information storage and processing. Recently, the genome of a new NCLDV, Cafeteria roenbergensis virus (CroV), was released. The present work aimed to determine if CroV supports the fourth domain of Life hypothesis.

Methods

A consensus phylogenetic tree of NCLDVs including CroV was generated from a concatenated alignment of four universal proteins of NCLDVs. Some features of the gene complement of CroV and its distribution along the genome were further analyzed. Phylogenetic and phyletic analyses were performed using the previously identified common set of informational genes present in Eukarya, Archaea, Bacteria, and NCLDVs, including CroV.

Findings

Phylogenetic reconstructions indicated that CroV is clearly related to the Mimiviridae family. The comparison between the gene repertoires of CroV and Mimivirus showed similarities regarding the gene contents and genome organization. In addition, the phyletic clustering based on the comparison of informational gene repertoire between Eukarya, Archaea, Bacteria, and NCLDVs unambiguously classified CroV with other NCLDVs and clearly included it in a fourth domain of Life. Taken together, these data suggest that Mimiviridae, including CroV, may have inherited a common gene content probably acquired from a common Mimiviridae ancestor.

Conclusions

This further analysis of the gene repertoire of CroV consolidated the fourth domain of Life hypothesis and contributed to outline a functional pan-genome for giant viruses infecting phagocytic protistan grazers.     

Exploring nucleo-cytoplasmic large DNA viruses in Tara Oceans microbial metagenomes

Nucleo-cytoplasmic large DNA viruses (NCLDVs) constitute a group of eukaryotic viruses that can have crucial ecological roles in the sea by accelerating the turnover of their unicellular hosts or by causing diseases in animals. To better characterize the diversity, abundance and biogeography of marine NCLDVs, we analyzed 17 metagenomes derived from microbial samples (0.2–1.6 μm size range) collected during the Tara Oceans Expedition. The sample set includes ecosystems under-represented in previous studies, such as the Arabian Sea oxygen minimum zone (OMZ) and Indian Ocean lagoons. By combining computationally derived relative abundance and direct prokaryote cell counts, the abundance of NCLDVs was found to be in the order of 10⁴–10⁵ genomes ml⁻¹ for the samples from the photic zone and 10²–10³ genomes ml⁻¹ for the OMZ. The Megaviridae and Phycodnaviridae dominated the NCLDV populations in the metagenomes, although most of the reads classified in these families showed large divergence from known viral genomes. Our taxon co-occurrence analysis revealed a potential association between viruses of the Megaviridae family and eukaryotes related to oomycetes. In support of this predicted association, we identified six cases of lateral gene transfer between Megaviridae and oomycetes. Our results suggest that marine NCLDVs probably outnumber eukaryotic organisms in the photic layer (per given water mass) and that metagenomic sequence analyses promise to shed new light on the biodiversity of marine viruses and their interactions with potential hosts.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Activation of tachykinin,neurokinin 3 receptors affects chromatin structure and gene expression by means of histone acetylation

The tachykinin, neurokinin 3 receptor (NK3R) is a g-protein coupled receptor that is broadly distributed in the nervous system and exerts its diverse physiological actions through multiple signaling pathways. Despite the role of the receptor system in a range of biological functions, the effects of NK3R activation on chromatin dynamics and gene expression have received limited attention. The present work determined the effects of senktide, a selective NK3R agonist, on chromatin organization, acetylation, and gene expression, using qRT-PCR, in a hypothalamic cell line (CLU 209) that expresses the NK3R. Senktide (1 nM, 10 nM) caused a relaxation of chromatin, an increase in global acetylation of histone H3 and H4, and an increase in the expression of a common set of genes involved in cell signaling, cell growth, and synaptic plasticity. Pretreatment with histone acetyltransferase (HAT) inhibitor (garcinol and 2-methylene y-butylactone), that inhibits p300, p300/CREB binding protein (CBP) associated factor (PCAF), and GCN 5, prevented the senktide-induced increase in expression of most, but not all, of the genes upregulated in response to 1 nM and 10 nM senktide. Treatment with 100 nM had the opposite effect: a reduction in chromatin relaxation and decreased acetylation. The expression of four genes was significantly decreased and the HAT inhibitor had a limited effect in blocking the upregulation of genes in response to 100 nM senktide. Activation of the NK3R appears to recruit multiple pathways, including acetylation, and possibly histone deactylases, histone methylases, or DNA methylases to affect chromatin structure and gene expression.     

Informational gene phylogenies do not support a fourth domain of life for nucleocytoplasmic large DNA viruses

Mimivirus is a nucleocytoplasmic large DNA virus (NCLDV) with a genome size (1.2 Mb) and coding capacity ( 1000 genes) comparable to that of some cellular organisms. Unlike other viruses, Mimivirus and its NCLDV relatives encode homologs of broadly conserved informational genes found in Bacteria, Archaea, and Eukaryotes, raising the possibility that they could be placed on the tree of life. A recent phylogenetic analysis of these genes showed the NCLDVs emerging as a monophyletic group branching between Eukaryotes and Archaea. These trees were interpreted as evidence for an independent "fourth domain" of life that may have contributed DNA processing genes to the ancestral eukaryote. However, the analysis of ancient evolutionary events is challenging, and tree reconstruction is susceptible to bias resulting from non-phylogenetic signals in the data. These include compositional heterogeneity and homoplasy, which can lead to the spurious grouping of compositionally-similar or fast-evolving sequences. Here, we show that these informational gene alignments contain both significant compositional heterogeneity and homoplasy, which were not adequately modelled in the original analysis. When we use more realistic evolutionary models that better fit the data, the resulting trees are unable to reject a simple null hypothesis in which these informational genes, like many other NCLDV genes, were acquired by horizontal transfer from eukaryotic hosts. Our results suggest that a fourth domain is not required to explain the available sequence data.     

I am what I eat and I eat what I am: acquisition of bacterial genes by giant viruses

Giant viruses are nucleocytoplasmic large DNA viruses (NCLDVs) that infect algae (phycodnaviruses) and amoebae (Mimivirus). We report an unexpected abundance in these giant viruses of islands of bacterial-type genes, including apparently intact prokaryotic mobile genetic elements, and hypothesize that NCLDV genomes undergo successive accretions of bacterial genes. The viruses could acquire bacterial genes within their bacteria-feeding eukaryotic hosts, and we suggest that such acquisition is driven by the intimate coupling of recombination and replication in NCLDVs.     

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (± standard deviation) is 197.8 ± 230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.     

The complete mitochondrial genome of Taeniogonalos taihorina (Bischoff) (Hymenoptera: Trigonalyidae) reveals a novel gene rearrangement pattern in the Hymenoptera

The family Trigonalyidae is considered to be one of the most basal lineages in the suborder Apocrita of Hymenoptera. Here, we determine the first complete mitochondrial genome of the Trigonalyidae, from the species Taeniogonalos taihorina (Bischoff, 1914). This mitochondrial genome is 15,927 bp long, with a high A + T-content of 84.60%. It contains all of the 37 typical animal mitochondrial genes and an A + T-rich region. The orders and directions of all genes are different from those of previously reported hymenopteran mitochondrial genomes. Eight tRNA genes, three protein-coding genes and the A + T-rich region were rearranged, with the dominant gene rearrangement events being translocation and local inversion. The arrangements of three tRNA clusters, trnY–trnM–trnI–trnQ, trnW–trnL2–trnC, and trnH–trnA–trnR–trnN–trnS–trnE–trnF, and the position of the cox1 gene, are novel to the Hymenoptera, even the insects. Six long intergenic spacers are present in the genome. The secondary structures of the RNA genes are normal, except for trnS2, in which the D-stem pairing is absent.     

The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

Key genes and proteins involved in CTCM-reducing microvascular endothelial cell permeability induced by SLT-IIv using gene chips and DIGE

An Affymetrix mouse genome array and differential in-gel electrophoresis (DIGE) techniques were used to investigate the pharmacological mechanisms of a mixture of herbs, designated CTCM, a compound of traditional Chinese medicine, for the treatment of increased permeability in mouse intestinal microvascular endothelial cells (MIMECs) induced by the Shiga-like toxin type II variant (SLT-IIv). MIMECs were challenged with 10 μg/ml SLT-IIv for 12 h and then treated with CTCM at a concentration of 200 μg/ml for 12 h. Total RNA and proteins from each treatment group were extracted from cultured MIMECs for analysis by the Affymetrix GeneChip® Mouse Genome 430 2.0 microarray and DIGE. The results obtained demonstrated that there were one genes downregulated and one genes upregulated, one protein downregulated and four proteins upregulated in the SLT-IIv group compared to the control group. In the CTCM group, four genes were upregulated, three genes were downregulated, a single protein was downregulated and a single protein was upregulated when compared to the control group. When the CTCM-treated group was compared to the SLT-IIv group, expression of one gene was found to be increased, and all other genes were decreased, with five proteins downregulated. Analysis of the data suggested that CTCM specifically and effectively reduced microvascular endothelial cell permeability to SLT-IIv in the treatment of pig edema disease. In the CTCM-treated group, hspa9 expression was increased in both gene chip and DIGE analysis, so it may be a key protein in reducing cell permeability and utilized in medical treatments.     

Phylogenetic and phyletic studies of informational genes in genomes highlight existence of a 4 domain of life including giant viruses

The discovery of Mimivirus, with its very large genome content, made it possible to identify genes common to the three domains of life (Eukarya, Bacteria and Archaea) and to generate controversial phylogenomic trees congruent with that of ribosomal genes, branching Mimivirus at its root. Here we used sequences from metagenomic databases, Marseillevirus and three new viruses extending the Mimiviridae family to generate the phylogenetic trees of eight proteins involved in different steps of DNA processing. Compared to the three ribosomal defined domains, we report a single common origin for Nucleocytoplasmic Large DNA Viruses (NCLDV), DNA processing genes rooted between Archaea and Eukarya, with a topology congruent with that of the ribosomal tree. As for translation, we found in our new viruses, together with Mimivirus, five proteins rooted deeply in the eukaryotic clade. In addition, comparison of informational genes repertoire based on phyletic pattern analysis supports existence of a clade containing NCLDVs clearly distinct from that of Eukarya, Bacteria and Archaea. We hypothesize that the core genome of NCLDV is as ancient as the three currently accepted domains of life.     

Distinct Gene Number-Genome Size Relationships for Eukaryotes and Non-Eukaryotes: Gene Content Estimation for Dinoflagellate Genomes

The ability to predict gene content is highly desirable for characterization of not-yet sequenced genomes like those of dinoflagellates. Using data from completely sequenced and annotated genomes from phylogenetically diverse lineages, we investigated the relationship between gene content and genome size using regression analyses. Distinct relationships between log₁₀-transformed protein-coding gene number (Y′) versus log₁₀-transformed genome size (X′, genome size in kbp) were found for eukaryotes and non-eukaryotes. Eukaryotes best fit a logarithmic model, Y′ = ln(-46.200+22.678X′, whereas non-eukaryotes a linear model, Y′ = 0.045+0.977X′, both with high significance (p<0.001, R²>0.91). Total gene number shows similar trends in both groups to their respective protein coding regressions. The distinct correlations reflect lower and decreasing gene-coding percentages as genome size increases in eukaryotes (82%–1%) compared to higher and relatively stable percentages in prokaryotes and viruses (97%–47%). The eukaryotic regression models project that the smallest dinoflagellate genome (3×10⁶ kbp) contains 38,188 protein-coding (40,086 total) genes and the largest (245×10⁶ kbp) 87,688 protein-coding (92,013 total) genes, corresponding to 1.8% and 0.05% gene-coding percentages. These estimates do not likely represent extraordinarily high functional diversity of the encoded proteome but rather highly redundant genomes as evidenced by high gene copy numbers documented for various dinoflagellate species.     

Infectious Bursal Disease Virus: Ribonucleoprotein Complexes of a Double-Stranded RNA Virus

Genome-binding proteins with scaffolding and/or regulatory functions are common in living organisms and include histones in eukaryotic cells, histone-like proteins in some double-stranded DNA (dsDNA) viruses, and the nucleocapsid proteins of single-stranded RNA viruses. dsRNA viruses nevertheless lack these ribonucleoprotein (RNP) complexes and are characterized by sharing an icosahedral T = 2 core involved in the metabolism and insulation of the dsRNA genome. The birnaviruses, with a bipartite dsRNA genome, constitute a well-established exception and have a single-shelled T = 13 capsid only. Moreover, as in many negative single-stranded RNA viruses, the genomic dsRNA is bound to a nucleocapsid protein (VP3) and the RNA-dependent RNA polymerase (VPg). We used electron microscopy and functional analysis to characterize these RNP complexes of infectious bursal disease virus, the best characterized member of the Birnaviridae family. Mild disruption of viral particles revealed that VP3, the most abundant core protein, present at ∼ 450 copies per virion, is found in filamentous material tightly associated with the dsRNA. We developed a method to purify RNP and VPg-dsRNA complexes. Analysis of these complexes showed that they are linear molecules containing a constant amount of protein. Sensitivity assays to nucleases indicated that VP3 renders the genomic dsRNA less accessible for RNase III without introducing genome compaction. Additionally, we found that these RNP complexes are functionally competent for RNA synthesis in a capsid-independent manner, in contrast to most dsRNA viruses.     

Revised karyotyping and gene mapping of the Biomphalaria glabrata embryonic (Bge) cell line

The fresh water snail Biomphalaria glabrata (2n = 36) belongs to the taxonomic class Gastropoda (family Planorbidae) and is integral to the spread of the human parasitic disease schistosomiasis. The importance of this mollusc is such that it has been selected as a model molluscan organism for whole genome sequencing. In order to understand the structure and organisation of the B. glabrata’s genome it is important that gene mapping studies are established. Thus, we have studied the genomes of two B. glabrata embryonic (Bge) cell line isolates 1 and 2 grown in separate laboratories, but both derived from Eder L. Hansen’s original culture from the 1970s. This cell line continues to be an important tool and model system for schistosomiasis and B. glabrata. Using these cell line isolates, we have investigated the genome content and established a revised karyotype based on chromosome size and centromere position for these cells. Unlike the original karyotype (2n = 36) established for the cell line, our investigations now show the existence of extensive aneuploidy in both cell line isolates to the extent that the total complement of chromosomes in both greatly exceeds the original cell line’s diploid number of 36 chromosomes. The isolates, designated Bge 1 and 2, had modal chromosome complements of 64 and 67, respectively (calculated from 50 metaphases). We found that the aneuploidy was most pronounced, for both isolates, amongst chromosomes of medium metacentric morphology. We also report, to our knowledge for the first time using Bge cells, the mapping of single-copy genes peroxiredoxin (BgPrx4) and P-element induced wimpy testis (piwi) onto Bge chromosomes. These B. glabrata genes were mapped onto pairs of homologous chromosomes using fluorescence in situ hybridization (FISH). Thus, we have now established a FISH mapping technique that can eventually be utilized for physical mapping of the snail genome.     

Construction of a bacterial artificial chromosome library for the Rongchang pig breed and its use for the identification of genes involved in intramuscular fat deposition

In a search for genes affecting intramuscular fat deposition, we constructed a bacterial artificial chromosome (BAC) library for the whole genome of Rongchang pig, a domestic Chinese swine breed. The library consisted of approximately 192,000 clones, with an averaged insert size of 116 kb. Frequency of non-insert clone of the BAC library was no higher than 1.8%, based on estimation of 220 BAC clones randomly selected. We estimated the coverage of the library to be more than seven porcine genome equivalents. Subsequent screening of the BAC library with a three-step PCR procedure resulted in identification of seven candidate genes that were potentially involved in intramuscular fat deposition. The number of positive BAC clones ranged from 2 to 4 for each of the seven genes. One positive clone, containing the lipin1 gene, was fully sequenced by shotgun method to generate 118,041 bp porcine genomic sequences. The BAC clone contained complete DNA sequence of porcine lipin1 gene including all the exons and introns. Our results indicate that this BAC library is a useful tool for gene identification and help to serve as an important resource for future porcine genomic study.