Genome structure and thymic expression of an endogenous retrovirus in zebrafish

Zebrafish-based platform for emerging bio-contaminants and virus inactivation research

Emerging bio-contaminants such as viruses have affected health and environment settings of every country. Viruses are the minuscule entities resulting in severe contagious diseases like SARS, MERS, Ebola, and avian influenza. Recent epidemic like the SARS-CoV-2, the virus has undergone mutations strengthen them and allowing to escape from the remedies. Comprehensive knowledge of viruses is essential for the development of targeted therapeutic and vaccination treatments. Animal models mimicking human biology like non-human primates, rats, mice, and rabbits offer competitive advantage to assess risk of viral infections, chemical toxins, nanoparticles, and microbes. However, their economic maintenance has always been an issue. Furthermore, the redundancy of experimental results due to aforementioned aspects is also in examine. Hence, exploration for the alternative animal models is crucial for risk assessments. The current review examines zebrafish traits and explores the possibilities to monitor emerging bio-contaminants. Additionally, a comprehensive picture of the bio contaminant and virus particle invasion and abatement mechanisms in zebrafish and human cells is presented. Moreover, a zebrafish model to investigate the emerging viruses such as coronaviridae and poxviridae has been suggested.

Susceptibility and Permissivity of Zebrafish (Danio rerio) Larvae to Cypriniviruses

The zebrafish (Danio rerio) represents an increasingly important model organism in virology. We evaluated its utility in the study of economically important viruses from the genus Cyprinivirus (anguillid herpesvirus 1, cyprinid herpesvirus 2 and cyprinid herpesvirus 3 (CyHV-3)). This revealed that zebrafish larvae were not susceptible to these viruses after immersion in contaminated water, but that infections could be established using artificial infection models in vitro (zebrafish cell lines) and in vivo (microinjection of larvae). However, infections were transient, with rapid viral clearance associated with apoptosis-like death of infected cells. Transcriptomic analysis of CyHV-3-infected larvae revealed upregulation of interferon-stimulated genes, in particular those encoding nucleic acid sensors, mediators of programmed cell death and related genes. It was notable that uncharacterized non-coding RNA genes and retrotransposons were also among those most upregulated. CRISPR/Cas9 knockout of the zebrafish gene encoding protein kinase R (PKR) and a related gene encoding a protein kinase containing Z-DNA binding domains (PKZ) had no impact on CyHV-3 clearance in larvae. Our study strongly supports the importance of innate immunity-virus interactions in the adaptation of cypriniviruses to their natural hosts. It also highlights the potential of the CyHV-3-zebrafish model, versus the CyHV-3-carp model, for study of these interactions.

Unique Structure and Distinctive Properties of the Ancient and Ubiquitous Gamma-Type Envelope Glycoprotein

After the onset of the AIDS pandemic, HIV-1 (genus Lentivirus) became the predominant model for studying retrovirus Env glycoproteins and their role in entry. However, HIV Env is an inadequate model for understanding entry of viruses in the Alpharetrovirus, Gammaretrovirus and Deltaretrovirus genera. For example, oncogenic model system viruses such as Rous sarcoma virus (RSV, Alpharetrovirus), murine leukemia virus (MLV, Gammaretrovirus) and human T-cell leukemia viruses (HTLV-I and HTLV-II, Deltaretrovirus) encode Envs that are structurally and functionally distinct from HIV Env. We refer to these as Gamma-type Envs. Gamma-type Envs are probably the most widespread retroviral Envs in nature. They are found in exogenous and endogenous retroviruses representing a broad spectrum of vertebrate hosts including amphibians, birds, reptiles, mammals and fish. In endogenous form, gamma-type Envs have been evolutionarily coopted numerous times, most notably as placental syncytins (e.g., human SYNC1 and SYNC2). Remarkably, gamma-type Envs are also found outside of the Retroviridae. Gp2 proteins of filoviruses (e.g., Ebolavirus) and snake arenaviruses in the genus Reptarenavirus are gamma-type Env homologs, products of ancient recombination events involving viruses of different Baltimore classes. Distinctive hallmarks of gamma-type Envs include a labile disulfide bond linking the surface and transmembrane subunits, a multi-stage attachment and fusion mechanism, a highly conserved (but poorly understood) "immunosuppressive domain", and activation by the viral protease during virion maturation. Here, we synthesize work from diverse retrovirus model systems to illustrate these distinctive properties and to highlight avenues for further exploration of gamma-type Env structure and function.

Epigenomic analysis reveals prevalent contribution of transposable elements to cis-regulatory elements, tissue-specific expression, and alternative promoters in zebrafish

Transposable elements (TEs) encode regulatory elements that impact gene expression in multiple species, yet a comprehensive analysis of zebrafish TEs in the context of gene regulation is lacking. Here, we systematically investigate the epigenomic and transcriptomic landscape of TEs across 11 adult zebrafish tissues using multidimensional sequencing data. We find that TEs contribute substantially to a diverse array of regulatory elements in the zebrafish genome and that 37% of TEs are positioned in active regulatory states in adult zebrafish tissues. We identify TE subfamilies enriched in highly specific regulatory elements among different tissues. We use transcript assembly to discover TE-derived transcriptional units expressed across tissues. Finally, we show that novel TE-derived promoters can initiate tissue-specific transcription of alternate gene isoforms. This work provides a comprehensive profile of TE activity across normal zebrafish tissues, shedding light on mechanisms underlying the regulation of gene expression in this widely used model organism.

Human Norovirus-Induced Gene Expression Biomarkers in Zebrafish

ABSTRACT

The challenges associated with development of an animal model system to replicate human norovirus (HuNoV) has hampered the study of the pathogenesis and therapeutic interventions for this virus. In this study, we replicated HuNoV GII.4 and evaluated virus gene expression in infected zebrafish. Three doses of inoculation resulted in successful virus replication. Genes for transmembrane transporters (tfa, cftr, slc26a3, and slc26a6), a heat shock chaperone (hspa8), and immune response cytokines (ifng1 and il1b) were highly expressed in HuNoV-infected zebrafish; however, expression levels of genes were reduced in zebrafish infected with thermally inactivated HuNoV. These results confirm HuNoV replication in juvenile zebrafish and will facilitate the investigation of biomarker gene expression during HuNoV infection.

