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

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

Insulin secretion deficits in a Prader-Willi syndrome β-cell model are associated with a concerted downregulation of multiple endoplasmic reticulum chaperones

Prader-Willi syndrome (PWS) is a multisystem disorder with neurobehavioral, metabolic, and hormonal phenotypes, caused by loss of expression of a paternally-expressed imprinted gene cluster. Prior evidence from a PWS mouse model identified abnormal pancreatic islet development with retention of aged insulin and deficient insulin secretion. To determine the collective roles of PWS genes in β-cell biology, we used genome-editing to generate isogenic, clonal INS-1 insulinoma lines having 3.16 Mb deletions of the silent, maternal- (control) and active, paternal-allele (PWS). PWS β-cells demonstrated a significant cell autonomous reduction in basal and glucose-stimulated insulin secretion. Further, proteomic analyses revealed reduced levels of cellular and secreted hormones, including all insulin peptides and amylin, concomitant with reduction of at least ten endoplasmic reticulum (ER) chaperones, including GRP78 and GRP94. Critically, differentially expressed genes identified by whole transcriptome studies included reductions in levels of mRNAs encoding these secreted peptides and the group of ER chaperones. In contrast to the dosage compensation previously seen for ER chaperones in Grp78 or Grp94 gene knockouts or knockdown, compensation is precluded by the stress-independent deficiency of ER chaperones in PWS β-cells. Consistent with reduced ER chaperones levels, PWS INS-1 β-cells are more sensitive to ER stress, leading to earlier activation of all three arms of the unfolded protein response. Combined, the findings suggest that a chronic shortage of ER chaperones in PWS β-cells leads to a deficiency of protein folding and/or delay in ER transit of insulin and other cargo. In summary, our results illuminate the pathophysiological basis of pancreatic β-cell hormone deficits in PWS, with evolutionary implications for the multigenic PWS-domain, and indicate that PWS-imprinted genes coordinate concerted regulation of ER chaperone biosynthesis and β-cell secretory pathway function.

Androgen-Binding Protein (Abp) Evolutionary History: Has Positive Selection Caused Fixation of Different Paralogs in Different Taxa of the Genus Mus?

Comparison of the androgen-binding protein (Abp) gene regions of six Mus genomes provides insights into the evolutionary history of this large murid rodent gene family. We identified 206 unique Abp sequences and mapped their physical relationships. At least 48 are duplicated and thus present in more than two identical copies. All six taxa have substantially elevated LINE1 densities in Abp regions compared with flanking regions, similar to levels in mouse and rat genomes, although nonallelic homologous recombination seems to have only occurred in Mus musculus domesticus. Phylogenetic and structural relationships support the hypothesis that the extensive Abp expansion began in an ancestor of the genus Mus. We also found duplicated Abpa27's in two taxa, suggesting that previously reported selection on a27 alleles may have actually detected selection on haplotypes wherein different paralogs were lost in each. Other studies reported that a27 gene and species trees were incongruent, likely because of homoplasy. However, L1MC3 phylogenies, supposed to be homoplasy-free compared with coding regions, support our paralog hypothesis because the L1MC3 phylogeny was congruent with the a27 topology. This paralog hypothesis provides an alternative explanation for the origin of the a27 gene that is suggested to be fixed in the three different subspecies of Mus musculus and to mediate sexual selection and incipient reinforcement between at least two of them. Finally, we ask why there are so many Abp genes, especially given the high frequency of pseudogenes and suggest that relaxed selection operates over a large part of the gene clusters.

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.

Consequence of Paradigm Shift with Repeat Landscapes in Reptiles: Powerful Facilitators of Chromosomal Rearrangements for Diversity and Evolution

Reptiles are notable for the extensive genomic diversity and species richness among amniote classes, but there is nevertheless a need for detailed genome-scale studies. Although the monophyletic amniotes have recently been a focus of attention through an increasing number of genome sequencing projects, the abundant repetitive portion of the genome, termed the "repeatome", remains poorly understood across different lineages. Consisting predominantly of transposable elements or mobile and satellite sequences, these repeat elements are considered crucial in causing chromosomal rearrangements that lead to genomic diversity and evolution. Here, we propose major repeat landscapes in representative reptilian species, highlighting their evolutionary dynamics and role in mediating chromosomal rearrangements. Distinct karyotype variability, which is typically a conspicuous feature of reptile genomes, is discussed, with a particular focus on rearrangements correlated with evolutionary reorganization of micro- and macrochromosomes and sex chromosomes. The exceptional karyotype variation and extreme genomic diversity of reptiles are used to test several hypotheses concerning genomic structure, function, and evolution.

Flying Around in the Genome: Characterization of LINE-1 in Chiroptera

