Intra- and interspecific variation of the CCR5 gene in higher primates

We have evaluated the molecular evolution of the chemokine receptor CCR5 in primates. The chemokine receptor CCR5 serves as a major co-receptor for human immunodeficiency virus/simian immunodeficiency virus (HIV/SIV) infection. Knowledge of evolution of the CCR5 molecule and selection on the CCR5 gene may shed light on its functional role. The comparison of differences between intraspecific polymorphisms and interspecific fixed substitutions provides useful information regarding modes of selection during the course of evolution. There is marked polymorphism in the CCR5 gene sequence within different primate species, whereas sequence divergence between different species is small. By using contingency tests, we compared synonymous (SS) and nonsynonymous (NS) CCR5 mutations occurring within and between a broad range of primates. Our results demonstrate that CCR5 evolution did not follow expectations of strict neutrality at the level of the whole gene. The proportion of NS to SS at the intraspecific level was significantly higher than that observed at the interspecific level. These results suggest that most CCR5 NS polymorphisms are slightly deleterious. However, at domains more closely correlated with its known biological functions, there was no obvious evidence to support deviation from neutrality.     

Evolutionary dynamics of duplicated microsatellites shared by sex chromosomes

Segmental duplications on sex chromosomes constitute an important proportion of recent duplications (approximately 30%). Among those, the evolution of duplicated noncoding DNA is still poorly investigated. We focus our work on repeated DNA sequences extensively used in population genetics and evolution: microsatellites. Six duplicated (CA), microsatellite loci, located on the homologous region of human sex chromosomes, were studied at the intraspecific level in Homo sapiens and by an orthologous comparison in eight primate species. At the intraspecific level, we evaluated the congruence in paralogous divergence between the flanking sequences of the six microsatellites and the approximately 2.2-kb surrounding sequences and observed that both phylogenies are congruent. At the interspecific level (8 species of primates: 54 individuals), we analyzed the sequence polymorphism and divergence of each orthologous locus for both the flanking sequence and the microsatellite. The results showed a lower divergence of flanking sequences than expected in noncoding DNA and a relative stability of the first nucleotides close to the microsatellite. The location of each CAIII locus in a Low Copy Repeated element containing duplicated VCX/Y genes (approximately 1 kb) suggested that direct or indirect selection could explain these results. Moreover, the substitution rates in the flanking sequences and in the microsatellites were correlated. Thus, the evolutionary dynamics of microsatellites seems closely linked to the variation of spontaneous mutations in the surrounding regions.     

A comprehensive analysis of recently integrated human Ta L1 elements

The Ta (transcribed, subset a) subfamily of L1 LINEs (long interspersed elements) is characterized by a 3-bp ACA sequence in the 3' untranslated region and contains approximately 520 members in the human genome. Here, we have extracted 468 Ta L1Hs (L1 human specific) elements from the draft human genomic sequence and screened individual elements using polymerase-chain-reaction (PCR) assays to determine their phylogenetic origin and levels of human genomic diversity. One hundred twenty-four of the elements amenable to complete sequence analysis were full length ( approximately 6 kb) and have apparently escaped any 5' truncation. Forty-four of these full-length elements have two intact open reading frames and may be capable of retrotransposition. Sequence analysis of the Ta L1 elements showed a low level of nucleotide divergence with an estimated age of 1.99 million years, suggesting that expansion of the L1 Ta subfamily occurred after the divergence of humans and African apes. A total of 262 Ta L1 elements were screened with PCR-based assays to determine their phylogenetic origin and the level of human genomic variation associated with each element. All of the Ta L1 elements analyzed by PCR were absent from the orthologous positions in nonhuman primate genomes, except for a single element (L1HS72) that was also present in the common (Pan troglodytes) and pygmy (P. paniscus) chimpanzee genomes. Sequence analysis revealed that this single exception is the product of a gene conversion event involving an older preexisting L1 element. One hundred fifteen (45%) of the Ta L1 elements were polymorphic with respect to insertion presence or absence and will serve as identical-by-descent markers for the study of human evolution.     

