Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Application of scaling to mouse spontaneous movement: Path curvature varies with speed and linear distance features isochrony

To conserve sequential behavior in relation to the topographic challenges of space, it is proposed that humans and nonhuman animals can organize behavior using different scaling principles. To deal with increases in linear distance, isochrony suggest that there is a corresponding increase in speed, whereas to deal with changes in curvature, speed is adjusted according to a power function. The present study investigates whether these principles provide a framework for describing the organization of mouse behavior in a variety of standard experimental tasks. The structure of movement was examined in ambulation during open field exploration; manipulation in a string-pulling task, in which a string is advanced hand over hand to retrieve food; and rung-walking, in which the limbs successively step from rung to rung on a horizontal ladder. Both principles were found to be conserved in the organization of mouse behavior across scales of movement. These principles provide novel measures of the temporal and geometric features of movement in the mouse and insights into how the temporal and geometric features of movement are conserved within different species.

Molecular evolution of Toll-like receptors in rodents

Toll-like receptors (TLRs), the key sensor molecules in vertebrates, trigger the innate immunity and prime the adaptive immune system. The TLR family of rodents, the largest order of mammals, typically contains 13 TLR genes. However, a clear picture of the evolution of the rodent TLR family has not yet emerged and the TLR evolutionary patterns are unclear in rodent clades. Here, we analyzed the natural variation and the evolutionary processes acting on the TLR family in rodents at both the interspecific and population levels. Our results showed that rodent TLRs were dominated by purifying selection, but a series of positively selected sites (PSSs) primarily located in the ligand-binding domain was also identified. The numbers of PSSs differed among TLRs, and nonviral-sensing TLRs had more PSSs than those in viral-sensing TLRs. Gene-conversion events were found between TLR1 and TLR6 in most rodent species. Population genetic analyses showed that TLR2, TLR8, and TLR12 were under positive selection in Rattus norvegicus and R. tanezumi, whereas positive selection also acted on TLR5 and TLR9 in the former species, as well as TLR1 and TLR7 in the latter species. Moreover, we found that the proportion of polymorphisms with potentially functional change was much lower in viral-sensing TLRs than in nonviral-sensing TLRs in both of these rat species. Our findings revealed the first thorough insight into the evolution of the rodent TLR genetic variability and provided important novel insights into the evolutionary history of TLRs over long and short timescales.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Reconciling Mouse and Human Immunology at the Altar of Genetics

Immunity to infection has been extensively studied in humans and mice bearing naturally occurring or experimentally introduced germline mutations. Mouse studies are sometimes neglected by human immunologists, on the basis that mice are not humans and the infections studied are experimental and not natural. Conversely, human studies are sometimes neglected by mouse immunologists, on the basis of the uncontrolled conditions of study and small numbers of patients. However, both sides would agree that the infectious phenotypes of patients with inborn errors of immunity often differ from those of the corresponding mutant mice. Why is that? We argue that this important question is best addressed by revisiting and reinterpreting the findings of both mouse and human studies from a genetic perspective. Greater caution is required for reverse-genetics studies than for forward-genetics studies, but genetic analysis is sufficiently strong to define the studies likely to stand the test of time. Genetically robust mouse and human studies can provide invaluable complementary insights into the mechanisms of immunity to infection common and specific to these two species.

Convergent and divergent neural circuit architectures that support acoustic communication

Vocal communication is used across extant vertebrates, is evolutionarily ancient, and been maintained, in many lineages. Here I review the neural circuit architectures that support intraspecific acoustic signaling in representative anuran, mammalian and avian species as well as two invertebrates, fruit flies and Hawaiian crickets. I focus on hindbrain motor control motifs and their ties to respiratory circuits, expression of receptors for gonadal steroids in motor, sensory, and limbic neurons as well as divergent modalities that evoke vocal responses. Hindbrain and limbic participants in acoustic communication are highly conserved, while forebrain participants have diverged between anurans and mammals, as well as songbirds and rodents. I discuss the roles of natural and sexual selection in driving speciation, as well as exaptation of circuit elements with ancestral roles in respiration, for producing sounds and driving rhythmic vocal features. Recent technical advances in whole brain fMRI across species will enable real time imaging of acoustic signaling partners, tying auditory perception to vocal production.

Hibernation is super complex: distribution, dynamics, and stability of electron transport system supercomplexes in Ictidomys tridecemlineatus

Complexes of the electron transport system can associate with each other to form supercomplexes (SCs) within mitochondrial membranes, perhaps increasing respiratory capacity or reducing reactive oxygen species production. In this study, we determined the abundance, composition, and stability of SCs in a mammalian hibernator, in which both whole-animal and mitochondria metabolism change greatly throughout winter. We isolated mitochondria from thirteen-lined ground squirrels (Ictidomys tridecemlineatus) in different hibernation states, as well as from rats (Rattus norvegicus). We extracted mitochondrial proteins using two non-ionic detergents of different strengths, and quantified SC abundance using two-dimensional gel electrophoresis and immunoblotting. Rat heart and liver had fewer SCs than ground squirrels. Within ground squirrels, SCs are dynamic, changing among hibernation states within a matter of hours. In brown adipose tissue, Complex III composition in different SCs differed between the torpid and interbout euthermic phase of a hibernation bout. In heart and liver, complex III composition changed between winter and summer. We also evaluated the stability of liver SCs using a stronger detergent and found that the stability of SCs differed: torpor SCs were more stable than the SCs of ground squirrels in other states and rats. This study is the first report of SC changes during hibernation, and the first to demonstrate their dynamics on a short timescale.