L1s are transposable elements that move by a copy-and-paste mechanism that continuously increases their copy number in the genome, such that each genome has a record of the L1 history in that host lineage. They make up about 20% of the genomes of eutherian mammals and have played a major role in shaping genome evolution. Chiroptera has the lowest average genome size among mammalian orders and the only documented case of L1 extinction affecting an entire mammalian family. Herein, L1 activity and extinction are characterized in all families of the order Chiroptera using a method that enriches for the youngest lineages of L1s in the genome. In addition to the previously reported L1 extinction in Pteropodidae, L1 extinction was documented to occur in Mormoops blainvilli, but this event did not affect all species of Mormoopidae. Further, there was no evidence of concordance between the evolution of L1s and their chiropteran host. There were two L1 lineages present before the divergence of all extant bats. Both lineages are extinct in the Pteropodidae. One or the other L1 lineage is extinct in almost all bat families, but Taphozous melanopogon maintains active members of both. Most intriguingly, some families within the Rhinolophoidea retain one active L1 lineage whereas other families retain the other, creating a deep discontinuity between L1 phylogeny and chiropteran phylogeny. These results indicate that there have been numerous losses of active L1 lineages over the history of chiropteran evolution, but that all chiropteran families except Pteropodidae have retained L1 activity.

Low diversity of Angiostrongylus cantonensis complete mitochondrial DNA sequences from Australia, Hawaii, French Polynesia and the Canary Islands revealed using whole genome next-generation sequencing

Background

Rats (Rattus spp.) invaded most of the world as stowaways including some that carried the rat lungworm, Angiostrongylus cantonensis, the cause of eosinophilic meningoencephalitis in humans and other warm-blooded animals. A high genetic diversity of A. cantonensis based on short mitochondrial DNA regions is reported from Southeast Asia. However, the identity of invasive A. cantonensis is known for only a minority of countries. The affordability of next-generation sequencing for characterisation of A. cantonensis genomes should enable new insights into rat lung worm invasion and parasite identification in experimental studies.

Methods

Genomic DNA from morphologically verified A. cantonensis (two laboratory-maintained strains and two field isolates) was sequenced using low coverage whole genome sequencing. The complete mitochondrial genome was assembled and compared to published A. cantonensis and Angiostrongylus malaysiensis sequences. To determine if the commonly sequenced partial cox1 can unequivocally identify A. cantonensis genetic lineages, the diversity of cox1 was re-evaluated in the context of the publicly available cox1 sequences and the entire mitochondrial genomes. Published experimental studies available in Web of Science were systematically reviewed to reveal published identities of A. cantonensis used in experimental studies.

Results

New A. cantonensis mitochondrial genomes from Sydney (Australia), Hawaii (USA), Canary Islands (Spain) and Fatu Hiva (French Polynesia), were assembled from next-generation sequencing data. Comparison of A. cantonensis mitochondrial genomes from outside of Southeast Asia showed low genetic diversity (0.02–1.03%) within a single lineage of A. cantonensis. Both cox1 and cox2 were considered the preferred markers for A. cantonensis haplotype identification. Systematic review revealed that unequivocal A. cantonensis identification of strains used in experimental studies is hindered by absence of their genetic and geographical identity.

Conclusions

Low coverage whole genome sequencing provides data enabling standardised identification of A. cantonensis laboratory strains and field isolates. The phenotype of invasive A. cantonensis, such as the capacity to establish in new territories, has a strong genetic component, as the A. cantonensis found outside of the original endemic area are genetically uniform. It is imperative that the genotype of A. cantonensis strains maintained in laboratories and used in experimental studies is unequivocally characterised.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3491-y) contains supplementary material, which is available to authorized users.

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.

Phylogenomics

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. The abundance of genomic data for an enormous variety of organisms has enabled phylogenomic inference of many groups, and this has motivated the development of many computer programs implementing the associated methods. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Population Genomics Reveals Speciation and Introgression between Brown Norway Rats and Their Sibling Species

Murine rodents are excellent models for study of adaptive radiations and speciation. Brown Norway rats (Rattus norvegicus) are successful global colonizers and the contributions of their domesticated laboratory strains to biomedical research are well established. To identify nucleotide-based speciation timing of the rat and genomic information contributing to its colonization capabilities, we analyzed 51 whole-genome sequences of wild-derived Brown Norway rats and their sibling species, R. nitidus, and identified over 20 million genetic variants in the wild Brown Norway rats that were absent in the laboratory strains, which substantially expand the reservoir of rat genetic diversity. We showed that divergence of the rat and its siblings coincided with drastic climatic changes that occurred during the Middle Pleistocene. Further, we revealed that there was a geographically widespread influx of genes between Brown Norway rats and the sibling species following the divergence, resulting in numerous introgressed regions in the genomes of admixed Brown Norway rats. Intriguing, genes related to chemical communications among these introgressed regions appeared to contribute to the population-specific adaptations of the admixed Brown Norway rats. Our data reveals evolutionary history of the Brown Norway rat, and offers new insights into the role of climatic changes in speciation of animals and the effect of interspecies introgression on animal adaptation.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