Molecular phylogeny of macaques: implications of nucleotide sequences from an 896-base pair region of mitochondrial DNA

We determined the nucleotide sequences of an 896-base pair region of mitochondrial DNA (mtDNA) from 20 primates representing 13 species of macaques, a baboon, and a patas. We compared these sequences and the homologous sequences from four macaques and a human against each other and deduced the phylogenetic relationships of macaques. The results from the phylogenetic analyses revealed five groups among the macaques: (1) Barbary macaque, (2) two species of Sulawesi macaques, (3) Japanese, rhesus, Taiwanese, crab-eating, and stump-tailed macaques, (4) toque, pig-tailed, and lion-tailed macaques, and (5) Assamese and bonnet macaques. The phylogenetic position of Tibetan macaque remains ambiguous as to whether it belongs to the fourth or fifth group. Phylogenetic trees revealed that Barbary macaque diverged first from the other Asian macaques. Subsequently, the four groups of Asian macaques diverged from one another in a relatively short period of time. Within each group, most of the species diverged in a relatively short period of time following the divergence of the groups. Assuming that the Asian macaques diverged from the outgroup Barbary macaque three million years ago (MYA), the divergence times among groups of Asian macaques were estimated at 2.1-2.5 MYA and within groups at 1.4- 2.2 MYA. The intraspecific nucleotide diversity observed among three rhesus macaques was so large that they did not form a monophyletic cluster in the phylogenetic trees. Instead, one of them formed a cluster with Japanese and Taiwanese macaques, whereas the other two formed a separate cluster. This implies that either polymorphisms of mtDNA sequences that existed before the divergence of these three species (ca. 700,000 years ago) have been retained in rhesus macaques or introgression has occurred among the three species.      

Internal Repetition and Intraindividual Variation in the rDNA ITS1 of the <Emphasis Type="Italic">Anopheles punctulatus</Emphasis> Group (Diptera: Culicidae): Multiple Units and Rates of Turnover

The rapid divergence of repetitive sequences makes them desirable markers for phylogenetic studies of closely related groups, provided that a high level of sequence homogeneity has been maintained within species. Intraspecific polymorphisms are found in an increasing number of studies now, and this highlights the need to determine why these occur. In this study we examined intraindividual variation present in the first ribosomal internal transcribed spacer (ITS1) from a group of cryptic mosquito species. Individuals of the Anopheles punctulatus group contained multiple ITS1 length variants that ranged from 1.2 to 8.0 kb. Nucleotide and copy number variation for several homologous internal repeats is common, yet the intraspecific sequence divergence of cloned PCR isolates is comparable to that of other mosquito species (~0.2–1.5%). Most of the length variation is comprised of a 5′-ITS1 repeat that was identified as a duplication of a conserved ITS2 region. Secondary structure conservation for this repeat is pronounced and several repeat types that are highly homogenized have formed. Significant interspecific divergence indicates a high rate of evolutionary change for this spacer. A maximum likelihood tree constructed here was congruent with previous phylogenetic hypotheses and suggests that concerted evolution is also accompanied by interpopulation divergence. The lack of interindividual differences and the presence of homogenized internal repeats suggest that a high rate of turnover has reduced the overall level of variation. However, the intraindividual variation also appears to be maintained by the absence of a single turnover rate and the complex dynamics of ongoing recombination within the spacer. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A 2.7-kb Deletion in the Human NLRP10 Gene Exon 2 Occurred After the Human–Chimpanzee Divergence