Conserved Signatures in Protein Sequences Reliably Demarcate Different Clades of Rodents/Glires Species and Consolidate Their Evolutionary Relationships

The grandorder Glires, consisting of the orders Rodentia and Lagomorpha, encompasses a significant portion of the extant mammalian species including Rat, Mouse, Squirrel, Guinea pig and Beaver. Glires species play an important role in the ecosystem and provide valuable animal models for genetic studies and animal testing. Thus, it is important to reliably determine their evolutionary relationships and identify molecular characteristics that are specific for different species groups within the Glires. In this work, we have constructed a phylogenetic tree for >30 genome sequenced Glires species based on concatenated sequences of 25 conserved proteins. In this tree, members of different orders, suborders, and families within Glires formed strongly supported clades, and their interrelationships were also generally reliably resolved. In parallel, we conducted comparative analyses on more than 1500 protein sequences from Glires species to identify highly conserved molecular markers. These markers were comprised of conserved signature indels (CSIs) in proteins, which are specific for different Rodentia/Glires clades. Of the 41 novel CSIs identified in this work, some are specific for the entire Glires, Rodentia, or Lagomorpha clades, whereas many others reliably demarcate different family/suborder level clades of Rodentia (viz. Myomorpha, Castorimorpha, Sciuromorpha, Hystricomorpha, and Muroidea). Additionally, some of the CSIs also provide information regarding the interrelationships among Rodentia subgroups. Our analysis has also identified one CSI that is commonly shared by the Glires and Scandentia species (tree shrew), however, its evolutionary significance is unclear. Several of the identifed rodents-specific CSIs are present in conserved disease-related proteins. Thus, they provide novel molecular markers for genetic and biochemical studies on the functions of these proteins.

OUP accepted manuscript

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 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim$\end{document}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.]

OUP accepted manuscript

We present Champagne, a whole-genome method for generating character matrices for phylogenomic analysis using large genomic indel events. By rigorously picking orthologous genes and locating large insertion and deletion events, Champagne delivers a character matrix that considerably reduces homoplasy compared with morphological and nucleotide-based matrices, on both established phylogenies and difficult-to-resolve nodes in the mammalian tree. Champagne provides ample evidence in the form of genomic structural variation to support incomplete lineage sorting and possible introgression in Paenungulata and human–chimp–gorilla which were previously inferred primarily through matrices composed of aligned single-nucleotide characters. Champagne also offers further evidence for Myomorpha as sister to Sciuridae and Hystricomorpha in the rodent tree. Champagne harbors distinct theoretical advantages as an automated method that produces nearly homoplasy-free character matrices on the whole-genome scale.

The existence and characteristics of rats and shrews in endemic leptospirosis areas and types of ectoparasites: a case study of West Jakarta, Indonesia

Background: This study aimed to determine the presence and species of the rats and shrews that can potentially cause leptospirosis in West Jakarta, Indonesia, and the species of ectoparasites found in them. Methods: The research was a descriptive study employing a cross-sectional approach. The study population was all species of rats and shrews in the region and the sample collection technique used was purposive sampling. The traps were installed in the homes of respondents who had suffered from leptospirosis and their closest neighbors, with a total of 521 traps. Leptospirosis data based on secondary data was obtained from West Jakarta Health Office (2016-August 2019). The technique for catching rats involved using humane live traps, while the identification of the rats and ectoparasites was done in the laboratory. Results: It was found that more rats were caught in Cengkareng Timur sub-district, Cengkareng District, with a percentage of 14.8%, while the least in Duri Kepa, Kapuk, Kedaung Kali Angke and Kedoya Utara with a percentage of 3.7%. The rats were mostly found in East Cengkareng Sub-District, with the most common type being Rattus rattus (74.1 %) and the least Suncus murinus (11.1%); more male rats were caught (66.7%) than female (33.3%). The type of ectoparasite found in the rats was fleas. Xenopsylla cheopis was the most common type, at 83.3% and more fleas were male, at 66.7%. The most common rat species was Rattus rattus. The ectoparasite most commonly found in them was the female flea Xenopsylla cheopis.Conclusions:  Rattus rattus and Xenopsylla cheopis were found in an East Cengkareng sub-district. Surveys, monitoring, and control of rats and ectoparasites are essential for the preparedness and development of an early warning system of possible diseases that they can cause.

High-altitude rodents have abundant collaterals that protect against tissue injury after cerebral, coronary and peripheral artery occlusion

Collateral number/density varies widely in brain and other tissues among strains of Mus musculus mice due to differences in genetic background. Recent studies have shown that prolonged exposure to reduced atmospheric oxygen induces additional collaterals to form, suggesting that natural selection may favor increased collaterals in populations native to high-altitude. High-altitude guinea pigs (Cavia) and deer mice (Peromyscus) were compared with lowland species of Peromyscus, Mus and Rattus (9 species/strains examined). Collateral density, diameter and other morphometrics were measured in brain where, importantly, collateral abundance reflects that in other tissues of the same individual. Guinea pigs and high-altitude deer mice had a greater density of pial collaterals than lowlanders. Consistent with this, guinea pigs and highlander mice evidenced complete and 80% protection against stroke, respectively. They also sustained significantly less ischemia in heart and lower extremities after arterial occlusion. Vessels of the circle of Willis, including the communicating collateral arteries, also exhibited unique features in the highland species. Our findings support the hypothesis that species native to high-altitude have undergone genetic selection for abundant collaterals, suggesting that besides providing protection in obstructive disease, collaterals serve a physiological function to optimize oxygen delivery to meet oxygen demand when oxygen is limiting.