High-throughput collagen fingerprinting of intact microfaunal remains; a low-cost method for distinguishing between murine rodent bones

RATIONALE

Microfaunal skeletal remains can be sensitive indicators of the contemporary ecosystem in which they are sampled and are often recovered in owl pellets in large numbers. Species identification of these remains can be obtained using a range of morphological criteria established for particular skeletal elements, but typically dominated by a reliance on cranial characters. However, this can induce biases under different environmental and taphonomic conditions. The aim of this research was to develop a high-throughput method of objectively identifying rodent remains from archaeological deposits using collagen fingerprinting, most notably the identification of rats from other myomorph rodents as a means to identify disturbances in the archaeofauna through the presence of invasive taxa not contemporary with the archaeological deposits.

METHODS

Collagen was extracted from complete microfaunal skeletal remains in such a manner as to leave the bones morphologically intact (i.e., weaker concentration of acid than previously used over shorter length of time). Acid-soluble collagen was then ultrafiltered into ammonium bicarbonate and digested with trypsin prior to dilution in the MALDI matrix and acquisition of peptide mass fingerprints using a matrix-assisted laser desorption/ionisation time-of-flight (MALDI-TOF) mass spectrometer.

RESULTS

Collagen fingerprinting was able to distinguish between Rattus, Mus, Apodemus and Micromys at the genus level; at the species level, R. rattus and R. norvegicus could be separated whereas A. flavicollis and A. sylvaticus could not. A total of 12,317 archaeological microvertebrate samples were screened for myomorph signatures but none were found to be invasive rats (Rattus) or mice (Mus). Of the contemporary murine fauna, no harvest mice (Micromys) were identified and only 24 field mouse (Apodemus) discovered.

CONCLUSIONS

As a result, no evidence of recent bioturbation could be inferred from the faunal remains of these archaeological deposits. More importantly this work presents a method for high-throughput screening of specific taxa and is the first application of collagen fingerprinting to microfaunal remains of archaeological specimens.

Proteomic study of the brackish water mussel Mytilopsis leucophaeata

BACKGROUND

We encountered the opportunity to study proteochemically a brackish water invertebrate animal, Mytilopsis leucophaeata, belonging to the bivalves which stem from the second half of the Cambrian Period (about 510 million years ago). This way, we were able to compare it with the vertebrate animal, the frilled shark (Chlamydoselachus anguineus) that stems from a much later period of geologic time (Permian: 245-286 MYA).

RESULTS

The mussel contains a well-adapted system of protein synthesis on the ER, protein folding on the ER, protein trafficking via COPI or clathrin-coated vesicles from endoplasmic reticulum (ER) to Golgi and plasmalemma, an equally well-developed system of actin filaments that with myosin forms the transport system for vesicular proteins and tubulin, which is also involved in ATP-driven vesicular protein transport via microtubules or transport of chromosomes in mitosis and meiosis. A few of the systems that we could not detect in M. leucophaeata in comparison with C. anguineus are the synaptic vesicle cycle components as synaptobrevin, cellubrevin (v-snare) and synaptosomal associated protein 25-A (t-snare), although one component: Ras-related protein (O-Rab1) could be involved in synaptic vesicle traffic. Another component that we did not find in M. leucophaeata was Rab11 that is involved in the tubulovesicular recycling process of H+/K+-ATPase in C. anguineus. We have not been able to trace the H+/K+-ATPase of M. leucophaeata, but Na+/K+-ATPase was present. Furthermore, we have studied the increase of percent protein expression between 1,070 MYA (the generation of the Amoeba Dictyostelium discoideum) and present (the generation of the mammal Sus scrofa = wild boar). In this time span, three proteomic uprises did occur: 600 to 500 MYA, 47.5 to 4.75 MYA, and 1.4 to 0 MYA. The first uprise covers the generation of bivalves, the second covers gold fish, chicken, brine shrimp, house mouse, rabbit, Japanese medaka and Rattus norvegicus, and the third covers cow, chimpanzee, Homo sapiens, dog, goat, Puccinia graminis and wild boar. We hypothesise that the latter two uprises are related to geological and climate changes and their compensation in protein function expression.

CONCLUSIONS

The proteomic and evolutionary data demonstrate that M. leucophaeata is a highly educatioanal animal to study.

Phylogenetic relationships in the Niviventer-Chiromyscus complex (Rodentia, Muridae) inferred from molecular data, with description of a new species

