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

Loss of LINE-1 activity in the megabats

LINE-1 (L1) retrotransposons are the most abundant type of mammalian retroelement. They have profound effects on genome plasticity and have been proposed to fulfill essential host functions, yet it remains unclear where they lie on the spectrum from parasitism to mutualism. Their ubiquity makes it difficult to determine the extent of their effects on genome evolution and gene expression because of the relative dearth of animal models lacking L1 activity. We have isolated L1 sequences from 11 megabat species by a method that enriches for recently inserted L1s and have done a bioinformatic examination of L1 sequences from a 12th species whose genome was recently shotgun sequenced. An L1 extinction event appears to have occurred at least 24 million years ago (MYA) in an ancestor of the megabats. The ancestor was unusual in having maintained two highly divergent long-term L1 lineages with different levels of activity, which appear, on an evolutionary scale, to have simultaneously lost that activity. These megabat species can serve as new animal models to ask what effect loss of L1 activity has on mammalian genome evolution and gene expression.

Assessment and disease comparisons of hybrid developmental defects

Rodents of the genus Peromyscus are among the most common North American mammals. Crosses between natural populations of two of these species, P. maniculatus (BW) and P. polionotus (PO), produce parent-of-origin effects on growth and development. BW females mated to PO males produce growth-retarded offspring. In contrast, PO females mated to BW males produce overgrown but dysmorphic conceptuses. Variation in imprinted loci and control of genomic imprinting appear to underlie the hybrid effects. Prior morphological and genetic analyses have focused on placental and post-natal growth. Here, we assess the frequency and scope of embryonic defects. The most frequent outcome of the PO x BW cross is death prior to embryonic day 13. Conceptuses lacking an embryo proper are also observed as in gestational trophoblast disease. Among the common embryonic phenotypes described and tabulated are edema, blood vessel enlargement/hemorrhaging, macroglossia, retention of nucleated erythrocytes, placentomegaly. We investigate expression of loci known to be mis-regulated in human growth/placental disorders and/or mouse knockouts with similar phenotypes. These loci are Igf2, Cdkn1c, Grb10, Gpc3, Phlda2 and Rb1. All exhibited significant differences in either placental or embryonic expression levels at one or more of the three timepoints examined. The data underscore the importance of placental gene expression on embryonic defects. We suggest that the hybrid defects offer a novel system to understand how natural allelic combinations interact to produce disease phenotypes. We propose that such interactions and their resulting epimutations may similarly underlie the phenotypic and causal heterogeneity seen in many human diseases.

Changes in cell cycle and extracellular matrix gene expression during placental development in deer mouse (Peromyscus) hybrids

Crosses between two species of the rodent genus Peromyscus produce defects in both growth and development. The defects are pronounced in the hybrid placentas. Peromyscuys maniculatus (strain BW) females mated to P. polionotus (strain PO) males produce placentas half the size of the parental species, as well as growth-retarded embryos. In contrast, PO females mated to BW males result in defective conceptuses that display embryonic and placental overgrowth. These ‘parent-of-origin’-dependent phenotypes are consistent with previous studies that demonstrated altered expression of imprinted genes and genetic linkage of the overgrowth phenotypes to imprinted domains. In the present study, we take a broader approach in assessing perturbations in hybrid placental gene expression through the use of Mus musculus cDNA microarrays. In verifying classes of genes identified in microarray screens differentially regulated during hybrid placental development, we focused on those influencing the cell cycle and extracellular matrix (ECM). Our work suggests that cell cycle regulators at the G1/S phase check-point are downregulated in the large hybrid placenta, whereas the small hybrid placenta is more variable. The ECM genes are typically downstream targets of cell cycle regulation and their misregulation is consistent with many of the dysmorphic phenotypes. Thus, these data suggest imbalances in proliferation and differentiation in hybrid placentation.

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.

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.

Deer Mice As Laboratory Animals

Although laboratory mice (Mus) and rats (Rattus) are the most widely used research rodents, deer mice (Peromyscus maniculatus) and their congeneric species are favored as nontraditional alternatives for some purposes. Mice of the native genus Peromyscus are the most abundant and widely distributed rodents in North America. They occur in a great diversity of habitats and play a significant role in natural ecosystems. Because of their abundance, peromyscines are commonly hosts for larva of ticks that transmit Lyme disease bacteria, and they are implicated in several other vector-borne diseases. Deer mice also are the principal carriers of the virus that causes hantaviral pulmonary syndrome, or "Four Corners disease." Deer mice are useful as laboratory models for a variety of other types of pure and applied research. They are easily maintained and bred in captivity using the husbandry protocols developed for other small laboratory rodent species. The Peromyscus Genetic Stock Center at the University of South Carolina maintains more than 50 laboratory-bred, well-characterized stocks of deer mice and other peromyscine species for research and educational use.

The Chironomus (Camptochironomus) tentans genome contains two non-LTR retrotransposons