NLRP10 is a member of the NLRP protein family, which is involved in inflammation and apoptosis. Genome sequence comparisons revealed that a 2.7-kb deletion occurred in the human NLRP10 gene exon 2 after the divergence of humans and chimpanzees, resulting in replacement of the entire 3′ untranslated region with the flanking LINE-1 element. The human NLRP10 protein lost 30 or more amino acids that are conserved in primates at its carboxy-terminus. The structural modification of the NLRP10 gene might have played a role in development or enhancement of human-specific traits during evolution.     

An in vitro-selected RNA-binding site for the KH domain protein PSI acts as a splicing inhibitor element

P element somatic inhibitor (PSI) is a 97-kDa RNA-binding protein with four KH motifs that is involved in the inhibition of splicing of the Drosophila P element third intron (IVS3) in somatic cells. PSI interacts with a negative regulatory element in the IVS3 5' exon. This element contains two pseudo-5' splice sites, termed F1 and F2. To identify high affinity binding sites for the PSI protein, in vitro selection (SELEX) was performed using a random RNA oligonucleotide pool. Alignment of high affinity PSI-binding RNAs revealed a degenerate consensus sequence consisting of a short core motif of CUU flanked by alternative purines and pyrimidines. Interestingly, this sequence resembles the F2 pseudo-5' splice site in the P element negative regulatory element. Additionally, a negative in vitro selection of PCR-mutagenized P element 5' exon regulatory element RNAs identified two U residues in the F1 and F2 pseudo-5' splice sites as important nucleotides for PSI binding and the U residue in the F2 region is a nearly invariant nucleotide in the consensus SELEX motif. The high affinity PSI SELEX sequence acted as a splicing inhibitor when placed in the context of a P element splicing pre-mRNA in vitro. Data from in vitro splicing assays, UV crosslinking and RNA-binding competition experiments indicates a strong correlation between the binding affinities of PSI for the SELEX sequences and their ability to modulate splicing of P element IVS3 in vitro.     

Molecular and cytological characterization of 5S rDNA in Oryza species: genomic organization and phylogenetic implications.

We investigated the molecular characteristics and chromosomal organization of 5S rDNA in the genus Oryza, including diploid and tetraploid species. A phylogenetic tree of Oryza species was constructed based on the non-transcribed spacer sequences of 5S rDNA, and some novel relationships were discovered. Specifically, comparative sequence analysis of 5S rDNA in several wild rice species showed unique characteristics inconsistent with the model of concerted evolution: (1) multiple distinct 5S rDNA types were detected within a species, leading to intraspecific divergence of 5S rDNA; (2) multiple identical 5S rDNA types were shared among species, resulting in interspecies clustering of 5S rDNA types; and (3) intraspecific nucleotide diversity was detected within a 5S rDNA class. Our results obtained by fluorescence in situ hybridization revealed that each rice species studied contained only one 5S rDNA locus with two hybridization sites, which were located on either chromosome 7 or chromosome 11. These results suggest that different 5S rDNA classes within the rice genome were arranged together and that one pair of 5S rDNA loci from a diploid progenitor of the tetraploid species might have been lost during evolution. Taken together, our data show that 5S rDNA in rice species is more informative at the gene level than at the chromosome level.     

Extensive intraspecific chloroplast DNA (cpDNA) variation in the alpine Draba aizoides L. (Brassicaceae): haplotype relationships and population structure

Chloroplast DNA (cpDNA) sequence variation is currently the most widely used tool for the inference of phylogenetic relationships among plants at all taxonomic levels. Generally, noncoding regions tend to evolve faster than coding sequences and have recently been applied to the study of phylogenetic relationships among closely related taxa. An implicit assumption of many of these studies is that intraspecific cpDNA variation is either absent or low and therefore will not interfere with the reconstruction of interspecific relationships. A survey of cpDNA sequence variation in the common alpine plant species Draba aizoides L. was undertaken to assess levels of intraspecific cpDNA sequence variation. These levels were compared to levels of interspecific sequence divergence between D. aizoides and related alpine Draba species. Intraspecific cpDNA sequence divergence was extensive in D. aizoides, and intraspecific differences were often larger than interspecific differences. cpDNA haplotype relationships were explored using a maximum parsimony approach and minimum-spanning networks. Results from both methods were largely congruent but comparisons provided interesting insights into the presumed evolutionary history of cpDNA haplotypes. A combined effect of cpDNA introgression and complex lineage sorting was inferred to explain the pattern of cpDNA variation found in D. aizoides. Our results suggest that intraspecific cpDNA variation can be extensive and that intraspecific variation needs to be taken into account when inferring phylogenetic relationships among closely related taxa.     