HIGHLIGHTS

Evolution and Genetic Diversity of the Retroviral Envelope in Anamniotes

Retroviruses are widely distributed in all vertebrates, as are their endogenous forms, endogenous retroviruses (ERV), which serve as "fossil" evidence to trace the ancient origins and history of virus-host interactions over millions of years. The retroviral envelope (Env) plays a significant role in host range determination, but major information on their genetic diversification and evolution in anamniotes is lacking. Here, by incorporating multiple-round in silico similarity search and phylogenomic analysis, more than 30,000 copies of ERV lineages with gamma-type Env (GTE), covalently associated Env, were discovered by searching against all fish and amphibian genomes and transcriptomic assemblies, but no beta-type Env (BTE), noncovalently associated Env, was found. Furthermore, a nine-type classification system of anamniote GTE was proposed by combining phylogenetic and domain/motif analyses. The elastic genomic organization and overall phylogenetic incongruence between anamniotic Env and its neighboring polymerase (Pol) implied that early retroviral diversification in anamniotic vertebrates was facilitated by frequent recombination. At last, host cellular opioid growth factor receptor (OGFr) gene capturing by anamniotic ERVs with GTE was reported for the first time. Overall, our findings overturn traditional Pol genotyping and reveal a complex evolutionary history of anamniotic retroviruses inferred by Env evolution. IMPORTANCE Although the retroviral envelope (Env) protein in amniotes has been well studied, its evolutionary history in anamniotic vertebrates is ambiguous. By analyzing more than 30,000 copies of ERV lineages with gamma-type Env (GTE) in anamniotes, several important evolutionary features were identified. First, GTE was found to be widely distributed among different amphibians and fish. Second, nine types of GTE were discovered and defined, revealing their great genetic diversity. Third, the incongruence between the Env and Pol phylogenies suggested that frequent recombination shaped the early evolution of anamniote retroviruses. Fourth, an ancient horizontal gene transfer event was discovered from anamniotes to ERVs with GTE. These findings reveal a complex evolution pattern for retroviral Env in anamniotes.

A zebrafish reporter line reveals immune and neuronal expression of endogenous retrovirus

Endogenous retroviruses (ERVs) are fossils left in our genome from retrovirus infections of the past. Their sequences are part of every vertebrate genome and their random integrations are thought to have contributed to evolution. Although ERVs are mainly silenced by the host genome, they have been found to be activated in multiple disease states, such as auto-inflammatory disorders and neurological diseases. However, the numerous copies in mammalian genomes and the lack of tools to study them make defining their role in health and diseases challenging. In this study, we identified eight copies of the zebrafish endogenous retrovirus zferv. We created and characterised the first in vivo ERV reporter line in any species. Using a combination of live imaging, flow cytometry and single-cell RNA sequencing, we mapped zferv expression to early T cells and neurons. Thus, this new tool identified tissues expressing ERV in zebrafish, highlighting a potential role of ERV during brain development and strengthening the hypothesis that ERV play a role in immunity and neurological diseases. This transgenic line is therefore a suitable tool to study the function of ERV in health and diseases.

Mapping the m1A, m5C, m6A and m7G methylation atlas in zebrafish brain under hypoxic conditions by MeRIP-seq

BACKGROUND

The epigenetic modifications play important regulatory roles in tissue development, maintenance of physiological functions and pathological process. RNA methylations, including newly identified m1A, m5C, m6A and m7G, are important epigenetic modifications. However, how these modifications are distributed in the transcriptome of vertebrate brains and whether their abundance is altered under pathological conditions are still poorly understood. In this study, we chose the model animal of zebrafish to conduct a systematic study to investigate the mRNA methylation atlas in the brain.

RESULTS

By performing unbiased analyses of the m1A, m5C, m6A and m7G methylation of mRNA, we found that within the whole brain transcriptome, with the increase of the gene expression levels, the overall level of each of these four modifications on the related genes was also progressively increased. Further bioinformatics analysis indicated that the zebrafish brain has an abundance of m1A modifications. In the hypoxia-treated zebrafish brains, the proportion of m1A is decreased, affecting the RNA splicing and zebrafish endogenous retroviruses.

CONCLUSIONS

Our study presents the first comprehensive atlas of m1A, m5C, m6A and m7G in the epitranscriptome of the zebrafish brain and reveals the distribution of these modifications in mRNA under hypoxic conditions. These data provide an invaluable resource for further research on the involvement of m1A, m5C, m6A and m7G in the regulation of miRNA and repeat elements in vertebrates, and provide new thoughts to study the brain hypoxic injury on the aspect of epitranscriptome.

Genome-Wide Characterization of Zebrafish Endogenous Retroviruses Reveals Unexpected Diversity in Genetic Organizations and Functional Potentials

Endogenous retroviruses (ERVs) occupy a substantial fraction of mammalian genomes. However, whether ERVs extensively exist in ancient vertebrates remains unexplored. Here, we performed a genome-wide characterization of ERVs in a zebrafish (Danio rerio) model. Approximately 3,315 ERV-like elements (DrERVs) were identified as Gypsy, Copia, Bel, and class I−III groups. DrERVs accounted for approximately 2.3% of zebrafish genome and were distributed in all 25 chromosomes, with a remarkable bias on chromosome 4. Gypsy and class I are the two most abundant groups with earlier insertion times. The vast majority of the DrERVs have varied structural defects. A total of 509 gag and 71 env genes with coding potentials were detected. The env-coding elements were well-characterized and classified into four subgroups. A ERV-E4.8.43-DanRer element shows high similarity with HERV9NC-int in humans and analogous sequences were detected in species spanning from fish to mammals. RNA-seq data showed that hundreds of DrERVs were expressed in embryos and tissues under physiological conditions, and most of them exhibited stage and tissue specificity. Additionally, 421 DrERVs showed strong responsiveness to virus infection. A unique group of DrERVs with immune-relevant genes, such as fga, ddx41, ftr35, igl1c3, and tbk1, instead of intrinsic viral genes were identified. These DrERVs are regulated by transcriptional factors binding at the long terminal repeats. This study provided a survey of the composition, phylogeny, and potential functions of ERVs in a fish model, which benefits the understanding of the evolutionary history of ERVs from fish to mammals.

