Tracking marsupial evolution using archaic genomic retroposon insertions

Unearthing the secrets of Australia's most enigmatic and cryptic mammal, the marsupial mole

The marsupial moles are arguably Australia's most enigmatic marsupials. Almost indistinguishable from placental (eutherian) moles, they provide a striking example of convergent evolution. Exploring the genome of the southern marsupial mole, we provide insights into its unusual biology. We show definitively by retrophylogenomic analysis that marsupial moles are most closely related to bandicoots and bilbies (order Peramelemorphia). We find evidence of a marked decline in marsupial mole effective population size, most likely preceding the arrival of humans in regions near its range, and potentially corresponding to periods of climatic change. Our analysis of loss of eye function-an adaptation to subterranean life-reveals a structured order of loss of gene function associated first with the lens, then cone, and finally rod cells. Last, we identify genetic changes suggestive of adaptation to an oxygen-poor environment and of its evolution of partially descended testes.

Expanding the evo-devo toolkit: generation of 3D mammary tissue from diverse mammals

The varying pathways of mammary gland development across species and evolutionary history are underexplored, largely due to a lack of model systems. Recent progress in organoid technology holds the promise of enabling in-depth studies of the developmental adaptations that have occurred throughout the evolution of different species, fostering beneficial phenotypes. The practical application of this technology for mammary glands has been mostly confined to rodents and humans. In the current study, we have successfully created next-generation 3D mammary gland organoids from eight eutherian mammals and the first branched organoid of a marsupial mammary gland. Using mammary organoids, we identified a role for ROCK protein in regulating branching morphogenesis, a role that manifests differently in organoids from different mammals. This finding demonstrates the utility of the 3D organoid model for understanding the evolution and adaptations of signaling pathways. These achievements highlight the potential for organoid models to expand our understanding of mammary gland biology and evolution, and their potential utility in studies of lactation or breast cancer.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Gross anatomy of the gluteal and posterior thigh muscles in koalas based on their innervations

Morphological and functional comparison of convergently-evolved traits in marsupials and eutherians is an important aspect of studying adaptive divergence in mammals. However, the anatomy of marsupials has been particularly difficult to evaluate for multiple reasons. First, previous studies on marsupial anatomy are often uniformly old and non-exhaustive. Second, muscle identification was historically based on muscle attachment sites, but attachment sites have since been declared insufficient for muscle identification due to extensive interspecific variation. For example, different names have been used for muscles that are now thought to be equivalent among several different species, which causes confusion. Therefore, descriptions of marsupial muscles have been inconsistent among previous studies, and their anatomical knowledge itself needs updating. In this study, the koala was selected as the representative marsupial, in part because koala locomotion may comprise primate (eutherian)-like and marsupial-like mechanics, making it an interesting phylogenetic group for studying adaptive divergence in mammals. Gross dissection of the lower limb muscles (the gluteal and the posterior thigh regions) was performed to permit precise muscle identification. We first resolved discrepancies among previous studies by identifying muscles according to their innervation; this recent, more reliable technique is based on the ontogenetic origin of the muscle, and it allows for comparison with other taxa (i.e., eutherians). We compared our findings with those of other marsupials and arboreal primates and identified traits common to both arboreal primates and marsupials as well as muscle morphological features unique to koalas.

Reunion of Australasian Possums by Shared SINE Insertions

Although first posited to be of a single origin, the two superfamilies of phalangeriform marsupial possums (Phalangeroidea: brushtail possums and cuscuses and Petauroidea: possums and gliders) have long been considered, based on multiple sequencing studies, to have evolved from two separate origins. However, previous data from these sequence analyses suggested a variety of conflicting trees. Therefore, we reinvestigated these relationships by screening $\sim$200,000 orthologous short interspersed element (SINE) loci across the newly available whole-genome sequences of phalangeriform species and their relatives. Compared to sequence data, SINE presence/absence patterns are evolutionarily almost neutral molecular markers of the phylogenetic history of species. Their random and highly complex genomic insertion ensures their virtually homoplasy-free nature and enables one to compare hundreds of shared unique orthologous events to determine the true species tree. Here, we identify 106 highly reliable phylogenetic SINE markers whose presence/absence patterns within multiple Australasian possum genomes unexpectedly provide the first significant evidence for the reunification of Australasian possums into one monophyletic group. Together, our findings indicate that nucleotide homoplasy and ancestral incomplete lineage sorting have most likely driven the conflicting signal distributions seen in previous sequence-based studies. [Ancestral incomplete lineage sorting; possum genomes; possum monophyly; retrophylogenomics; SINE presence/absence.].

SINEs as Credible Signs to Prove Common Ancestry in the Tree of Life: A Brief Review of Pioneering Case Studies in Retroposon Systematics

Currently, the insertions of SINEs (and other retrotransposed elements) are regarded as one of the most reliable synapomorphies in molecular systematics. The methodological mainstream of molecular systematics is the calculation of nucleotide (or amino acid) sequence divergences under a suitable substitution model. In contrast, SINE insertion analysis does not require any complex model because SINE insertions are unidirectional and irreversible. This straightforward methodology was named the “SINE method,” which resolved various taxonomic issues that could not be settled by sequence comparison alone. The SINE method has challenged several traditional hypotheses proposed based on the fossil record and anatomy, prompting constructive discussions in the Evo/Devo era. Here, we review our pioneering SINE studies on salmon, cichlids, cetaceans, Afrotherian mammals, and birds. We emphasize the power of the SINE method in detecting incomplete lineage sorting by tracing the genealogy of specific genomic loci with minimal noise. Finally, in the context of the whole-genome era, we discuss how the SINE method can be applied to further our understanding of the tree of life.

Incomplete lineage sorting and phenotypic evolution in marsupials

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

The Reaction Specificity of Mammalian ALOX15 Orthologs is Changed During Late Primate Evolution and These Alterations Might Offer Evolutionary Advantages for Hominidae

