Identification and classification of endogenous retroviruses in cattle

The aim of this study was to identify the endogenous retrovirus (ERV) sequences in a bovine genome. We subjected bovine genomic DNA to PCR with degenerate or ovine ERV (OERV) family-specific primers that aimed to amplify the retroviral pro/pol region. Sequence analysis of 113 clones obtained by PCR revealed that 69 were of retroviral origin. On the basis of the OERV classification system, these clones from degenerate PCR could be divided into the beta3, gamma4, and gamma9 families. PCR with OERV family-specific primers revealed an additional ERV that was classified into the bovine endogenous retrovirus (BERV) gamma7 family. In conclusion, here we report the results of a genome scale study of the BERV. Our study shows that the ERV family expansion in cattle may be somewhat limited, while more diverse family members of ERVs have been reported from other artiodactyls, such as pigs and sheep.     

Characterization of porcine endogenous retrovirus gamma pro-pol nucleotide sequences

Endogenous retroviral sequences in the pig genome (PERV) represent a potential infectious risk in xenotransplantation. All known infectious PERV have been asssigned to the PERV gamma1 family, consisting of the subfamilies A, B, and C. The aim of the study was the concise examination of PERV gamma by the analysis of the retroviral pro-pol sequences. The analysis of 52 pro-pol clones amplified in this study revealed eight PERV gamma families. In addition to four already-described families (gamma1, gamma4, gamma5, gamma6), four novel families (gamma7, gamma8, gamma9, gamma10) were identified. Quantitative analysis of the novel PERV gamma sequences in selected breeds revealed variations in the endogenous retroviral load. Open reading frames (ORF) in the amplified proviral fragment were only found for PERV gamma1. In addition, novel ORF-containing PERV gamma1 clones consisting of hybrid sequences were revealed. Sequence comparison from published full-length PERV gamma1 clones of the PERV subfamilies A, B, and C resulted in a lack of strict correlation of the classification of pro-pol and env. The results indicated the occurrence of causative recombination events between retroviral genomes. Thus, our study on PERV gamma provides new data for the evaluation and selection of pigs intended to be used in xenotransplantation.     

Characterization of the bovine endogenous retrovirus beta3 genome

We recently used degenerate PCR and locus-specific PCR methods to identify the endogenous retroviruses (ERV) in the bovine genome. Using the ovine ERV classification system, the bovine ERVs (BERVs) could be classified into four families. Here, we searched the most recently released bovine genome database with the partial nucleotide sequence of the pro/pol region of the BERV beta3 family. This allowed us to obtain and analyze the complete genome of BERV beta3. The BERV beta3 genome is 7666 nucleotides long and has the typical retroviral organization, namely, 5'-long terminal repeat (LTR)-gag-pro-pol-env-LTR-3'. The deduced open reading frames for gag, pro, pol and env of BERV Beta en- code 507, 271, 879 and 603 amino acids, respectively. BERV beta3 showed little amino acid similarity to other betaretroviruses. Phylogenetic analysis showed that it clusters with HERV-K. This is the first report describing the genetic structure and sequence of an entire BERV.     

Structural characterization of the genome of BERV gamma4, the most abundant endogenous retrovirus family in cattle

The genome of replication-competent BERV gamma4 provirus, which is the most abundant ERV family in the bovine genome, was characterized in detail. The BERV gamma4 genome showed that BERV gamma4 harbors 8576 nucleotides and has the typical 5'-long terminal repeat (LTR)-gag-pro-pol-env-LTR-3' retroviral organization with a long leader region positioned before the gag open reading frame. Multiple sequences analysis showed that the nucleotide difference between 5' and 3' LTRs was 4.2% (mean value 0.042) in average, suggesting that the provirus formed at most 13.3 million years ago. Gag separated by a stop codon from pro-pol in the same reading frame, while env resides in another reading frame lacking of a functional surface domain. According to the current bovine genome sequence assembly, the full-length BERV gamma4 provirus sequences were only found in the chromosomes 1, 2, 6, 10, 15, 23, 26, 28, X, and unassigned, although the partial sequences almost evenly distributed in the entire bovine genome. This is the first detailed study describing the genome structure of BERV gamma4, the most abundant ERV family present in bovine genome. Combined with our recent reports on characterization of ERVs in bovine, this study will contribute to illuminate ERVs in the cattle of which no information was previously available.     

