Characterisation of retroviruses in the horse genome and their transcriptional activity via transcriptome sequencing

Differential Expression Pattern of Retroviral Envelope Gene in the Equine Placenta

Endogenous retroviruses (ERVs) are proviral phases of exogenous retroviruses, which have coevolved with vertebrate genomes for millions of years. The conservation of ERV genes throughout evolution suggests their beneficial effects on their hosts' survival. An example of such positive selection is demonstrated by the syncytin gene, which encodes a protein with affinity for various mammalian placentas that is involved in the formation of syncytiotrophoblasts. Although the horse has an epitheliochorial placenta, in which the fetal trophoblasts are simply apposed to the intact uterine epithelium, we have previously demonstrated that the equine ERV (EqERV) env RNA is unexpectedly expressed in placental tissue. In the present study, we investigated the mRNA expression pattern of the EqERV env gene in different parts of the equine placenta, to gain more insight into its putative role in the fetal–maternal relationship. To this end, we used reverse transcription–quantitative PCR (RT–qPCR) and in situ hybridization assays to analyze different target areas of the equine placenta. The retroviral env gene is expressed in the equine placenta, even though there is no syncytium or erosion of the uterine endometrium. The gene is also expressed in all the sampled areas, although with some quantitative differences. We suggest that these differences are attributable to variations in the density, height, and degree of morphological complexity of the chorionic villi forming the microcotyledons. The involvement of the EqERV env gene in different functional pathways affecting the fetus–mother relationship can be hypothesized.

Distribution, Diversity, and Evolution of Endogenous Retroviruses in Perissodactyl Genomes

The evolution of mammalian genomes has been shaped by interactions with endogenous retroviruses (ERVs). In this study, we investigated the distribution and diversity of ERVs in the mammalian order Perissodactyla, with a view to understanding their impact on the evolution of modern equids (family Equidae). We characterize the major ERV lineages in the horse genome in terms of their genomic distribution, ancestral genome organization, and time of activity. Our results show that subsequent to their ancestral divergence from rhinoceroses and tapirs, equids acquired four novel ERV lineages. We show that two of these ERV lineages proliferated extensively in the lineage leading to modern horses, and one contains loci that are actively transcribed in specific tissues. In addition, we show that the white rhinoceros has resisted germ line colonization by retroviruses for more than 54 million years-longer than any other extant mammalian species. The map of equine ERVs that we provide here will be of great utility to future studies aiming to investigate the potential functional roles of equine ERVs and their impact on equine evolution.IMPORTANCE ERVs in the host genome are highly informative about the long-term interactions of retroviruses and hosts. They are also interesting because they have influenced the evolution of mammalian genomes in various ways. In this study, we derive a calibrated timeline describing the process through which ERV diversity has been generated in the equine germ line. We determined the distribution and diversity of perissodactyl ERV lineages and inferred their retrotranspositional activity during evolution, thereby gaining insight into the long-term coevolutionary history of retroviruses and mammals. Our study provides a platform for future investigations to identify equine ERV loci involved in physiological processes and/or pathological conditions.

Identification and expression analysis of a novel miRNA derived from ERV-E1 LTR in Equus caballus

Horses (Equus caballus) have been domesticated and bred to enhance speed, strength, and agility. Members of the Equus caballus Endogenous Retrovirus (EqERV) family affect several of these abilities in horses. EqERV elements have been integrated in the horse genome during evolution and generate repeat elements such as long terminal repeats (LTRs). LTR sequences are involved in retrovirus replication and play an essential function in post-transcriptional control mechanisms, such as by providing binding sites for microRNAs (miRNAs) or generating miRNA precursors. In this study, we identified a novel miRNA derived from EqERV-E1 LTR using various bioinformatics tools. To examine the relationship between EqERV-E1 LTR and similar elements, we used BLAST2seq and phylogenetic analysis. LTR sequences were located in the untranslated region (UTR) of mRNAs and also formed the stem-loop secondary structure. The sequence was registered in the DDBJ database as LTR derived miRNA under the accession number corresponding to LC383797 (referred to eca-miR-1804). Quantitative polymerase chain reaction (qPCR) to confirm the expression of eca-miR-1804 and the similar miR-1255a, showed an almost identical expression pattern in eight different equine tissues. Therefore, these data imply that the LTR could function as an miRNA, which is expressed in the examined equine tissues. In addition, the current study provides inputs for additional functional studies concerning the LTR of other EqERV families.

