Amplification of the ancient murine Lx family of long interspersed repeated DNA occurred during the murine radiation

Subfamily-specific quantification of endogenous mouse L1 retrotransposons by droplet digital PCR

Long interspersed element type 1 (LINE-1; L1) mobilizes during early embryogenesis, neurogenesis, and germ cell development, accounting for 25% of disease-causing heritable insertions and 98% of somatic insertions in cancer. To better understand the regulation and impact of L1 mobilization in the genome, reliable methods for measuring L1 copy number variation (CNV) are needed. Here we present a comprehensive analysis of a droplet digital PCR (ddPCR) based method for quantifying endogenous mouse L1. We provide experimental evidence that ddPCR assays can be designed to target specific L1 subfamilies using diagnostic single nucleotide polymorphisms (SNPs). The target and off-target L1 subfamilies form distinct droplet clusters, which were experimentally verified using both synthetic gene fragments and endogenous L1 derived plasmid clones. We further provide a roadmap for in silico assay design and evaluation of target specificity, ddPCR testing, and optimization for L1 CNV quantification. The assay can achieve a sensitivity of 5% CNV with 8 technical replicates. With 24 technical replicates, it can detect 2% CNV because of the increased precision. The same approach will serve as a guide for the development of ddPCR based assays for quantifying human L1 copy number and any other high copy genomic target sequences.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0117-4) contains supplementary material, which is available to authorized users.

Transposable elements and small RNAs: Genomic fuel for species diversity

While transposable elements (TE) have long been suspected of involvement in species diversification, identifying specific roles has been difficult. We recently found evidence of TE-derived regulatory RNAs in a species-rich family of bats. The TE-derived small RNAs are temporally associated with the burst of species diversification, suggesting that they may have been involved in the processes that led to the diversification. In this commentary, we expand on the ideas that were briefly touched upon in that manuscript. Specifically, we suggest avenues of research that may help to identify the roles that TEs may play in perturbing regulatory pathways. Such research endeavors may serve to inform evolutionary biologists of the ways that TEs have influenced the genomic and taxonomic diversity around us.

Jumping the fine LINE between species: horizontal transfer of transposable elements in animals catalyses genome evolution

Horizontal transfer (HT) is the transmission of genetic material between non-mating species, a phenomenon thought to occur rarely in multicellular eukaryotes. However, many transposable elements (TEs) are not only capable of HT, but have frequently jumped between widely divergent species. Here we review and integrate reported cases of HT in retrotransposons of the BovB family, and DNA transposons, over a broad range of animals spanning all continents. Our conclusions challenge the paradigm that HT in vertebrates is restricted to infective long terminal repeat (LTR) retrotransposons or retroviruses. This raises the possibility that other non-LTR retrotransposons, such as L1 or CR1 elements, believed to be only vertically transmitted, can horizontally transfer between species. Growing evidence indicates that the process of HT is much more general across different TEs and species than previously believed, and that it likely shapes eukaryotic genomes and catalyses genome evolution.

Phylogenetic relationships of extant zokors (Myospalacinae) (Rodentia, Spalacidae) inferred from mitochondrial DNA sequences

In this study, we use three mitochondrial markers, cytochrome b gene (Cyt b), NADH dehydrogenase subunit 4 (ND4) and control region (D-loop) to investigate the phylogenetic relationships of extant zokor species in Mysopalacinae. The phylogenetic tree constructed based on Cyt b strongly supports the monophyly genera Eospalax and Myospalax with E. fontanierii being the most ancient species in Eospalax. Further phylogenetic analyses of four species of Eospalax based on ND4 and D-loop sequences revealed two clades that correspond to two geographical distributions. The basal clade includes E. cansus which is mainly found on Loess Plateau (LP) and another clade including E. baileyi, E. smithii and E. rufescens that inhabits areas above 2000 m on Qinghai-Tibetan Plateau (QTP) and Qinling Mountains. Geographical events of QTP and LP may have played a major role in the diversification and evolution of Mysopalacinae.