A major difference between the divergence patterns within the lines-1 families in mice and voles

L1 retroposons are represented in mice by subfamilies of interspersed sequences of varied abundance. Previous analyses have indicated that subfamilies are generated by duplicative transposition of a small number of members of the L1 family, the progeny of which then become a major component of the murine L1 population, and are not due to any active processes generating homology within preexisting groups of elements in a particular species. In mice, more than a third of the L1 elements belong to a clade that became active approximately 5 Mya and whose elements are > or = 95% identical. We have collected sequence information from 13 L1 elements isolated from two species of voles (Rodentia: Microtinae: Microtus and Arvicola) and have found that divergence within the vole L1 population is quite different from that in mice, in that there is no abundant subfamily of homologous elements. Individual L1 elements from voles are very divergent from one another and belong to a clade that began a period of elevated duplicative transposition approximately 13 Mya. Sequence analyses of portions of these divergent L1 elements (approximately 250 bp each) gave no evidence for concerted evolution having acted on the vole L1 elements since the split of the two vole lineages approximately 3.5 Mya; that is, the observed interspecific divergence (6.7%-24.7%) is not larger than the intraspecific divergence (7.9%-27.2%), and phylogenetic analyses showed no clustering into Arvicola and Microtus clades.      

Phylogeny of hoplocercine lizards (Squamata: Iguania) with estimates of relative divergence times

Hoplocercine lizards form a clade of 11 currently recognized species traditionally placed in three genera (Enyalioides, Hoplocercus, and Morunasaurus) that occur in the lowlands on both sides of the Andes between Panama and the Brazilian Cerrado. We analyze 11 mitochondrial and two nuclear loci using probabilistic methods and different partitioning strategies to (1) infer the phylogenetic relationships among species of Hoplocercinae, (2) examine amounts of inter- and intraspecific sequence divergence, (3) address monophyly of four species, (4) test previous phylogenetic hypotheses, and (5) estimate divergence times. Our preferred hypothesis places H. spinosus as the sister taxon to all other species of hoplocercines, with M. annularis nested within Enyalioides. Species with multiple samples are monophyletic except for Enyalioides oshaughnessyi, which is paraphyletic relative to an undescribed species of Enyalioides. All previously published phylogenetic hypotheses for hoplocercines are rejected. Monophyly of Enyalioides cannot be rejected and, consequently, the position of Morunasaurus remains unclear. The most recent common ancestor of Hoplocercinae probably occurred east of the Andes; western taxa included in our analyses originated from at least two separate colonizations whether pre- or post-dating vicariance resulting from uplift of the Andes.     

Evolution of a repeat sequence in the parathyroid hormone-related peptide gene in primates