Arachidonic acid lipoxygenases (ALOXs) have been implicated in the immune response of mammals. The reaction specificity of these enzymes is decisive for their biological functions and ALOX classification is based on this enzyme property. Comparing the amino acid sequences and the functional properties of selected mammalian ALOX15 orthologs we previously hypothesized that the reaction specificity of these enzymes can be predicted based on their amino acid sequences (Triad Concept) and that mammals, which are ranked in evolution below gibbons, express arachidonic acid 12-lipoxygenating ALOX15 orthologs. In contrast, Hominidae involving the great apes and humans possess 15-lipoxygenating enzymes (Evolutionary Hypothesis). These two hypotheses were based on sequence data of some 60 mammalian ALOX15 orthologs and about half of them were functionally characterized. Here, we compared the ALOX15 sequences of 152 mammals representing all major mammalian subclades expressed 44 novel ALOX15 orthologs and performed extensive mutagenesis studies of their triad determinants. We found that ALOX15 genes are absent in extant Prototheria but that corresponding enzymes frequently occur in Metatheria and Eutheria. More than 90% of them catalyze arachidonic acid 12-lipoxygenation and the Triad Concept is applicable to all of them. Mammals ranked in evolution above gibbons express arachidonic acid 15-lipoxygenating ALOX15 orthologs but enzymes with similar specificity are only present in less than 5% of mammals ranked below gibbons. This data suggests that ALOX15 orthologs have been introduced during Prototheria-Metatheria transition and put the Triad Concept and the Evolutionary Hypothesis on a much broader and more reliable experimental basis.

A chromosome-level genome of Antechinus flavipes provides a reference for an Australian marsupial genus with male death after mating

The 15 species of small carnivorous marsupials that comprise the genus Antechinus exhibit semelparity, a rare life-history strategy in mammals where synchronized death occurs after one breeding season. Antechinus males, but not females, age rapidly (demonstrate organismal senescence) during the breeding season and show promise as new animal models of ageing. Some antechinus species are also threatened or endangered. Here, we report a chromosome-level genome of a male yellow-footed antechinus Antechinus flavipes. The genome assembly has a total length of 3.2 Gb with a contig N50 of 51.8 Mb and a scaffold N50 of 636.7 Mb. We anchored and oriented 99.7% of the assembly on seven pseudochromosomes and found that repetitive DNA sequences occupy 51.8% of the genome. Draft genome assemblies of three related species in the subfamily Phascogalinae, two additional antechinus species (Antechinus argentus and A. arktos) and the iteroparous sister species Murexia melanurus, were also generated. Preliminary demographic analysis supports the hypothesis that climate change during the Pleistocene isolated species in Phascogalinae and shaped their population size. A transcriptomic profile across the A. flavipes breeding season allowed us to identify genes associated with aspects of the male die-off. The chromosome-level A. flavipes genome provides a steppingstone to understanding an enigmatic life-history strategy and a resource to assist the conservation of antechinuses.

The ecology and evolution of the monito del monte, a relict species from the southern South America temperate forests

The arboreal marsupial monito del monte (genus Dromiciops, with two recognized species) is a paradigmatic mammal. It is the sole living representative of the order Microbiotheria, the ancestor lineage of Australian marsupials. Also, this marsupial is the unique frugivorous mammal in the temperate rainforest, being the main seed disperser of several endemic plants of this ecosystem, thus acting as keystone species. Dromiciops is also one of the few hibernating mammals in South America, spending half of the year in a physiological dormancy where metabolism is reduced to 10% of normal levels. This capacity to reduce energy expenditure in winter contrasts with the enormous energy turnover rate they experience in spring and summer. The unique life history strategies of this living Microbiotheria, characterized by an alternation of life in the slow and fast lanes, putatively represent ancestral traits that permitted these cold-adapted mammals to survive in this environment. Here, we describe the ecological role of this emblematic marsupial, summarizing the ecophysiology of hibernation and sociality, updated phylogeographic relationships, reproductive cycle, trophic relationships, mutualisms, conservation, and threats. This marsupial shows high densities, despite presenting slow reproductive rates, a paradox explained by the unique characteristics of its three-dimensional habitat. We finally suggest immediate actions to protect these species that may be threatened in the near future due to habitat destruction and climate change.

The Coevolution of Placentation and Cancer

Analogies between placentation, in particular the behavior of trophoblast cells, and cancer have been noted since the beginning of the twentieth century. To what degree these can be explained as a consequence of the evolution of placentation has been unclear. In this review, we conclude that many similarities between trophoblast and cancer cells are shared with other, phylogenetically older processes than placentation. The best candidates for cancer hallmarks that can be explained by the evolution of eutherian placenta are mechanisms of immune evasion. Another dimension of the maternal accommodation of the placenta with an impact on cancer malignancy is the evolution of endometrial invasibility. Species with lower degrees of placental invasion tend to have lower vulnerability to cancer malignancy. We finally identify several areas in which one could expect to see coevolutionary changes in placental and cancer biology but that, to our knowledge, have not been explored. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Dispersal, not vicariance, explains the biogeographic origin of iguanas on Madagascar

Lizards of the clade Iguanidae (sensu lato) are primarily a New World group. Thus, the remarkable presence of an endemic lineage of iguanas (family Opluridae) on the isolated Indian Ocean island of Madagascar has long been considered a biogeographic anomaly. Previous work attributed this disjunct extant distribution to: (1) vicariance at about 140-165 Ma, caused by the breakup of Gondwana and the separation of South America, Africa, and Madagascar (with subsequent extinction of iguanas on Africa, and potentially other Gondwanan landmasses), (2) vicariance at about 80-90 Ma, caused by the sundering of hypothesized land-bridge connections between South America, Antarctica, India, and Madagascar, or (3) long-distance overwater dispersal from South America to Madagascar. Each hypothesis has been supported with molecular divergence dating analyses, and thus the biogeographic origin of the Opluridae is not yet well resolved. Here we utilize genetic sequences of ultraconserved elements for all Iguania families and the majority of Iguanidae (s.l.) genera, and morphological data for extant and fossil taxa (used for divergence dating analyses), to produce the most comprehensive dataset applied to date to test these origin hypotheses. We find strong support for a sister relationship between the Opluridae (Madagascar) and Leiosauridae (South America). Divergence of the Opluridae from Leiosauridae is dated to between the late Cretaceous and mid-Paleogene, at a time when Madagascar was already an island and was isolated from all other Gondwanan landmasses. Consequently, our results support a hypothesis of long-distance overwater dispersal of the Opluridae lineage, either directly from South America to Madagascar or potentially via Antarctica or Africa, leading to this radiation of iguanas in the Indian Ocean.

