Identification of a functional envelope protein from the HERV-K family of human endogenous retroviruses

Genome-wide screening of sequence databases for human endogenous retroviruses (HERVs) has led to the identification of 18 coding env genes, among which two-the syncytin genes-encode fusogenic ENV proteins possibly involved in placenta physiology. Here we show that a third ENV, originating from the most "recent" HERV-K(HML2) family, is functional. Immunofluorescence analysis of env-transduced cells demonstrates expression of the protein at the cell surface, and we show that the protein confers infectivity to simian immunodeficiency virus pseudotypes. Western blot analysis of the pseudotyped virions further discloses the expected specific cleavage of the ENV precursor protein. This functional ENV could play a role in the amplification--via infection of the germ line--of the HERV-K genomic copies, all the more as coding HERV-K gag and pol genes can similarly be found in the human genome, which could therefore generate infectious virions of a fully endogenous origin.     

Identification, characterization, and comparative genomic distribution of the HERV-K (HML-2) group of human endogenous retroviruses

Background

Integration of retroviral DNA into a germ cell may lead to a provirus that is transmitted vertically to that host's offspring as an endogenous retrovirus (ERV). In humans, ERVs (HERVs) comprise about 8% of the genome, the vast majority of which are truncated and/or highly mutated and no longer encode functional genes. The most recently active retroviruses that integrated into the human germ line are members of the Betaretrovirus-like HERV-K (HML-2) group, many of which contain intact open reading frames (ORFs) in some or all genes, sometimes encoding functional proteins that are expressed in various tissues. Interestingly, this expression is upregulated in many tumors ranging from breast and ovarian tissues to lymphomas and melanomas, as well as schizophrenia, rheumatoid arthritis, and other disorders.

Results

No study to date has characterized all HML-2 elements in the genome, an essential step towards determining a possible functional role of HML-2 expression in disease. We present here the most comprehensive and accurate catalog of all full-length and partial HML-2 proviruses, as well as solo LTR elements, within the published human genome to date. Furthermore, we provide evidence for preferential maintenance of proviruses and solo LTR elements on gene-rich chromosomes of the human genome and in proximity to gene regions.

Conclusions

Our analysis has found and corrected several errors in the annotation of HML-2 elements in the human genome, including mislabeling of a newly identified group called HML-11. HML-elements have been implicated in a wide array of diseases, and characterization of these elements will play a fundamental role to understand the relationship between endogenous retrovirus expression and disease.     

Structural characterization of the fusion core in syncytin, envelope protein of human endogenous retrovirus family W

Syncytin is a captive retroviral envelope protein, possibly involved in the formation of the placental syncytiotrophoblast layer generated by trophoblast cell fusion at the maternal-fetal interface. We found that syncytin and type I viral envelope proteins shared similar structural profiling, especially in the regions of N- and C-terminal heptad repeats (NHR and CHR). We expressed the predicted regions of NHR (41aa) and CHR (34aa) in syncytin as a native single chain (named 2-helix protein) to characterize it. 2-Helix protein exists as a trimer and is highly alpha-helix, thermo-stable, and denatured by low pH. NHR and CHR could form a protease-resistant complex. The complex structure built by the molecular docking demonstrated that NHR and CHR associated in an antiparallel manner. Overall, the 2-helix protein could form a thermo-stable coiled coil trimer. The fusion core structure of syncytin was first demonstrated in endogenous retrovirus. These results support the explanation how syncytin mediates cytotrophoblast cell fusion involved in placental morphogenesis.     

Comparative sequence analysis and predictions for the envelope glycoproteins of insect endogenous retroviruses

Insect endogenous retroviruses (IERVs) are present in the genome of several species. Previous studies have shown a relationship between the envelope glycoproteins (Envs) and fusion proteins (FPs) of several baculoviruses. We used this sequence similarity to predict fusion domains in the Envs of IERVs. We suggest that FPs and Envs share several specific sequence and structural motifs with other RNA viruses in the viral transmembrane protein superfamily.     