IMPORTANCE Endogenous retroviruses (ERVs) are relics of past infection that constitute up to 8% of the human genome. Understanding the genetic evolution of the ERV family and the interplay of ERVs and encoded RNAs and proteins with host function has become a new frontier in biology. Fish, as the most primitive vertebrate host for retroviruses, is an indispensable integral part for such investigations. In the present study, we report the genome-wide characterization of ERVs in zebrafish, an attractive model organism of ancient vertebrates from multiple perspectives, including composition, genomic organization, chromosome distribution, classification, phylogeny, insertion time, characterization of gag and env genes, and expression profiles in embryos and tissues. The result helps uncover the evolutionarily conserved and fish-specific ERVs, as well as the immune-relevant ERVs in response to virus infection. This study demonstrates the previously unrecognized abundance, diversification, and extensive activity of ERVs at the early stage of ERV evolution.

A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis

BACKGROUND

Retroviruses exist as exogenous infectious agents and as endogenous retroviruses (ERVs) integrated into host chromosomes. Such endogenous retroviruses (ERVs) are grouped into three classes roughly corresponding to the seven genera of infectious retroviruses: class I (gamma-, epsilonretroviruses), class II (alpha-, beta-, delta-, lentiretroviruses) and class III (spumaretroviruses). Some ERVs have counterparts among the known infectious retroviruses, while others represent paleovirological relics of extinct or undiscovered retroviruses.

RESULTS

Here we identify an intact ERV in the Anuran amphibian, Xenopus tropicalis. XtERV-S has open reading frames (ORFs) for gag, pol (polymerase) and env (envelope) genes, with a small additional ORF in pol and a serine tRNA primer binding site. It has unusual features and domain relationships to known retroviruses. Analyses based on phylogeny and functional motifs establish that XtERV-S gag and pol genes are related to the ancient env-less class III ERV-L family but the surface subunit of env is unrelated to known retroviruses while its transmembrane subunit is class I-like. LTR constructs show transcriptional activity, and XtERV-S transcripts are detected in embryos after the maternal to zygotic mid-blastula transition and before the late tailbud stage. Tagged Gag protein shows typical subcellular localization. The presence of ORFs in all three protein-coding regions along with identical 5' and 3' LTRs (long terminal repeats) indicate this is a very recent germline acquisition. There are older, full-length, nonorthologous, defective copies in Xenopus laevis and the distantly related African bullfrog, Pyxicephalus adspersus. Additional older, internally deleted copies in X. tropicalis carry a 300 bp LTR substitution.

CONCLUSIONS

XtERV-S represents a genera-spanning member of the largely env-less class III ERV that has ancient and modern copies in Anurans. This provirus has an env ORF with a surface subunit unrelated to known retroviruses and a transmembrane subunit related to class I gammaretroviruses in sequence and organization, and is expressed in early embryogenesis. Additional XtERV-S-related but defective copies are present in X. tropicalis and other African frog taxa. XtERV-S is an unusual class III ERV variant, and it may represent an important transitional retroviral form that has been spreading in African frogs for tens of millions of years.

Modeling Virus-Induced Inflammation in Zebrafish: A Balance Between Infection Control and Excessive Inflammation

The inflammatory response to viral infection in humans is a dynamic process with complex cell interactions that are governed by the immune system and influenced by both host and viral factors. Due to this complexity, the relative contributions of the virus and host factors are best studied in vivo using animal models. In this review, we describe how the zebrafish (Danio rerio) has been used as a powerful model to study host-virus interactions and inflammation by combining robust forward and reverse genetic tools with in vivo imaging of transparent embryos and larvae. The innate immune system has an essential role in the initial inflammatory response to viral infection. Focused studies of the innate immune response to viral infection are possible using the zebrafish model as there is a 4-6 week timeframe during development where they have a functional innate immune system dominated by neutrophils and macrophages. During this timeframe, zebrafish lack a functional adaptive immune system, so it is possible to study the innate immune response in isolation. Sequencing of the zebrafish genome has revealed significant genetic conservation with the human genome, and multiple studies have revealed both functional conservation of genes, including those critical to host cell infection and host cell inflammatory response. In addition to studying several fish viruses, zebrafish infection models have been developed for several human viruses, including influenza A, noroviruses, chikungunya, Zika, dengue, herpes simplex virus type 1, Sindbis, and hepatitis C virus. The development of these diverse viral infection models, coupled with the inherent strengths of the zebrafish model, particularly as it relates to our understanding of macrophage and neutrophil biology, offers opportunities for far more intensive studies aimed at understanding conserved host responses to viral infection. In this context, we review aspects relating to the evolution of innate immunity, including the evolution of viral pattern recognition receptors, interferons and interferon receptors, and non-coding RNAs.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Chromosome Mapping of H1 and H4 Histones in Parodontidae (Actinopterygii: Characiformes): Dispersed and/or Co-Opted Transposable Elements?

The karyotypes of the family Parodontidae consist of 2n = 54 chromosomes. The main chromosomal evolutionary changes of its species are attributed to chromosome rearrangements in repetitive DNA regions in their genomes. Physical mapping of the H1 and H4 histones was performed in 7 Parodontidae species to analyze the chromosome rearrangements involved in karyotype diversification in the group. In parallel, the observation of a partial sequence of an endogenous retrovirus (ERV) retrotransposon in the H1 histone sequence was evaluated to verify molecular co-option of the transposable elements (TEs) and to assess paralogous sequence dispersion in the karyotypes. Six of the studied species had an interstitial histone gene cluster in the short arm of the autosomal pair 13. Besides this interstitial cluster, in Apareiodon davisi, a probable further site was detected in the terminal region of the long arm in the same chromosome pair. The H1/H4 clusters in Parodon cf. pongoensis were located in the smallest chromosomes (pair 20). In addition, scattered H1 signals were observed on the chromosomes in all species. The H1 sequence showed an ERV in the open reading frame (ORF), and the scattered H1 signals on the chromosomes were attributed to the ERV's location. The H4 sequence had no similarity to the TEs and displayed no dispersed signals. Furthermore, the degeneration of the inner ERV in the H1 sequence (which overlapped a stretch of the H1 ORF) was discussed regarding the likelihood of molecular co-option of this retroelement in histone gene function in Parodontidae.

