Amplification of an ancestral mammalian L1 family of long interspersed repeated DNA occurred just before the murine radiation

The retroelement Lx9 puts a brake on the immune response to virus infection

The notion that mobile units of nucleic acid known as transposable elements can operate as genomic controlling elements was put forward over six decades ago^1,2. However, it was not until the advancement of genomic sequencing technologies that the abundance and repertoire of transposable elements were revealed, and they are now known to constitute up to two-thirds of mammalian genomes^3,4. The presence of DNA regulatory regions including promoters, enhancers and transcription-factor-binding sites within transposable elements^5-8 has led to the hypothesis that transposable elements have been co-opted to regulate mammalian gene expression and cell phenotype^8-14. Mammalian transposable elements include recent acquisitions and ancient transposable elements that have been maintained in the genome over evolutionary time. The presence of ancient conserved transposable elements correlates positively with the likelihood of a regulatory function, but functional validation remains an essential step to identify transposable element insertions that have a positive effect on fitness. Here we show that CRISPR-Cas9-mediated deletion of a transposable element-namely the LINE-1 retrotransposon Lx9c11-in mice results in an exaggerated and lethal immune response to virus infection. Lx9c11 is critical for the neogenesis of a non-coding RNA (Lx9c11-RegoS) that regulates genes of the Schlafen family, reduces the hyperinflammatory phenotype and rescues lethality in virus-infected Lx9c11^-/- mice. These findings provide evidence that a transposable element can control the immune system to favour host survival during virus infection.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Multiple and Independent Phases of Transposable Element Amplification in the Genomes of Piciformes (Woodpeckers and Allies)

The small and conserved genomes of birds are likely a result of flight-related metabolic constraints. Recombination-driven deletions and minimal transposable element (TE) expansions have led to continually shrinking genomes during evolution of many lineages of volant birds. Despite constraints of genome size in birds, we identified multiple waves of amplification of TEs in Piciformes (woodpeckers, honeyguides, toucans, and barbets). Relative to other bird species’ genomic TE abundance (< 10% of genome), we found ∼17–30% TE content in multiple clades within Piciformes. Several families of the retrotransposon superfamily chicken repeat 1 (CR1) expanded in at least three different waves of activity. The most recent CR1 expansions (∼4–7% of genome) preceded bursts of diversification in the woodpecker clade and in the American barbets + toucans clade. Additionally, we identified several thousand polymorphic CR1 insertions (hundreds per individual) in three closely related woodpecker species. Woodpecker CR1 insertion polymorphisms are maintained at lower frequencies than single nucleotide polymorphisms indicating that purifying selection is acting against additional CR1 copies and that these elements impose a fitness cost on their host. These findings provide evidence of large scale and ongoing TE activity in avian genomes despite continual constraint on genome size.

Contrasting Rates of LINE-1 Amplification among New World Primates of the Atelidae Family

LINE-1 (L1) retrotransposons constitute the dominant category of transposons in mammalian genomes. L1 elements are active in the vast majority of mammals, and only a few cases of L1 extinction have been documented. The only possible case of extinction in primates was suggested for South American spider monkeys. However, these previous studies were based on a single species. We revisited this question with a larger phylogenetic sample, covering all 4 genera of Atelidae and 3 species of spider monkeys. We used an enrichment method to clone recently inserted L1 elements and performed an evolutionary analysis of the sequences. We were able to identify young L1 elements in all taxa, suggesting that L1 is probably still active in all Atelidae examined. However, we also detected considerable variations in the proportion of recent elements indicating that the rate of L1 amplification varies among Atelidae by a 3-fold factor. The extent of L1 amplification in Atelidae remains overall lower than in other New World monkeys. Multiple factors can affect the amplification of L1, such as the demography of the host and the control of transposition. These factors are discussed in the context of host life history.

Differential Effect of Selection against LINE Retrotransposons among Vertebrates Inferred from Whole-Genome Data and Demographic Modeling