Theoretical and practical considerations when using retroelement insertions to estimate species trees in the anomaly zone

A potential shortcoming of concatenation methods for species tree estimation is their failure to account for incomplete lineage sorting. Coalescent methods address this problem but make various assumptions that, if violated, can result in worse performance than concatenation. Given the challenges of analyzing DNA sequences with both concatenation and coalescent methods, retroelement insertions (RIs) have emerged as powerful phylogenomic markers for species tree estimation. Here, we show that two recently proposed quartet-based methods, SDPquartets and ASTRAL_BP, are statistically consistent estimators of the unrooted species tree topology under the coalescent when RIs follow a neutral infinite-sites model of mutation and the expected number of new RIs per generation is constant across the species tree. The accuracy of these (and other) methods for inferring species trees from RIs has yet to be assessed on simulated data sets, where the true species tree topology is known. Therefore, we evaluated eight methods given RIs simulated from four model species trees, all of which have short branches and at least three of which are in the anomaly zone. In our simulation study, ASTRAL_BP and SDPquartets always recovered the correct species tree topology when given a sufficiently large number of RIs, as predicted. A distance-based method (ASTRID_BP) and Dollo parsimony also performed well in recovering the species tree topology. In contrast, unordered, polymorphism, and Camin-Sokal parsimony typically fail to recover the correct species tree topology in anomaly zone situations with more than four ingroup taxa. Of the methods studied, only ASTRAL_BP automatically estimates internal branch lengths (in coalescent units) and support values (i.e. local posterior probabilities). We examined the accuracy of branch length estimation, finding that estimated lengths were accurate for short branches but upwardly biased otherwise. This led us to derive the maximum likelihood (branch length) estimate for when RIs are given as input instead of binary gene trees; this corrected formula produced accurate estimates of branch lengths in our simulation study, provided that a sufficiently large number of RIs were given as input. Lastly, we evaluated the impact of data quantity on species tree estimation by repeating the above experiments with input sizes varying from 100 to 100 000 parsimony-informative RIs. We found that, when given just 1 000 parsimony-informative RIs as input, ASTRAL_BP successfully reconstructed major clades (i.e clades separated by branches > 0.3 CUs) with high support and identified rapid radiations (i.e. shorter connected branches), although not their precise branching order. The local posterior probability was effective for controlling false positive branches in these scenarios.

A Chromosome-Level Genome of the Agile Gracile Mouse Opossum (Gracilinanus agilis)

There are more than 100 species of American didelphid marsupials (opossums and mouse opossums). Limited genomic resources for didelphids exists, with only two publicly available genome assemblies compared with dozens in the case of their Australasian counterparts. This discrepancy impedes evolutionary and ecological research. To address this gap, we assembled a high-quality chromosome-level genome of the agile gracile mouse opossum (Gracilinanus agilis) using a combination of stLFR sequencing, polishing with mate-pair data, and anchoring onto pseudochromosomes using Hi-C. This species employs a rare life-history strategy, semelparity, and all G. agilis males and most females die at the end of their first breeding season after succumbing to stress and exhaustion. The 3.7-Gb chromosome-level assembly, with 92.6% anchored onto pseudochromosomes, has a scaffold N50 of 683.5 Mb and a contig N50 of 56.9 kb. The genome assembly shows high completeness, with a mammalian BUSCO score of 88.1%. Around 49.7% of the genome contains repetitive elements. Gene annotation yielded 24,425 genes, of which 83.9% were functionally annotated. The G. agilis genome is an important resource for future studies of marsupial biology, evolution, and conservation.

Heterothermy as the Norm, Homeothermy as the Exception: Variable Torpor Patterns in the South American Marsupial Monito del Monte (Dromiciops gliroides)

Hibernation (i.e., multiday torpor) is considered an adaptive strategy of mammals to face seasonal environmental challenges such as food, cold, and/or water shortage. It has been considered functionally different from daily torpor, a physiological strategy to cope with unpredictable environments. However, recent studies have shown large variability in patterns of hibernation and daily torpor ("heterothermic responses"), especially in species from tropical and subtropical regions. The arboreal marsupial "monito del monte" (Dromiciops gliroides) is the last living representative of the order Microbiotheria and is known to express both short torpor episodes and also multiday torpor depending on environmental conditions. However, only limited laboratory experiments have documented these patterns in D. gliroides. Here, we combined laboratory and field experiments to characterize the heterothermic responses in this marsupial at extreme temperatures. We used intraperitoneal data loggers and simultaneous measurement of ambient and body temperatures (T A and T B, respectively) for analyzing variations in the thermal differential, in active and torpid animals. We also explored how this differential was affected by environmental variables (T A, natural photoperiod changes, food availability, and body mass changes), using mixed-effects generalized linear models. Our results suggest that: (1) individuals express short bouts of torpor, independently of T A and even during the reproductive period; (2) seasonal torpor also occurs in D. gliroides, with a maximum bout duration of 5 days and a mean defended T B of 3.6 ± 0.9°C (one individual controlled T B at 0.09°C, at sub-freezing T A); (3) the best model explaining torpor occurrence (Akaike information criteria weight = 0.59) discarded all predictor variables except for photoperiod and a photoperiod by food interaction. Altogether, these results confirm that this marsupial expresses a dynamic form of torpor that progresses from short torpor to hibernation as daylength shortens. These data add to a growing body of evidence characterizing tropical and sub-tropical heterothermy as a form of opportunistic torpor, expressed as daily or seasonal torpor depending on environmental conditions.

Mobility of the forearm skeleton in the Asiatic black (Ursus thibetanus), brown (U. arctos) and polar (U. maritimus) bears