Vitamin A1/A2 chromophore exchange: Its role in spectral tuning and visual plasticity

Vertebrate rod and cone photoreceptors detect light via a specialized organelle called the outer segment. This structure is packed with light-sensitive molecules known as visual pigments that consist of a G-protein-coupled, seven-transmembrane protein known as opsin, and a chromophore prosthetic group, either 11-cis retinal ('A₁') or 11-cis 3,4-didehydroretinal ('A₂'). The enzyme cyp27c1 converts A₁ into A₂ in the retinal pigment epithelium. Replacing A₁ with A₂ in a visual pigment red-shifts its spectral sensitivity and broadens its bandwidth of absorption at the expense of decreased photosensitivity and increased thermal noise. The use of vitamin A₂-based visual pigments is strongly associated with the occupation of aquatic habitats in which the ambient light is red-shifted. By modulating the A₁/A₂ ratio in the retina, an organism can dynamically tune the spectral sensitivity of the visual system to better match the predominant wavelengths of light in its environment. As many as a quarter of all vertebrate species utilize A₂, at least during a part of their life cycle or under certain environmental conditions. A₂ utilization therefore represents an important and widespread mechanism of sensory plasticity. This review provides an up-to-date account of the A₁/A₂ chromophore exchange system.

Enhancing mitogenomic phylogeny and resolving the relationships of extinct megafaunal placental mammals

Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.

Evolution of arboreality and fossoriality in squirrels and aplodontid rodents: Insights from the semicircular canals of fossil rodents

Reconstructing locomotor behaviour for fossil animals is typically done with postcranial elements. However, for species only known from cranial material, locomotor behaviour is difficult to reconstruct. The semicircular canals (SCCs) in the inner ear provide insight into an animal's locomotor agility. A relationship exists between the size of the SCCs relative to body mass and the jerkiness of an animal's locomotion. Additionally, studies have also demonstrated a relationship between SCC orthogonality and angular head velocity. Here, we employ two metrics for reconstructing locomotor agility, radius of curvature dimensions and SCC orthogonality, in a sample of twelve fossil rodents from the families Ischyromyidae, Sciuridae and Aplodontidae. The method utilizing radius of curvature dimensions provided a reconstruction of fossil rodent locomotor behaviour that is more consistent with previous studies assessing fossil rodent locomotor behaviour compared to the method based on SCC orthogonality. Previous work on ischyromyids suggests that this group displayed a variety of locomotor modes. Members of Paramyinae and Ischyromyinae have relatively smaller SCCs and are reconstructed to be relatively slower compared to members of Reithroparamyinae. Early members of the Sciuroidea clade including the sciurid Cedromus wilsoni and the aplodontid Prosciurus relictus are reconstructed to be more agile than ischyromyids, in the range of extant arboreal squirrels. This reconstruction supports previous inferences that arboreality was likely an ancestral trait for this group. Derived members of Sciuridae and Aplodontidae vary in agility scores. The fossil squirrel Protosciurus cf. rachelae is inferred from postcranial material as arboreal, which is in agreement with its high agility, in the range of extant arboreal squirrels. In contrast, the fossil aplodontid Mesogaulus paniensis has a relatively low agility score, similar to the fossorial Aplodontia rufa, the only living aplodontid rodent. This result is in agreement with its postcranial reconstruction as fossorial and with previous indications that early aplodontids were more arboreal than their burrowing descendants.

Orientation Preference Maps in Microcebus murinus Reveal Size-Invariant Design Principles in Primate Visual Cortex

Orientation preference maps (OPMs) are a prominent feature of primary visual cortex (V1) organization in many primates and carnivores. In rodents, neurons are not organized in OPMs but are instead interspersed in a “salt and pepper” fashion, although clusters of orientation-selective neurons have been reported. Does this fundamental difference reflect the existence of a lower size limit for orientation columns (OCs) below which they cannot be scaled down with decreasing V1 size? To address this question, we examined V1 of one of the smallest living primates, the 60-g prosimian mouse lemur (Microcebus murinus). Using chronic intrinsic signal imaging, we found that mouse lemur V1 contains robust OCs, which are arranged in a pinwheel-like fashion. OC size in mouse lemurs was found to be only marginally smaller compared to the macaque, suggesting that these circuit elements are nearly incompressible. The spatial arrangement of pinwheels is well described by a common mathematical design of primate V1 circuit organization. In order to accommodate OPMs, we found that the mouse lemur V1 covers one-fifth of the cortical surface, which is one of the largest V1-to-cortex ratios found in primates. These results indicate that the primate-type visual cortical circuit organization is constrained by a size limitation and raises the possibility that its emergence might have evolved by disruptive innovation rather than gradual change.