Based on molecular data for mitochondrial (Cyt b, COI) and nuclear (IRBP, GHR) genes, and morphological examinations of museum specimens, we examined diversity, species boundaries, and relationships within and between the murine genera Chiromyscus and Niviventer. Phylogenetic patterns recovered demonstrate that Niviventer sensu lato is not monophyletic but instead includes Chiromyscuschiropus, the only previously recognized species of Chiropus. To maintain the genera Niviventer and Chiropus as monophyletic lineages, the scope and definition of the genus Chiromyscus is revised to include at least three distinct species: Chiromyscuschiropus (the type species of Chiromyscus), Chiromyscuslangbianis (previously regarded as a species of Niviventer), and a new species, described in this paper under the name Chiromyscusthomasi sp. n.

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.

Phylogenetic species identification in Rattus highlights rapid radiation and morphological similarity of New Guinean species

The genus Rattus is highly speciose, the taxonomy is complex, and individuals are often difficult to identify to the species level. Previous studies have demonstrated the usefulness of phylogenetic approaches to identification in Rattus but some species, especially among the endemics of the New Guinean region, showed poor resolution. Possible reasons for this are simple misidentification, incomplete gene lineage sorting, hybridization, and phylogenetically distinct lineages that are unrecognised taxonomically. To assess these explanations we analysed 217 samples, representing nominally 25 Rattus species, collected in New Guinea, Asia, Australia and the Pacific. To reduce misidentification problems we sequenced museum specimens from earlier morphological studies and recently collected tissues from samples with associated voucher specimens. We also reassessed vouchers from previously sequenced specimens. We inferred combined and separate phylogenies from two mitochondrial DNA regions comprising 550 base pair D-loop sequences and both long (655 base pair) and short (150 base pair) cytochrome oxidase I sequences. Our phylogenetic species identification for 17 species was consistent with morphological designations and current taxonomy thus reinforcing the usefulness of this approach. We reduced misidentifications and consequently the number of polyphyletic species in our phylogenies but the New Guinean Rattus clades still exhibited considerable complexity. Only three of our eight New Guinean species were monophyletic. We found good evidence for either incomplete mitochondrial lineage sorting or hybridization between species within two pairs, R. leucopus/R. cf. verecundus and R. steini/R. praetor. Additionally, our results showed that R. praetor, R. niobe and R. verecundus each likely encompass more than one species. Our study clearly points to the need for a revised taxonomy of the rats of New Guinea, based on broader sampling and informed by both morphology and phylogenetics. The remaining taxonomic complexity highlights the recent and rapid radiation of Rattus in the Australo-Papuan region.

Mitochondrial DNA phylogeography of the Norway rat

Central Eastern Asia, foremost the area bordering northern China and Mongolia, has been thought to be the geographic region where Norway rats (Rattus norvegicus) have originated. However recent fossil analyses pointed to their origin in southern China. Moreover, whereas analyses of fossils dated the species' origin as ∼1.2–1.6 million years ago (Mya), molecular analyses yielded ∼0.5–2.9 Mya. Here, to study the geographic origin of the Norway rat and its spread across the globe we analyzed new and all published mitochondrial DNA cytochrome-b (cyt-b; N = 156) and D-loop (N = 212) sequences representing wild rats from four continents and select inbred strains. Our results are consistent with an origin of the Norway rat in southern China ∼1.3 Mya, subsequent prehistoric differentiation and spread in China and Asia from an initially weakly structured ancestral population, followed by further spread and differentiation across the globe during historic times. The recent spreading occurred mostly from derived European populations rather than from archaic Asian populations. We trace laboratory strains to wild lineages from Europe and North America and these represent a subset of the diversity of the rat; leaving Asian lineages largely untapped as a resource for biomedical models. By studying rats from Europe we made the observation that mtDNA diversity cannot be interpreted without consideration of pest control and, possibly, the evolution of rodenticide resistance. However, demographic models explored by forward-time simulations cannot fully explain the low mtDNA diversity of European rats and lack of haplotype sharing with their source from Asia. Comprehensive nuclear marker analyses of a larger sample of Norway rats representing the world are needed to better resolve the evolutionary history of wild rats and of laboratory rats, as well as to better understand the evolution of anticoagulant resistance.

Host cell/Orientia tsutsugamushi interactions: evolution and expression of syndecan-4 in Asian rodents (Rodentia, Muridae)