In several primates and carnivores, pronation/supination angles of the forearm skeleton were examined, and it is thought that a larger angle is useful to acquire dexterous behaviors in feeding and/or life style, including climbing. In this study, the pronation/supination angles in Asiatic black, brown and polar bears were nondestructively examined. These specimens were classified as adult or non-adult. Three or four carcasses of each group of Asiatic black and brown bears were used for CT analysis, whereas only one adult polar bear was used. The forearms were positioned within the gantry of a CT scanner in both maximally supinated and pronated states. Extracted cross-sectional CT images of two positions were superimposed by overlapping the outlines of each ulna. The centroids of the radii were detected, and then the centroid of each radius and the midpoint of a line which connects between both ends of the surface of each radius facing the ulna, were connected by lines to measure the angle of rotation as an index of pronation/supination. In adult brown and polar bears, the angles were smaller as compared with the other groups (Asiatic black and non-adult brown bears). Asiatic black and non-adult brown bears can climb trees, whereas adult brown bears and polar bears cannot. This suggests that the pronation/supination angle is related to arboreal activity in Ursidae.

Unraveling Genome Evolution Throughout Visual Analysis: The XCout Portal

Due to major breakthroughs in sequencing technologies throughout the last decades, the time and cost per sequencing experiment have reduced drastically, overcoming the data generation barrier during the early genomic era. Such a shift has encouraged the scientific community to develop new computational methods that are able to compare large genomic sequences, thus enabling large-scale studies of genome evolution. The field of comparative genomics has proven itself invaluable for studying the evolutionary mechanisms and the forces driving genome evolution. In this line, a full genome comparison study between 2 species requires a quadratic number of comparisons in terms of the number of sequences (around 400 chromosome comparisons in the case of mammalian genomes); however, when studying conserved syntenies or evolutionary rearrangements, many sequence comparisons can be skipped for not all will contain significant signals. Subsequently, the scientific community has developed fast heuristics to perform multiple pairwise comparisons between large sequences to determine whether significant sets of conserved similarities exist. The data generation problem is no longer an issue, yet the limitations have shifted toward the analysis of such massive data. Therefore, we present XCout, a Web-based visual analytics application for multiple genome comparisons designed to improve the analysis of large-scale evolutionary studies using novel techniques in Web visualization. XCout enables to work on hundreds of comparisons at once, thus reducing the time of the analysis by identifying significant signals between chromosomes across multiple species. Among others, XCout introduces several techniques to aid in the analysis of large-scale genome rearrangements, particularly (1) an interactive heatmap interface to display comparisons using automatic color scales based on similarity thresholds to ease detection at first sight, (2) an overlay system to detect individual signal contributions between chromosomes, (3) a tracking tool to trace conserved blocks across different species to perform evolutionary studies, and (4) a search engine to search annotations throughout different species.

100-My history of bornavirus infections hidden in vertebrate genomes

Although viruses have threatened our ancestors for millions of years, prehistoric epidemics of viruses are largely unknown. Endogenous bornavirus-like elements (EBLs) are ancient bornavirus sequences derived from the viral messenger RNAs that were reverse transcribed and inserted into animal genomes, most likely by retrotransposons. These elements can be used as molecular fossil records to trace past bornaviral infections. In this study, we systematically identified EBLs in vertebrate genomes and revealed the history of bornavirus infections over nearly 100 My. We confirmed that ancient bornaviral infections have occurred in diverse vertebrate lineages, especially in primate ancestors. Phylogenetic analyses indicated that primate ancestors were infected with various bornaviral lineages during evolution. EBLs in primate genomes formed clades according to their integration ages, suggesting that bornavirus lineages infected with primate ancestors had changed chronologically. However, some bornaviral lineages may have coexisted with primate ancestors and underwent repeated endogenizations for tens of millions of years. Moreover, a bornaviral lineage that coexisted with primate ancestors also endogenized in the genomes of some ancestral bats. The habitats of these bat ancestors have been reported to overlap with the migration route of primate ancestors. These results suggest that long-term virus-host coexistence expanded the geographic distributions of the bornaviral lineage along with primate migration and may have spread their infections to these bat ancestors. Our findings provide insight into the history of bornavirus infections over geological timescales that cannot be deduced from research using extant viruses alone, thus broadening our perspective on virus-host coevolution.

Molecular analysis suggests that Namibian cheetahs (Acinonyx jubatus) are definitive hosts of a so far undescribed Besnoitia species

Background

Besnoitia darlingi, B. neotomofelis and B. oryctofelisi are closely related coccidian parasites with felids as definitive hosts. These parasites use a variety of animal species as intermediate hosts. North American opossums (Didelphis virginiana), North American southern plains woodrats (Neotoma micropus) and South American domestic rabbits (Oryctolagus cuniculus) are intermediate hosts of B. darlingi, B. neotomofelis and B. oryctofelisi, respectively. Based on conserved regions in the internal transcribed spacer-1 (ITS1) sequence of the ribosomal DNA (rDNA), a real-time PCR for a sensitive detection of these Besnoitia spp. in tissues of intermediate hosts and faeces of definitive hosts has recently been established. Available sequence data suggest that species such as B. akodoni and B. jellisoni are also covered by this real-time PCR. It has been hypothesised that additional Besnoitia spp. exist worldwide that are closely related to B. darlingi or B. darlingi-like parasites (B. neotomofelis, B. oryctofelisi, B. akodoni or B. jellisoni). Also related, but not as closely, is B. besnoiti, the cause of bovine besnoitiosis.

Methods

Faecal samples from two free-ranging cheetahs (Acinonyx jubatus) from Namibia that had previously tested positive for coccidian parasites by coproscopy were used for this study. A conventional PCR verified the presence of coccidian parasite DNA. To clarify the identity of these coccidia, the faecal DNA samples were further characterised by species-specific PCRs and Sanger sequencing.

Results

One of the samples tested positive for B. darlingi or B. darlingi-like parasites by real-time PCR, while no other coccidian parasites, including Toxoplasma gondii, Hammondia hammondi, H. heydorni, B. besnoiti and Neospora caninum, were detected in the two samples. The rDNA of the B. darlingi-like parasite was amplified and partially sequenced. Comparison with existing sequences in GenBank revealed a close relationship to other Besnoitia spp., but also showed clear divergences.

Conclusions

Our results suggest that a so far unknown Besnoitia species exists in Namibian wildlife, which is closely related to B. darlingi, B. neotomofelis, B. oryctofelisi, B. akodoni or B. jellisoni. The cheetah appears to be the definitive host of this newly discovered parasite, while prey species of the cheetah may act as intermediate hosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04697-3.