De novo assembly, characterization, functional annotation and expression patterns of the black tiger shrimp (Penaeus monodon) transcriptome

The black tiger shrimp (Penaeus monodon) remains the second most widely cultured shrimp species globally; however, issues with disease and domestication have seen production levels stagnate over the past two decades. To help identify innovative solutions needed to resolve bottlenecks hampering the culture of this species, it is important to generate genetic and genomic resources. Towards this aim, we have produced the most complete publicly available P. monodon transcriptome database to date based on nine adult tissues and eight early life-history stages (BUSCO - Complete: 98.2% [Duplicated: 51.3%], Fragmented: 0.8%, Missing: 1.0%). The assembly resulted in 236,388 contigs, which were then further segregated into 99,203 adult tissue specific and 58,678 early life-history stage specific clusters. While annotation rates were low (approximately 30%), as is typical for a non-model organisms, annotated transcript clusters were successfully mapped to several hundred functional KEGG pathways. Transcripts were clustered into groups within tissues and early life-history stages, providing initial evidence for their roles in specific tissue functions, or developmental transitions. We expect the transcriptome to provide an essential resource to investigate the molecular basis of commercially relevant-significant traits in P. monodon and other shrimp species.

Loss of DNA methylation in zebrafish embryos activates retrotransposons to trigger antiviral signaling

Complex cytoplasmic nucleotide-sensing mechanisms can recognize foreign DNA based on a lack of methylation and initiate an immune response to clear the infection. Zebrafish embryos with global DNA hypomethylation caused by mutations in the ubiquitin-like with PHD and ring finger domains 1 (uhrf1) or DNA methyltransferase 1 (dnmt1) genes exhibit a robust interferon induction characteristic of the first line of defense against viral infection. We found that this interferon induction occurred in non-immune cells and examined whether intracellular viral sensing pathways in these cells were the trigger. RNA-seq analysis of uhrf1 and dnmt1 mutants revealed widespread induction of Class I retrotransposons and activation of cytoplasmic DNA viral sensors. Attenuating Sting, phosphorylated Tbk1 and, importantly, blocking reverse transcriptase activity suppressed the expression of interferon genes in uhrf1 mutants. Thus, activation of transposons in cells with global DNA hypomethylation mimics a viral infection by activating cytoplasmic DNA sensors. This suggests that antiviral pathways serve as surveillance of cells that have derepressed intragenomic parasites due to DNA hypomethylation.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Endogenous Retroviruses in Fish Genomes: From Relics of Past Infections to Evolutionary Innovations?

The increasing availability of fish genome sequences has allowed to gain new insights into the diversity and host distribution of retroviruses in fish and other vertebrates. This distribution can be assessed through the identification and analysis of endogenous retroviruses, which are proviral remnants of past infections integrated in genomes. Retroviral sequences are probably important for evolution through their ability to induce rearrangements and to contribute regulatory and coding sequences; they may also protect their host against new infections. We argue that the current mass of genome sequences will soon strongly improve our understanding of retrovirus diversity and evolution in aquatic animals, with the identification of new/re-emerging elements and host resistance genes that restrict their infectivity.

Recommendations for Health Monitoring and Reporting for Zebrafish Research Facilities

The presence of subclinical infection or clinical disease in laboratory zebrafish may have a significant impact on research results, animal health and welfare, and transfer of animals between institutions. As use of zebrafish as a model of disease increases, a harmonized method for monitoring and reporting the health status of animals will facilitate the transfer of animals, allow institutions to exclude diseases that may negatively impact their research programs, and improve animal health and welfare. All zebrafish facilities should implement a health monitoring program. In this study, we review important aspects of a health monitoring program, including choice of agents, samples for testing, available testing methodologies, housing and husbandry, cost, test subjects, and a harmonized method for reporting results. Facilities may use these recommendations to implement their own health monitoring program.

Viral Outbreak in Corals Associated with an In Situ Bleaching Event: Atypical Herpes-Like Viruses and a New Megavirus Infecting Symbiodinium

Previous studies of coral viruses have employed either microscopy or metagenomics, but few have attempted to comprehensively link the presence of a virus-like particle (VLP) to a genomic sequence. We conducted transmission electron microscopy imaging and virome analysis in tandem to characterize the most conspicuous viral types found within the dominant Pacific reef-building coral genus Acropora. Collections for this study inadvertently captured what we interpret as a natural outbreak of viral infection driven by aerial exposure of the reef flat coincident with heavy rainfall and concomitant mass bleaching. All experimental corals in this study had high titers of viral particles. Three of the dominant VLPs identified were observed in all tissue layers and budding out from the epidermis, including viruses that were ∼70, ∼120, and ∼150 nm in diameter; these VLPs all contained electron dense cores. These morphological traits are reminiscent of retroviruses, herpesviruses, and nucleocytoplasmic large DNA viruses (NCLDVs), respectively. Some 300–500 nm megavirus-like VLPs also were observed within and associated with dinoflagellate algal endosymbiont (Symbiodinium) cells. Abundant sequence similarities to a gammaretrovirus, herpesviruses, and members of the NCLDVs, based on a virome generated from five Acropora aspera colonies, corroborated these morphology-based identifications. Additionally sequence similarities to two diagnostic genes, a MutS and (based on re-annotation of sequences from another study) a DNA polymerase B gene, most closely resembled Pyramimonas orientalis virus, demonstrating the association of a cosmopolitan megavirus with Symbiodinium. We also identified several other virus-like particles in host tissues, along with sequences phylogenetically similar to circoviruses, phages, and filamentous viruses. This study suggests that viral outbreaks may be a common but previously undocumented component of natural bleaching events, particularly following repeated episodes of multiple environmental stressors.