Identification of a novel family of human endogenous retroviruses and characterization of one family member, HERV-K(C4), located in the complement C4 gene cluster.

We have identified a novel family of about 10-50 human endogenous retrovirus elements (HERVs) and have characterized one family member (HERV-KC4). This retrovirus element is integrated within intron 9 of and complement C4A genes and also in some C4B genes, and is a principal contribution to interlocus and interallelic length heterogeneity of C4 genes. The HERV-K(C4) sequence has a typical retrovirus structure with elements of gag, pol and env domains, flanked by two long terminal repeats (LTRs) and is similar to type A, B and D retroviruses. Multiple termination codons preclude the existence of long open reading frames, suggesting that the HERV-K(C4) sequence is no longer functional. Zoo blot hybridization reveals that New World monkeys appear to lack sequences similar to HERV-K(C4), suggesting that integration has occurred after the divergence of Old and New World monkeys. Retrotransposition of prototype viruses is presumed to have led to the amplification and integration of the members of the family in different loci, which in humans, appear to be dispersed over several chromosomes. The absence of the HERV-K(C4) element in some C4B genes in both humans and orangutangs indicate that the retrovirus inserted into the C4A gene after the duplication of the cluster. Subsequent spread of the HERV-K(C4) sequence to C4B genes presumably occurred by interlocus sequence exchange mechanisms, such as unequal crossover and gene conversion-like mechanisms.     

Envelope gene of the human endogenous retrovirus HERV-W encodes a functional retrovirus envelope

A member of the human endogenous retrovirus (HERV) family termed HERV-W encodes a highly fusogenic membrane glycoprotein that appears to be expressed specifically in the placenta. It is unclear whether the glycoproteins of the HERVs can serve as functional retrovirus envelope proteins to confer infectivity on retrovirus particles. We found that the HERV-W envelope glycoprotein can form pseudotypes with human immunodeficiency virus type 1 virions and confers tropism for CD4-negative cells. Thus, the HERV-W env gene represents the first HERV env gene demonstrated to encode the functional properties of a retrovirus envelope glycoprotein.     

A large number of the human microRNAs target lentiviruses, retroviruses, and endogenous retroviruses

Retroelements (including transposons, retrotransposons, retroviruses, and lentiviruses) make up a significant portion of eukaryotic genomes. Given their ability to mutate genes these mobile elements always present a threat to the integrity of the host genomes. Recent studies have revealed complex molecular mechanisms that silence the mutagenic ability of these RE as well strategically express the pieces of the incorporated RE that are utilized to silence human endogenous retroviruses (HERVs) or invading exogenous retroviruses (IERV). We have hypothesized that small endogenous RNA originally evolved to quell “foreign” IERV-genes and subsequently emerged into elaborate silencing systems that include RNA interference, miRNA-based gene regulation and other gene silencing mechanisms. Here, we present evidence that the replication of complex RE are most likely silenced or regulated by homologous miRNA that are found as a part of the cellular repertoire. We analyzed Homo sapiens miRNAs for possible target genetic sequences in selected HERVs and IERV found in humans and other large primates. We identified several miRNAs that have >80% sequence homology with human HERVs; -L, -W, and -K, and IERV like SIVcpz, HTLV-1, and HTLV-2. We found an inverse correlation between the numbers and relative degree of homology of miRNAs to the relative replication capacity of a specific RE. Therefore, larger numbers of miRNAs with greater degree of homology are found against the least active RE and the least numbers of miRNAs with smaller degree of homology are found against the most active RE (i.e. HERV-K). Implications of these observations in RE disease and therapy are discussed.     

The human immunodeficiency virus (HIV) type 2 envelope protein is a functional complement to HIV type 1 Vpu that enhances particle release of heterologous retroviruses.