The complete mitochondrial genome of the vulnerable Australian crest-tailed mulgara (Dasycercus cristicauda)

In this announcement, we report the complete mitogenome of the vulnerable Crest-tailed Mulgara (Dasycercus cristicauda) (Krefft, 1867). The mitogenome was 17,085 bp in length and contained 13 protein-coding genes, two rRNA genes, 22 tRNAs and a 1583 bp variable control region (D-loop). The features of the D. cristicauda mitogenome are consistent with other vertebrate mitogenomes but, in contrast to other marsupials, appears to contain a functional tRNA-Lysine with a UUU anticodon. Phylogenetic analysis of available entire mitogenomes reveals it forms a cluster with other marsupials in the Dasyuromorphia order within the Australidelphian clade, being most closely related to the Northern Quoll and the Tasmanian Devil.

Intra-skeletal vascular density in a bipedal hopping macropod with implications for analyses of rib histology

Human ribs are thought to be less affected by mechanical strain at the microscopic level than limb bones, implying that rib remodelling better reflects bone physiological homeostasis. Here, we test the hypothesis that rib tissue will be well vascularized and thus enhance susceptibility to metabolic influence. An intra-skeletal comparison of bone vascular canal density was conducted using a macropod animal model adapted to bipedal habitual hopping. The right humerus, ulna, radius, femur, tibia, fibula, a mid-thoracic and upper-thoracic rib of an eastern grey kangaroo (Macropus giganteus) were sectioned at the midshaft, from which histological sections were prepared. Bone vascularity from a maximum of 12 mm² of sub-periosteal parallel-fibred and lamellar bone was recorded, resulting in a total of 2047 counted vessels. Vascular canal density data were corrected by cortical width, maximum length, and midshaft circumference robusticity indices computed for each bone. The fibula consistently had the highest vascular canal density, even when corrected for maximum length, cortical width and midshaft circumference robusticities. This was followed by the mid- and upper-thoracic ribs. Vascularity differences between bones were relatively consistent whether vascular canal density was controlled for by cortical width or midshaft circumference robusticities. Vascular canal density and robusticity indices were also positively and negatively correlated (p < 0.05). Results confirm that the ribs are well vascularized, which facilitates bone metabolic processes such as remodelling, but the fibula also appears to be a well vascularized bone. Future research investigating human bone metabolism will benefit from examining thoracic rib or fibula samples.

Endometrial receptivity and embryo implantation in carnivores-commonalities and differences with other mammalian species

Endometrial receptivity and embryo implantation processes are a major point of pregnancy failure in many mammalian species, including humans. Although reproductive biology in many carnivore species remains enigmatic, the few that have been studied so far are invaluable comparative models. The goals of this review are to (1) summarize current data on the mechanisms involved in uterine receptivity and embryo implantation in carnivores, including commonalities and differences with other mammalian species and (2) identify research priorities to better understand a key phenomenon in a critical group of mammals. Besides unique reproductive traits in some carnivores (induced vs. spontaneous ovulation in cats, ovulation at the germinal vesicle stage in dogs), preimplantation embryo development is comparable with other orders. However, the timing of implantation varies, especially in species having an embryonic diapause. Mechanisms involved in endometrial receptivity and decidualization still remain to be fully understood, but specific markers have already been identified. Importantly, the use of endogenous hormones to control the ovarian activity may impact endometrial receptivity and subsequent embryo implantation. Next, research efforts should take advantage of advanced technologies to further study embryo implantation in carnivores and to provide more relevant models to reproductive medicine or for the conservation of rare and endangered species.

Boehmiella wilsoni (Nematoda, Heligmosomoidea, Boehmiellidae fam. nov.), found in Amazonian rodents

The genus Boehmiella was initially described as a member of the family Trichostrongylidae. Subsequently, it was assigned to the subfamily Haemonchinae in the family Haemonchidae. We analyzed parasites of spiny tree-rats, Mesomys hispidus, collected in the Amazon rainforest, which were identified as B. wilsoni based on integrative taxonomy. Using morphology, morphometry, and scanning electron microscopy (SEM), we added new data to the original description of the species. We also inferred phylogenetic hypotheses for its relationships within the Trichostrongylina, based on partial nuclear 18S and 28S rRNA genes, through Maximum Likelihood and Bayesian analyses. In conclusion, B. wilsoni does not belong to the family Haemonchidae, nor is it closely related to any other trichostrongylin family, and therefore, we propose the establishment of a new family, Boehmiellidae fam. nov., to which the genus Boehmiella is allocated.

•

Helminth found in arboreal rodent of the Amazon region, with description of new host, new geographical distribution, new taxonomic data, and proposal of a new family to allocate the genus according to molecular and morphological analyses.

Variability and variation in Rhyncholestes raphanurus Osgood (Paucituberculata, Caenolestidae)

Background

Caenolestids are a group of poorly known South American marsupials with a restricted distribution in Páramo and Subpáramo environments of the Andes from Colombia and western Venezuela to Bolivia (represented by the genera Caenolestes and Lestoros), and Valdivian rainforest in southern Chile (including a separate population in Chiloé Island) and Argentina, where a single species lives: the Long-nosed shrew opossum (Rhyncholestes raphanurus). The objectives of this work were to analyze the intraspecific variability of R. raphanurus, which includes an anatomical description of the skull and dentition, describe its geographic variation, test for sexual dimorphism, and assess potential differences between continental and Island populations.

Methods

Linear Mossimann-transformed variables were used to assess sexual differences within a large population (La Picada), compare sexes within other continental populations, and in a separate analysis, compare continental from Island samples. A full model Principal Components Analysis was performed to assess differences between males and females of the continental and Island populations. A thorough description of the skull and teeth of the species and comparisons with other living Caenolestidae is presented.

Results

Rhyncholestes raphanurus presents little geographic variation, even between Island and continental populations. Similarly, we found no significant difference between sexes of this species in cranial and dental measurements. We provide a detailed description of cranial morphology and its variation, and also, the first description of the occipital bones, which haven’t been previously described for any Paucituberculata.

Conclusions

Comparative studies of continental and Chiloé Island specimens support the treatment of R. raphanurus as a single valid species, especially since morphologic and morphometric differences fall within the extremes of continental populations. The morphology of R. raphanurus clearly separates this genus from other extant Caenolestidae, and in a much greater degree than the differences found between Lestoros and Caenolestes.