Phylogeny and biogeography of African Murinae based on mitochondrial and nuclear gene sequences, with a new tribal classification of the subfamily

Background

Within the subfamily Murinae, African murines represent 25% of species biodiversity, making this group ideal for detailed studies of the patterns and timing of diversification of the African endemic fauna and its relationships with Asia. Here we report the results of phylogenetic analyses of the endemic African murines through a broad sampling of murine diversity from all their distribution area, based on the mitochondrial cytochrome b gene and the two nuclear gene fragments (IRBP exon 1 and GHR).

Results

A combined analysis of one mitochondrial and two nuclear gene sequences consistently identified and robustly supported ten primary lineages within Murinae. We propose to formalize a new tribal arrangement within the Murinae that reflects this phylogeny. The diverse African murine assemblage includes members of five of the ten tribes and clearly derives from multiple faunal exchanges between Africa and Eurasia. Molecular dating analyses using a relaxed Bayesian molecular clock put the first colonization of Africa around 11 Mya, which is consistent with the fossil record. The main period of African murine diversification occurred later following disruption of the migration route between Africa and Asia about 7–9 Mya. A second period of interchange, dating to around 5–6.5 Mya, saw the arrival in Africa of Mus (leading to the speciose endemic Nannomys), and explains the appearance of several distinctive African lineages in the late Miocene and Pliocene fossil record of Eurasia.

Conclusion

Our molecular survey of Murinae, which includes the most complete sampling so far of African taxa, indicates that there were at least four separate radiations within the African region, as well as several phases of dispersal between Asia and Africa during the last 12 My. We also reconstruct the phylogenetic structure of the Murinae, and propose a new classification at tribal level for this traditionally problematic group.

Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates

We investigated the evolution of the families of LINE-1 (L1) retrotransposons that have amplified in the human lineage since the origin of primates. We identified two phases in the evolution of L1. From approximately 70 million years ago (Mya) until approximately 40 Mya, three distinct L1 lineages were simultaneously active in the genome of ancestral primates. In contrast, during the last 40 million years (Myr), i.e., during the evolution of anthropoid primates, a single lineage of families has evolved and amplified. We found that novel (i.e., unrelated) regulatory regions (5'UTR) have been frequently recruited during the evolution of L1, whereas the two open-reading frames (ORF1 and ORF2) have remained relatively conserved. We found that L1 families coexisted and formed independently evolving L1 lineages only when they had different 5'UTRs. We propose that L1 families with different 5'UTR can coexist because they don't rely on the same host-encoded factors for their transcription and therefore do not compete with each other. The most prolific L1 families (families L1PA8 to L1PA3) amplified between 40 and 12 Mya. This period of high activity corresponds to an episode of adaptive evolution in a segment of ORF1. The correlation between the high activity of L1 families and adaptive evolution could result from the coevolution of L1 and a host-encoded repressor of L1 activity.

SINE extinction preceded LINE extinction in sigmodontine rodents: implications for retrotranspositional dynamics and mechanisms

Short Interspersed Nuclear Elements, or SINEs, retrotranspose despite lacking protein-coding capability. It has been proposed that SINEs utilize enzymes produced in trans by Long Interspersed Nuclear Elements, or LINEs. Strong support for this hypothesis is found in LINE and SINE pairs that share sequence homology; however, LINEs and SINEs in primates and rodents are only linked by an insertion site motif. We have now profiled L1 LINE and B1 SINE activity in 24 rodent species including candidate taxa for the first documented L1 extinction. As expected, there was no evidence for recent activity of B1s in species that also lack L1 activity. However, B1 silencing appears to have preceded L1 extinction, since B1 activity is also lacking in the genus most closely related to those lacking active L1s despite the presence of active L1s in this genus. A second genus with active L1s but inactive B1s was also identified.

Extinction of LINE-1 activity coincident with a major mammalian radiation in rodents

LINE-1 transposable elements (L1s) are ubiquitous in mammals and are thought to have remained active since before the mammalian radiation. Only one L1 extinction event, in South American rodents in the genus Oryzomys, has been convincingly demonstrated. Here we examine the phylogenetic limits and evolutionary tempo of that extinction event by characterizing L1s in related rodents. Fourteen genera from five tribes within the Sigmodontinae subfamily were examined. Only the Sigmodontini, the most basal tribe in this group, demonstrate recent L1 activity. The Oryzomyini, Akodontini, Phyllotini, and Thomasomyini contain only L1s that appear to have inserted long ago; their L1s lack open reading frames, have mutations at conserved amino acid residues, and show numerous private mutations. They also lack restriction site-defined L1 subfamilies specific to any species, genus or tribe examined, and fail to form monophyletic species, genus or tribal L1 clusters. We determine here that this L1 extinction event occurred roughly 8.8 million years ago, near the divergence of Sigmodon from the remaining Sigmodontinae species. These species appear to be ideal model organisms for studying the impact of L1 inactivity on mammalian genomes.

The phylogenetic position of the zokors (Myospalacinae) and comments on the families of muroids (Rodentia)

Recent molecular studies have concluded that the genus Myospalax evolved from within the rodent subfamily Cricetinae. This conclusion was tested using the complete sequences from the mitochondrial 12S rRNA and cytochrome b genes. Based on our analyses, Myospalax appears to be sister to a clade containing the subfamilies Spalacinae and Rhizomyinae, and all three of these lineages appear to be basal to the superfamily Muroidea. Based on the position of these three lineages, we suggest that they be placed in a distinct family, the Spalacidae, rather than subsumed as subfamilies in the family Muridae. Finally, our analyses suggest that the earlier placement of Myospalax as a member of the Cricetinae is the result of a single misidentified specimen, which was not a Myospalax.

Different rates of LINE-1 (L1) retrotransposon amplification and evolution in New World monkeys

LINE-1 (L1) elements constitute the major family of retrotransposons in mammalian genomes. Here we report the first investigation of L1 evolution in New World monkeys (NWM). Two regions of the second open-reading frame were analyzed by two methods in three NWM species, the squirrel monkey (Saimiri sciureus), the tamarin (Saguinus oedipus), and the spider monkey (Ateles paniscus). Since these three species diverged, L1 has amplified in the Saimiri and Saguinus lineages but L1 activity seems to have been strongly reduced in the Ateles lineage. In addition, the active L1 lineage has evolved rapidly in Saimiri and Saguinus, generating species-specific subfamilies. In contrast, we found no evidence for a species-specific subfamily in Ateles, a result consistent with the low L1 activity in this species for the last approximately 25 My.

The spermatozoon of Eurasian murine rodents: Its morphological diversity and evolution

The murine rodents are the most speciose subfamily of mammals. Here the morphology of the spermatozoon, as determined by scanning and transmission electron microscopy of representative species from four Eurasian clades, is described. Much interspecific variability in all components of the spermatozoon was found to occur, although most species have a bilaterally flattened sperm head with a single apical hook of variable length and orientation. Ultrastructural observations indicate that this apical hook invariably contains a nuclear projection as well as a large extension of the subacrosomal cytoskeleton, as a perforatorium rostrally, and a complex asymmetrical acrosomal extension. These spermatozoa also have relatively long tails that are attached to the lower concave surface of the sperm head. Uniquely, in species in the Apodemus clade, the apical hook is orientated caudally. In a few species a highly derived sperm head morphotype that does not contain an apical hook is present. These sperm heads vary in morphology from being globular in two species of Bandicota, to bilaterally flattened and paddle-shaped in Tokudaia and Micromys. In spermatozoa of the latter two genera the subacrosomal cytoskeleton, which is less extensive than in species with a hooked sperm head, forms an apical extension, but that is not the case in Bandicota. In all species where the sperm head lacks an apical hook the acrosome is more symmetrical. The sperm tail is much shorter in these species, with attachment to the head occurring on the ventral surface in Tokudaia and basal in Micromys and the two species of Bandicota. As the sperm head morphotype with a complex apical hook is present in all the major clades of murine rodents, it is likely to be a plesiomorphic character within each of these clades, with the nonhooked sperm heads, which vary greatly in structure between species of the different lineages, probably being independently derived. The ultrastructural organization of the sperm head of Bandicota, but not those of Micromys or Tokudaia, suggest divergence in some of the morphological events associated with sperm-egg interaction at the time of fertilization.

Phylogeny of the genus Apodemus with a special emphasis on the subgenus Sylvaemus using the nuclear IRBP gene and two mitochondrial markers: cytochrome b and 12S rRNA

Phylogenetic relationships among 17 extant species of Murinae, with special reference to the genus Apodemus, were investigated using sequence data from the nuclear protein-coding gene IRBP (15 species) and the two mitochondrial genes cytochrome b and 12S rRNA (17 species). The analysis of the three genes does not resolve the relationships between Mus, Apodemus, and Rattus but separates Micromys from these three genera. The analysis of the two mitochondrial regions supported an association between Apodemus and Tokudaia and indicated that these two genera are more closely related to Mus than to Rattus or Micromys. Within Apodemus, the mitochondrial data sets indicated that 8 of the 9 species analyzed can be sorted into two main groups: an Apodemus group, with A. agrarius, semotus, and peninsulae, and a Sylvaemus group, with uralensis, flavicollis, alpicola, sylvaticus, and hermonensis. The position of Apodemus mystacinus is ambiguous and might be either included in Sylvaemus or considered a distinct subgenus, Karstomys, more closely related to Sylvaemus than to Apodemus. Estimation of the divergence time for these taxa suggests a separation between 7 and 8 My ago for the three groups (mystacinus and the two subgenera Apodemus and Sylvaemus). Within each subgenus, divergence times are between 5.4 and 6 My for Apodemus and between 2.2 and 3.5 My for Sylvaemus and mystacinus.

Black rat ( Rattus rattus) genomic variability characterized by chromosome painting

Black rats are of outstanding interest in parasitology and infective disease analysis. We used chromosome paints from both the mouse ( Mus musculus) and the Norway rat ( Rattus norvegicus) to characterize the genome of two Black rat subspecies from Italy. Both subspecies have two large metacentrics (n. 1, 4) not present in the Norway rat (2n = 42). Rattus rattus rattus has a diploid number of 2n = 38, while Rattus rattus frugivorous has two small metacentric "supernumerary" or B chromosomes for a diploid number of 2n = 38 + 2B. The 21 mouse paints gave 38 signals on the R. r. rattus karyotype and 39 signals in the R. r. frugivorous karyotype. The two metacentrics, not present in R. norvegicus, were hybridized by mouse 16/1/17 and mouse 4/10/15. These chromosomes are homologous to: RRA1 = RNO 5/7, and RRA4 = RNO 9/11 and not "4/7" and "11/12" as previously reported. Furthermore, the synteny of Chr 13 of the R. r. frugivorous with R. norvegicus Chr 16 and mouse Chrs 8/14 is not complete, because there is a small pericentromeric insertion of RNO Chr 18 (mouse Chr 18). If we consider only the two metacentrics, RRA1 and RRA4, the principal differences between R. norvegicus and R. rattus, then we can propose the derived synteny of 124 genes in the black rat. A comparison of the Z index between rats and mice shows an acceleration of genomic evolution among genus, species, and subspecies. The chromosomal differences between R. r. rattus x R. r. frugivorous suggest that they may be classified as different species because hybrids would produce 50% unbalanced gametes.

Sequence analysis of mouse vomeronasal receptor gene clusters reveals common promoter motifs and a history of recent expansion

We have analyzed the organization and sequence of 73 V1R genes encoding putative pheromone receptors to identify regulatory features and characterize the evolutionary history of the V1R family. The 73 V1Rs arose from seven ancestral genes around the time of mouse-rat speciation through large local duplications, and this expansion may contribute to speciation events. Orthologous V1R genes appear to have been lost during primate evolution. Exceptional noncoding homology is observed across four V1R subfamilies at one cluster and thus may be important for locus-specific transcriptional regulation.

The bushlike radiation of muroid rodents is exemplified by the molecular phylogeny of the LCAT nuclear gene

Phylogenetic relationships among 40 extant species of rodents, with an emphasis on the taxonomic sampling of Muridae and Dipodidae, were studied using sequences of the nuclear protein-coding gene LCAT (lecithin cholesterol acyl transferase). Analysis of 804 bp from the exonic regions of LCAT confirmed many traditional groupings in and around Muridae. A strong support was found for the families Muridae (represented by 29 species) and Dipodidae (5 species). Compared with Sciuridae, Gliridae, and Caviomorpha, the Dipodidae family appeared the closest relative of Muridae, confirming the suprafamilial Myodonta concept. Within the speciose family Muridae, the first branching leads to the fossorial Spalacinae and semifossorial Rhyzomyinae. The remaining components of Muridae appear as a polytomy from which are issued Sigmodontinae, Calomyscinae, Arvicolinae, Cricetinae, Mystromyinae, Nesomyinae, and some Dendromurinae (Steatomys and Dendromus). This phylogeny is interpreted as the result of a bushlike radiation at the end of the early Miocene, leading to emergence of most living subfamilies. The separation between three additional taxa, Murinae, Gerbillinae, and "Acomyinae" (which comprises the genera Acomys, Deomys, Uranomys, and Lophuromys), has occurred more recently from a common ancestor issued from the main basal radiation. As previously shown by other molecular studies, the vlei rats, Otomyinae, are nested within Old World Murinae. In the same way, the zokors, Myospalacinae, appear strongly nested within the hamsters, Cricetinae. Finally, we propose a sister group relationship between Malagasy Nesomyinae and south African Mystromyinae.

Molecular phylogeny of European muroid rodents based on complete cytochrome b sequences

Phylogenetic relationships among 18 species of mainly European muroid rodents that belong to three subfamilies were estimated using complete sequences of the mitochondrial cytochrome b gene. The inferred monophyly of the subfamilies Murinae (mice and rats) and Arvicolinae (voles, lemmings, and muskrats) is in agreement with previous studies. Within the Murinae, the morphology-based division of the genus Apodemus into three subgenera is supported by these DNA sequence data. The relationships among the different genera of the Murinae were generally poorly resolved, and the relationships of Micromys and Acomys to the other murine genera remained unresolved. Within the subfamily Arvicolinae, the relations of the genera Arvicola, Clethrionomys, and Microtus remained tentative with our data. However, within the Microtus group, there is a good molecular support for the phylogenetic relationships. These findings suggest that the origin of the different murine and arvicoline lineages was rapid, indicating an adaptive radiation with fast speciation.

L1 (LINE-1) retrotransposon evolution and amplification in recent human history

L1 (LINE-1) elements constitute a large family of mammalian retrotransposons that have been replicating and evolving in mammals for more than 100 Myr and now compose 20% or more of the DNA of some mammals. Here, we investigated the evolutionary dynamics of the active human Ta L1 family and found that it arose approximately 4 MYA and subsequently differentiated into two major subfamilies, Ta-0 and Ta-1, each of which contain additional subsets. Ta-1, which has not heretofore been described, is younger than Ta-0 and now accounts for at least 50% of the Ta family. Although Ta-0 contains some active elements, the Ta-1 subfamily has replaced it as the replicatively dominant subfamily in humans; 69% of the loci that contain Ta-1 inserts are polymorphic for the presence or absence of the insert in human populations, as compared with 29% of the loci that contain Ta-0 inserts. This value is 90% for loci that contain Ta-1d inserts, which are the youngest subset of Ta-1 and now account for about two thirds of the Ta-1 subfamily. The successive emergence and amplification of distinct Ta L1 subfamilies shows that L1 evolution has been as active in recent human history as it has been found to be for rodent L1 families. In addition, Ta-1 elements have been accumulating in humans at about the same rate per generation as recently evolved active rodent L1 subfamilies.

The biological properties and evolutionary dynamics of mammalian LINE-1 retrotransposons

Mammalian LINE-1 (L1) elements belong to the superfamily of autonomously replicating retrotransposable elements that lack the long terminal repeated (LTR) sequences typical of retroviruses and retroviral-like retrotransposons. The non-LTR superfamily is very ancient and L1-like elements are ubiquitous in nature, having been found in plants, fungi, invertebrates, and various vertebrate classes from fish to mammals. L1 elements have been replicating and evolving in mammals for at least the past 100 million years and now constitute 20% or more of some mammalian genomes. Therefore, L1 elements presumably have had a profound, perhaps defining, effect on the evolution, structure, and function of mammalian genomes. L1 elements contain regulatory signals and encode two proteins: one is an RNA-binding protein and the second one presumably functions as an integrase-replicase, because it has both endonuclease and reverse transcriptase activities. This work reviews the structure and biological properties of L1 elements, including their regulation, replication, evolution, and interaction with their mammalian hosts. Although each of these processes is incompletely understood, what is known indicates that they represent challenging and fascinating biological phenomena, the resolution of which will be essential for fully understanding the biology of mammals.

LINE-1 (L1) lineages in the mouse

Recently, a rapidly amplifying family of mouse LINE-1 (L1) has been identified and named T(F). The evolutionary context surrounding the derivation of the T(F) family was examined through phylogenetic analysis of sequences in the 3' portion of the repeat. The Mus musculus domesticus T(F) family was found to be the terminal subfamily of the previously identified L1Md4 lineage. The L1Md4 lineage joins the other prototypical mouse LINE-1 lineage (the L1MdA2 lineage) approximately 1 MYA at about the time of the common ancestor of M. m. domesticus, Mus spicilegus, and Mus spretus. However, the T(F) family from M. m. domesticus was found to join to the previously reported M. spretus Ms475 and Ms7024 LINE-1 families at just 0.5 MYA, indicating horizontal transfer. The T(F) family from M. m. domesticus was then found to be even more recently related to LINE-1's from another species, M. spicilegus. A separate spretus A2 lineage was found through a directed search of a PCR library. This lineage, in contrast to the spretus T(F) lineage, does join domesticus at about 1 MYA, as would be expected in the absence of horizontal transfer. A third major family was also found that splits off from the L1Md4 lineage shortly after its departure from the L1MdA2 lineage. The new family, named the Z family, was found to contain the de novo LINE-1 inserts causing the beige and med mutations. Whether the split with the Z family was before or after the recombination that introduced the F-type promoters and defined the inception of T(F) as a lineage is unclear. In enumerating copies of the various LINE-1 families, we found that T(F) 3' ends were not much more numerous than the reported number of 5' ends, suggesting that T(F) may not be subjected to the 90% truncation pattern typical of LINE-1 as a whole.

The phylogenetic position of "Acomyinae" (Rodentia, Mammalia) as sister group of a Murinae + Gerbillinae clade: evidence from the nuclear ribonuclease gene

The phylogenetic relationships of Acomys and Uranomys within Muridae were investigated using nuclear pancreatic ribonuclease A gene sequences. The various kinds of substitutions in the data matrix (15 taxa x 375 nucleotides) were examined for saturation, in order to apply a weighted parsimony approach. Phylogenies were derived by maximum parsimony (weighted and unweighted) and maximum likelihood procedures, using a dormouse (Gliridae) as outgroup. Maximum likelihood gave the most robust results. All analyses cluster some traditional taxa with a strong robustness, such as three species of the genus Mus, two South-East Asian rats, and two genera in each of the gerbil and vole families. When analyzed with those of other murid rodents representing Murinae, Gerbillinae, Arvicolinae, Cricetinae, and Sigmodontinae, sequences of the ribonuclease gene suggest that Acomys and Uranomys constitute a monophyletic clade at the subfamily level, denoted "Acomyinae." The relationships between the six subfamilies of Muridae appear poorly resolved, except for a clade uniting Murinae, Acomyinae, and Gerbillinae. Within this clade, the sister group of Acomyinae could not be identified, as the branch length defining a Gerbillinae + Murinae cluster is extremely short. The poor resolution of our phylogenetic inferences is probably the result of two confounding factors, namely the limited size of the pancreatic ribonuclease sequence and the probable short time intervals during the radiation of the six murid subfamilies involved in this study.

Determining and dating recent rodent speciation events by using L1 (LINE-1) retrotransposons

Phylogenies based on the inheritance of shared derived characters will be ambiguous when the shared characters are not the result of common ancestry. Such characters are called homoplasies. Phylogenetic analysis also can be problematic if the characters have not changed sufficiently, as might be the case for rapid or recent speciations. The latter are of particular interest because evolutionary processes may be more accessible the more recent the speciation. The repeated DNA subfamilies generated by the mammalian L1 (LINE-1) retrotransposon are apparently homoplasy-free phylogenetic characters. L1 retrotransposons are transmitted only by inheritance and rapidly generate novel variants that produce distinct subfamilies of mostly defective copies, which then "age" as they diverge. Here we show that the L1 character can both resolve and date recent speciation events within the large group of very closely related rats known as Rattus sensu stricto. This lineage arose 5-6 million years ago (Mya) and subsequently underwent two episodes of speciation: an intense one, approximately 2.7 Mya, produced at least five lineages in <0.3 My; a second began approximately 1.2 Mya and may still be continuing.

Molecular phylogeny of rodents, with special emphasis on murids: evidence from nuclear gene LCAT

Phylogenetic relationships among 19 extant species of rodents, with special emphasis on rats, mice, and allied Muroidea, were studied using sequences of the nuclear protein-coding gene LCAT (lecithin:cholesterol acyltransferase), an enzyme of cholesterol metabolism. Analysis of 705 base pairs from the exonic regions of LCAT confirmed known groupings in and around Muroidea. Strong support was found for the families Sciuridae (squirrel and marmot) and Gliridae (dormice) and for suprafamilial taxa Muroidea and Caviomorpha (guinea pig and allies). Within Muroidea, the first branching leads to the fossorial mole rats Spalacinae and bamboo rats Rhizomyinae. The other Muroidea appear as a polytomy from which are issued Gerbillinae (gerbils), Murinae (rats and mice), Sigmodontinae (New World cricetids), Cricetinae (hamsters), and Arvicolinae (voles). Evidence from LCAT sequences agrees with that from a number of previous molecular and morphological studies, both concerning branching orders inside Muroidea and the bush-like radiation of rodent suprafamilial taxa (caviomorphs, sciurids, glirids, muroids), thus suggesting that this nuclear gene is an appropriate candidate for addressing questions of rodents relationships.

Determination of the evolutionary relationships in Rattus sensu lato (Rodentia : Muridae) using L1 (LINE-1) amplification events

We determined approximately 215 bp of DNA sequence from the 3'-untranslated region (UTR) of 240 cloned L1 (LINE-1) elements isolated from 22 species of Rattus sensu lato and Rattus sensu stricto murine rodents. The sequences were sorted into different L1 subfamilies, and oligonucleotides cognate to them were hybridized to genomic DNA of various taxa. From the distribution of the L1 subfamilies in the various species, we inferred the partial phylogeny of Rattus sensu lato. The four Maxomys species comprise a well-defined clade separate from a monophyletic cluster that contains the two Leopoldamys and four Niviventer species. The Niviventer/Leopoldamys clade, in turn, shares a node with the clade that contains Berylmys, Sundamys, Bandicota, and Rattus sensu stricto. The evolutionary relationships that we deduced agree with and significantly extend the phylogeny of Rattus sensu lato established by other molecular criteria. Furthermore, the L1 amplification events scored here produced a unique phylogenetic tree, that is, in no case did a character (a given L1 amplification event) appear on more than one branch. The lack of homoplasy found in this study supports the robustness of L1 amplification events as phylogenetic markers for the study of mammalian evolution.

Retrotransposon Mys was active during evolution of the Peromyscus leucopus-maniculatus complex

Mys is a retrovirus-like transposable element found throughout the genus Peromyscus. Several mys subfamilies identified on the basis of restriction site variation occur in more than one species. The distribution of these subfamilies is consistent with the accepted species phylogeny, suggesting that mys was present in the ancestor of Peromyscus and has been active through much of the evolution of this genus. Quantitative Southern blot analysis was used to examine the variability of subfamilies in P. leucopus and maniculatus. We found that subfamilies with phylogenetically narrow distributions were more variable in copy number both within and between species than subfamilies with a broader distribution. Taken together, our data suggest that mys has undergone multiple rounds of transposition since the peromyscine radiation, and that five subfamilies have been amplified during the evolution of the leucopus-maniculatus species complex.

Sporadic amplification of ID elements in rodents

ID sequences are members of a short interspersed element (SINE) repetitive DNA family within the rodent genome. The copy number of individual ID elements varies by up to three orders of magnitude between species. This amplification has been highly sporadic in the order Rodentia and does not follow any phylogenetic trend. Using library screening and dot-blot analysis, we estimate there are 25,000 copies of ID elements in the deer mouse, 1,500 copies in the gerbil (both cricetid rodents), and 60,000 copies of either ID or ID-like elements in a sciurid rodent (squirrel). By dot-blot analysis, we estimate there are 150,000, 4,000, 1,000, and 200 copies of ID elements in the rat, mouse, hamster, and guinea pig, respectively (which is consistent with previous reports) and 200 copies in the hystricognath rodent, nutria. Therefore, a rapid amplification took place not only after the divergence of rat and mouse but also following the deer mouse (Peromyscus) and hamster split, with no evidence of increased amplifications in hystricognath rodents. No notable variations of sequences from the BC1 genes of several myomorphic rodents were observed that would possibly explain the varied levels of ID amplification. We did observe subgenera and species-group-specific variation in the ID core sequence of the BC1 gene within the genus Peromyscus. Sequence analysis of cloned ID elements in Peromyscus show most ID elements in this genus arose prior to Peromyscus subgenus divergence. Correspondence of the consensus sequence of individual ID elements in gerbil and deer mouse further confirms BC1 as a master gene in ID amplification. Several possible mechanisms responsible for the quantitative variations are explored.

Computer simulation of transposable element evolution: random template and strict master models

It has been proposed that the most extensively studied mammalian retrotransposons replicate by some form of a master template model. This conclusion has been drawn largely from DNA sequence analysis and is based on phylogenetic tree topology, the presence and ordering of shared variants, the degree of divergence between elements within a subfamily, and the shape of the distribution of pairwise differences between elements. To investigate how robust these parameters are as predictors of the model of transposition, computer simulations of the two most extreme transposition models, the Random Template Model and the Strict Master Model, were carried out. A prototype of a computer simulator for studying retrotransposition is presented. The simulator is a versatile digital workbench that maintains DNA sequence data and allows manipulation of a range of factors including reverse transcriptase and in situ mutation rates, transposition template, and transposition rate. All parameters previously used as predictors of the model of transposition were markedly different for the two extreme models when evaluated using large sample sizes of sequences from experiments simulating up to 15 million years of evolution.