A polymorphism of the variable number of tandem repeat (VNTR) type is located 97 bp downstream of exon VI of the parathyroid hormone-related peptide (PTHrP) gene in humans. The repeat unit has the general sequence G(TA)_nC, where n equals 4–11. In order to characterize the evolutionary history of this VNTR, we initially tested for its presence in 13 different species representing four main groups of living primates. The sequence is present in the human, great apes, and Old World monkeys, but not in New World monkeys; and this region failed to PCR amplify in the Loris group. Thus, the evolution of the sequence as part of the PTHrP gene started at least 25–35 millions years ago, after divergence of the Old World and New World monkeys, but before divergence of Old World monkeys and great apes and humans. The structural changes occurring during evolution are characterized by a relatively high degree of sequence divergence. In general, the tandem repeat region tends to be longer and more complex in higher primates with the repeat unit motifs all being based on a TA-dinucleotide repeat sequence. Intra-species variability of the locus was demonstrated only in humans and gorilla. The divergence of the TA-dinucleotide repeat sequence and the variable mutation rates observed in different primate species are in contrast to the relative conservation of the flanking sequences during primate evolution. This suggests that the nature of the TA-dinucleotide repeat sequence, rather than its flanking sequences, is responsible for generating variability. Particular features of the sequence may allow it to form stable secondary structures during DNA replication, and this, in turn, could promote slipped-strand mispairing to occur.     

The ITS1-5.8S-ITS2 sequence region in the Musaceae: structure, diversity and use in molecular phylogeny

Genes coding for 45S ribosomal RNA are organized in tandem arrays of up to several thousand copies and contain 18S, 5.8S and 26S rRNA units separated by internal transcribed spacers ITS1 and ITS2. While the rRNA units are evolutionary conserved, ITS show high level of interspecific divergence and have been used frequently in genetic diversity and phylogenetic studies. In this work we report on the structure and diversity of the ITS region in 87 representatives of the family Musaceae. We provide the first detailed information on ITS sequence diversity in the genus Musa and describe the presence of more than one type of ITS sequence within individual species. Both Sanger sequencing of amplified ITS regions and whole genome 454 sequencing lead to similar phylogenetic inferences. We show that it is necessary to identify putative pseudogenic ITS sequences, which may have negative effect on phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstruction based on ITS sequence showed that the genus Musa is divided into two distinct clades--Callimusa and Australimusa and Eumusa and Rhodochlamys. Most of the intraspecific banana hybrids analyzed contain conserved parental ITS sequences, indicating incomplete concerted evolution of rDNA loci. Independent evolution of parental rDNA in hybrids enables determination of genomic constitution of hybrids using ITS. The observation of only one type of ITS sequence in some of the presumed interspecific hybrid clones warrants further study to confirm their hybrid origin and to unravel processes leading to evolution of their genomes.     

Full-Length L1 Elements Have Arisen Recently in the Same 1-kb Region of the Gorilla and Human Genomes

New copies of the mammalian retrotransposon L1 arise in the germline at an undetermined rate. Each new L1 copy appears at a specific evolutionary time point that can be estimated by phylogenetic analysis. In humans, the active L1 sequence L1.2 resides at the genomic locus LRE1. Here we analyzed the region surrounding the LRE1 locus in humans and gorillas to determine the evolutionary history of the region and to estimate the age of L1.2. We found that the region was composed of an ancient L1, L1Hs-Lrg, which was significantly divergent from all other L1 sequences available in the databases. We also determined that L1.2 was absent from the gorilla genome and arose in humans after the divergence of gorilla and human lineages. In the gorilla LRE1 region, we discovered a different full-length L1 element, L1Gg-1, which was allelic and present at a high gene frequency in gorillas but absent from other primates. We determined the nucleotide sequence of L1Gg-1 and found that it was 98% identical to L1.2, suggesting a close relationship between active L1s in gorillas and humans. Received: 28 December 1997 / Accepted: 20 March 1998     

Intraspecific Concerted Evolution of the rDNA ITS1 in Anopheles farauti Sensu Stricto (Diptera: Culicidae) Reveals Recent Patterns of Population Structure