Marsupial and monotreme milk-a review of its nutrient and immune properties

All mammals are characterized by the ability of females to produce milk. Marsupial (metatherian) and monotreme (prototherian) young are born in a highly altricial state and rely on their mother’s milk for the first part of their life. Here we review the role and importance of milk in marsupial and monotreme development. Milk is the primary source of sustenance for young marsupials and monotremes and its composition varies at different stages of development. We applied nutritional geometry techniques to a limited number of species with values available to analyze changes in macronutrient composition of milk at different stages. Macronutrient energy composition of marsupial milk varies between species and changes concentration during the course of lactation. As well as nourishment, marsupial and monotreme milk supplies growth and immune factors. Neonates are unable to mount a specific immune response shortly after birth and therefore rely on immunoglobulins, immunological cells and other immunologically important molecules transferred through milk. Milk is also essential to the development of the maternal-young bond and is achieved through feedback systems and odor preferences in eutherian mammals. However, we have much to learn about the role of milk in marsupial and monotreme mother-young bonding. Further research is warranted in gaining a better understanding of the role of milk as a source of nutrition, developmental factors and immunity, in a broader range of marsupial species, and monotremes.

First fossil frog from Antarctica: implications for Eocene high latitude climate conditions and Gondwanan cosmopolitanism of Australobatrachia

Cenozoic ectothermic continental tetrapods (amphibians and reptiles) have not been documented previously from Antarctica, in contrast to all other continents. Here we report a fossil ilium and an ornamented skull bone that can be attributed to the Recent, South American, anuran family Calyptocephalellidae or helmeted frogs, representing the first modern amphibian found in Antarctica. The two bone fragments were recovered in Eocene, approximately 40 million years old, sediments on Seymour Island, Antarctic Peninsula. The record of hyperossified calyptocephalellid frogs outside South America supports Gondwanan cosmopolitanism of the anuran clade Australobatrachia. Our results demonstrate that Eocene freshwater ecosystems in Antarctica provided habitats favourable for ectothermic vertebrates (with mean annual precipitation ≥900 mm, coldest month mean temperature ≥3.75 °C, and warmest month mean temperature ≥13.79 °C), at a time when there were at least ephemeral ice sheets existing on the highlands within the interior of the continent.

Regional-scale variation on Dromiciops gliroides occurrence, abundance, and activity patterns along a habitat disturbance gradient

Habitat structure may have a significant influence on the occurrence, abundance, and activity patterns of forest mammals. However, anthropogenic habitat disturbance changes habitat structure, which may alter those patterns of activity. We assessed occurrence, relative abundance, and activity patterns of Dromiciops gliroides, an arboreal marsupial endemic to the temperate rainforests of southern South America, contrasting four forest conditions at a regional scale: old-growth, second-growth, and logged forests, and abandoned exotic plantations. We conducted a camera-trap assessment in two consecutive austral summers across most of the Chilean range of D. gliroides, and compared habitat structure along a disturbance gradient. All structural features assessed differed among forest conditions. Dromiciops gliroides was present in all forest conditions, but its abundance decreased and activity got narrower as disturbance increased, being significantly lower in the exotic plantations. Activity patterns were variable among forest conditions and months, and were significantly more restricted temporally at exotic plantations. Although D. gliroides is tolerant to habitat disturbance, we show that structural alteration results in lower abundances and narrower activity patterns.

Uterine epithelial remodelling during pregnancy in the marsupial Monodelphis domestica (Didelphidae): Implications for mammalian placental evolution

Mammalian pregnancy involves remodelling of the uterine epithelium to enable placentation. In marsupials, such remodelling has probably played a key role in the transition from ancestral invasive placentation to non-invasive placentation. Identifying uterine alterations that are unique to marsupials with non-invasive placentation can thus elucidate mechanisms of marsupial placental evolution. We identified apical alterations to uterine epithelial cells prior to implantation in Monodelphis domestica, a member of the least derived living marsupial clade (Didelphidae) with invasive (endotheliochorial) placentation. We then compared these traits with those of Macropus eugenii (Macropodidae) and Trichosurus vulpecula (Phalangeridae), both with non-invasive placentation, to identify which alterations to the uterine epithelium are ancestral and which facilitate secondarily evolved non-invasive placentation. In M. domestica, remodelling of the uterine epithelium involves reduced cellular heterogeneity and development of uterodome-like cells, suggesting that similar alterations may also have occurred in the marsupial common ancestor. These alterations also overlap with those of both T. vulpecula and Ma. eugenii, suggesting that the placental shift from invasive to non-invasive placentation in marsupials involves essential, conserved characteristics, irrespective of placental mode. However, unique apical alterations of both T. vulpecula and Ma. eugenii, relative to M. domestica, imply that lineage-specific alterations underpin the evolutionary shift to non-invasive placentation in marsupials.

ILS-Aware Analysis of Low-Homoplasy Retroelement Insertions: Inference of Species Trees and Introgression Using Quartets

DNA sequence alignments have provided the majority of data for inferring phylogenetic relationships with both concatenation and coalescent methods. However, DNA sequences are susceptible to extensive homoplasy, especially for deep divergences in the Tree of Life. Retroelement insertions have emerged as a powerful alternative to sequences for deciphering evolutionary relationships because these data are nearly homoplasy-free. In addition, retroelement insertions satisfy the “no intralocus-recombination” assumption of summary coalescent methods because they are singular events and better approximate neutrality relative to DNA loci commonly sampled in phylogenomic studies. Retroelements have traditionally been analyzed with parsimony, distance, and network methods. Here, we analyze retroelement data sets for vertebrate clades (Placentalia, Laurasiatheria, Balaenopteroidea, Palaeognathae) with 2 ILS-aware methods that operate by extracting, weighting, and then assembling unrooted quartets into a species tree. The first approach constructs a species tree from retroelement bipartitions with ASTRAL, and the second method is based on split-decomposition with parsimony. We also develop a Quartet-Asymmetry test to detect hybridization using retroelements. Both ILS-aware methods recovered the same species-tree topology for each data set. The ASTRAL species trees for Laurasiatheria have consecutive short branch lengths in the anomaly zone whereas Palaeognathae is outside of this zone. For the Balaenopteroidea data set, which includes rorquals (Balaenopteridae) and gray whale (Eschrichtiidae), both ILS-aware methods resolved balaeonopterids as paraphyletic. Application of the Quartet-Asymmetry test to this data set detected 19 different quartets of species for which historical introgression may be inferred. Evidence for introgression was not detected in the other data sets.