We have recently shown that the envelope glycoprotein of the ROD10 isolate of human immunodeficiency virus type 2 (HIV-2) has the ability to positively regulate HIV-2 viral particle release. The activity provided by the ROD10 Env was remarkably similar to that of the HIV-1 Vpu protein, thus raising the possibility that the two proteins act in a related fashion. We now show that the ROD10 Env can functionally replace Vpu to enhance the rate of HIV-1 particle release. When provided in trans, both Vpu and the ROD10 Env restored wild-type levels of particle release in a Vpu-deficient mutant of the NL4-3 molecular clone with indistinguishable efficiencies. This effect was independent of the presence of the HIV-1 envelope protein. The ROD10 Env also enhanced HIV-1 particle release in the context of HIV-2 chimeric viruses containing the HIV-1 gag-pol, indicating a lack of need for additional HIV-1 products in this process. In addition, we show for the first time that HIV-1 Vpu, as well as ROD10 Env, has the ability to enhance simian immunodeficiency virus (SIV) particle release. The effects of Vpu and ROD10 Env on SIV particle release were indistinguishable and were observed in the context of full-length SIVmac239 and simian-human immunodeficiency virus chimeras. These results further demonstrate that ROD10 Env can functionally complement Vpu with respect to virus release. In contrast, we found no evidence of a destabilizing activity of ROD10 Env on the CD4 molecule. HIV-1 and HIV-2 thus appear to have evolved genetically distinct but functionally similar strategies to resolve the common problem of efficient release of progeny virus from infected cells.     

Phylogeny, sequence conservation, and functional complementation of the SBDS protein family

The Shwachman-Bodian-Diamond syndrome (SBDS) protein family occurs widely in nature, although its function has not been determined. Comprehensive database searches revealed SBDS homologues from 159 species, including examples from all sequenced archaeal and eukaryotic genomes and all eukaryotic kingdoms. Sequence alignment with ClustalX and MUSCLE algorithms led to the identification of conserved residues that occurred predominantly in the amino-terminal FYSH domain where they appeared to contribute to protein folding or stability. Only SBDS residue Gly91 was invariant in all species. Four distantly related protists were found to have two divergent SBDS genes in their genomes. In each case, phylogenetic analyses and the identification of shared sequence features suggested that one gene was derived from lateral gene transfer. We also identified a shared C-terminal zinc finger domain fusion in flowering plants and chromalveolates that may shed light on the function of the protein family and the evolutionary histories of these kingdoms. To assess the extent of SBDS functional conservation, we carried out complementation studies of SBDS homologues and interspecies chimeras in Saccharomyces cerevisiae. We determined that the FYSH domain was widely interchangeable among eukaryotes, while domain 2 imparted species specificity to protein function. Domain 3 was largely dispensable for function in our yeast complementation assay. Overall, the phylogeny of SBDS was shared with a group of proteins that were markedly enriched for RNA metabolism and/or ribosome-associated functions. These findings link Shwachman-Diamond syndrome to other bone marrow failure syndromes with defects in nucleolus-associated processes, including Diamond-Blackfan anemia, cartilage-hair hypoplasia, and dyskeratosis congenita.     

Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches

Based on integrative computational hybrid approaches that combined statistical coupling analysis (SCA), molecular dynamics (MD), and normal mode analysis (NMA), evolutionarily coupled residues involved in functionally relevant motion in the adenylate kinase protein family were identified. The hybrids identified four top-ranking site pairs that belong to a conserved hydrogen bond network that is involved in the enzyme's flexibility. A second group of top-ranking site pairs was identified in critical regions for functional dynamics, such as those related to enzymatic turnover. The high consistency of the results obtained by SCA with NMA (SCA.NMA) and by SCA.MD hybrid analyses suggests that suitable replacement of the matrix of cross-correlation analysis of atomic fluctuations (derived by using NMA) with those based on MD contributes to the identification of such sites by means of a fast computational calculation. The analysis presented here strongly supports the hypothesis that evolutionary forces, such as coevolution at the sequence level, have promoted functional dynamic properties of the adenylate kinase protein family. Finally, these hybrid approaches can be used to identify, at the residue level, protein motion coordination patterns not previously observed, such as in hinge regions.     