A cDNA library from salivary gland cells of Chironomus tentans was screened with a probe containing the NLRCth1 non-LTR (long terminal repeat) retrotransposon from Chironomus thummi. Several positive clones were obtained and one of them, p62, was characterized by in situ hybridization and sequencing. The sequencing analysis showed that this clone contained a 4607 bp nucleotide sequence of a new transposable element that hybridized in situ to more than 100 sites over all four C. tentans chromosomes. The detailed analysis of this sequence revealed the presence of the 3'-end of open reading frame 1 (ORF1), a complete ORF2, and a 1.3-kb 3'-end untranslated region (UTR). The new element has been designated NLRCt2 (non-LTR retrotransposon 2 from C. tentans). A comparison of the nucleotide sequences of NLRCth1 and NLRCt2 showed 30% similarity in the region of ORF1 and 70% similarity in the region of ORF2. Based on the results of Southern blot analysis, two transposable elements have been found in the C. tentans genome, one of which is identical to NLRCth1 from C. thummi. This may be explained by horizontal transmission. The second element, NLRCt2, has been found in two different forms in the C. tentans genome. These can be distinguished by the presence of the 1.3-kb 3'-end UTR in one of the forms. Since the cDNA clone investigated was isolated from a tissue-specific cDNA library, the data showed that NRLCt2 is expressed in somatic cells.

DNA sequences of RAPD fragments in artiodactyls

A bovine RAPD profile, generated by a 10-mer primer, was analysed by sequencing the major fragments. Three of four different fragments showed homologies to previously characterized mammalian sequences. One was 61-66% identical to LINE sequences and another was 78.5% identical to a human chromosome 2 sequence tagged site. The third fragment was 93.1% identical to the human type 2 inositol 1,4,5-trisphosphate receptor gene. This fragment had counterparts in white-tailed deer and reindeer; fragments of slightly different size in these species showed high sequence similarity and the size differences were due to varying numbers of dinucleotide microsatellite repeats inside the fragment.

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.

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.

Rapidly evolving repetitive DNAs in a conservative genome: a test of factors that affect chromosomal evolution

The hypothesis that tandemly repeated DNA sequences may facilitate chromosomal rearrangements was tested by comparing a conservatively evolving karyotype of a bat species (Macrotus waterhousii) with data published for a rapidly evolving karyotype of an equid species (Equus zebra). Empirical data generated from the phylogenetic screening of rapidly evolving repetitive DNAs from approximately 0.1% of the M. waterhousii genome showed only one sequence that was repetitive in M. waterhousii but low in copy number or absent from the outgroup Artibeus jamaicensis. This compares to 34 such clones containing sequences which were repetitive in E. zebra but were low in copy number or absent from the outgroup Ceratotherium simum. The bat sequence represents a single family of repeated sequences, whereas six families of sequences were identified in E. zebra. Southern blot analysis suggested that the sequence from M. waterhousii is interspersed rather than tandemly repeated, as are the sequences in E. zebra. These data support the above hypothesis and suggest that species with conservatively evolving karyotypes have fewer numbers and families of rapidly evolving DNA sequences than do species such as the equids that possess a karyotype that is considered to have undergone rapid karyotypic evolution.

Comparison of chromosomal distribution of a retroposon (LINE) and a retrovirus-like element mys in Peromyscus maniculatus and P. leucopus

Chromosomal distribution for two interspersed elements (LINEs and mys) that are thought to have established their chromosomal position primarily by transposition was compared between two species of deer mice (Peromyscus leucopus and P. maniculatus). Both LINEs and mys generally produced an autosomal banding pattern reflective of G-bands and both hybridized preferentially to the sex chromosomes. The pattern on the long arm of the X was unique for each, with mys reflecting the G-bands (four bands with the telomeric most prominent) and LINE producing five equally spaced bands of equal intensity. LINE also preferentially hybridized to the short arm of the longest autosomal pair. Some aspects of these patterns are explained adequately with proposed mechanisms that would produce a non-random pattern of chromosomal distribution (i.e. both reflect autosomal G-bands and both preferentially insert into AT-rich regions characteristic of G-bands). However, other aspects such as the differences observed on the long arm of the X do not appear to fit any predictions of proposed mechanisms.

Genome organization of repetitive elements in the rodent, Peromyscus leucopus

To document the frequency and distribution of repetitive elements in Peromyscus leucopus, the white-footed mouse, a cosmid genomic library was examined. Two thousand thirteen randomly chosen recombinants, with an average insert size of 35 kb and representing 2.35% of the haploid genome of P. leucopus, were screened with probes representing microsatellites, tandem repeats, and transposable elements. Of the four dinucleotides, (GT)n was present in 87% of the clones, (CT)n was present in 59% of the clones, and (AT)n and (GC)n each was represented in our sample by a single clone (0.05%). (TCC)n was present in 8% of the clones. Of the tandem repeats, the 28S ribosomal probe and the (TTAGGG)n telomere probe were not represented in the library, whereas a heterochromatic fragment was present in 9% of the clones. A transposable element, mys, was estimated to occur in 4700 copies, whereas a long interspersed element (LINE) was estimated to occur in about 41,000 copies per haploid genome. LINE and mys occurred together in the same clones more frequently than expected on the basis of chance. Hybridizing the library to genomic DNA from P. leucopus, Reithrodontomys fulvescens, Mus musculus, and human produced general agreement between phylogenetic relatedness and intensity of hybridization. However, dinucleotide repeats appeared to account for a disproportionately high number of positive clones in the more distantly related taxa.