Molecular characterization of the HERV-W env gene in humans and primates: expression, FISH, phylogeny, and evolution

We characterized the human endogenous retrovirus (HERV-W) family in humans and primates. In silico expression data indicated that 22 complete HERV-W families from human chromosomes 1-3, 5-8, 10-12, 15, 19, and X are randomly expressed in various tissues. Quantitative real-time RT-PCR analysis of the HERV-W env gene derived from human chromosome 7q21.2 indicated predominant expression in the human placenta. Several copies of repeat sequences (SINE, LINE, LTR, simple repeat) were detected within the complete or processed pseudo HERV-W of the human, chimpanzee, and rhesus monkey. Compared to other regions (5'LTR, Gag, Gag-Pol, Env, 3'LTR), the repeat family has been mainly integrated into the region spanning the 5'LTRs of Gag (1398 bp) and Pol (3242 bp). FISH detected the HERV-W probe (fosWE1) derived from a gorilla fosmid library in the metaphase chromosomes of all primates (five hominoids, three Old World monkeys, two New World monkeys, and one prosimian), but not in Tupaia. This finding was supported by molecular clock and phylogeny data using the divergence values of the complete HERV-W LTR elements. The data suggested that the HERV-W family was integrated into the primate genome approximately 63 million years (Myr) ago, and evolved independently during the course of primate radiation.     

Molecular evolution of a primate-specific microRNA family

Lineage-specific microRNA (miRNA) families may contribute to developmental novelties during evolution. However, little is known about the origin and evolution of new miRNA families. We report evidence of an Alu-mediated rapid expansion of miRNA genes in a previously identified primate-specific miRNA family, drawn from sequencing and comparative analysis of 9 diverse primate species. Evolutionary analysis reveals similar divergence among miRNA copies whether they are within or between species, lineage-specific gain and loss of miRNAs, and gene pseudolization in multiple species. These observations support a birth-and-death process of miRNA genes in this family, implicating functional diversification during primate evolution. In addition, both secondary structure conservation and reduced single nucleotide polymorphisms density attest to functional constraint of this family in primates. Finally, we observed preferential expression of miRNAs in human placenta and fetal brain, suggesting a functional importance of this family for primate development.     

Human endogenous retroviral elements as indicators of ectopic recombination events in the primate genome

HERV elements make up a significant fraction of the human genome and, as interspersed repetitive elements, have the capacity to provide substrates for ectopic recombination and gene conversion events. To understand the extent to which these events occur and gain further insight into the complex evolutionary history of these elements in our genome, we undertook a phylogenetic study of the long terminal repeat sequences of 15 HERV-K(HML-2) elements in various primate species. This family of human endogenous retroviruses first entered the primate genome between 35 and 45 million years ago. Throughout primate evolution, these elements have undergone bursts of amplification. From this analysis, which is the largest-scale study of HERV sequence dynamics during primate evolution to date, we were able to detect intraelement gene conversion and recombination at five HERV-K loci. We also found evidence for replacement of an ancient element by another HERV-K provirus, apparently reflecting an occurrence of retroviral integration by homologous recombination. The high frequency of these events casts doubt on the accuracy of integration time estimates based only on divergence between retroelement LTRs.     

The monosaccharide transporter gene family in Arabidopsis and rice: a history of duplications, adaptive evolution, and functional divergence

Current hypotheses of gene duplicate divergence propose that surviving members of a gene duplicate pair may evolve, under conditions of purifying or nearly neutral selection, in one of two ways: with new function arising in one duplicate while the other retains original function (neofunctionalization [NF]) or partitioning of the original function between the 2 paralogs (subfunctionalization [SF]). More recent studies propose that SF followed by NF (subneofunctionalization [SNF]) explains the divergence of many duplicate genes. In this analysis, we evaluate these hypotheses in the context of the large monosaccharide transporter (MST) gene families in Arabidopsis and rice. MSTs have an ancient origin, predating plants, and have evolved in the seed plant lineage to comprise 7 subfamilies. In Arabidopsis, 53 putative MST genes have been identified, with one subfamily greatly expanded by tandem gene duplications. We searched the rice genome for members of the MST gene family and compared them with the MST gene family in Arabidopsis to determine subfamily expansion patterns and estimate gene duplicate divergence times. We tested hypotheses of gene duplicate divergence in 24 paralog pairs by comparing protein sequence divergence rates, estimating positive selection on codon sites, and analyzing tissue expression patterns. Results reveal the MST gene family to be significantly larger (65) in rice with 2 subfamilies greatly expanded by tandem duplications. Gene duplicate divergence time estimates indicate that early diversification of most subfamilies occurred in the Proterozoic (2500-540 Myr) and that expansion of large subfamilies continued through the Cenozoic (65-0 Myr). Two-thirds of paralog pairs show statistically symmetric rates of sequence evolution, most consistent with the SF model, with half of those showing evidence for positive selection in one or both genes. Among 8 paralog pairs showing asymmetric divergence rates, most consistent with the NF model, nearly half show evidence of positive selection. Positive selection does not appear in any duplicate pairs younger than approximately 34 Myr. Our data suggest that the NF, SF, and SNF models describe different outcomes along a continuum of divergence resulting from initial conditions of relaxed constraint after duplication.     