Orientation preference maps are a hallmark of V1 organization in all primates studied thus far, yet they are absent in rodents. It is uncertain whether these structures scale with body or brain size. Using intrinsic signal imaging, Ho et al. reveal the presence of such maps in the V1 of the world’s smallest primate, the mouse lemur (Microcebus murinus).

MicroRNA retrocopies generated via L1-mediated retrotransposition in placental mammals help to reveal how their parental genes were transcribed

In mammalian genomes, most retrocopies emerged via the L1 retrotransposition machinery. The hallmarks of an L1-mediated retrocopy, i.e., the intronlessness, the presence of a 3′ poly-A tail, and the TSDs at both ends, were frequently used to identify retrotransposition events. However, most previous studies only focused on protein-coding genes as their possible parental sources and thus only a few retrocopies derived from non-coding genes were reported. Remarkably, none of them was from microRNAs. Here in this study, we found several retrocopies generated from the mir-302–367 cluster gene (MIR302CHG), and identified a novel alternatively spliced exon encoding mir-302a. The other recognized microRNA retrotransposition events are primate-specific with mir-373 and mir-498 as their parental genes. The 3′ poly-A tracts of these two retrocopy groups were directly attached to the end of the microRNA precursor homologous regions, which suggests that their parental transcripts might alternatively terminate at the end of mir-373 and mir-498. All the three parental microRNAs are highly expressed in specific tissues with elevated retrotransposon activity, such as the embryonic stem cells and the placenta. This might be the reason that our first microRNA retrocopy findings were derived from these three microRNA genes.

Whole mitochondrial genome sequence and phylogenetic relationships of Williams's jerboa (Scarturus williamsi) from Turkey

Williams’s jerboa (Scarturus williamsi), a medium-sized jerboa distributed in Anatolia and its adjacent regions, is a member of the four- and five-toed jerboas found mostly in Asia. Disagreements about the taxonomy of this taxon at the genus/species level continue to exist. Here, we report the first effort to sequence and assemble the mitochondrial genome of Williams’s jerboa from Turkey. The mitochondrial genome of S. williamsi was 16,653 bp in total length and contained 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and two non-coding regions (the D-loop and O_L region) with intergenic spacer. All of the genes, except ND6 and eight tRNAs, were encoded on the heavy chain strand, similar to the features of mitogenomes of other rodents. When compared with all available rodent mitochondrial genomes, Williams’s jerboa showed (1) a serine deletion at the 3′-end of the ATP8 gene, (2) the ND5 gene terminated with a TAG codon and (3) a tandem repeat cluster (273 bp in length) in the control region. Williams’s jerboa and Siberian jerboa grouped as sister taxa despite the high genetic distance (17.6%) between them, belonging to Allactaginae. This result is consistent with the latest pre-revision, which suggests that Williams’s jerboa and the Siberian jerboa may belong to separate genera, as Scarturus and Orientallactaga, respectively. The present study provides a reference mitochondrial genome for Williams’s jerboa for further molecular studies of other species of Dipodoidea and Rodentia.

Retrotransposons spread potential cis-regulatory elements during mammary gland evolution

Acquisition of cis-elements is a major driving force for rewiring a gene regulatory network. Several kinds of transposable elements (TEs), mostly retrotransposons that propagate via a copy-and-paste mechanism, are known to possess transcription factor binding motifs and have provided source sequences for enhancers/promoters. However, it remains largely unknown whether retrotransposons have spread the binding sites of master regulators of morphogenesis and accelerated cis-regulatory expansion involved in common mammalian morphological features during evolution. Here, I demonstrate that thousands of binding sites for estrogen receptor α (ERα) and three related pioneer factors (FoxA1, GATA3 and AP2γ) that are essential regulators of mammary gland development arose from a spreading of the binding motifs by retrotransposons. The TE-derived functional elements serve primarily as distal enhancers and are enriched around genes associated with mammary gland morphogenesis. The source TEs occurred via a two-phased expansion consisting of mainly L2/MIR in a eutherian ancestor and endogenous retrovirus 1 (ERV1) in simian primates and murines. Thus the build-up of potential sources for cis-elements by retrotransposons followed by their frequent utilization by the host (co-option/exaptation) may have a general accelerating effect on both establishing and diversifying a gene regulatory network, leading to morphological innovation.

A phylogenomic rodent tree reveals the repeated evolution of masseter architectures

Understanding the number of times a trait has evolved is a necessary foundation for comprehending its potential relationships with selective regimes, developmental constraints and evolutionary diversification. Rodents make up over 40% of extant mammalian species, and their ecological and evolutionary success has been partially attributed to the increase in biting efficiency that resulted from a forward shift of one or two portions of the masseter muscle from the zygomatic arch onto the rostrum. This forward shift has occurred in three discrete ways, but the number of times it has occurred has never been explicitly quantified. We estimated an ultrametric phylogeny, the first to include all rodent families, using thousands of ultraconserved elements. We examined support for evolutionary relationships among the five rodent suborders and then incorporated relevant fossils, fitted models of character evolution, and used stochastic character mapping to determine that a portion of the masseter muscle has moved forward onto the rostrum at least seven times (with one reversal) during the approximately 70 Myr history of rodents. Combined, the repeated evolution of this key innovation, its increasing prevalence through time, and the species diversity of clades with this character underscores the adaptive value of improved biting efficiency and the relative ease with which some advantageous traits arise.