Scrub typhus is an acute febrile zoonotic disease and worldwide more than a billion people may be at risk for infection. Orientia tsutsugamushi, the causative agent of scrub typhus, is an obligate intracellular bacterium. Rodents are reported to be the primary reservoir hosts of the disease and according to the most recent surveys, all species within the Rattus sensu lato complex of the tribe Rattini are carriers of scrub typhus. There is no evidence that any of mouse (Mus) species serves as the primary reservoir of the bacterium even when occurring in sympatry with wild infected rats. This contrast in the host/syndecan-4 interactions between Rattini and Asian Murini may be due to intrinsic (i.e., genetic) differences. Herein we compare the sequence and expression levels of syndecan-4 (the putative cell receptor of O. tsutsugamushi) between Rattini species that are known to be natural reservoirs for the typhus agents, and Murini species that are not. Although it was not possible to conclusively link the structural variations detected in syndecan-4 with carrier status in either Rattini and Murini, our findings indicate the absence of a strong Orientia-mediated selective regime acting on gene structure. In contrast, variable spleen-specific syndecan-4 expression levels show a strong correlation between under-expression of syndecan-4 in Murini and seropositive Rattini, compared to seronegative Rattini rodents. We postulate that two divergent responses may be at work in Murini and Rattini, both linked with differential expression of syndecan-4: (i) reduced syndecan-4 transcription in Murini decreases the likelihood that the host cells will become infected by the Orientia bacterium, while (ii) reduced syndecan-4 expression in seropositive Rattini limits the pathogenicity of Orientia and consequently improves the longevity of the rat hosts. These patterns may underpin the poor carrier status of wild mice on the one hand, and the effective role of wild rats as reservoir hosts on the other.

Chromosomal evolution in Rattini (Muridae, Rodentia)

The Rattini (Muridae, Murinae) includes the biologically important model species Rattus norvegicus (RNO) and represents a group of rodents that are of clinical, agricultural and epidemiological importance. We present a comparative molecular cytogenetic investigation of ten Rattini species representative of the genera Maxomys, Leopoldamys, Niviventer, Berylmys, Bandicota and Rattus using chromosome banding, cross-species painting (Zoo-fluorescent in situ hybridization or FISH) and BAC-FISH mapping. Our results show that these taxa are characterised by slow to moderate rates of chromosome evolution that contrasts with the extensive chromosome restructuring identified in most other murid rodents, particularly the mouse lineage. This extends to genomic features such as NOR location (for example, NORs on RNO 3 are present on the corresponding chromosomes in all species except Bandicota savilei and Niviventer fulvescens, and the NORs on RNO 10 are conserved in all Rattini with the exception of Rattus). The satellite I DNA family detected and characterised herein appears to be taxon (Rattus) specific, and of recent origin (consistent with a feedback model of satellite evolution). BAC-mapping using clones that span regions responsible for the morphological variability exhibited by RNO 1, 12 and 13 (acrocentric/submetacentric) and their orthologues in Rattus species, demonstrated that the differences are most likely due to pericentric inversions as exemplified by data on Rattus tanezumi. Chromosomal characters detected using R. norvegicus and Maxomys surifer whole chromosome painting probes were mapped to a consensus sequence-based phylogenetic tree thus allowing an objective assessment of ancestral states for the reconstruction of the putative Rattini ancestral karyotype. This is thought to have comprised 46 chromosomes that, with the exception of a single pair of metacentric autosomes, were acrocentric in morphology.

Evolutionary relationships and divergence times among the native rats of Australia

BACKGROUND

The genus Rattus is highly speciose and has a complex taxonomy that is not fully resolved. As shown previously there are two major groups within the genus, an Asian and an Australo-Papuan group. This study focuses on the Australo-Papuan group and particularly on the Australian rats. There are uncertainties regarding the number of species within the group and the relationships among them. We analysed 16 mitochondrial genomes, including seven novel genomes from six species, to help elucidate the evolutionary history of the Australian rats. We also demonstrate, from a larger dataset, the usefulness of short regions of the mitochondrial genome in identifying these rats at the species level.

RESULTS

Analyses of 16 mitochondrial genomes representing species sampled from Australo-Papuan and Asian clades of Rattus indicate divergence of these two groups ~2.7 million years ago (Mya). Subsequent diversification of at least 4 lineages within the Australo-Papuan clade was rapid and occurred over the period from ~ 0.9-1.7 Mya, a finding that explains the difficulty in resolving some relationships within this clade. Phylogenetic analyses of our 126 taxon, but shorter sequence (1952 nucleotides long), Rattus database generally give well supported species clades.

CONCLUSIONS

Our whole mitochondrial genome analyses are concordant with a taxonomic division that places the native Australian rats into the Rattus fuscipes species group. We suggest the following order of divergence of the Australian species. R. fuscipes is the oldest lineage among the Australian rats and is not part of a New Guinean radiation. R. lutreolus is also within this Australian clade and shallower than R. tunneyi while the R. sordidus group is the shallowest lineage in the clade. The divergences within the R. sordidus and R. leucopus lineages occurring about half a million years ago support the hypotheses of more recent interchanges of rats between Australia and New Guinea. While problematic for inference of deeper divergences, we report that the analysis of shorter mitochondrial sequences is very useful for species identification in rats.

Revisiting the taxonomy of the Rattini tribe: a phylogeny-based delimitation of species boundaries

Background