Concerted evolution of the tandem array encoding primate U2 snRNA (the RNU2 locus) is accompanied by dramatic remodeling of the junctions with flanking chromosomal sequences.

The genes encoding primate U2 snRNA are organized as a nearly perfect tandem array (the RNU2 locus) that has been evolving concertedly for >35 Myr since the divergence of baboons and humans. Thus the repeat units of the tandem array are essentially identical within each species, but differ between species. Homogeneity is maintained because any change in one repeat unit is purged from the array or fixed in all other repeats. Intriguingly, the cytological location of RNU2 has remained unchanged despite concerted evolution of the tandem array. We had found previously that junction sequences between the U2 tandem array and flanking DNA were subject to remodeling over a region of 200-300 bp during the past 5 Myr in the hominid lineage. Here we show that the junctions between the U2 tandem array and flanking DNA have undergone dramatic rearrangements over a region of 1 to >10 kbp in the 35 Myr since divergence of the Old World Monkey and hominid lineages. We argue that these rearrangements reflect the high level of genetic activity required to sustain concerted evolution, and propose a model to explain why maintenance of homogeneity within a tandemly repeated multigene family would lead to junctional diversity.     

Monophyletic Origin of Alu Elements in Primates

To get insight into the early evolution of the primate Alu elements, we characterized sequences of these repeats from the Malagasy prosimians, lemurs (Lemuridae) and sifakas (Indriidae), as well as from galagos (Lorisidae). These sequences were compared with the oldest Alu species known from the human genome: dimeric Alu J and S and free Alu monomers. Our analysis indicates that about 60 Myr ago, before the prosimian divergence, free left and right monomers formed an Alu heterodimer connected by a 19-nucleotide-long A-rich linker. The resulting elements successfully propagated in diverging primate lineages until about ∼20 Myr ago, conserving similar sequence features and essentially the same Alu RNA secondary structure. We suggest that until that time the same ``retropositional niche', molecular machinery making possible the proliferation by retroposition, constrained the evolution of Alu elements in extant primate species. These constraints became subsequently relaxed. In the Malagasy prosimians the dimeric Alu continued to amplify after acquiring a 34- to 36-nucleotide extension of their linker segment, whereas in the galago genome the ``retropositional niche' was occupied by novel short elements. Received: 1 December 1997 / Accepted: 30 January 1998     

Structure, distribution, and expression of an ancient murine endogenous retroviruslike DNA family.

An endogenous retroviruslike DNA, B-26, was cloned from a BALB/c mouse embryo gene library by using a generalized murine leukemia virus DNA probe. Southern blot hybridization and nucleotide sequence analyses indicated that B-26 DNA might be a novel member of the GLN DNA family (A. Itin and E. Keshet, J. Virol. 59:301-307, 1986) which contains murine leukemia virus-related pol and env sequences. Northern analysis indicated that B-26-related RNAs of 8.4 and 3.0 kilobases were transcribed in thymus, spleen, brain, and liver tissues of 6-week-old BALB/c mice.     

The Molecular Clock Runs at Different Rates Among Closely Related Members of a Gene Family