The genomic locus of the human hemopoietic-specific cell protein tyrosine kinase (PTK)-encoding gene (HCK) confirms conservation of exon-intron structure among human PTKs of the src family.

Protein tyrosine kinases (PTKs) are implicated in the control of cell growth by virtue of their frequent appearance as products of retroviral oncogenes, as intracellular signal transducers, and as growth factor receptors or components thereof. The knowledge of the structure and sequence of family genes encoding PTKs is still limited. To date, the complete genomic structure of human src family members is only available for the C-FGR gene (encoding p55 Fgr, PTK). Sequence analysis and characterization of the intron/exon organization of the human HCK gene, encoding a hemopoietic-specific cell PTK of the src-related family, revealed a length of over 16 kb for the seven 3'-exons. All intron/exon splice junctions agree with the GT/AG rule. In each case where a boundary occurs at a Gly codon, GGG or GGA, the triplet is split between the first and second nucleotide (nt). A total of eight complete and one partial Alu repeats were identified within the introns. The nt sequence of the genomic clones resolves existing discrepancies among two published sequences of HCK cDNAs. Human HCK, C-SRC (encoding p60 Src PTK), C-FGR and LCK (encoding p56 Lck, PTK) genes thus share very similar exon/intron structures for the conserved exons. These results provide additional evidence that the different PTKs of the src-like family most likely arose by duplication of an ancestral src-like gene.     

Identification and functional characterization of a new member of the human Mcm protein family: hMcm8

The six minichromosome maintenance proteins (Mcm2–7) are required for both the initiation and elongation of chromosomal DNA, ensuring that DNA replication takes place once, and only once, during the S phase. Here we report on the cloning of a new human Mcm gene (hMcm8) and on characterisation of its protein product. The hMcm8 gene contains the central Mcm domain conserved in the Mcm2–7 gene family, and is expressed in a range of cell lines and human tissues. hMcm8 mRNA accumulates during G₁/S phase, while hMcm8 protein is detectable throughout the cell cycle. Immunoprecipitation-based studies did not reveal any participation of hMcm8 in the Mcm3/5 and Mcm2/4/6/7 subcomplexes. hMcm8 localises to the nucleus, although it is devoid of a nuclear localisation signal, suggesting that it binds to a nuclear protein. In the nucleus, the hMcm8 structure-bound fraction is detectable in S, but not in G₂/M, phase, as for hMcm3. However, unlike hMcm3, the hMcm8 structure-bound fraction is not detectable in G₁ phase. Overall, our data identify a new Mcm protein, which does not form part of the Mcm2–7 complex and which is only structure-bound during S phase, thus suggesting its specific role in DNA replication.     

Identification and functional characterization of the BAG protein family in Arabidopsis thaliana

The genes that control mammalian programmed cell death are conserved across wide evolutionary distances. Although plant cells can undergo apoptosis-like cell death, plant homologs of mammalian regulators of apoptosis have, in general, not been found. This is in part due to the lack of primary sequence conservation between animal and putative plant regulators of apoptosis. Thus, alternative approaches beyond sequence similarities are required to find functional plant homologs of apoptosis regulators. Here, we present the results of using advanced bioinformatic tools to uncover the Arabidopsis family of BAG proteins. The mammalian BAG (Bcl-2-associated athanogene) proteins are a family of chaperone regulators that modulate a number of diverse processes ranging from proliferation to growth arrest and cell death. Such proteins are distinguished by a conserved BAG domain that directly interacts with Hsp70 and Hsc70 proteins to regulate their activity. Our searches of the Arabidopsis thaliana genome sequence revealed seven homologs of the BAG protein family. We further show that plant BAG family members are also multifunctional and remarkably similar to their animal counterparts, as they regulate apoptosis-like processes ranging from pathogen attack to abiotic stress and development.