Identification and Expression Analyses of Equine Endogenous Retroviruses in Horses

Endogenous retroviruses (ERVs) have been integrated into vertebrate genomes and have momentously affected host organisms. Horses (Equus caballus) have been domesticated and selected for elite racing ability over centuries. ERVs played an important role in the evolutionary diversification of the horse genome. In the present study, we identified six equine ERV families (EqERVs-E1, I1, M2, P1, S1, and Y4), their full-length viral open reading frames (ORFs), and elucidated their phylogenetic relationships. The divergence time of EqERV families assuming an evolutionary rate of 0.2%/Myr indicated that EqERV-S3 (75.4 million years ago; mya) on chromosome 10 is an old EqERV family and EqERV-P5 (1.2 Mya) on chromosome 12 is a young member. During the evolutionary diversification of horses, the EqERV-I family diverged 1.7 Mya to 38.7 Mya. Reverse transcription quantitative real-time PCR (RT-qPCR) amplification of EqERV pol genes showed greater expression in the cerebellum of the Jeju horse than the Thoroughbred horse. These results could contribute further dynamic studies for horse genome in relation to EqERV gene function.

RNA expression of TLR10 in normal equine tissues

Background

Toll like receptors are one of the major innate immune system pathogen recognition systems. There is little data on the expression of the TLR10 member of this family in the horse.

Results

This paper describes the genetic structure of the Equine TLR10 gene and its RNA expression in a range of horse tissues. It describes the phylogenetic analysis of the Equine TLR1,6,10,2 annotations in the horse genome, firmly identifying them in their corresponding gene clades compared to other species and firmly placing the horse gene with other TLR10 genes from odd-toed ungulates. Additional 3’ transcript extensions to that annotated for TLR10 in the horse genome have been identified by analysis of RNAseq data. RNA expression of the equine TLR10 gene was highest in peripheral blood mononucleocytes and lymphoid tissue (lymph nodes and spleen), however some expression was detected in all tissues tested (jejunum, caudal mesenteric lymph nodes, bronchial lymph node, spleen, lung, colon, kidney and liver). Additional data on RNAseq expression of all equine TLR genes (1–4 and 6–10) demonstrate higher expression of TLR4 than other equine TLRs in all tissues.

Conclusion

The equine TLR10 gene displays significant homology to other mammalian TLR10 genes and could be reasonably assumed to have similar fuctions. Its RNA level expression is higher in resting state PBMCs in horses than in other tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-2161-9) contains supplementary material, which is available to authorized users.

High Expression of Endogenous Retroviral Envelope Gene in the Equine Fetal Part of the Placenta

Endogenous retroviruses (ERVs) are proviral phases of exogenous retroviruses that have co-evolved with vertebrate genomes for millions of years. Previous studies have identified the envelope (env) protein genes of retroviral origin preferentially expressed in the placenta which suggests a role in placentation based on their membrane fusogenic capacity and therefore they have been named syncytins. Until now, all the characterized syncytins have been associated with three invasive placentation types: the endotheliochorial (Carnivora), the synepitheliochorial (Ruminantia), and the hemochorial placentation (human, mouse) where they play a role in the syncytiotrophoblast formation. The purpose of the present study was to evaluate whether EqERV env RNA is expressed in horse tissues as well and investigate if the horse, possessing an epitheliochorial placenta, has "captured" a common retroviral env gene with syncytin-like properties in placental tissues. Interestingly, although in the equine placenta there is no syncytiotrophoblast layer at the maternal-fetal interface, our results showed that EqERV env RNA is highly expressed at that level, as expected for a candidate syncytin-like gene but with reduced abundance in the other somatic tissues (nearly 30-fold lower) thus suggesting a possible role in the placental tissue. Although the horse is one of the few domestic animals with a sequenced genome, few studies have been conducted about the EqERV and their expression in placental tissue has never been investigated.

Bone marrow transcriptome and epigenome profiles of equine common variable immunodeficiency patients unveil block of B lymphocyte differentiation