Sex steroids influence the plasma membrane transformation in the uterus of the fat-tailed dunnart (Sminthopsis crassicaudata, Marsupialia)

The uterine epithelium undergoes remodelling to become receptive to blastocyst implantation during pregnancy in a process known as the plasma membrane transformation. There are commonalities in ultrastructural changes to the epithelium, which, in eutherian, pregnancies are controlled by maternal hormones, progesterone and oestrogens. The aim of this study was to determine the effects that sex steroids have on the uterine epithelium in the fat-tailed dunnart Sminthopsis crassicaudata, the first such study in a marsupial. Females were exposed to exogenous hormones while they were reproductively quiescent, thus not producing physiological concentrations of ovarian hormones. We found that changes to the protein E-cadherin, which forms part of the adherens junction, are controlled by progesterone and that changes to the desmoglein-2 protein, which forms part of desmosomes, are controlled by 17β-oestradiol. Exposure to a combination of progesterone and 17β-oestradiol causes changes to the microvilli on the apical surface and to the ultrastructure of the uterine epithelium. There is a decrease in lateral adhesion when the uterus is exposed to progesterone and 17β-oestradiol that mimics the hormone environment of uterine receptivity. We conclude that uterine receptivity and the plasma membrane transformation in marsupial and eutherian pregnancies are under the same endocrine control and may be an ancestral feature of therian mammals.

Physiological and anatomical investigation of the auditory brainstem in the Fat-tailed dunnart (Sminthopsis crassicaudata)

The fat-tailed dunnart (Sminthopsis crassicaudata) is a small (10–20 g) native marsupial endemic to the south west of Western Australia. Currently little is known about the auditory capabilities of the dunnart, and of marsupials in general. Consequently, this study sought to investigate several electrophysiological and anatomical properties of the dunnart auditory system. Auditory brainstem responses (ABR) were recorded to brief (5 ms) tone pips at a range of frequencies (4–47.5 kHz) and intensities to determine auditory brainstem thresholds. The dunnart ABR displayed multiple distinct peaks at all test frequencies, similar to other mammalian species. ABR showed the dunnart is most sensitive to higher frequencies increasing up to 47.5 kHz. Morphological observations (Nissl stain) revealed that the auditory structures thought to contribute to the first peaks of the ABR were all distinguishable in the dunnart. Structures identified include the dorsal and ventral subdivisions of the cochlear nucleus, including a cochlear nerve root nucleus as well as several distinct nuclei in the superior olivary complex, such as the medial nucleus of the trapezoid body, lateral superior olive and medial superior olive. This study is the first to show functional and anatomical aspects of the lower part of the auditory system in the Fat-tailed dunnart.

Composition of marsupial zona pellucida: a molecular and phylogenetic approach