Variation in LINE composition is one of the major determinants for the substantial size and structural differences among vertebrate genomes. In particular, the larger genomes of mammals are characterized by hundreds of thousands of copies from a single LINE clade, L1, whereas nonmammalian vertebrates possess a much greater diversity of LINEs, yet with orders of magnitude less in copy number. It has been proposed that such variation in copy number among vertebrates is due to differential effect of LINE insertions on host fitness. To investigate LINE selection, we deployed a framework of demographic modeling, coalescent simulations, and probabilistic inference against population-level whole-genome data sets for four model species: one population each of threespine stickleback, green anole, and house mouse, as well as three human populations. Specifically, we inferred a null demographic background utilizing SNP data, which was then exploited to simulate a putative null distribution of summary statistics that was compared with LINE data. Subsequently, we applied the inferred null demographic model with an additional exponential size change parameter, coupled with model selection, to test for neutrality as well as estimate the strength of either negative or positive selection. We found a robust signal for purifying selection in anole and mouse, but a lack of clear evidence for selection in stickleback and human. Overall, we demonstrated LINE insertion dynamics that are not in accordance to a mammalian versus nonmammalian dichotomy, and instead the degree of existing LINE activity together with host-specific demographic history may be the main determinants of LINE abundance.

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.

Two new species of fossil Leggadina (Rodentia: Muridae) from Northwestern Queensland

Only three species of fossil murine have been described to date in Australia even though they are often found in fossil deposits and can be highly useful in understanding environmental change over time. Until now the genus Leggadina, a group of short-tailed mice that is particularly well adapted to an arid environment, was only known from two extant species: L. forresti and L. lakedownensis. Here two new fossil species of the genus are described from sites in northwestern Queensland. Leggadina gregoriensis sp. nov. comes from the Early Pleistocene Rackham’s Roost Site in the Riversleigh World Heritage Area and Leggadina macrodonta sp. nov. is from the Plio-Pleistocene Site 5C at Floraville Station. The evolution of the genus Leggadina and the lineage’s response to palaeoecological factors is considered. Taphonomy of the two fossil deposits is examined and shows marked differences in both faunal composition of the assemblages and preservation. Description of L. gregoriensis and L. macrodonta extends the known temporal range of the Leggadina lineage by over 2 million years.

LINE-1 retrotransposons: from 'parasite' sequences to functional elements

Long interspersed nuclear elements-1 (LINE-1) are the most abundant and active retrotransposons in the mammalian genomes. Traditionally, the occurrence of LINE-1 sequences in the genome of mammals has been explained by the selfish DNA hypothesis. Nevertheless, recently, it has also been argued that these sequences could play important roles in these genomes, as in the regulation of gene expression, genome modelling and X-chromosome inactivation. The non-random chromosomal distribution is a striking feature of these retroelements that somehow reflects its functionality. In the present study, we have isolated and analysed a fraction of the open reading frame 2 (ORF2) LINE-1 sequence from three rodent species, Cricetus cricetus, Peromyscus eremicus and Praomys tullbergi. Physical mapping of the isolated sequences revealed an interspersed longitudinal AT pattern of distribution along all the chromosomes of the complement in the three genomes. A detailed analysis shows that these sequences are preferentially located in the euchromatic regions, although some signals could be detected in the heterochromatin. In addition, a coincidence between the location of imprinted gene regions (as Xist and Tsix gene regions) and the LINE-1 retroelements was also observed. According to these results, we propose an involvement of LINE-1 sequences in different genomic events as gene imprinting, X-chromosome inactivation and evolution of repetitive sequences located at the heterochromatic regions (e.g. satellite DNA sequences) of the rodents' genomes analysed.

Evolutionary conservation of the functional modularity of primate and murine LINE-1 elements

LINE-1 (L1) retroelements emerged in mammalian genomes over 80 million years ago with a few dominant subfamilies amplifying over discrete time periods that led to distinct human and mouse L1 lineages. We evaluated the functional conservation of L1 sequences by comparing retrotransposition rates of chimeric human-rodent L1 constructs to their parental L1 counterparts. Although amino acid conservation varies from ∼35% to 63% for the L1 ORF1p and ORF2p, most human and mouse L1 sequences can be functionally exchanged. Replacing either ORF1 or ORF2 to create chimeric human-mouse L1 elements did not adversely affect retrotransposition. The mouse ORF2p retains retrotransposition-competency to support both Alu and L1 mobilization when any of the domain sequences we evaluated were substituted with human counterparts. However, the substitution of portions of the mouse cys-domain into the human ORF2p reduces both L1 retrotransposition and Alu trans-mobilization by 200-1000 fold. The observed loss of ORF2p function is independent of the endonuclease or reverse transcriptase activities of ORF2p and RNA interaction required for reverse transcription. In addition, the loss of function is physically separate from the cysteine-rich motif sequence previously shown to be required for RNP formation. Our data suggest an additional role of the less characterized carboxy-terminus of the L1 ORF2 protein by demonstrating that this domain, in addition to mediating RNP interaction(s), provides an independent and required function for the retroelement amplification process. Our experiments show a functional modularity of most of the LINE sequences. However, divergent evolution of interactions within L1 has led to non-reciprocal incompatibilities between human and mouse ORF2 cys-domain sequences.