We examined the intraindividual variation present in the first ribosomal internal transcribed spacer (ITS1) of Anopheles farauti to determine the level of divergence among populations for this important malarial vector. We isolated 187 clones from 70 individuals and found regional variation among four internal tandem repeats. The data were partitioned prior to analysis given the presence of a paralogous ITS2 sequence, called the 5'-subrepeat, inserted in the ITS1 of most clones. A high level of homogenization and population differentiation was observed for this repeat, which indicates a higher rate of turnover relative to the adjacent 'core' region. Bayesian analysis was performed using several substitutional models on both a combined and a partitioned data set. On the whole, the ITS1 phylogeny and geographic origin of the samples appear to be congruent. Some interesting exceptions indicate the spread of variant repeats between populations and the retention of ancestral polymorphism. Our data clearly demonstrate concerted evolution at the intraspecific level despite intraindividual variation and a complex internal repeat structure from a species that occupies a continuous coastal distribution. A high rate of genomic turnover in combination with a high level of sequence divergence appears to be a major factor leading to its concerted evolution within these populations.     

mtDNA diversity in rhesus monkeys reveals overestimates of divergence time and paraphyly with neighboring species

Reconstructions of the human-African great ape phylogeny by using mitochondrial DNA (mtDNA) have been subject to considerable debate. One confounding factor may be the lack of data on intraspecific variation. To test this hypothesis, we examined the effect of intraspecific mtDNA diversity on the phylogenetic reconstruction of another Plio- Pleistocene radiation of higher primates, the fascicularis group of macaque (Macaca) monkey species. Fifteen endonucleases were used to identify 10 haplotypes of 40-47 restriction sites in M. mulatta, which were compared with similar data for the other members of this species group. Interpopulational, intraspecific mtDNA diversity was large (0.5%- 4.5%), and estimates of divergence time and branching order incorporating this variation were substantially different from those based on single representatives of each species. We conclude that intraspecific mtDNA diversity is substantial in at least some primate species. Consequently, without prior information on the extent of genetic diversity within a particular species, intraspecific variation must be assessed and accounted for when reconstructing primate phylogenies. Further, we question the reliability of hominoid mtDNA phylogenies, based as they are on one or a few representatives of each species, in an already depauperate superfamily of primates.      

Alternative splicing of type II procollagen exon 2 is regulated by the combination of a weak 5' splice site and an adjacent intronic stem-loop cis element

Alternative splicing of the type II procollagen gene (COL2A1) is developmentally regulated during chondrogenesis. Chondroprogenitor cells produce the type IIA procollagen isoform by splicing (including) exon 2 during pre-mRNA processing, whereas differentiated chondrocytes synthesize the type IIB procollagen isoform by exon 2 skipping (exclusion). Using a COL2A1 mini-gene and chondrocytes at various stages of differentiation, we identified a non-classical consensus splicing sequence in intron 2 adjacent to the 5' splice site, which is essential in regulating exon 2 splicing. RNA mapping confirmed this region contains secondary structure in the form of a stem-loop. Mutational analysis identified three cis elements within the conserved double-stranded stem region that are functional only in the context of the natural weak 5' splice site of exon 2; they are 1) a uridine-rich enhancer element in all cell types tested except differentiated chondrocytes; 2) an adenine-rich silencer element, and 3) an enhancer cis element functional in the context of secondary structure. This is the first report identifying key cis elements in the COL2A1 gene that modulate the cell type-specific alternative splicing switch of exon 2 during cartilage development.     

Origins,Diversity and Relationships of Lemurs