A Novel Test for Absolute Fit of Evolutionary Models Provides a Means to Correctly Identify the Substitution Model and the Model Tree

A novel test is described that visualizes the absolute model-data fit of the substitution and tree components of an evolutionary model. The test utilizes statistics based on counts of character state matches and mismatches in alignments of observed and simulated sequences. This comparison is used to assess model-data fit. In simulations conducted to evaluate the performance of the test, the test estimator was able to identify both the correct tree topology and substitution model under conditions where the Goldman-Cox test-which tests the fit of a substitution model to sequence data and is also based on comparing simulated replicates with observed data-showed high error rates. The novel test was found to identify the correct tree topology within a wide range of DNA substitution model misspecifications, indicating the high discriminatory power of the test. Use of this test provides a practical approach for assessing absolute model-data fit when testing phylogenetic hypotheses.

Generalist Eimeria species in rodents: Multilocus analyses indicate inadequate resolution of established markers

Intracellular parasites of the genus Eimeria are described as tissue/host-specific. Phylogenetic classification of rodent Eimeria suggested that some species have a broader host range than previously assumed. We explore whether Eimeria spp. infecting house mice are misclassified by the most widely used molecular markers due to a lack of resolution, or whether, instead, these parasite species are indeed infecting multiple host species.With the commonly used markers (18S/COI), we recovered monophyletic clades of E. falciformis and E. vermiformis from Mus that included E. apionodes identified in other rodent host species (Apodemus spp., Myodes glareolus, and Microtus arvalis). A lack of internal resolution in these clades could suggest the existence of a species complex with a wide host range infecting murid and cricetid rodents. We question, however, the power of COI and 18S markers to provide adequate resolution for assessing host specificity. In addition to the rarely used marker ORF470 from the apicoplast genome, we present multilocus genotyping as an alternative approach. Phylogenetic analysis of 35 nuclear markers differentiated E. falciformis from house mice from isolates from Apodemus hosts. Isolates of E. vermiformis from Mus are still found in clusters interspersed with non-Mus isolates, even with this high-resolution data.In conclusion, we show that species-level resolution should not be assumed for COI and 18S markers in coccidia. Host-parasite cospeciation at shallow phylogenetic nodes, as well as contemporary coccidian host ranges more generally, is still open questions that need to be addressed using novel genetic markers with higher resolution.

Necroptosis and Caspase-2-Mediated Apoptosis of Astrocytes and Neurons, but Not Microglia, of Rat Hippocampus and Parenchyma Caused by Angiostrongylus cantonensis Infection

Infection with the roundworm Angiostrongylus cantonensis is the main cause of eosinophilic meningitis worldwide. The underlying molecular basis of the various pathological outcomes in permissive and non-permissive hosts infected with A. cantonensis remains poorly defined. In the present study, the histology of neurological disorders in the central nervous system (CNS) of infected rats was assessed by using hematoxylin and eosin staining. Quantitative reverse transcription polymerase chain reaction (RT-qPCR), western blot and immunofluorescence (IF) were used in evolutions of the transcription and translation levels of the apoptosis-, necroptosis-, autophagy-, and pyroptosis-related genes. The distribution of apoptotic and necroptotic cells in the rat hippocampus and parenchyma was further detected using flow cytometry, and the features of the ultrastructure of the cells were examined by transmission electron microscopy (TEM). The inflammatory response upon CNS infection with A. cantonensis evolved, as characterized by the accumulation of a small number of inflammatory cells under the thickened meninges, which peaked at 21 days post-infection (dpi) and returned to normal by 35 dpi. The transcription levels and translation of caspase-2, caspase-8, RIP1 and RIP3 increased significantly at 21 and 28 dpi but decreased sharply at 35 dpi compared to those in the normal control group. However, the changes in the expression of caspase-1, caspase-3, caspase-11, Beclin-1 and LC3B were not obvious, suggesting that apoptosis and necroptosis but not autophagy or pyroptosis occurred in the brains of infected animals at 21 and 28 dpi. The results of RT-qPCR, western blot analysis, IF, flow cytometry and TEM further illustrated that necroptosis and caspase-2-mediated apoptosis occurred in astrocytes and neurons but not in microglia in the parenchyma and hippocampus of infected animals. This study provides the first evidence that neuronal and astrocytic necroptosis and caspase-2-mediated apoptosis are induced by A. cantonensis infection in the parenchymal and hippocampal regions of rats at 21 and 28 dpi but these processes are negligible at 35 dpi. These findings enhance our understanding of the pathogenesis of A. cantonensis infection and provide new insights into therapeutic approaches targeting the occurrence of cell death in astrocytes and neurons in infected patients.

Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation

Big, time-scaled phylogenies are fundamental to connecting evolutionary processes to modern biodiversity patterns. Yet inferring reliable phylogenetic trees for thousands of species involves numerous trade-offs that have limited their utility to comparative biologists. To establish a robust evolutionary timescale for all approximately 6,000 living species of mammals, we developed credible sets of trees that capture root-to-tip uncertainty in topology and divergence times. Our "backbone-and-patch" approach to tree building applies a newly assembled 31-gene supermatrix to two levels of Bayesian inference: (1) backbone relationships and ages among major lineages, using fossil node or tip dating, and (2) species-level "patch" phylogenies with nonoverlapping in-groups that each correspond to one representative lineage in the backbone. Species unsampled for DNA are either excluded ("DNA-only" trees) or imputed within taxonomic constraints using branch lengths drawn from local birth-death models ("completed" trees). Joining time-scaled patches to backbones results in species-level trees of extant Mammalia with all branches estimated under the same modeling framework, thereby facilitating rate comparisons among lineages as disparate as marsupials and placentals. We compare our phylogenetic trees to previous estimates of mammal-wide phylogeny and divergence times, finding that (1) node ages are broadly concordant among studies, and (2) recent (tip-level) rates of speciation are estimated more accurately in our study than in previous "supertree" approaches, in which unresolved nodes led to branch-length artifacts. Credible sets of mammalian phylogenetic history are now available for download at http://vertlife.org/phylosubsets, enabling investigations of long-standing questions in comparative biology.

Two novel cricetine mitogenomes: Insight into the mitogenomic characteristics and phylogeny in Cricetinae (Rodentia: Cricetidae)

Both Cricetus cricetus and Phodopus sungorus mitochondrial genomes (mitogenomes) were sequenced and elaborated for the first time in the present study. Their mitogenomes contained 37 genes and showed typical characteristics of the vertebrate mitogenome. Comparative analysis of 10 cricetine mitogenomes indicated that they shared similar characteristics with those of other cricetines in terms of genes arrangement, nucleotide composition, codon usage, tRNA structure, nucleotide skew and the origin of replication of light strand. Phylogenetic relationship of the subfamily Cricetinae was reconstructed using mitogenomes data with the methods of Bayesian Inference and Maximum Likelihood. Phylogenetic analysis indicated that Cricetulus kamensis was at basal position and phylogenetically distant from all other Cricetulus species but had a close relationship with the group of Phodopus, and supported that the genus Urocricetus deserved as a separate genus rank. The phylogenetic status of Tscherskia triton represented a separate clade corresponding to a diversified cricetine lineage (Cricetulus, Allocricetulus, and Cricetus).

Emergence of hystricognathous rodents: Palaeogene fossil record, phylogeny, dental evolution and historical biogeography

Although phylogenetic trees imply Asia as the ancestral homeland of the Hystricognathi clade (Rodentia: Ctenohystrica), curiously the oldest known fossil occurrences of hystricognathous rodents are not from Asia, but from Africa and South America, where they appear suddenly in the fossil record of both landmasses by the Late Middle Eocene. Here we performed cladistic and Bayesian (standard and tip-dating analyses) assessments of the dental evidence documenting early ctenohystricans, including several Asian ‘ctenodactyloids’, virtually all Palaeogene Asian and African hystricognaths known thus far and two representatives of the earliest known South American hystricognaths. Our results provide a phylogenetic context of early hystricognaths (with implications on systematics) and suggest that some Eocene Asian ‘ctenodactyloids’ could be considered as stem hystricognaths and pre-hystricognaths, although they were not recognized as such originally. However, this view does not fill the gap of the Eocene Asian hystricognath record, as the proposed results imply many ghost lineages extending back to the Middle Eocene for several Asian and African taxa. They also imply a complex early historical biogeography of the group, involving multiple dispersal events from Asia to Africa (and possibly from Africa back to Asia) and then to South America sometime during the Middle Eocene. Based on these phylogenetic considerations, we discuss the emergence of hystricognathous rodents from a morpho-anatomical perspective by analysing the differentiation of their masticatory apparatus and chewing movements, notably through the evolution of their dental patterns.

Structural organization, GABAergic and tyrosine hydroxylase expression in the striatum and globus pallidus of the South American plains vizcacha, Lagostomus maximus (Rodentia, Caviomorpha)

The striatum is an essential component of the basal ganglia that regulatessensory processing, motor, cognition, and behavior. Depending on the species, the striatum shows a unique structure called caudate-putamen as in mice, or its separation into two regions called caudate and lenticular nuclei, the latter formed by putamen and globus pallidus areas, as in primates. These structures have two compartments, striosome and matrix. We investigated the structural organization, GABAergic and tyrosine hydroxylase (TH) expression in the striatum and globus pallidus of the South American plains vizcacha, Lagostomus maximus. Its striatum showed regionalization arising from the presence of an internal capsule, and a similar organization to a striosome-matrix compartmentalization. GABAergic neurons in the matrix of caudate exhibited parvalbumin, calretinin, calbindin, GAD65, and NADPH-d-immunoreactivity. These were also expressed in cells of the putamen with the exception of calretinin showing neurofibers localization. Globus pallidus showed parvalbumin- and GAD65-immunoreactive cells, and calretinin- and calbindin-immunoreactive neuropil, plus GABA-A-immunoreactive neurofibers. NADPH-d-, GAD65- and GABA-A-immunoreactive neurons were larger than parvalbumin-, calretinin-, and calbindin-immunoreactive cells, whereas calbindin-immunoreactive cells were the most abundant. In addition, TH-immunoreactive neuropil was observed in the matrix of the striatum. A significant larger TH-immunoreactive area and neuron number was found in females compared to males. The presence of an internal capsule suggests an adaptive advantage concerning motor and cognitive abilities favoring reaction time in response to predators. In an anatomy-evolutive perspective, the striatum of vizcacha seems to be closer to that of humans than to that of laboratory traditional models such as mouse.