Inactivation dates of the human and guinea pig vitamin C genes

The capacity to biosynthesize ascorbic acid has been lost in a number of species including primates, guinea pigs, teleost fishes, bats, and birds. This inability results from mutations in the GLO gene coding for L-gulono-γ-lactone oxidase, the enzyme responsible for catalyzing the last step in the vitamin C biosynthetic pathway. We analyzed available primate and rodent GLO gene sequences to determine their evolutionary history. We used a method based on sequence comparisons of lineages with and without functional GLO genes to calculate inactivation dates of 61 and 14 MYA for the primate and guinea pig genes, respectively. These estimates are consistent with previous phylogeny-based estimates. An analysis of transposable element distribution in the primate and rodent GLO sequences did not reveal conclusive evidence that illegitimate recombination between repeats has contributed to the loss of exons in the primate and guinea pig genes.

X accumulation of LINE-1 retrotransposons in Tokudaia osimensis, a spiny rat with the karyotype XO

The observation that LINE-1 transposable elements are enriched on the X in comparison to the autosomes led to the hypothesis that LINE-1s play a role in X chromosome inactivation. If this hypothesis is correct, loss of LINE-1 activity would be expected to result in species extinction or in an alternate pathway of dosage compensation. One such alternative pathway would be to evolve a karyotype that does not require dosage compensation between the sexes. Two of the three extant species of the Ryukyu spiny rat Tokudaia have such a karyotype; both males and females are XO. We asked whether this karyotype arose due to loss of LINE-1 activity and thus the loss of a putative component in the X inactivation pathway. Although XO Tokudaia has no need for dosage compensation, LINE-1s have been recently active in Tokudaia osimensis and show higher density on the lone X than on the autosomes.

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.

LINE-1 amplification accompanies explosive genome repatterning in rodents

Transposable elements (TEs) sometimes induce karyotypic changes following recombination, breakage and rearrangement. We used FISH and Southern blot analyses to investigate the amount and distribution of LINE-1 retrotransposons in rodents (genus Taterillus, Muridae, Gerbillinae) that have recently undergone an important genome repatterning. Our results were interpreted in a known phylogenetic framework and clearly showed that LINE-1 elements were greatly amplified and non-randomly distributed in the most rearranged karyotypes. A comparison between FISH and conventional banding patterns provided evidence that LINE-1 insertion sites and chromosome breakpoints were not strongly correlated, thus suggesting that LINE-1 amplification subsequently accompanied Taterillus chromosome evolution. Similar patterns are observed in some cases of genomic stresses (hybrid genomes, cancer and DNA-damaged cells) and usually associated with DNA hypomethylation. We propose that intensively repatterned genomes face transient stress phases during which some epigenetic features, such as DNA methylation, are relaxed, thus allowing TE amplification.

LINE-1 distribution in Afrotheria and Xenarthra: implications for understanding the evolution of LINE-1 in eutherian genomes

Long interspersed nuclear elements (LINEs) comprise about 21% of the human genome (of which L1 is most abundant) and are preferentially accumulated in AT-rich regions, as well as the X and Y chromosomes. Most knowledge of L1 distribution in mammals is restricted to human and mouse. Here we report the first investigation of L1 distribution in the genomes of a wide variety of eutherian mammals, including species in the two basal clades, Afrotheria and Xenarthra. Our results show L1 accumulation on the X of all eutherian mammals, an observation consistent with an ancestral involvement of these elements in the X-inactivation process (the Lyon repeat hypothesis). Surprisingly, conspicuous accumulation of L1 in AT-rich regions of the genome was not observed in any species outside of Euarchontoglires (represented by human, mouse and rabbit). Although several features were common to most species investigated, our comprehensive survey shows that the patterns observed in human and mouse are, in many aspects, far from typical for all mammals. We discuss these findings with reference to models that have previously been proposed to explain the AT distribution bias of L1 in human and mouse, and how this relates to the evolution of these elements in other eutherian genomes.

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.

Adaptive evolution in LINE-1 retrotransposons