I review new evidence on origins and adaptive radiation of Malagasy lemurs, a remarkably diverse group containing 13% of living primate species. The number of recognized lemur species has increased significantly, partly due to research revealing specific subdivisions within known populations but mainly because of discovery of new populations through fieldwork. Some species feared to be extinct have also been rediscovered. Specific numbers have increased particularly in small-bodied, cryptic genera for which continued research will surely reveal even more species.Adaptative radiation of lemurs has been essentially confined to Madagascar. The high density of lemur species on that island, associated with very small geographical ranges, has major implications both for their evolutionary divergence and for conservation. Reconstructions of phylogenetic relationships among primates have been considerably enhanced by DNA sequence data. Sufficient data are now available from both nuclear and mitochondrial sequences to examine relationships among and within the major groups of living primates. Most studies have confirmed that lemurs constitute a monophyletic sister-group of the lorisiform clade and all exclude a specific relationship between cheirogaleids and lorisiforms repeatedly inferred from morphological evidence. However, some analyses indicate that the aye-aye may have branched away before the divergence between other lemurs and lorisiforms. DNA sequence analyses have also yielded a broad consensus for relationships between Eulemur, Hapalemur, Lemur and Varecia: Varecia branched away first, while Lemur is more closely related to Hapalemur than to Eulemur. As debate about phylogenetic relationships among lemurs and other primates seems to have been settled in favor of lemur monophyly (possibly excluding the aye-aye), only a single invasion of Madagascar is required; but it must still be explained how ancestral lemurs could have migrated there at an appropriate time. Separation between Madagascar and Africa was apparently complete by about 120 Ma, too far in the past for direct overland migration. A recent hypothesis suggested that uplifted land in the Mozambique Channel assisted colonization of Madagascar 26-45 Ma, seemingly agreeing with an estimated date of about 40 Ma for divergence of lemurs from other primates. However, mounting evidence suggests that divergence occurred significantly earlier. Because the earliest known fossil representatives of several modern orders of placental mammals (including primates) are dated no earlier than the early Tertiary, it is widely accepted that their divergence took place after the Cretaceous/Tertiary mass extinction. Yet the known fossil record can only yield minimum divergence times; if sampling is poor and/or biased there may be a considerable discrepancy between minimum and actual dates. There is, for example, virtually no known fossil record for lemurs in Madagascar and the earliest known representatives are subfossil lemurs, so in this case a direct reading of the fossil record would indicate that the lemurs first originated just a few thousand years ago! Examination of underestimation of times of origin because of poor sampling in the fossil record has confirmed previous suggestions that primates originated considerably earlier than generally believed. Several recent phylogenetic reconstructions based on DNA sequence data and using calibration dates derived from groups other than primates provide independent support for this inference. Overall, it now seems that primates originated at around 90 Ma rather than the 55 Ma indicated by direct reading of the known fossil record. Hence, colonization of Madagascar by lemurs would have taken place at about 80 Ma, double the date usually accepted, and should be interpreted in terms of contemporary continental relationships.     

The evolution of coexisting highly divergent line-1 subfamilies within the rodent genus Peromyscus

Summary Two distinct members of the LINE-1 (L1) family in Peromyscus were characterized. The two clones, denoted L1Pm55 and L1Pm62, were 1.5 kb and 1.8 kb in length, respectively, and align to the identical region of the L1 sequence of Mus domesticus. Sequence similarity was on the order of 70% between L1Pm55 and L1Pm62, which approximates that between either Peromyscus sequence and Mus Ll. L1Pm62 represents a more prevalent subfamily than L1Pm55. L1Pm62 exists in about 500 copies per haploid genome, while L1Pm55 exists in about 100 copies. The existence of major and minor subpopulations of L1 within Peromyscus is in contrast to murine rodents and higher primates, where L1 copy number is on the order of 20,000 to 100,000, and where levels of intraspecific divergence among L1 elements are typically less than 15–20%. Additional Peromyscus clones are similarly divergent from both L1Pm62 and LIPm55, implying the existence of more than two distinct L1 subfamilies. The highly divergent L1 subfamilies in Peromyscus apparently have been evolving independently for more than 25 million years, preceding the divergence of cricetine and murine rodents. Investigations of the evolution of L1 within Peromyscus by restriction and Southern analysis was performed using species groups represented by the partially interfertile species pairs P. maniculatus-P. polionotus, P. leucopus-P. gossypinus, and P. truei-P. difficilis of the nominate subgenus and P. californicus of the Haplomylomys subgenus. Changes in L1 and species group taxonomic boundaries frequently coincided. The implications for phylogeny are discussed.