Characterizing novel endogenous retroviruses from genetic variation inferred from short sequence reads

From Illumina sequencing of DNA from brain and liver tissue from the lion, Panthera leo, and tumor samples from the pike-perch, Sander lucioperca, we obtained two assembled sequence contigs with similarity to known retroviruses. Phylogenetic analyses suggest that the pike-perch retrovirus belongs to the epsilonretroviruses, and the lion retrovirus to the gammaretroviruses. To determine if these novel retroviral sequences originate from an endogenous retrovirus or from a recently integrated exogenous retrovirus, we assessed the genetic diversity of the parental sequences from which the short Illumina reads are derived. First, we showed by simulations that we can robustly infer the level of genetic diversity from short sequence reads. Second, we find that the measures of nucleotide diversity inferred from our retroviral sequences significantly exceed the level observed from Human Immunodeficiency Virus infections, prompting us to conclude that the novel retroviruses are both of endogenous origin. Through further simulations, we rule out the possibility that the observed elevated levels of nucleotide diversity are the result of co-infection with two closely related exogenous retroviruses.

Characterization of the fusion core in zebrafish endogenous retroviral envelope protein

Zebrafish endogenous retrovirus (ZFERV) is the unique endogenous retrovirus in zebrafish, as yet, containing intact open reading frames of its envelope protein gene in zebrafish genome. Similarly, several envelope proteins of endogenous retroviruses in human and other mammalian animal genomes (such as syncytin-1 and 2 in human, syncytin-A and B in mouse) were identified and shown to be functional in induction of cell-cell fusion involved in placental development. ZFERV envelope protein (Env) gene appears to be also functional in vivo because it is expressible. After sequence alignment, we found ZFERV Env shares similar structural profiles with syncytin and other type I viral envelopes, especially in the regions of N- and C-terminal heptad repeats (NHR and CHR) which were crucial for membrane fusion. We expressed the regions of N + C protein in the ZFERV Env (residues 459-567, including predicted NHR and CHR) to characterize the fusion core structure. We found N + C protein could form a stable coiled-coil trimer that consists of three helical NHR regions forming a central trimeric core, and three helical CHR regions packing into the grooves on the surface of the central core. The structural characterization of the fusion core revealed the possible mechanism of fusion mediated by ZFERV Env. These results gave comprehensive explanation of how the ancient virus infects the zebrafish and integrates into the genome million years ago, and showed a rational clue for discovery of physiological significance (e.g., medicate cell-cell fusion).

Embryonic expression of endogenous retroviral RNAs in somatic tissues adjacent to the Oikopleura germline

Selective pressure to maintain small genome size implies control of transposable elements, and most old classes of retrotransposons are indeed absent from the very compact genome of the tunicate Oikopleura dioica. Nonetheless, two families of retrotransposons are present, including the Tor elements. The gene organization within Tor elements is similar to that of LTR retrotransposons and retroviruses. In addition to gag and pol, many Tor elements carry a third gene encoding viral envelope-like proteins (Env) that may mediate infection. We show that the Tor family contains distinct classes of elements. In some classes, env mRNA is transcribed from the 5′LTR as in retroviruses. In others, env is transcribed from an additional promoter located downstream of the 5′LTR. Tor Env proteins are membrane-associated glycoproteins which exhibit some features of viral membrane fusion proteins. Whereas some elements are expressed in the adult testis, many others are specifically expressed in embryonic somatic cells adjacent to primordial germ cells. Such embryonic expression depends on determinants present in the Tor elements and not on their surrounding genomic environment. Our study shows that unusual modes of transcription and expression close to the germline may contribute to the proliferation of Tor elements.

Multiple groups of endogenous epsilon-like retroviruses conserved across primates

Several types of cancer in fish are caused by retroviruses, including those responsible for major outbreaks of disease, such as walleye dermal sarcoma virus and salmon swim bladder sarcoma virus. These viruses form a phylogenetic group often described as the epsilonretrovirus genus. Epsilon-like retroviruses have become endogenous retroviruses (ERVs) on several occasions, integrating into germ line cells to become part of the host genome, and sections of fish and amphibian genomes are derived from epsilon-like retroviruses. However, epsilon-like ERVs have been identified in very few mammals. We have developed a pipeline to screen full genomes for ERVs, and using this pipeline, we have located over 800 endogenous epsilon-like ERV fragments in primate genomes. Genomes from 32 species of mammals and birds were screened, and epsilon-like ERV fragments were found in all primate and tree shrew genomes but no others. These viruses appear to have entered the genome of a common ancestor of Old and New World monkeys between 42 million and 65 million years ago. Based on these results, there is an ancient evolutionary relationship between epsilon-like retroviruses and primates. Clearly, these viruses had the potential to infect the ancestors of primates and were at some point a common pathogen in these hosts. Therefore, this result raises questions about the potential of epsilonretroviruses to infect humans and other primates and about the evolutionary history of these retroviruses.

IMPORTANCE

Epsilonretroviruses are a group of retroviruses that cause several important diseases in fish. Retroviruses have the ability to become a permanent part of the DNA of their host by entering the germ line as endogenous retroviruses (ERVs), where they lose their infectivity over time but can be recognized as retroviruses for millions of years. Very few mammals are known to have epsilon-like ERVs; however, we have identified over 800 fragments of endogenous epsilon-like ERVs in the genomes of all major groups of primates, including humans. These viruses seem to have circulated and infected primate ancestors 42 to 65 million years ago. We are now interested in how these viruses have evolved and whether they have the potential to infect modern humans or other primates.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

Through the looking glass: witnessing host-virus interplay in zebrafish

Host-pathogen interactions can be very complex at all scales; understanding organ- or organism-level events require in vivo approaches. Besides traditional host models such as mice, the zebrafish offers an attractive cocktail of optical accessibility and genetic tractability, blended with a vertebrate-type immunity, where innate responses can easily be separated from adaptive ones. Applied to viral infections, this model has revealed unexpected idiosyncrasies among organs, which we believe may apply to the human situation. We also argue that the dynamic analysis of virus spread and immune response in zebrafish make this model particularly well suited to the exploration of the concept of infection tolerance and resistance in relation to viral diseases.