Common variable immunodeficiency (CVID) is a late-onset humoral deficiency characterized by B lymphocyte dysfunction or loss, decreased immunoglobulin production, and recurrent bacterial infections. CVID is the most frequent human primary immunodeficiency but still presents challenges in the understanding of its etiology and treatment. CVID in equine patients manifests with a natural impairment of B lymphocyte differentiation, and is a unique model to identify genetic and epigenetic mechanisms of disease. Bone marrow transcriptome analyses revealed decreased expression of genes indicative of the pro-B cell differentiation stage, importantly PAX5 (p≤0.023). We hypothesized that aberrant epigenetic regulation caused PAX5 gene silencing, resulting in the late-onset and non-familial manifestation of CVID. A significant increase in PAX5 enhancer region methylation was identified in equine CVID patients by genome-wide reduced-representation bisulfite sequencing and bisulfite PCR sequencing (p=0.000). Thus, we demonstrate that integrating transcriptomics and epigenetics in CVID enlightens potential mechanisms of dysfunctional B lymphopoiesis or function.

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.

Endogenous retroviruses in domestic animals

Endogenous retroviruses (ERVs) are genomic elements that are present in a wide range of vertebrates. Although the study of ERVs has been carried out mainly in humans and model organisms, recently, domestic animals have become important, and some species have begun to be analyzed to gain further insight into ERVs. Due to the availability of complete genomes and the development of new computer tools, ERVs can now be analyzed from a genome-wide viewpoint. In addition, more experimental work is being carried out to analyze the distribution, expression and interplay of ERVs within a host genome. Cats, cattle, chicken, dogs, horses, pigs and sheep have been scrutinized in this manner, all of which are interesting species in health and economic terms. Furthermore, several studies have noted differences in the number of endogenous retroviruses and in the variability of these elements among different breeds, as well as their expression in different tissues and the effects of their locations, which, in some cases, are near genes. These findings suggest a complex, intriguing relationship between ERVs and host genomes. In this review, we summarize the most important in silico and experimental findings, discuss their implications and attempt to predict future directions for the study of these genomic elements.

Characterisation of the horse transcriptome from immunologically active tissues

The immune system of the horse has not been well studied, despite the fact that the horse displays several features such as sensitivity to bacterial lipopolysaccharide that make them in many ways a more suitable model of some human disorders than the current rodent models. The difficulty of working with large animal models has however limited characterisation of gene expression in the horse immune system with current annotations for the equine genome restricted to predictions from other mammals and the few described horse proteins. This paper outlines sequencing of 184 million transcriptome short reads from immunologically active tissues of three horses including the genome reference “Twilight”. In a comparison with the Ensembl horse genome annotation, we found 8,763 potentially novel isoforms.

Genome-wide characterization of endogenous retroviruses in the bat Myotis lucifugus reveals recent and diverse infections

Bats are increasingly recognized as reservoir species for a variety of zoonotic viruses that pose severe threats to human health. While many RNA viruses have been identified in bats, little is known about bat retroviruses. Endogenous retroviruses (ERVs) represent genomic fossils of past retroviral infections and, thus, can inform us on the diversity and history of retroviruses that have infected a species lineage. Here, we took advantage of the availability of a high-quality genome assembly for the little brown bat, Myotis lucifugus, to systematically identify and analyze ERVs in this species. We mined an initial set of 362 potentially complete proviruses from the three main classes of ERVs, which were further resolved into 13 major families and 86 subfamilies by phylogenetic analysis. Consensus or representative sequences for each of the 86 subfamilies were then merged to the Repbase collection of known ERV/long terminal repeat (LTR) elements to annotate the retroviral complement of the bat genome. The results show that nearly 5% of the genome assembly is occupied by ERV-derived sequences, a quantity comparable to findings for other eutherian mammals. About one-fourth of these sequences belong to subfamilies newly identified in this study. Using two independent methods, intraelement LTR divergence and analysis of orthologous loci in two other bat species, we found that the vast majority of the potentially complete proviruses identified in M. lucifugus were integrated in the last ~25 million years. All three major ERV classes include recently integrated proviruses, suggesting that a wide diversity of retroviruses is still circulating in Myotis bats.