Congruence, fossils and the evolutionary tree of rodents and lagomorphs

Given an evolutionary process, we expect distinct categories of heritable data, sampled in ever larger amounts, to converge on a single tree of historical relationships. We tested this assertion by undertaking phylogenetic analyses of a new morphology-DNA dataset for mammals, focusing on Glires and including the oldest known skeletons of geomyoid and Ischyromys rodents. Our results support geomyoids in the mouse-related clade (Myomorpha) and a ricochetal locomotor pattern for the common ancestor of geomyoid rodents. They also support Ischyromys in the squirrel-related clade (Sciuromorpha) and the evolution of sciurids and Aplodontia from extinct, 'protrogomorph'-grade rodents. Moreover, ever larger samples of characters from our dataset increased congruence with an independent, well-corroborated tree. Addition of morphology from fossils increased congruence to a greater extent than addition of morphology from extant taxa, consistent with fossils' temporal proximity to the common ancestors of living species, reflecting the historical, phylogenetic signal present in our data, particularly in morphological characters from fossils. Our results support the widely held but poorly tested intuition that fossils resemble the common ancestors shared by living species, and that fossilizable hard tissues (i.e. bones and teeth) help to reconstruct the evolutionary tree of life.

Adaptive evolutionary expansion of the ribonuclease 6 in Rodentia

Ribonuclease 6 (RNase6 or RNase K6) is a protein that belongs to a superfamily thought to be the sole verte-brate-specific enzyme known for a wide range of physiological functions, including digestion, cytotoxicity, angiogenesis, male reproduction and host defense. In our study, 51 functional genes and 11 pseudogenes were identified from 27 Rodentia species. Intriguingly, in the 3 main lineages of rodents there were multiple RNase6s identified in all species of Ctenohystrica, whereas only a single RNase6 was observed in other Rodentia species examined except for 2 species in the mouse-related clade. The evolutionary scenario of "birth (gene duplication) and death (gene deactivation)" and gene sorting have been demonstrated in Ctenohystrica. In addition, bursts of positive selection, diversification of isoelectric point and positive net charge have been identified in Ctenohystrica, especially at two key sites that are involved in antimicrobial function. Site Trp30 has undergone positive selection and Ile45 has changed into other residues in Group B and Group C of the Ctenohystrica. Our results demonstrated a complex and intriguing evolutionary pattern of rodent RNase6, and indicated that functional modification may have occurred, which establishes an important theoretical foundation for future functional assays in rodent RNase6.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Dynamic temperature-sensitive A-to-I RNA editing in the brain of a heterothermic mammal during hibernation

RNA editing diversifies genomically encoded information to expand the complexity of the transcriptome. In ectothermic organisms, including Drosophila and Cephalopoda, where body temperature mirrors ambient temperature, decreases in environmental temperature lead to increases in A-to-I RNA editing and cause amino acid recoding events that are thought to be adaptive responses to temperature fluctuations. In contrast, endothermic mammals, including humans and mice, typically maintain a constant body temperature despite environmental changes. Here, A-to-I editing primarily targets repeat elements, rarely results in the recoding of amino acids, and plays a critical role in innate immune tolerance. Hibernating ground squirrels provide a unique opportunity to examine RNA editing in a heterothermic mammal whose body temperature varies over 30°C and can be maintained at 5°C for many days during torpor. We profiled the transcriptome in three brain regions at six physiological states to quantify RNA editing and determine whether cold-induced RNA editing modifies the transcriptome as a potential mechanism for neuroprotection at low temperature during hibernation. We identified 5165 A-to-I editing sites in 1205 genes with dynamically increased editing after prolonged cold exposure. The majority (99.6%) of the cold-increased editing sites are outside of previously annotated coding regions, 82.7% lie in SINE-derived repeats, and 12 sites are predicted to recode amino acids. Additionally, A-to-I editing frequencies increase with increasing cold-exposure, demonstrating that ADAR remains active during torpor. Our findings suggest that dynamic A-to-I editing at low body temperature may provide a neuroprotective mechanism to limit aberrant dsRNA accumulation during torpor in the mammalian hibernator.

Mammalian transposable elements and their impacts on genome evolution

Transposable elements (TEs) are genetic elements with the ability to mobilize and replicate themselves in a genome. Mammalian genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. As a result, TEs have had significant impacts on mammalian evolution. Here we summarize the current understanding of TE content in mammal genomes and find that, with a few exceptions, most fall within a predictable range of observations. First, one third to one half of the genome is derived from TEs. Second, most mammalian genomes are dominated by LINE and SINE retrotransposons, more limited LTR retrotransposons, and minimal DNA transposon accumulation. Third, most mammal genome contains at least one family of actively accumulating retrotransposon. Finally, horizontal transfer of TEs among lineages is rare. TE exaptation events are being recognized with increasing frequency. Despite these beneficial aspects of TE content and activity, the majority of TE insertions are neutral or deleterious. To limit the deleterious effects of TE proliferation, the genome has evolved several defense mechanisms that act at the epigenetic, transcriptional, and post-transcriptional levels. The interaction between TEs and these defense mechanisms has led to an evolutionary arms race where TEs are suppressed, evolve to escape suppression, then are suppressed again as the defense mechanisms undergo compensatory change. The result is complex and constantly evolving interactions between TEs and host genomes.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