Viral diseases in zebrafish: what is known and unknown

Naturally occurring viral infections have the potential to introduce confounding variability that leads to invalid and misinterpreted data. Whereas the viral diseases of research rodents are well characterized and closely monitored, no naturally occurring viral infections have been characterized for the laboratory zebrafish (Danio rerio), an increasingly important biomedical research model. Despite the ignorance about naturally occurring zebrafish viruses, zebrafish models are rapidly expanding in areas of biomedical research where the confounding effects of unknown infectious agents present a serious concern. In addition, many zebrafish research colonies remain linked to the ornamental (pet) zebrafish trade, which can contribute to the introduction of new pathogens into research colonies, whereas mice used for research are purpose bred, with no introduction of new mice from the pet industry. Identification, characterization, and monitoring of naturally occurring viruses in zebrafish are crucial to the improvement of zebrafish health, the reduction of unwanted variability, and the continued development of the zebrafish as a model organism. This article addresses the importance of identifying and characterizing the viral diseases of zebrafish as the scope of zebrafish models expands into new research areas and also briefly addresses zebrafish susceptibility to experimental viral infection and the utility of the zebrafish as an infection and immunology model.

A retrospective study of the prevalence and classification of intestinal neoplasia in zebrafish (Danio rerio)

For over a decade, spontaneous intestinal neoplasia has been observed in zebrafish (Danio rerio) submitted to the ZIRC (Zebrafish International Resource Center) diagnostic service. In addition, zebrafish displayed preneoplastic intestinal changes including hyperplasia, dysplasia, and enteritis. A total of 195 zebrafish, representing 2% of the total fish submitted to the service, were diagnosed with these lesions. Neoplastic changes were classified either as adenocarcinoma or small cell carcinoma, with a few exceptions (carcinoma not otherwise specified, tubular adenoma, and tubulovillous adenoma). Tumor prevalence appeared similarly distributed between sexes and generally occurred in zebrafish greater than 1 year of age, although neoplastic changes were observed in fish 6 months of age. Eleven lines displayed these preneoplastic and neoplastic changes, including wild-types and mutants. Affected zebrafish originated from 18 facilities, but the majority of fish were from a single zebrafish research facility (hereafter referred to as the primary facility) that has submitted numerous samples to the ZIRC diagnostic service. Zebrafish from the primary facility submitted as normal sentinel fish demonstrate that these lesions are most often subclinical. Fish fed the diet from the primary facility and held at another location did not develop intestinal lesions, indicating that diet is not the etiologic agent.

Replacement of the human cytomegalovirus promoter with fish enhancer and core elements to control the expression of the G gene of viral haemorrhagic septicemia virus (VHSV)

This work explores some of the possibilities to replace human cytomegalovirus (CMV) core and/or enhancer promoter control elements to create new expression vectors for use with fish. The work is relevant to fish vaccination, since DNA vaccines use eukaryotic expression plasmids controlled by the human cytomegalovirus (CMV) promoter to be effective against novirhabdoviruses, such as viral haemorrhagic septicemia virus (VHSV), one of the most devastating fish viral European diseases. To reduce possible homologous recombination with fish genome, core and enhancer sequences from fish origin, such as trout interferon-inducible myxovirus protein (Mx), zebrafish retrovirus long terminal repeat (LTR) and carp β-actin (AE6), were combined with those of CMV to design alternative hybrid promoters. The substitution of CMV core and/or enhancer with the corresponding elements of Mx or the LTR core maintained a similar in vitro protein G expression level than that obtained by using the CMV promoter. Vectors using the dsRNA-inducible Mx enhancer followed either by the LTR or the AE6 cores showed the highest in vitro protein G expression levels. Furthermore, synthetic constructs using the Mx enhancer maintained their polyI:C induction capabilities despite the core used. Some of these hybrid promoters might contribute to the development of all-fish-vectors for DNA vaccines while others might be useful for more basic studies.

Genomic amplification of an endogenous retrovirus in zebrafish T-cell malignancies

Genomic instability plays a crucial role in oncogenesis. Somatically acquired mutations can disable some genes and inappropriately activate others. In addition, chromosomal rearrangements can amplify, delete, or even fuse genes, altering their functions and contributing to malignant phenotypes. Using array comparative genomic hybridization (aCGH), a technique to detect numeric variations between different DNA samples, we examined genomes from zebrafish (Danio rerio) T-cell leukemias of three cancer-prone lines. In all malignancies tested, we identified recurring amplifications of a zebrafish endogenous retrovirus. This retrovirus, ZFERV, was first identified due to high expression of proviral transcripts in thymic tissue from larval and adult fish. We confirmed ZFERV amplifications by quantitative PCR analyses of DNA from wild-type fish tissue and normal and malignant D. rerio T cells. We also quantified ZFERV RNA expression and found that normal and neoplastic T cells both produce retrovirally encoded transcripts, but most cancers show dramatically increased transcription. In aggregate, these data imply that ZFERV amplification and transcription may be related to T-cell leukemogenesis. Based on these data and ZFERV's phylogenetic relation to viruses of the murine-leukemia-related virus class of gammaretroviridae, we posit that ZFERV may be oncogenic via an insertional mutagenesis mechanism.

EZR1: a novel family of highly expressed retroelements induced by TCDD and regulated by a NF-κB-like factor in embryos of zebrafish (Danio rerio)