We traced the sequence evolution of the active lineage of LINE-1 (L1) retrotransposons over the last approximately 25 Myr of human evolution. Five major families (L1PA5, L1PA4, L1PA3B, L1PA2, and L1PA1) of elements have succeeded each other as a single lineage. We found that part of the first open-reading frame (ORFI) had a higher rate of nonsynonymous (amino acid replacement) substitution than synonymous substitution during the evolution of the ancestral L1PA5 through the L1PA3B families. This segment encodes the coiled coil region of the protein-protein interaction domain of the ORFI protein (ORFIp). Statistical analysis of these changes indicates that positive selection had been acting on this region. In contrast, the coiled coil segment hardly changed during the evolution of the L1PA3B to the present L1PA1 family. Therefore, selective pressure on the coiled coil segment has changed over time. We suggest that the fast rate of amino acid replacement in the coiled coil segment reflects the adaptation of L1 either to a changing genomic environment or to host repression factors. In contrast, the second open-reading frame and the nucleic acid-binding domain of the first open-reading frame are extremely well conserved, attesting to the strong purifying selection acting on these regions.

Identification of a deletion hotspot on distal mouse chromosome 4 by YAC fingerprinting

Using repetitive elements as probes, genomic DNA fingerprints of four randomly selected yeast artificial chromosome (YAC) clones (two human and two mouse-derived YAC) were analyzed to determine the mutation level following X-ray exposure. Because the repetitive probes were derived from the mammalian host DNA, most of the fingerprint bands originated from the artificial chromosomes and not from the yeast genome. For none of the YAC clones was the mutation frequency elevated following X-ray exposure. However, for one mouse-derived YAC, the mutation level was unusually high (7%; 42 mutants of 607 clones analyzed), whereas for the other three YACs, the mutation level was nearly 0%. Surprisingly, 40 of the 42 mutations were deletions occurring only at three of the 20 mouse specific fingerprint bands. One of the frequently deleted fragments was cloned, sequenced and mapped to distal mouse chromosome 4, which has been repeatedly reported to be the most unstable region of the whole mouse genome, associated with various tumors. Deletion mapping of six YAC mutants revealed this fragment to be completely deleted in four YACs. In the other two mutants, recombination occurred within the fragment, in each case initiated at the same LINE-1 element. In conclusion, the presented YAC fingerprint is a useful tool for detecting and characterizing unstable regions in mammalian genomes.

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.

The end of the LINE?: lack of recent L1 activity in a group of South American rodents

L1s (LINE-1: Long Interspersed Nuclear Element 1) are present in all mammals examined to date. They occur in both placental mammals and marsupials and thus are thought to have been present in the genome prior to the mammalian radiation. This unusual conservation of a transposable element family for over 100 million years has led to speculation that these elements provide an advantage to the genomes they inhabit. We have recently identified a group of South American rodents, including rice rats (Oryzomys), in which L1s appear to be quiescent or extinct. Several observations support this conclusion. First, genomic Southern blot analysis fails to reveal genus-specific bands in Oryzomys. Second, we were unable to find recently inserted elements. Procedures to enrich for young elements did not yield any with an intact open reading frame for reverse transcriptase; all elements isolated had numerous insertions, deletions, and stop codons. Phylogenetic analysis failed to yield species-specific clusters among the L1 elements isolated, and all Oryzomys sequences had numerous private mutations. Finally, in situ hybridization of L1 to Oryzomys chromosomes failed to reveal the characteristic L1 distribution in Oryzomys with either a homologous or heterologous probe. Thus, Oryzomys is a viable candidate for L1 extinction from a mammalian host.

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.

T cells are necessary and critical for xenograft rejection in new concordant cardiac xenotransplant model

BACKGROUND

A new vascularized concordant xenotransplant model using the Chinese hamster as donor and mouse as recipient species is reported. This model takes advantage of the wealth of informative immune reagents and knockout and transgenic backgrounds available for the mouse.

METHODS

Heterotopic auxillary cardiac transplantation was performed. The mean survival time was assessed by daily palpation. Xenoreactive antibody production was measured by flow cytometry, and cardiac xenografts were examined by light microscopy.

RESULTS

The tempo of xenograft rejection in this model is consistent with concordant species combination. IgM and IgG3 responses were not critical for the concordant xenograft rejection. Long-term survival (>100 days) of the concordant cardiac xenografts was observed without any immunosuppression in nude mice. Reconstitution of nude mice with CD3+ T cells induced the xenograft rejection in 5.7 days (P<0.01).

CONCLUSION

This new concordant cardiac xenotransplant model demonstrates that T-dependent xenogeneic immune response is necessary and critical for the xenograft rejection.