ERV3 and related sequences in humans: structure and RNA expression

The ERV3 locus at chromosome 7q11 is a much studied human endogenous retroviral (HERV) sequence, owing to an env open reading frame (ORF) and placental RNA and protein expression. An analysis of the human genome demonstrated that ERV3 is one of a group of 41 highly related elements (ERV3-like HERVs) which use proline, isoleucine, or arginine tRNA in their primer binding sites. In addition to elements closely related to ERV3, the group included the previously known retinoic acid-inducible element, RRHERVI, also referred to as HERV15, but was separate from the related HERV-E elements. The ERV3-like elements are defective. The only element with an ORF among gag, pro, pol, and env genes was the env ORF of the original ERV3 locus. A search in dbEST revealed ERV3 RNA expression in placenta, skin, carcinoid tumor, and adrenal glands. Expression was also studied with newly developed real-time quantitative PCRs (QPCR) of ERV3 and HERV-E(4-1) env sequences. Results from a novel histone 3.3 RNA QPCR result served as the expression control. QPCR results for ERV3 were compatible with previously published results, with a stronger expression in adrenal gland and placenta than in 15 other human tissues. The expression of the envelope (env) of ERV3 at chromosome 7q11 was also studied by using stringent in situ hybridization. Expression was found in corpus luteum, testis, adrenal gland, Hassal's bodies in thymus, brown fat, pituitary gland, and epithelium of the lung. We conclude that ERV3 env is most strongly expressed in adrenal and sebaceous glands as well as in placenta.     

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.     

Chromosomal distribution, localization and expression of the human endogenous retrovirus ERV9

ERV9 is a class I family of human endogenous retroviral sequences. Somatic cell hybrid genomic hybridization experiments using a mono-chromosomal panel indicate the presence of approximately 120 ERV9 loci in the human genome distributed on most chromosomes. Fluorescence in situ hybridization (FISH) using an ERV9 cDNA probe containing gag, pol and env sequences, verified this observation and a consistent signal was found at the chromosome region 11q13.3-->q13.5. By analysis of a panel of radiation hybrids, an ERV9 locus was mapped to within a 300-kbp region at the chromosome site 11q13. The marker cCLGW567 and the locus MAP3K11/D11S546 centromeric and telomeric flanked it, respectively. Northern blot analysis, using an ERV9 LTR probe, indicated that most normal tissues examined expressed low abundant ERV9 LTR driven mRNAs of various sizes. The most prominent expression was found in adrenal glands and testis. However, the level of expression varied in the same tissues among different individuals indicating that ERV9 mRNA expression probably is inducible in certain tissues or at various cell stages.     

Multiple groups of novel retroviral genomes in pigs and related species

In view of the concern over potential infection hazards in the use of porcine tissues and organs for xenotransplantation to humans, we investigated the diversity of porcine endogenous retrovirus (PERV) genomes in the DNA of domestic pigs and related species. In addition to the three known envelope subgroups of infectious gamma retroviruses (PERV-A, -B, and -C), classed together here as PERV group gamma 1, four novel groups of gamma retrovirus (gamma 2 to gamma 5) and four novel groups of beta retrovirus (beta 1 to beta 4) genomes were detected in pig DNA using generic and specific PCR primers. PCR quantification indicated that the retroviral genome copy number in the Landrace x Duroc F(1) hybrid pig ranged from 2 (beta 2 and gamma 5) to approximately 50 (gamma 1). The gamma 1, gamma 2, and beta 4 genomes were transcribed into RNA in adult kidney tissue. Apart from gamma 1, the retroviral genomes are not known to be infectious, and sequencing of a small number of amplified genome fragments revealed stop codons in putative open reading frames in several cases. Analysis of DNA from wild boar and other species of Old World pigs (Suidae) and New World peccaries (Tayassuidae) showed that one retrovirus group, beta 2, was common to all species tested, while the others were present among all Old World species but absent from New World species. The PERV-C subgroup of gamma1 genomes segregated among domestic pigs and were absent from two African species (red river hog and warthog). Thus domestic swine and their phylogenetic relatives harbor multiple groups of hitherto undescribed PERV genomes.     