Rodents are recognized as hosts for at least 60 zoonotic diseases and may represent a serious threat for human health. In the context of global environmental changes and increasing mobility of humans and animals, contacts between pathogens and potential animal hosts and vectors are modified, amplifying the risk of disease emergence. An accurate identification of each rodent at a specific level is needed in order to understand their implications in the transmission of diseases. Among the Muridae, the Rattini tribe encompasses 167 species inhabiting South East Asia, a hotspot of both biodiversity and emerging and re-emerging diseases. The region faces growing economical development that affects habitats, biodiversity and health. Rat species have been demonstrated as significant hosts of pathogens but are still difficult to recognize at a specific level using morphological criteria. DNA-barcoding methods appear as accurate tools for rat species identification but their use is hampered by the need of reliable identification of reference specimens. In this study, we explore and highlight the limits of the current taxonomy of the Rattini tribe.

Results

We used the DNA sequence information itself as the primary information source to establish group membership and estimate putative species boundaries. We sequenced two mitochondrial and one nuclear genes from 122 rat samples to perform phylogenetic reconstructions. The method of Pons and colleagues (2006) that determines, with no prior expectations, the locations of ancestral nodes defining putative species was then applied to our dataset. To give an appropriate name to each cluster recognized as a putative species, we reviewed information from the literature and obtained sequences from a museum holotype specimen following the ancient DNA criteria.

Conclusions

Using a recently developed methodology, this study succeeds in refining the taxonomy of one of the most difficult groups of mammals. Most of the species expected within the area were retrieved but new putative species limits were also indicated, in particular within Berylmys and Rattus genera, where future taxonomic studies should be directed. Our study lays the foundations to better investigate rodent-born diseases in South East Asia and illustrates the relevance of evolutionary studies for health and medical sciences.

Analysis of vkorc1 polymorphisms in Norway rats using the roof rat as outgroup

Background

Certain mutations in the vitamin K epoxide reductase subcomponent 1 gene (vkorc1) mediate rodent resistance to warfarin and other anticoagulants. Testing for resistance often involves analysis of the vkorc1. However, a genetic test for the roof rat (Rattus rattus) has yet to be developed. Moreover, an available roof rat vkorc1 sequence would enable species identification based on vkorc1 sequence and the evaluation of natural selection on particular vkorc1 polymorphisms in the Norway rat (R. norvegicus).

Results

We report the coding sequence, introns and 5' and 3' termini for the vkorc1 gene of roof rats (R. r. alexandrinus and R. r. frugivorus) from Uganda, Africa. Newly designed PCR primers now enable genetic testing of the roof rat and Norway rat. Only synonymous and noncoding polymorphisms were found in roof rats from Uganda. Both nominal subspecies of roof rats were indistinguishable from each other but were distinct from R. losea and R. flavipectus; however, the roof rat also shares at least three coding sequence polymorphisms with R. losea and R. flavipectus. Many of recently published vkorc1 synonymous and non-synonymous single nucleotide polymorphisms (SNPs) in Norway rats are likely SNPs from roof rats and/or other Rattus species. Tests applied to presumably genuine Norway rat vkorc1 SNPs are consistent with a role for selection in two populations carrying the derived Phe63Cys and Tyr139Cys mutations.

Conclusion

Geographic mapping of vkorc1 SNPs in roof rats should be facilitated by our report. Our assay should be applicable to most species of Rattus, which are intermediate in genetic distance from roof and Norway rats. Vkorc1-mediated resistance due to non-synonymous coding SNPs is not segregating in roof rats from Uganda. By using the roof rat sequence as a reference vkorc1, SNPs now can be assigned to the correct rat species with more confidence. Sampling designs and genotyping strategies employed so far have helped detect candidate mutations underlying vkorc1-mediated resistance, but generally provided unsuitable data to test for selection. We propose that our understanding of vkorc1-mediated evolution of resistance in rodents would benefit from the adoption of sampling and genotyping designs that enable tests for selection on vkorc1.

Dating of divergences within the Rattus genus phylogeny using whole mitochondrial genomes

The timing and order of divergences within the genus Rattus have, to date, been quite speculative. In order to address these important issues we sequenced six new whole mitochondrial genomes from wild-caught specimens from four species, Rattus exulans, Rattus praetor, Rattus rattus and Rattus tanezumi. The only rat whole mitochondrial genomes available previously were all from Rattus norvegicus specimens. Our phylogenetic and dating analyses place the deepest divergence within Rattus at approximately 3.5 million years ago (Mya). This divergence separates the New Guinean endemic R. praetor lineage from the Asian lineages. Within the Asian/Island Southeast Asian clade R. norvegicus diverged earliest at approximately 2.9Mya. R. exulans and the ancestor of the sister species R. rattus and R. tanezumi subsequently diverged at approximately 2.2Mya, with R. rattus and R. tanezumi separating as recently as approximately 0.4Mya. Our results give both a better resolved species divergence order and diversification dates within Rattus than previous studies.

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.

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.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