Molecular evolution of two lineages of L1 (LINE-1) retrotransposons in the california mouse, Peromyscus californicus

The large number of L1 [long interspersed elements (LINE)-1] sequences found in the genome is due to the insertion of copies of the retrotransposon over evolutionary time. The majority of copies appear to be replicates of a few active, or "master" templates. A continual replacement of master templates over time gives rise to lineages distinguishable by their own unique set of shared-sequence variants. A previous analysis of L1 sequences in deer mice, Peromyscus maniculatus and P. leucopus, revealed two active L1 lineages, marked by different rates of evolution, whose most recent common ancestor predates the expansion of the Peromyscus species. Here we exploit lineage-specific, shared-sequence variants to reveal a paucity of Lineage 2 sequences in at least one species, P. californicus. The dearth of Lineage 2 copies in P. californicus suggests that Lineage 2 may have been unproductive until after the most recent common ancestor of P. californicus and P. maniculatus. We also show that Lineage 1 appears to have a higher rate of evolution in P. maniculatus relative to either P. californicus or P. leucopus. As a phylogenetic tool, L1 lineage-specific variants support a close affinity between P. californicus and P. eremicus relative to the other species examined.

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.

Rapid evolution of a young L1 (LINE-1) clade in recently speciated Rattus taxa

L1 elements are retrotransposons that have been replicating and evolving in mammalian genomes since before the mammalian radiation. Rattus norvegicus shares the young L1mlvi2 clade only with its sister taxon, Rattus cf moluccarius. Here we compared the L1mlvi2 clade in these recently diverged species and found that it evolved rapidly into closely related but distinct clades: the L1mlvi2-rm clade (or subfamily), characterized here from R. cf moluccarius, and the L1mlvi2-rn clade, originally described in R. norvegicus. In addition to other differences, these clades are distinguished by a cluster of amino acid replacement substitutions in ORF I. Both rat species contain the L1mlvi2-rm clade, but the L1mlvi2-rn clade is restricted to R. norvegicus. Therefore, the L1mlvi2-rm clade arose prior to the divergence of R. norvegicus and R. cf moluccarius, and the L1mlvi2-rn clade amplified after their divergence. The total number of L1mlvi2-rm elements in R. cf moluccarius is about the same as the sum of the L1mlvi2-rm and L1mlvi2-rn elements in R. norvegicus. The possibility that L1 amplification is in some way limited so that the two clades compete for replicative supremacy as well as the implications of the other distinguishing characteristic of the L1mlvi2-rn and L1mlvi2-rm clades are discussed.

Recombination creates novel L1 (LINE-1) elements in Rattus norvegicus

Mammalian L1 (long interspersed repeated DNA. LINE-1) retrotransposons consist of a 5' untranslated region (UTR) with regulatory properties, two protein encoding regions (ORF I, ORF II, which encodes a reverse transcriptase) and a 3' UTR. L1 elements have been evolving in mammals for > 100 million years and this process continues to generate novel L1 subfamilies in modern species. Here we characterized the youngest known subfamily in Rattus norvegicus, L1mlvi2, and unexpectedly found that this element has a dual ancestry. While its 3' UTR shares the same lineage as its nearest chronologically antecedent subfamilies, L13 and L14, its ORF I sequence does not. The L1mlvi2 ORF I was derived from an ancestral ORF I sequence that was the evolutionary precursor of the L13 and L14 ORF I. We suggest that an ancestral ORF I sequence was recruited into the modern L1mlvi2 subfamily by recombination that possibly could have resulted from template strand switching by the reverse transcriptase during L1 replication. This mechanism could also account for some of the structural features of rodent L1 5' UTR and ORF I sequences including one of the more dramatic features of L1 evolution in mammals, namely the repeated acquisition of novel 5' UTRs.

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

We identified and characterized the relics of an ancient rodent L1 family, referred to as Lx, which was extensively amplified at the time of the murine radiation about 12 million years ago, and which we showed was ancestral to the modern L1 families in rat and mouse. Here we have extended our analysis of the Lx amplification by examining more murine and nonmurine species for Lx sequences using both blot hybridization and the polymerase chain reaction for a total of 36 species. In addition we have determined the relative copy number and sequence divergence, or age, of Lx elements in representative murine genera. Our results show that while Lx sequences are confined to murine genera, the extent of the amplification was different in the different murine lineages, indicating that the amplification of Lx did not precede, but was coincident with, the murine radiation. The implications of our findings for the evolutionary dynamics of L1 families and the utility of ancestral amplification events for systematics are discussed.