Transcript profiling using a zebrafish heart cDNA library previously revealed abundant expressed sequence tags (ESTs) upregulated in zebrafish embryos treated with the aryl hydrocarbon receptor (AHR) agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Here, we identify those ESTs as LTR-containing retroelements termed EZR1 (Expressed-Zebrafish-Retroelement group 1). EZR1 is highly redundant in the genome and includes canonical long terminal repeats (LTRs) flanking an integrase-like open reading frame and a region similar to retroviral envelope protein genes. EZR1 sequences lack reverse transcriptase, RNase H, or protease, indicating retrotransposition would be nonautonomous. No AHR binding motifs were found in the EZR1 promoter region. A putative NF-κB-binding site was found, and TCDD-treated zebrafish embryos had significantly increased levels of nuclear protein(s) binding to this sequence. Protein-EZR1 DNA complex formation was partially competed by a mammalian consensus κB sequence, consistent with NF-κB-like activation contributing to increased protein binding to this site. Mobility of the TCDD-induced protein-EZR1 complex differed from that of authentic NF-κB protein bound to the consensus κB site. The results suggest that EZR1 is regulated by interaction with NF-κB or NF-κB-like protein(s) different from the NF-κB protein binding to the consensus κB site. The nature of the NF-κB-like protein and the relationship between EZR1 induction and cardiovascular toxicity caused by TCDD warrant further investigation.

Characterization of a Xenopus tropicalis endogenous retrovirus with developmental and stress-dependent expression

We report on the identification and characterization of XTERV1, a full-length endogenous retrovirus (ERV) within the genome of the western clawed frog (Xenopus tropicalis). XTERV1 contains all the basic genetic elements common to ERVs, including the classical 5'-long terminal repeat (LTR)-gag-pol-env-3'-LTR architecture, as well as conserved functional motifs inherent to each retroviral protein. Using phylogenetic analysis, we show that XTERV1 is related to the Epsilonretrovirus genus. The X. tropicalis genome harbors a single full-length copy with intact gag and pol open reading frames that localizes to the centromeric region of chromosome 5. About 10 full-length defective copies of XTERV1 are found interspersed in the genome, and 2 of them could be assigned to chromosomes 1 and 3. We find that XTERV1 genes are zygotically transcribed in a regulated spatiotemporal manner during frog development, including metamorphosis. Moreover, XTERV1 transcription is upregulated under certain cellular stress conditions, including cytotoxic and metabolic stresses. Interestingly, XTERV1 Env is found to be homologous to FR47, a protein upregulated following cold exposure in the freeze-tolerant wood frog (Rana sylvatica). In addition, we find that R. sylvatica FR47 mRNA originated from a retroviral element. We discuss the potential role(s) of ERVs in physiological processes in vertebrates.

The source of self: genetic parasites and the origin of adaptive immunity

Stable colonization of the host by viruses (genetic parasites) can alter the systems of host identity and provide immunity against related viruses. To attain the needed stability, some viruses of prokaryotes (P1 phage) use a strategy called an addiction module. The linked protective and destructive gene functions of an addiction module insures both virus persistence but will also destroy cells that interrupt this module and thereby prevent infection by competitors. Previously, I have generalized this concept to also include persistent and lytic states of virus infection, which can be considered as a virus addiction module. Such states often involve defective viruses. In this report, I examine the origin of the adaptive immune system from the perspective of a virus addiction module. The likely role of both endogenous and exogenous retroviruses, DNA viruses, and their defective elements is considered in the origin of all the basal components of adaptive immunity (T-cell receptor, RAG-mediated gene rearrangement, clonal lymphocyte proliferation, antigen surface presentation, apoptosis, and education of immune cells). It is concluded that colonization by viruses and their defectives provides a more coherent explanation for the origin of adaptive immunity.

OVEX1, a novel chicken endogenous retrovirus with sex-specific and left-right asymmetrical expression in gonads

BACKGROUND

In chickens, as in most birds, female gonad morphogenesis is asymmetrical. Gonads appear first rather similarly, but only the left one undergoes full differentiation and gives rise to a functional ovary. The right gonad, in which the cortex does not develop, remains restricted to the medulla and finally regresses. Opportunity was taken of this left-right asymmetry to perform a suppression subtractive hybridization screening to select for transcripts preferentially expressed in the developing left ovary as compared to the right one, and thus identify genes that are potentially involved in the process of ovarian differentiation.

RESULTS

One of these transcripts, named Ovex1 according to its expression profile, corresponds to an endogenous retrovirus that has not been previously characterized. It is transcribed as full-length and singly spliced mRNAs and contains three uninterrupted open reading frames coding potentially for proteins with homology to Gag and Pro-Pol retroviral polyproteins and a third protein showing only a weak similarity with Env glycoproteins. Ovex1 is severely degenerated; it is devoid of typical long terminal repeats and displays some evidence of recombination. An orthologous Ovex1 locus was identified in the genome of zebra finch, a member of a different bird order, and similar sequences were detected in turkey, guinea fowl, and duck DNA. The relationship between these sequences follows the bird phylogeny, suggesting vertical transmission of the endogenous retrovirus for more than 100 million years. Ovex1 is transcribed in chicken gonads with a sex-dependent and left-right asymmetrical pattern. It is first expressed in the cortex of the left indifferent gonads of both sexes. Expression is transient in the left testis and absent in the right one. In developing ovaries, Ovex1 transcription increases sharply in the left cortex and is weakly detected in the medulla. After folliculogenesis, Ovex1-expressing cells constitute the follicular granulosa cell layer. Ovex1 expression highlights a striking desquamation process that leads to profound cortical remodeling associated with follicle morphogenesis.

CONCLUSION

Evidence for a selection pressure at the protein level suggests that this endogenous retrovirus, expressed in the ovarian supporting cell lineage, might play an active role in bird ovarian physiology.

In vitro search for alternative promoters to the human immediate early cytomegalovirus (IE-cMV) to express the G gene of viral haemorrhagic septicemia virus (VHSV) in fish epithelial cells