The rise and fall of the CR1 subfamily in the lineage leading to penguins

The evolution of penguins has been investigated extensively, although inconclusively, by morphologists, biogeographers and molecular phylogeneticists. We investigated this issue using retroposon analysis of insertions of CR1, which is a member of the LINE (long interspersed element) family, in the genomes of penguins and penguin relatives. The retroposon method is a powerful tool for identifying monophyletic groups. Because retroposons often show different relative frequencies of retroposition during evolution, it is first necessary to identify a certain subgroup that was specifically active during the period when the species in question diverged. Hence, we systematically analyzed many CR1 members isolated from penguin and penguin-related genomes. These CR1s are divided into at least three distinct subgroups that share diagnostic nucleotide insertions and/or deletions, namely, penguin CR1 Sph I, Sph II type A and Sph II type B. The analysis of the inserted retroposons by PCR revealed that different CR1 subfamilies or types had amplified at different rates among different periods during penguin evolution. Namely, the penguin CR1 Sph I subfamily had higher rates of retroposition in a common ancestor of all orders examined in this study or at least in a common ancestor of all extant penguins, and the subfamily Sph II type A also had the same tendency. Therefore, these CR1 members can be used to elucidate the phylogenetic relationships of Sphenisciformes (penguins) among different avian orders. In contrast, the penguin CR1 Sph II type B subfamily had higher rates of retroposition just before and after the emergence of the extant genera in Spheniscidae, suggesting that they are useful for elucidating the intra-relationships among extant penguins. This is the first report for the characterization among the members of CR1 family in avian genomes excluding those of chickens. Hence, this work will be a cornerstone for elucidating the phylogenetic relationships in penguin evolution using the retroposon method.

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.

Primate phylogeny: molecular evidence from retroposons

In these postgenomic times where aspects of functional genetics and character evolution form a focal point of human-mouse comparative research, primate phylogenetic research gained a widespread interest in evolutionary biology. Nevertheless, it also remains a controversial subject. Despite the surge in available primate sequences and corresponding phylogenetic interpretations, primate origins as well as several branching events in primate divergence are far from settled. The analysis of SINEs - short interspersed elements - as molecular cladistic markers represents a particularly interesting complement to sequence data. The following summarizes and discusses potential applications of this new approach in molecular phylogeny and outlines main results obtained with SINEs in the context of primate evolutionary research. Another molecular cladistic marker linking the tarsier with the anthropoid primates is also presented. This eliminates any possibility of confounding phylogenetic interpretations through lineage sorting phenomena and makes use of a new point of view in settling the phylogenetic relationships of the primate infraorders.

Phylogeny and Divergence-Date Estimates of Rapid Radiations in Muroid Rodents Based on Multiple Nuclear Genes

The muroid rodents are the largest superfamily of mammals, containing nearly one third of all mammal species. We report on a phylogenetic study comprising 53 genera sequenced for four nuclear genes, GHR, BRCA1, RAG1, and c-myc, totaling up to 6400 nucleotides. Most relationships among the subfamilies are resolved. All four genes yield nearly identical phylogenies, differing only in five key regions, four of which may represent particularly rapid radiations. Support is very strong for a fundamental division of the mole rats of the subfamilies Spalacinae and Rhizomyinae from all other muroids. Among the other "core" muroids, a rapid radiation led to at least four distinct lineages: Asian Calomyscus, an African clade of at least four endemic subfamilies, including the diverse Nesomyinae of Madagascar, a hamster clade with maximum diversity in the New World, and an Old World clade including gerbils and the diverse Old World mice and rats (Murinae). The Deomyinae, recently removed from the Murinae, is well supported as the sister group to the gerbils (Gerbillinae). Four key regions appear to represent rapid radiations and, despite a large amount of sequence data, remain poorly resolved: the base of the "core" muroids, among the five cricetid (hamster) subfamilies, within a large clade of Sigmodontinae endemic to South America, and among major geographic lineages of Old World Murinae. Because of the detailed taxon sampling within the Murinae, we are able to refine the fossil calibration of a rate-smoothed molecular clock and apply this clock to date key events in muroid evolution. We calculate rate differences among the gene regions and relate those differences to relative contribution of each gene to the support for various nodes. The among-gene variance in support is greatest for the shortest branches. We present a revised classification for this largest but most unsettled mammalian superfamily.

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.

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.

B1 and related SINEs in mammalian genomes