Identification of a new, abundant superfamily of mammalian LTR-transposons

A new superfamily of mammalian transposable genetic elements is described with an estimated 40,000 to 100,000 members in both primate and rodent genomes. Sequences known before as MT, ORR-1, MstII, MER15 and MER18 are shown to represent (part of) the long terminal repeats of retrotransposon-like elements related to THE1 in humans. These transposons have structural similarities to retroviruses. However, the putative product of a 1350 base pair open reading frame detected in the consensus internal sequence of THE1 does not resemble retroviral proteins. The elements are named 'Mammalian apparent LTR-retrotransposons' (MaLRs). The internal sequence is usually found to be excised. Their presence in rodents, artiodactyls, lagomorphs, and primates, the divergence of the individual elements from their consensus, and the existence of a probably orthologous element in mouse and man suggest that the first MaLRs were distributed before the radiation of eutherian mammals 80-100 million years ago. MaLRs may prove to be very helpful in determining the evolutionary branching pattern of mammalian orders and suborders.

Molecular evidence that the spiny mouse (Acomys) is more closely related to gerbils (Gerbillinae) than to true mice (Murinae)

Spiny mice of the genus Acomys traditionally have been classified as members of the Murinae, a subfamily of rodents that also includes rats and mice with which spiny mice share a complex set of morphological characters, including a unique molar pattern. The origin and evolution of this molar pattern, documented by many fossils from Southern Asia, support the hypothesis of the monophyly of Acomys and all other Murinae. This view has been challenged by immunological studies that have suggested that Acomys is as distantly related to mice (Mus) as are other subfamilies (e.g., hamsters: Cricetinae) of the muroid rodents. We present molecular evidence derived from DNA.DNA hybridization data that indicate that the spiny mouse Acomys and two African genera of Murinae, Uranomys and Lophuromys, constitute a monophyletic clade, a view that was recently suggested on the basis of dental characters. However, our DNA.DNA hybridization data also indicate that the spiny mice (Acomys) are more closely related to gerbils (Gerbillinae) than to the true mice and rats (Murinae) with which they have been classified. Because Acomys and the brush-furred mice Uranomys and Lophuromys share no derived morphological characters with the Gerbillinae, their murine morphology must have evolved by convergence, including the molar pattern previously considered to support the monophyly of the Murinae.

The evolution of long interspersed repeated DNA (L1, LINE 1) as revealed by the analysis of an ancient rodent L1 DNA family

All modern mammals contain a distinctive, highly repeated (> or = 50,000 members) family of long interspersed repeated DNA called the L1 (LINE 1) family. While the modern L1 families were derived from a common ancestor that predated the mammalian radiation approximately 80 million years ago, most of the members of these families were generated within the last 5 million years. However, recently we demonstrated that modern murine (Old World rats and mice) genomes share an older long interspersed repeated DNA family that we called Lx. Here we report our analysis of the DNA sequence of Lx family members and the relationship of this family to the modern L1 families in mouse and rat. The extent of DNA sequence divergence between Lx members indicates that the Lx amplification occurred about 12 million years ago, around the time of the murine radiation. Parsimony analysis revealed that Lx elements were ancestral to both the modern rat and mouse L1 families. However, we found that few if any of the evolutionary intermediates between the Lx and the modern L1 families were extensively amplified. Because the modern L1 families have evolved under selective pressure, the evolutionary intermediates must have been capable of replication. Therefore, replication-competent L1 elements can reside in genomes without undergoing extensive amplification. We discuss the bearing of our findings on the evolution of L1 DNA elements and the mammalian genome.

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

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 L1. 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 L1Pm55, 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.

Software tools for analyzing pairwise alignments of long sequences

Pairwise comparison of long stretches of genomic DNA sequence can identify regions conserved across species, which often indicate functional significance. However, the novel insights frequently must be windowed from a flood of information; for instance, running an alignment program on two 50-kilobase sequences might yield over a hundred pages of alignments. Direct inspection of such a volume of printed output is infeasible, or at best highly undesirable, and computer tools are needed to summarize the information, to assist in its analysis, and to report the findings. This paper describes two such software tools. One tool prepares publication-quality pictorial representations of alignments, while another facilitates interactive browsing of pairwise alignment data. Their effectiveness is illustrated by comparing the beta-like globin gene clusters between humans and rabbits. A second example compares the chloroplast genomes of tobacco and liverwort.