The Small and the Dead: A Review of Ancient DNA Studies Analysing Micromammal Species

The field of ancient DNA (aDNA) has recently been in a state of exponential growth, largely driven by the uptake of Next Generation Sequencing (NGS) techniques. Much of this work has focused on the mammalian megafauna and ancient humans, with comparatively less studies looking at micromammal fauna, despite the potential of these species in testing evolutionary, environmental and taxonomic theories. Several factors make micromammal fauna ideally suited for aDNA extraction and sequencing. Micromammal subfossil assemblages often include the large number of individuals appropriate for population level analyses, and, furthermore, the assemblages are frequently found in cave sites where the constant temperature and sheltered environment provide favourable conditions for DNA preservation. This review looks at studies that include the use of aDNA in molecular analysis of micromammal fauna, in order to examine the wide array of questions that can be answered in the study of small mammals using new palaeogenetic techniques. This study highlights the bias in current aDNA studies and assesses the future use of aDNA as a tool for the study of micromammal fauna.

Early embryonic development and spatiotemporal localization of mammalian primordial germ cell-associated proteins in the basal rodent Lagostomus maximus

The gene network controlling primordial germ cell (PGC) specification in eutherian mammals has been exhaustively investigated in mice. The egg-cylinder morphology of the mouse embryo is the key event enabling inductive signals from the extra-embryonic ectoderm (ExE) to specify epiblast cells as PGCs early on. We investigated the embryonic development and the spatiotemporal localization of PGC-associated proteins in the basal Hystricognathi rodent Lagostomus maximus. L. maximus develops through a flat-disc epiblast far apart from the ExE. In the primitive streak stage, OCT4-positive cells are detected in the posterior pole of the embryo disc in the mesoderm of the proximal epiblast. In the neural plate stage, a reduced 8 to 12 OCT4-positive cell population transiently expresses FRAGILIS, STELLA and SOX17 in the posterior streak. Soon after translocation to the hindgut, pluripotent OCT4 cells start expressing VASA, and then, STELLA and FRAGILIS are turned on during migration toward the genital ridge. L. maximus shows a spatiotemporal pattern of PGC-associated markers divergent from the early PGC restriction model seen in mice. This pattern conforms to alternative models that are based on a pluripotent population in the embryonic axis, where PGCs are specified later during development.

The Beaver's Phylogenetic Lineage Illuminated by Retroposon Reads

Solving problematic phylogenetic relationships often requires high quality genome data. However, for many organisms such data are still not available. Among rodents, the phylogenetic position of the beaver has always attracted special interest. The arrangement of the beaver's masseter (jaw-closer) muscle once suggested a strong affinity to some sciurid rodents (e.g., squirrels), placing them in the Sciuromorpha suborder. Modern molecular data, however, suggested a closer relationship of beaver to the representatives of the mouse-related clade, but significant data from virtually homoplasy-free markers (for example retroposon insertions) for the exact position of the beaver have not been available. We derived a gross genome assembly from deposited genomic Illumina paired-end reads and extracted thousands of potential phylogenetically informative retroposon markers using the new bioinformatics coordinate extractor fastCOEX, enabling us to evaluate different hypotheses for the phylogenetic position of the beaver. Comparative results provided significant support for a clear relationship between beavers (Castoridae) and kangaroo rat-related species (Geomyoidea) (p < 0.0015, six markers, no conflicting data) within a significantly supported mouse-related clade (including Myodonta, Anomaluromorpha, and Castorimorpha) (p < 0.0015, six markers, no conflicting data).

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Identification and characterisation of Short Interspersed Nuclear Elements in the olive tree (Olea europaea L.) genome

Short Interspersed Nuclear Elements (SINEs) are nonautonomous retrotransposons in the genome of most eukaryotic species. While SINEs have been intensively investigated in humans and other animal systems, SINE identification has been carried out only in a limited number of plant species. This lack of information is apparent especially in non-model plants whose genome has not been sequenced yet. The aim of this work was to produce a specific bioinformatics pipeline for analysing second generation sequence reads of a non-model species and identifying SINEs. We have identified, for the first time, 227 putative SINEs of the olive tree (Olea europaea), that constitute one of the few sets of such sequences in dicotyledonous species. The identified SINEs ranged from 140 to 362 bp in length and were characterised with regard to the occurrence of the tRNA domain in their sequence. The majority of identified elements resulted in single copy or very lowly repeated, often in association with genic sequences. Analysis of sequence similarity allowed us to identify two major groups of SINEs showing different abundances in the olive tree genome, the former with sequence similarity to SINEs of Scrophulariaceae and Solanaceae and the latter to SINEs of Salicaceae. A comparison of sequence conservation between olive SINEs and LTR retrotransposon families suggested that SINE expansion in the genome occurred especially in very ancient times, before LTR retrotransposon expansion, and presumably before the separation of the rosids (to which Oleaceae belong) from the Asterids. Besides providing data on olive SINEs, our results demonstrate the suitability of the pipeline employed for SINE identification. Applying this pipeline will favour further structural and functional analyses on these relatively unknown elements to be performed also in other plant species, even in the absence of a reference genome, and will allow establishing general evolutionary patterns for this kind of repeats in plants.