Although B1 and Alu were the first discovered Short Interspersed Elements (SINEs), the studies of these genomic repeats were mostly limited to mice and humans and little data on their presence in other animals were available. Here we report the presence of these SINEs in a wide range of rodents (in all 15 tested families) as well as primates and tree-shrews and their absence in other mammals. Distribution pattern of these SINEs in mammals supports close relationship between rodents and primates as well as tree-shrews. Sequence analysis of these elements, apparently descending from cellular 7SL RNA indicates their rearrangements such as dimerization (Alu), quasi-dimerization (B1), acquiring a tRNA-related unit (B1-dID), extended deletions, etc., preceding their active expansion in the genomes. The revealed common pattern of microenvironment of some rearrangement hot spots in SINEs (internal duplications and deletions) suggests involvement of short direct repeats in the mechanism of such rearrangements. This hypothesis allows us to explain short rearrangements in these and other short retroposons.

Growth and decline of introns

To gauge the processes that might direct the length of introns, I studied the balance of indels (insertions or deletions, determined using Alu and LINE1 retroposon repeats) and the density of these repeats in the introns of the human genome. The indel balance is biased in favour of deletions and correlated with the divergence of repeats. At fixed repeat divergence, the indel bias correlated with the intron size: the shorter the intron, the more deletions were favoured over insertions. This correlation with the intron size was stronger than with the gene-wide or isochore-wide parameters. The density of repeats (the number of repeats in a unit of intron length) correlated positively with the intron size. Thus, quite different mechanisms, the indel bias and the integration and/or persistence of retroposons, act in the same direction in regards to intron size, which suggests selection for the size of individual introns.

Species identity: conserved inverted LINE repeat clusters (ILRC) in the vertebrate genome as indicators of population boundaries

The avian LINE CR1 generates multiply-superimposed insertions, resulting in apparently fortuitous inverted LINE repeat clusters (ILRCs). These loci display size micro-heterogeneity within populations, with few or no presence/absence polymorphisms, and yet only very closely related species share loci. Whilst the CR1 sequences that the ILRC loci are derived from are not species-specific, the loci themselves appear to be species-characteristic if not totally species specific. The mammalian LINE L1 is shown to act similarly to CR1 and also forms ILRCs. It is proposed that whilst the formation of these loci may be from a near-random process of super-insertion certain of them are in some way functional, explaining their conservation and rapid spread to population boundaries, whilst non-functional or inactive loci are quickly lost. ILRCs appear to decay from the element as formed by the accumulation of point mutations. ILRCs appear to an unusual example of non-polymorphic sequences being younger than polymorphic sequences with no obvious selective reason.

Comparative genomic aspects of rat, mouse and human MHC class I gene regions

In this review a particular aspect of the genomic structure of the major histocompatibility complex (MHC), the organization of MHC class I regions, will be discussed for the rat in comparison to mouse and human.

Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

R1 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of arthropods. One aspect of R1 evolution that has been difficult to explain is the presence of divergent lineages of R1 in the rDNA loci of the same species. Multiple lineages should compete for a limited number of insertion sites, in addition to being subject to the concerted evolution processes homogenizing the rRNA genes. The presence of multiple lineages suggests either the ability of the elements to overcome these factors and diverge within rDNA loci, or the introduction of new lineages by horizontal transmission. To address this issue, we attempted to characterize the complete set of R1 elements in the rDNA locus from five Drosophila species groups (melanogaster, obscura, testacea, quinaria, and repleta). Two major R1 lineages, A and B, that diverged about 100 MYA were found to exist in Drosophila. Elements of the A lineage were found in all 35 Drosophila species tested, while elements of the B lineage were found in only 11 species from three species groups. Phylogenetic analysis of the R1 elements, supported by comparison of their rates of nucleotide sequence substitution, revealed that both the A and the B lineages have been maintained by vertical descent. The B lineage was less stable and has undergone numerous, independent elimination events, while the A lineage has diverged into three sublineages, which were, in turn, differentially stable. We conclude that while the differential retention of multiple lineages greatly complicates its phylogenetic history, the available R1 data continue to be consistent with the strict vertical descent of these elements.

[Strain-, sex- and species-related differences of acetohexamide reductase and 20 beta-hydroxysteroid dehydrogenase activities in liver microsomes of experimental animals]

We examined physiological and genetic factors affecting acetohexamide reductase (AHR) and 20 beta-hydroxysteroid dehydrogenase (20 beta-HSD) activities in liver microsomes of experimental animals. Pronounced strain-related differences were found in both activities of AHR and 20 beta-HSD present in liver microsomes of male rats. Among rat strains tested in this study, even though a Wistar-Imamichi (WIM) rat strain was taken to lack AHR activity, it exhibited a significant 20 beta-HSD activity. These findings appeared to be in conflict with our conclusion reported so far, which AHR and 20 beta-HSD present in liver microsomes of male rats are identical enzymes. Thus the reason for this discrepancy was discussed. Furthermore, AHR and 20 beta-HSD activities were little or not observed in liver microsomes of female rats or male experimental animals other than the rat, indicating the existence of sex- and species-related differences in these two enzyme activities.