The zona pellucida (ZP) is an extracellular matrix that surrounds mammalian oocytes. In eutherians it is formed from three or four proteins (ZP1, ZP2, ZP3, ZP4). In the few marsupials that have been studied, however, only three of these have been characterised (ZP2, ZP3, ZP4). Nevertheless, the composition in marsupials may be more complex, since a duplication of the ZP3 gene was recently described in one species. The aim of this work was to elucidate the ZP composition in marsupials and relate it to the evolution of the ZP gene family. For that, an in silico and molecular analysis was undertaken, focusing on two South American species (gray short-tailed opossum and common opossum) and five Australian species (brushtail possum, koala, Bennett's wallaby, Tammar wallaby and Tasmanian devil). This analysis identified the presence of ZP1 mRNA and mRNA from two or three paralogues of ZP3 in marsupials. Furthermore, evidence for ZP1 and ZP4 pseudogenes in the South American subfamily Didelphinae and for ZP3 pseudogenes in two marsupials is provided. In conclusion, two different composition models are proposed for marsupials: a model with four proteins (ZP1, ZP2 and ZP3 (two copies)) for the South American species and a model with six proteins (ZP1, ZP2, ZP3 (three copies) and ZP4) for the Australasian species.

Functional Conservation of a Developmental Switch in Mammals since the Jurassic Age

ThPOK is a "master regulator" of T lymphocyte lineage choice, whose presence or absence is sufficient to dictate development to the CD4 or CD8 lineages, respectively. Induction of ThPOK is transcriptionally regulated, via a lineage-specific silencer element, SilThPOK. Here, we take advantage of the available genome sequence data as well as site-specific gene targeting technology, to evaluate the functional conservation of ThPOK regulation across mammalian evolution, and assess the importance of motif grammar (order and orientation of TF binding sites) on SilThPOK function in vivo. We make three important points: First, the SilThPOK is present in marsupial and placental mammals, but is not found in available genome assemblies of nonmammalian vertebrates, indicating that it arose after divergence of mammals from other vertebrates. Secondly, by replacing the murine SilThPOK in situ with its marsupial equivalent using a knockin approach, we demonstrate that the marsupial SilThPOK supports correct CD4 T lymphocyte lineage-specification in mice. To our knowledge, this is the first in vivo demonstration of functional equivalency for a silencer element between marsupial and placental mammals using a definitive knockin approach. Finally, we show that alteration of the position/orientation of a highly conserved region within the murine SilThPOK is sufficient to destroy silencer activity in vivo, demonstrating that motif grammar of this "solid" synteny block is critical for silencer function. Dependence of SilThPOK function on motif grammar conserved since the mid-Jurassic age, 165 Ma, suggests that the SilThPOK operates as a silenceosome, by analogy with the previously proposed enhanceosome model.

NGS-Indel Coder: A pipeline to code indel characters in phylogenomic data with an example of its application in milkweeds (Asclepias)

Targeted genome sequencing approaches allow characterization of evolutionary relationships using a considerable number of nuclear genes and informative characters. However, most phylogenomic analyses only utilize single nucleotide polymorphisms (SNPs). Studies at the species level, especially in groups that have recently radiated, often recover low amounts of phylogenetically informative variation in coding regions, and require non-coding sequences, which are richer in indels, to resolve gene trees. Here, NGS-Indel Coder, a pipeline to detect and omit false positive indels inferred from assemblies of short read sequence data, was developed to resolve the relationships among and within major clades of the American milkweeds (Asclepias), which are the result of a rapid and recent evolutionary radiation, and whose phylogeny has been difficult to resolve. This pipeline was applied to a Hyb-Seq data set of 768 loci including targeted exons and flanking intron regions from 33 milkweed species. Robust species tree inference was improved by excluding small alignment partitions (<100 bp) that increased gene tree ambiguity and incongruence. To further investigate the robustness of indel coding, data sets that included small and large indels were explored, and species trees derived from concatenated loci versus coalescent methods based on gene trees were compared. The phylogeny of Asclepias obtained using nuclear data was well resolved, and phylogenetic information from indels improved resolution of specific nodes. The Temperate North American, Mexican Highland, and Incarnatae clades were well supported as monophyletic. Asclepias coulteri, which has been considered part of the Sonoran Desert clade based on plastome analyses, was placed as sister to all the other milkweed species studied here, rather than as a member of that clade. Two groups within the Temperate North American and Mexican clades were not resolved, and the inferred relationships strongly conflicted when comparing results based on data sets that did or did not include indel characters. This new pipeline represents a step forward in making maximal use of the information content in phylogenomic data sets.

Retrophylogenomics in rorquals indicate large ancestral population sizes and a rapid radiation

Background

Baleen whales (Mysticeti) are the largest animals on earth and their evolutionary history has been studied in detail, but some relationships still remain contentious. In particular, reconstructing the phylogenetic position of the gray whales (Eschrichtiidae) has been complicated by evolutionary processes such as gene flow and incomplete lineage sorting (ILS). Here, whole-genome sequencing data of the extant baleen whale radiation allowed us to identify transposable element (TE) insertions in order to perform phylogenomic analyses and measure germline insertion rates of TEs. Baleen whales exhibit the slowest nucleotide substitution rate among mammals, hence we additionally examined the evolutionary insertion rates of TE insertions across the genomes.

Results

In eleven whole-genome sequences representing the extant radiation of baleen whales, we identified 91,859 CHR-SINE insertions that were used to reconstruct the phylogeny with different approaches as well as perform evolutionary network analyses and a quantification of conflicting phylogenetic signals. Our results indicate that the radiation of rorquals and gray whales might not be bifurcating. The morphologically derived gray whales are placed inside the rorqual group, as the sister-species to humpback and fin whales. Detailed investigation of TE insertion rates confirm that a mutational slow down in the whale lineage is present but less pronounced for TEs than for nucleotide substitutions.

Conclusions

Whole genome sequencing based detection of TE insertions showed that the speciation processes in baleen whales represent a rapid radiation. Large genome-scale TE data sets in addition allow to understand retrotransposition rates in non-model organisms and show the potential for TE calling methods to study the evolutionary history of species.

Electronic supplementary material

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

Molecular characterisation of Interleukin-2 in two Australian marsupials (the tammar wallaby, Notamacropus eugenii, and the Tasmanian devil, Sarcophilus harrisii) facilitates the development of marsupial-specific immunological reagents

Interleukin-2 (IL-2) is an important regulator of cellular immunity in mammals. For many years, our inability to identify the expression of this cytokine in marsupials hindered our capacity to progress studies in metatherian immunology. Here, we report the use of molecular techniques to characterise the IL-2 gene for the tammar wallaby (Notamacropus eugenii) and the Tasmanian devil (Sarcophilus harrisii), which allowed the prediction of the structure and probable functions of the IL-2 proteins of these species. Deduced marsupial IL-2 proteins show considerable sequence identity to each other and to common brushtail possum (Trichosurus vulpecula) IL-2 (≥65%) but shared only 35% (tammar wallaby) and 32% (Tasmanian devil) identity with human IL-2. This difference means that reagents used to study IL-2 in human and other eutherians are unlikely to cross-react with marsupials. As a key step in furthering our ability to study cellular immune responses in marsupials and, more specifically, the susceptibility of macropodoid marsupials to intracellular pathogens, a polyclonal antibody was designed for the detection and future investigation of tammar wallaby IL-2 protein expression. The molecular data and polyclonal antibody described herein will support our development of gene probes and immunological reagents that will aid studies of infection and disease in marsupials.

Transposable Elements: Classification, Identification, and Their Use As a Tool For Comparative Genomics

Most genomes are populated by hundreds of thousands of sequences originated from mobile elements. On the one hand, these sequences present a real challenge in the process of genome analysis and annotation. On the other hand, they are very interesting biological subjects involved in many cellular processes. Here we present an overview of transposable elements biodiversity, and we discuss different approaches to transposable elements detection and analyses.

LINEs Contribute to the Origins of Middle Bodies of SINEs besides 3' Tails

Short interspersed elements (SINEs), which are nonautonomous transposable elements, require the transposition machinery of long interspersed elements (LINEs) to mobilize. SINEs are composed of two or more independently originating parts. The 5′ region is called the “head” and is derived mainly from small RNAs, and the 3′ region (“tail”) originates from the 3′ region of LINEs and is responsible for being recognized by counterpart LINE proteins. The origin of the middle “body” of SINEs is enigmatic, although significant sequence similarities among SINEs from very diverse species have been observed. Here, a systematic analysis of the similarities among SINEs and LINEs deposited on Repbase, a comprehensive database of eukaryotic repeat sequences was performed. Three primary findings are described: 1) The 5′ regions of only two clades of LINEs, RTE and Vingi, were revealed to have contributed to the middle parts of SINEs; 2) The linkage of the 5′ and 3′ parts of LINEs can be lost due to occasional tail exchange of SINEs; and 3) The previously proposed Ceph-domain was revealed to be a fusion of a CORE-domain and a 5′ part of RTE clade of LINE. Based on these findings, a hypothesis that the 5′ parts of bipartite nonautonomous LINEs, which possess only the 5′ and 3′ regions of the original LINEs, can contribute to the undefined middle part of SINEs is proposed.