Genome-wide assessments reveal extremely high levels of polymorphism of two active families of mouse endogenous retroviral elements

Endogenous retroviral elements (ERVs) in mice are significant genomic mutagens, causing ~10% of all reported spontaneous germ line mutations in laboratory strains. The majority of these mutations are due to insertions of two high copy ERV families, the IAP and ETn/MusD elements. This significant level of ongoing retrotranspositional activity suggests that inbred mice are highly variable in content of these two ERV groups. However, no comprehensive genome-wide studies have been performed to assess their level of polymorphism. Here we compared three test strains, for which sufficient genomic sequence is available, to each other and to the reference C57BL/6J genome and detected very high levels of insertional polymorphism for both ERV families, with an estimated false discovery rate of only 0.4%. Specifically, we found that at least 60% of IAP and 25% of ETn/MusD elements detected in any strain are absent in one or more of the other three strains. The polymorphic nature of a set of 40 ETn/MusD elements found within gene introns was confirmed using genomic PCR on DNA from a panel of mouse strains. For some cases, we detected gene-splicing abnormalities involving the ERV and obtained additional evidence for decreased gene expression in strains carrying the insertion. In total, we identified nearly 700 polymorphic IAP or ETn/MusD ERVs or solitary LTRs that reside in gene introns, providing potential candidates that may contribute to gene expression differences among strains. These extreme levels of polymorphism suggest that ERV insertions play a significant role in genetic drift of mouse lines.     

A Genetic Linkage Map of the Male Goat Genome

This paper presents a first genetic linkage map of the goat genome. Primers derived from the flanking sequences of 612 bovine, ovine and goat microsatellite markers were gathered and tested for amplification with goat DNA under standardized PCR conditions. This screen made it possible to choose a set of 55 polymorphic markers that can be used in the three species and to define a panel of 223 microsatellites suitable for the goat. Twelve half-sib paternal goat families were then used to build a linkage map of the goat genome. The linkage analysis made it possible to construct a meiotic map covering 2300 cM, i.e., >80% of the total estimated length of the goat genome. Moreover, eight cosmids containing microsatellites were mapped by fluorescence in situ hybridization in goat and sheep. Together with 11 microsatellite-containing cosmids previously mapped in cattle (and supposing conservation of the banding pattern between this species and the goat) and data from the sheep map, these results made the orientation of 15 linkage groups possible. Furthermore, 12 coding sequences were mapped either genetically or physically, providing useful data for comparative mapping.     

Identification of Allelic Polymorphism in the Ovine Leptin Gene

Leptin is an adipocyte-derived hormone/cytokine that influences the physiological control of numerous biological functions and links nutritional status with both neuroendocrine and immune functions. In livestock, variation in the leptin (LEP) gene has been characterized in cattle and pig, but it has not been reported in sheep. In this study, variation in the exon 3 coding sequence of the ovine LEP gene was investigated by polymerase chain reaction–single-strand conformational polymorphism (PCR–SSCP) analysis and DNA sequencing. Five novel SSCP patterns, representing five different sequences, were identified under a combination of two different electrophoresis conditions. Either one or two different sequences were detected in individual sheep and all the sequences identified shared high homology with the LEP sequences from a variety of species, suggesting that these sequences represent alleles of the ovine LEP gene. Four single nucleotide polymorphisms (SNPs) were detected, and three of these resulted in amino acid changes. Variation detected here might have an impact on leptin activity and function.