The serum albumin gene family is composed of four members that have arisen by a series of duplications from a common ancestor. From sequence differences between members of the gene family, we infer that a gene duplication some 580 Myr ago gave rise to the vitamin D–binding protein (DBP) gene and a second lineage, which reduplicated about 295 Myr ago to give the albumin (ALB) gene and a common precursor to α-fetoprotein (AFP) and α-albumin (ALF). This precursor itself duplicated about 250 Myr ago, giving rise to the youngest family members, AFP and ALF. It should be possible to correlate these dates with the phylogenetic distribution of members of the gene family among different species. All four genes are found in mammals, but AFP and ALF are not found in amphibia, which diverged from reptiles about 360 Myr ago, before the divergence of the AFP-ALF progenitor from albumin. Although individual family members display an approximate clock-like evolution, there are significant deviations—the rates of divergence for AFP differ by a factor of 7, the rates for ALB differ by a factor of 2.1. Since the progenitor of this gene family itself arose by triplication of a smaller gene, the rates of evolution of individual domains were also calculated and were shown to vary within and between family members. The great variation in the rates of the molecular clock raises questions concerning whether it can be used to infer evolutionary time from contemporary sequence differences. Received: 28 February 1995 / Accepted: 6 October 1997     

Comparative structure and evolution of murine CR2. The homolog of the human C3d/EBV receptor (CD21)

The complete nucleotide sequence of murine complement receptor type 2 (CR2) was determined from two overlapping cDNA clones derived from a lambda gt11 library of late pre-B cell origin. Comparison of the predicted sequence of the 1014 amino acid murine homolog with that of human CR2 revealed marked evolutionary conservation. The murine molecule was 65% identical to human CR2 overall, lacking a single repetitive sequence variably present in man. The 15 approximately 60-75 amino acid short consensus repeats (SCR) that constitute the entire extracellular domain of murine CR2 were 53 to 81% identical to and could be directly aligned with the human protein. As reported, the cytoplasmic tail shared 79% amino acid identity with human CR2, whereas that of the transmembrane was only 33%. Murine CR2 contained 16 potential N-linked glycosylation sites of which 6 were conserved, 4 altered, and 6 lost during human evolution. The hydropathicity profile of the two molecules was nearly colinear with some variation in the N-terminal region of the first repeat, as well as within the sixth and twelfth repeats. RNA blot analysis revealed a approximately 4.0 to 5.0 kb message in murine B lymphocytes, which was absent in T lymphocytes (thymus and spleen), liver, brain, lung, kidney, and heart. A method was devised to more precisely compare the repeat structures. An identity matrix analysis suggests that human ancestral CR2 evolved before divergence of the rodent and primate branches of the evolutionary tree through a series of predictable gene duplications, possibly giving rise to the precursor of human CR1 and murine CRY. The marked structural similarity between the human and murine receptors suggests functional conservation as well.     

Rabbit globin pseudogene psi beta 2 is a hybrid of delta- and beta-globin gene sequences

The evolutionary history of the rabbit globin pseudogene psi beta 2 was studied by completing its nucleotide sequence and aligning the sequence with that of the rabbit adult globin gene beta 1 and the human minor adult globin gene delta. The 5' flanking region and exon 1 of psi beta 2 were most similar to rabbit beta 1, but the large intervening sequence and the 3' untranslated region were most similar to human delta. Intron 1 and exon 2 were equally similar to both delta and beta 1. This pattern indicates that psi beta 2 was originally a delta-like gene that acquired the 5' portion of gene beta 1 by intrachromosomal gene conversion. The presence of a delta-globin gene sequence in both rabbits and humans shows that it is an ancient gene, predating the mammalian radiation that occurred over 85 Myr ago. Delta has shown a pronounced tendency to be altered in its 5' end during the course of mammalian evolution. Quantitative divergence analysis shows that the ancestor to rabbit psi beta 2 was active until 20-30 Myr ago, during which time the lagomorph beta-globin gene family apparently functioned without a pseudogene.      

Expression and identification of HERV-W family in Japanese monkeys (Macaca fuscata)

We investigated structural genes (gag, pol, env) of HERV-W family in the Macaca fuscata (Japanese monkey). Those genes are expressed in various tissues (testis, prostate, kidney, cerebellum, thymus, pancreas, intestine, stomach, ovary) of the Japanese monkey in RT-PCR and sequencing analyses. Nine clones for gag, thirty-one clones for pol and thirty-four clones for env fragments of the HERV-W family in monkey tissues were identified and analyzed. These clones showed a high degree of sequence similarity, 82.2-84.7% for gag, 88.4-91.7% for pol, and 90.8-95.4% for env, to those of HERV-W family. Translation to amino acids in all clones derived from the monkey indicated that they showed multiple interruptions of frameshifts and termination codons by deletion/insertion or point mutation. Identical sequences from different tissues of the monkey were found in env and pol clones of the HERV-W family.     