Present DNA vaccines against fish rhabdoviruses require intramuscular injection (fish-to-fish vaccination) of their G-protein gene under the control of the human immediate early cytomegalovirus (IE-CMV) promoter, while immersion delivery (mass DNA vaccination), for instance, by using fish epithelial-specific promoters, would be more practical for aquaculture. To find fish epithelial-specific promoters alternative to the IE-CMV, a comparative study of the effectiveness of different fish promoters constitutively expressing the G gene of the viral haemorrhagic septicemia virus (VHSV) in the epithelial papulosum cyprini (EPC) cell line was performed. The study included MCV1.4 (an alternative IE-CMV promoter version), AE6 (a version of the carp beta-actin promoter), long terminal repeats (LTR) of zebrafish or walleye retroviruses, trout Mx1, carp myosin-heavy-chain and flatfish pleurocidin promoters and salmonid sleeping beauty (SB)/medaka Tol2 transposon repeats. The G-protein expression in transfected EPC cells was studied by estimating the number of cells expressing the G-protein in their membrane and the average expression level per cell. In addition, in an attempt to reduce their sizes, some regions of the MCV1.4 and AE6 promoters were deleted and expression levels compared to those observed for full-length promoters. Since both zebrafish LTR and carp AE6 promoters were the most effective regulatory sequences for expressing the VHSV G-protein in EPC cells, these sequences might be candidates for new DNA vaccine vectors for fish epithelial tissues avoiding the IE-CMV promoter. Furthermore, known transcription factor binding sites (TFBS) common to most of the fish G-expressing promoters, might enable the future design of fully synthetic or hybrid promoters with improved efficacy of VHSV G-protein expression in epithelial fish cells.

Powerful mutators lurking in the genome

The human genome encodes numerous enzymes capable of deaminating polynucleotides. While they are capable of exquisite specificity, occasionally they result in hypermutation where up to 90 per cent of cytidine or adenosine residues may be edited. As such, they constitute a formidable anti-viral barrier, for no virus can survive such a high mutation rate. As the APOBEC3 group of cytidine deaminases edit single-stranded viral DNA, the crucial question is can they hyperedit chromosomal DNA? Everything points to a positive answer. Nonetheless, hypermutants per se have not yet been described, probably being countered by highly efficient mismatch repair. For the APOBEC3 genes, not only is their physiological function unknown, but also their role in the induction of cancer remains to be determined. Yet given the pace of research, all this is certain to change in the next few years.

Identification and characterization of an exogenous retrovirus from atlantic salmon swim bladder sarcomas

A novel piscine retrovirus has been identified in association with an outbreak of leiomyosarcoma in the swim bladders of Atlantic salmon. The complete nucleotide sequence of the Atlantic salmon swim bladder sarcoma virus (SSSV) provirus is 10.9 kb in length and shares a structure and transcriptional profile similar to those of murine leukemia virus-like simple retroviruses. SSSV appears unique to simple retroviruses by not harboring sequences in the Atlantic salmon genome. Additionally, SSSV differs from other retroviruses in potentially utilizing a methionine tRNA primer binding site. SSSV-associated tumors contain high proviral copy numbers (greater than 30 per cell) and a polyclonal integration pattern. Phylogenetic analysis based on reverse transcriptase places SSSV with zebrafish endogenous retrovirus (ZFERV) between the Gammaretrovirus and Epsilonretrovirus genera. Large regions of continuous homology between SSSV and ZFERV Gag, Pol, and Env suggest that these viruses represent a new group of related piscine retroviruses.

Fish transposons and their potential use in aquaculture

A large part of repetitive DNA of vertebrate genomes have been identified as transposon elements (TEs) or mobile sequences. Although TEs detected to date in most vertebrates are inactivated, active TEs have been found in fish and a salmonid TE has been successfully reactivated by molecular genetic manipulation from inactive genomic copies (Sleeping Beauty, SB). Progress in the understanding of the dynamics, control and evolution of fish TEs will allow the insertion of selected sequences into the fish genomes of germ cells to obtain transgenics or to identify genes important for growth and/or of somatic cells to improve DNA vaccination. Expectations are high for new possible applications to fish of this well developed technology for mammals. Here, we review the present state of knowledge of inactive and active fish TEs and briefly discuss how their possible future applications might be used to improve fish production in aquaculture.

Functional characterization of two newly identified Human Endogenous Retrovirus coding envelope genes

A recent in silico search for coding sequences of retroviral origin present in the human genome has unraveled two new envelope genes that add to the 16 genes previously identified. A systematic search among the latter for a fusogenic activity had led to the identification of two bona fide genes, named syncytin-1 and syncytin-2, most probably co-opted by primate genomes for a placental function related to the formation of the syncytiotrophoblast by cell-cell fusion. Here, we show that one of the newly identified envelope gene, named envP(b), is fusogenic in an ex vivo assay, but that its expression – as quantified by real-time RT-PCR on a large panel of human tissues – is ubiquitous, albeit with a rather low value in most tissues. Conversely, the second envelope gene, named envV, discloses a placenta-specific expression, but is not fusogenic in any of the cells tested. Altogether, these results suggest that at least one of these env genes may play a role in placentation, but most probably through a process different from that of the two previously identified syncytins.

Identification of endogenous retroviral reading frames in the human genome

BACKGROUND

Human endogenous retroviruses (HERVs) comprise a large class of repetitive retroelements. Most HERVs are ancient and invaded our genome at least 25 million years ago, except for the evolutionary young HERV-K group. The far majority of the encoded genes are degenerate due to mutational decay and only a few non-HERV-K loci are known to retain intact reading frames. Additional intact HERV genes may exist, since retroviral reading frames have not been systematically annotated on a genome-wide scale.

RESULTS

By clustering of hits from multiple BLAST searches using known retroviral sequences we have mapped 1.1% of the human genome as retrovirus related. The coding potential of all identified HERV regions were analyzed by annotating viral open reading frames (vORFs) and we report 7836 loci as verified by protein homology criteria. Among 59 intact or almost-intact viral polyproteins scattered around the human genome we have found 29 envelope genes including two novel gammaretroviral types. One encodes a protein similar to a recently discovered zebrafish retrovirus (ZFERV) while another shows partial, C-terminal, homology to Syncytin (HERV-W/FRD).

CONCLUSIONS

This compilation of HERV sequences and their coding potential provide a useful tool for pursuing functional analysis such as RNA expression profiling and effects of viral proteins, which may, in turn, reveal a role for HERVs in human health and disease. All data are publicly available through a database at http://www.retrosearch.dk.