The approximately 30-million-year-old ERVPb1 envelope gene is evolutionarily conserved among hominoids and Old World monkeys

Most human endogenous retroviruses (HERVs) are ancient and their genes are rendered nonfunctional by debilitating mutations. One exception is a recently discovered envelope gene located on chromosome 14. This envelope protein was also recently shown to be expressed in various human tissues and to mediate cell-cell fusion ex vivo. In this study, we demonstrate that this locus (designated ERVPb1) is preserved in Old World monkeys and that the reading frame is maintained. This is congruent with the entry of the HERV-P(b) group between 27 and 36 million years ago as suggested by long terminal repeat divergence. Although the coding capacity is generally lost in the HERV-IP supergroup, the analysis of nucleotide substitutions, lack of stop codons, and single-nucleotide polymorephisms strongly indicates a selective advantage of the ERVPb1 envelope genes during primate evolution. The purifying selection and tissue-specific expression of the human ERVPb1 envelope gene provide strong evidence of a beneficial role for the host.     

The human endogenous retrovirus family HERV-K(HML-3)

A substantial amount of the human genome is composed of human endogenous retroviruses (HERVs). Manifold HERV families have been identified, among them several so-called HERV-K(HML) families. Although the HERV-K(HML-2) family has been studied in detail, other HERV-K families are not as well characterized. We describe here the HERV-K HML-3 family in more detail. We estimate that there are about 140 proviral loci or remains of such per haploid genome. Most loci are severely mutated. Proviruses displaying larger deletions in gag and pol are common. A multiple alignment of 73 HERV-K(HML-3) sequences displays several potentially important differences compared with the HERVK9I sequence in Repbase. A consensus sequence with open reading frames for all retroviral genes was generated, for which intact dUTPase motifs and env gene variants with different coding capacities are observed. Phylogenetic analysis shows near-monophyly with distinction of two closely related subgroups. Proviruses formed about 36 million years ago. However, no continuous activity through primate evolution is indicated.     

Human leukocyte formin: a novel protein expressed in lymphoid malignancies and associated with Akt

The very large family of Formin proteins is involved in processes such as morphogenesis, embryonic differentiation, cell polarity, and cytokinesis. A novel human gene from the Formin family, denominated human leukocyte formin gene, was cloned. The cDNA of the gene was determined to be 3959bp long with an open reading frame of 3302bp and computational analysis located this gene on chromosome 17, suggesting that it is composed of 27 exons. Northern blot analysis revealed a restricted expression of mRNA in the thymus, spleen, and peripheral blood leukocytes in normal human tissues. Western blot analysis demonstrated that the protein encoded by this gene is overexpressed in lymphoid malignancies; cancer cell lines and peripheral blood leukocyte from chronic lymphocytic leukemia (CLL) patients. Furthermore, the human leukocyte formin protein was observed to associate with Akt, a critical survival regulator in many different cell types.     

Evolutionary dynamics of the T-cell receptor VB gene family as inferred from the human and mouse genomic sequences

The diversity of T-cell receptors is generated primarily by the variable-region gene families, each of which is composed of a large number of member genes. The entire genomic sequence of the variable region (VB) of the T- cell receptor beta chain from humans and mice has become available. To understand the evolutionary dynamics of the VB gene family, we conducted a phylogenetic analysis of all VB genes from humans and mice, as well as a detailed analysis of internal DNA duplications in the human genomic VB region. The phylogenetic tree obtained shows that human and mouse VB genes intermingle extensively rather than forming two separate clusters and that many gene duplications occurred both before and after the divergence between primates and rodents. Analyzing the genomic maps of transposable elements (e.g., LINEs and SINEs) and relic VB genes in the VB gene region, we present evidence that a 20-kb VB region duplicated tandemly four times in the human lineage during the last 32 Myr, and 6 out of the 15 VB genes in this region have become nonfunctional during this period. Our results show that the VB gene family is subject to evolution by a birth-and-death process rather than to concerted evolution.