Where do rodents fit? Evidence from the complete mitochondrial genome of Sciurus vulgaris

Complete mitogenomes of Turkish tree squirrels, Sciurus anomalus and S. vulgaris, (Sciuridae: Rodentia: Mammalia) and their phylogenetic status within the tribe Sciurini

The genus Sciurus, a member of the family Sciuridae, is widely distributed in the Holarctic region. To better understand mitogenomic characteristics and to reveal internal phylogenetic relationships of the genus, 20 complete mitogenomes of Turkish tree squirrels were successfully sequenced for the first time, including 19 for S. anomalus (from 16,505 bp to 16,510 bp) and one for S. vulgaris (16,511 bp). The mitogenomes of two species were AT-biased. All tRNAs for two species displayed a typical clover-leaf structure, except for tRNA^Ser(AGY). The tRNA Serine1 (S1)-GCT structure lacked the dihydrouridine (DHU) loop and stem. Based on mitogenomic dataset for phylogeny of Sciurinae, phylogenetic analyses (Bayesian Inference and Maximum Likelihood) did not support monophyly of Sciurus and proposed that S. anomalus, the most basal taxa in the Sciurini tribe, had at least five mitogenome lineages, which were also supported by network analysis. The dissimilarities among the five linegaes of S. anomalus ranged from 0.0042 (0.42%) to 0.0062 (0.62%) using K2P sequence pairwise distances. In addition to this mitogenomic analysis result, phylogenetic analyses using the CYTB + D-loop dataset proposed the existence of at least nine lineages for S. anomalus, which was different than those of the previous studies. The current study proposed that the use of mitogenomic data for reconstructing the phylogeny of Turkey' Sciurus holds an important value for revealing evolutionary relationships.

Phylogenies from unaligned proteomes using sequence environments of amino acid residues

Alignment-free methods for sequence comparison and phylogeny inference have attracted a great deal of attention in recent years. Several algorithms have been implemented in diverse software packages. Despite the great number of existing methods, most of them are based on word statistics. Although they propose different filtering and weighting strategies and explore different metrics, their performance may be limited by the phylogenetic signal preserved in these words. Herein, we present a different approach based on the species-specific amino acid neighborhood preferences. These differential preferences can be assessed in the context of vector spaces. In this way, a distance-based method to build phylogenies has been developed and implemented into an easy-to-use R package. Tests run on real-world datasets show that this method can reconstruct phylogenetic relationships with high accuracy, and often outperforms other alignment-free approaches. Furthermore, we present evidence that the new method can perform reliably on datasets formed by non-orthologous protein sequences, that is, the method not only does not require the identification of orthologous proteins, but also does not require their presence in the analyzed dataset. These results suggest that the neighborhood preference of amino acids conveys a phylogenetic signal that may be of great utility in phylogenomics.

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.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

Nonlinearity of intracochlear motion and local cochlear microphonic: Comparison between guinea pig and gerbil

Studying the in-vivo mechanical and electrophysiological cochlear responses in several species helps us to have a comprehensive view of the sensitivity and frequency selectivity of the cochlea. Different species might use different mechanisms to achieve the sharp frequency-place map. The outer hair cells (OHC) play an important role in mediating frequency tuning. In the present work, we measured the OHC-generated local cochlear microphonic (LCM) and the motion of different layers in the organ of Corti using optical coherence tomography (OCT) in the first turn of the cochlea in guinea pig. In the best frequency (BF) band, our observations were similar to our previous measurements in gerbil: a nonlinear peak in LCM responses and in the basilar membrane (BM) and OHC-region displacements, and higher motion in the OHC region than the BM. Sub-BF the responses in the two species were different. In both species the sub-BF displacement of the BM was linear and LCM was nonlinear. Sub-BF in the OHC-region, nonlinearity was only observed in a subset of healthy guinea pig cochleae while in gerbil, robust nonlinearity was observed in all healthy cochleae. The differences suggest that gerbils and guinea pigs employ different mechanisms for filtering sub-BF OHC activity from BM responses. However, it cannot be ruled out that the differences are due to technical measurement differences across the species.

Phylogenetic relationships and taxonomic position of genus Hyperacrius (Rodentia: Arvicolinae) from Kashmir based on evidences from analysis of mitochondrial genome and study of skull morphology

In this article, we present the nearly complete mitochondrial genome of the Subalpine Kashmir vole Hyperacrius fertilis (Arvicolinae, Cricetidae, Rodentia), assembled using data from Illumina next-generation sequencing (NGS) of the DNA from a century-old museum specimen. De novo assembly consisted of 16,341 bp and included all mitogenome protein-coding genes as well as 12S and 16S RNAs, tRNAs and D-loop. Using the alignment of protein-coding genes of 14 previously published Arvicolini tribe mitogenomes, seven Clethrionomyini mitogenomes, and also Ondatra and Dicrostonyx outgroups, we conducted phylogenetic reconstructions based on a dataset of 13 protein-coding genes (PCGs) under maximum likelihood and Bayesian inference. Phylogenetic analyses robustly supported the phylogenetic position of this species within the tribe Arvicolini. Among the Arvicolini, Hyperacrius represents one of the early-diverged lineages. This result of phylogenetic analysis altered the conventional view on phylogenetic relatedness between Hyperacrius and Alticola and prompted the revision of morphological characters underlying the former assumption. Morphological analysis performed here confirmed molecular data and provided additional evidence for taxonomic replacement of the genus Hyperacrius from the tribe Clethrionomyini to the tribe Arvicolini.

Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary

Deciphering the timing of the placental mammal radiation is a longstanding problem in evolutionary biology, but consensus on the tempo and mode of placental diversification remains elusive. Nevertheless, an accurate timetree is essential for understanding the role of important events in Earth history (e.g., Cretaceous Terrestrial Revolution, KPg mass extinction) in promoting the taxonomic and ecomorphological diversification of Placentalia. Archibald and Deutschman described three competing models for the diversification of placental mammals, which are the Explosive, Long Fuse, and Short Fuse Models. More recently, the Soft Explosive Model and Trans-KPg Model have emerged as additional hypotheses for the placental radiation. Here, we review molecular and paleontological evidence for each of these five models including the identification of general problems that can negatively impact divergence time estimates. The Long Fuse Model has received more support from relaxed clock studies than any of the other models, but this model is not supported by morphological cladistic studies that position Cretaceous eutherians outside of crown Placentalia. At the same time, morphological cladistics has a poor track record of reconstructing higher-level relationships among the orders of placental mammals including the results of new pseudoextinction analyses that we performed on the largest available morphological data set for mammals (4,541 characters). We also examine the strengths and weaknesses of different timetree methods (node dating, tip dating, and fossilized birth-death dating) that may now be applied to estimate the timing of the placental radiation. While new methods such as tip dating are promising, they also have problems that must be addressed if these methods are to effectively discriminate among competing hypotheses for placental diversification. Finally, we discuss the complexities of timetree estimation when the signal of speciation times is impacted by incomplete lineage sorting (ILS) and hybridization. Not accounting for ILS results in dates that are older than speciation events. Hybridization, in turn, can result in dates than are younger or older than speciation dates. Disregarding this potential variation in "gene" history across the genome can distort phylogenetic branch lengths and divergence estimates when multiple unlinked genomic loci are combined together in a timetree analysis.

K2 and K2*: efficient alignment-free sequence similarity measurement based on Kendall statistics

Motivation

Alignment-free sequence comparison methods can compute the pairwise similarity between a huge number of sequences much faster than sequence-alignment based methods.

Results

We propose a new non-parametric alignment-free sequence comparison method, called K2, based on the Kendall statistics. Comparing to the other state-of-the-art alignment-free comparison methods, K2 demonstrates competitive performance in generating the phylogenetic tree, in evaluating functionally related regulatory sequences, and in computing the edit distance (similarity/dissimilarity) between sequences. Furthermore, the K2 approach is much faster than the other methods. An improved method, K2*, is also proposed, which is able to determine the appropriate algorithmic parameter (length) automatically, without first considering different values. Comparative analysis with the state-of-the-art alignment-free sequence similarity methods demonstrates the superiority of the proposed approaches, especially with increasing sequence length, or increasing dataset sizes.

Availability and implementation

The K2 and K2* approaches are implemented in the R language as a package and is freely available for open access (http://community.wvu.edu/daadjeroh/projects/K2/K2_1.0.tar.gz).

Contact

yueljiang@163.com.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias)

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

Genetic characterization of the Turkish gray hamster (Cricetulus migratorius) based on mitochondrial cytochrome b and 12S rRNA sequences

Although genetic diversity and phylogenetic status of the gray hamster (Cricetulus migratorius) have been investigated from different regions in previous studies, genetic data on this species from Turkey are still lacking, since previous data have been based on a limited number of gray hamsters sampled across the Anatolian part of Turkey. The aim of this study was to determine the genetic diversity of the Anatolian population and to reveal the phylogenetic relationships among the Anatolian population and conspecific populations of the gray hamster. The complete and partial fragments of mitochondrial Cyt b and 12S rRNA from the 20 Turkish samples were amplified and sequenced. Ten 12S rRNA (901 bp) and 15 Cyt b (1140 bp) haplotypes found in this work were not previously reported. Based on Bayesian, Maximum Likelihood, Neighbour-Joining and Median-Joining network analyses by using mitochondrial data under the name Cricetulus, the results of phylogenetic and network analyses indicated that there was a deep separation among the distinct lineages within the genus Cricetulus. When considering the position of the Turkish haplotypes in median joining network, the Anatolian part of Turkey may have hosted a source population of the gray hamster for expansion to adjacent regions in the past period. Additionally, the Anatolian population of gray hamster had relatively high haplotype diversity and the present study propounded the importance of data obtained from the Anatolian population of gray hamster to reveal the phylogenetic relationships among conspecific populations of the gray hamster.

Clustering DNA sequences using the out-of-place measure with reduced n-grams

The alignment-free n-gram based method with the out-of-place measures as the distance has been successfully applied to automatic text or natural languages categorization in real time. However, it is not clear about its performance and the selection of n for comparing genome sequences. Here we propose a symmetric version of the out-of-place measure and a new approach for finding the optimal range of n to construct a phylogenetic tree with the symmetric out-of-place measures. Our method is then applied to real genome sequence datasets. The resulting phylogenetic trees are matching with the standard biological classification. It shows that our proposed method is a very powerful tool for phylogenetic analysis in terms of both classification accuracy and computation efficiency.

Primate phylogenetic relationships and divergence dates inferred from complete mitochondrial genomes

The origins and the divergence times of the most basal lineages within primates have been difficult to resolve mainly due to the incomplete sampling of early fossil taxa. The main source of contention is related to the discordance between molecular and fossil estimates: while there are no crown primate fossils older than 56Ma, most molecule-based estimates extend the origins of crown primates into the Cretaceous. Here we present a comprehensive mitogenomic study of primates. We assembled 87 mammalian mitochondrial genomes, including 62 primate species representing all the families of the order. We newly sequenced eleven mitochondrial genomes, including eight Old World monkeys and three strepsirrhines. Phylogenetic analyses support a strong topology, confirming the monophyly for all the major primate clades. In contrast to previous mitogenomic studies, the positions of tarsiers and colugos relative to strepsirrhines and anthropoids are well resolved. In order to improve our understanding of how fossil calibrations affect age estimates within primates, we explore the effect of seventeen fossil calibrations across primates and other mammalian groups and we select a subset of calibrations to date our mitogenomic tree. The divergence date estimates of the Strepsirrhine/Haplorhine split support an origin of crown primates in the Late Cretaceous, at around 74Ma. This result supports a short-fuse model of primate origins, whereby relatively little time passed between the origin of the order and the diversification of its major clades. It also suggests that the early primate fossil record is likely poorly sampled.

Phylogenetic analysis of genome rearrangements among five mammalian orders

Evolutionary relationships among placental mammalian orders have been controversial. Whole genome sequencing and new computational methods offer opportunities to resolve the relationships among 10 genomes belonging to the mammalian orders Primates, Rodentia, Carnivora, Perissodactyla and Artiodactyla. By application of the double cut and join distance metric, where gene order is the phylogenetic character, we computed genomic distances among the sampled mammalian genomes. With a marsupial outgroup, the gene order tree supported a topology in which Rodentia fell outside the cluster of Primates, Carnivora, Perissodactyla, and Artiodactyla. Results of breakpoint reuse rate and synteny block length analyses were consistent with the prediction of random breakage model, which provided a diagnostic test to support use of gene order as an appropriate phylogenetic character in this study. We discussed the influence of rate differences among lineages and other factors that may contribute to different resolutions of mammalian ordinal relationships by different methods of phylogenetic reconstruction.

Complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Sciuromorpha, Sciuridae) and revision of rodent phylogeny

In this study, the complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Rodentia, Sciuromorpha, Sciuridae) was sequenced and characterized in detail. The entire mitochondrial genome of P. volans consisted of 16,513 bp and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. Its gene arrangement pattern was consistent with the mammalian ground pattern. The overall base composition and AT contents were similar to those of other rodent mitochondrial genomes. The light-strand origin generally identified between tRNA ( Asn ) and tRNA ( Cys ) consisted of a secondary structure with an 11-bp stem and an 11-bp loop. The large control region was constructed of three characteristic domains, ETAS, CD, and CSB without any repeat sequences. Each domain contained ETAS1, subsequences A, B, and C, and CSB1, respectively. In order to examine phylogenetic contentious issues of the monophyly of rodents and phylogenetic relationships among five rodent suborders, here, phylogenetic analyses based on nucleotide sequence data of the 35 rodent and 3 lagomorph mitochondrial genomes were performed using the Bayesian inference and maximum likelihood method. The result strongly supported the rodent monophyly with high node confidence values (BP 100 % in ML and BPP 1.00 in BI) and also monophylies of four rodent suborders (BP 85-100 % in ML and BPP 1.00 in BI), except for Anomalumorpha in which only one species was examined here. Also, phylogenetic relationships among the five rodent suborders were suggested and discussed in detail.

Generation of trophoblast stem cells from rabbit embryonic stem cells with BMP4

Trophoblast stem (TS) cells are ideal models to investigate trophectoderm differentiation and placental development. Herein, we describe the derivation of rabbit trophoblast stem cells from embryonic stem (ES) cells. Rabbit ES cells generated in our laboratory were induced to differentiate in the presence of BMP4 and TS-like cell colonies were isolated and expanded. These cells expressed the molecular markers of mouse TS cells, were able to invade, give rise to derivatives of TS cells, and chimerize placental tissues when injected into blastocysts. The rabbit TS-like cells maintained self-renewal in culture medium with serum but without growth factors or feeder cells, whilst their proliferation and identity were compromised by inhibitors of FGFs and TGF-β receptors. Taken together, our study demonstrated the derivation of rabbit TS cells and suggested the essential roles of FGF and TGF-β signalings in maintenance of rabbit TS cell self-renewal.

An efficient binomial model-based measure for sequence comparison and its application

Sequence comparison is one of the major tasks in bioinformatics, which could serve as evidence of structural and functional conservation, as well as of evolutionary relations. There are several similarity/dissimilarity measures for sequence comparison, but challenges remains. This paper presented a binomial model-based measure to analyze biological sequences. With help of a random indicator, the occurrence of a word at any position of sequence can be regarded as a random Bernoulli variable, and the distribution of a sum of the word occurrence is well known to be a binomial one. By using a recursive formula, we computed the binomial probability of the word count and proposed a binomial model-based measure based on the relative entropy. The proposed measure was tested by extensive experiments including classification of HEV genotypes and phylogenetic analysis, and further compared with alignment-based and alignment-free measures. The results demonstrate that the proposed measure based on binomial model is more efficient.

Dynamic refolding of IFN-gamma mRNA enables it to function as PKR activator and translation template

Interferon-gamma mRNA activates the RNA-dependent protein kinase PKR, which in turn strongly attenuates translation of interferon-gamma mRNA. Unlike riboswitches restricted to noncoding regions, the interferon-gamma RNA domain that activates PKR comprises the 5' UTR and 26 translated codons. Extensive interferon-gamma coding sequence is thus dedicated to activating PKR and blocking interferon-gamma synthesis. This implies that the PKR activator is disrupted by ribosomes during translation initiation and must refold promptly to restore PKR activation. The activator structure harbors an essential kink-turn, probably to allow formation of a pseudoknot that is critical for PKR activation. Three indispensable short helices, bordered by orientation-sensitive base pairs, align with the pseudoknot stem, generating RNA helix of sufficient length to activate PKR. Through gain-of-function mutations, we show that the RNA activator can adopt alternative conformations that activate PKR. This flexibility promotes efficient refolding of interferon-gamma mRNA, which is necessary for its dual function as translation template and activator of PKR, and which thus prevents overexpression of this inflammatory cytokine.

Posttranscriptional regulation of pineal melatonin synthesis in Octodon degus

Small laboratory animals have provided significant information about melatonin regulation, yet most of these organisms are nocturnal and regulate melatonin synthesis by mechanisms that diverge from those of humans. For example, in all rodents examined, melatonin secretion occurs with a time lag of several hours after the onset of darkness; in addition, arylalkylamine N-acetyltransferase (AANAT), the key enzyme in melatonin synthesis, displays dynamic transcriptional activation specifically at night in all rodents studied to date. In ungulates and primates including humans, on the other hand, melatonin secretion occurs immediately during the early night and is controlled by circadian posttranscriptional regulation of AANAT. We hypothesize that the diurnal Octodon degus (an Hystricognath rodent) could serve as an improved experimental model for studies of human melatonin regulation. To test this, we monitored melatonin production in degus using pineal microdialysis and characterized the regulation of melatonin synthesis by analyzing degu Aanat. Degu pineal melatonin rises with little latency at night, as in ungulates and primates. In addition, degu Aanat mRNA expression displays no detectable diurnal variation, suggesting that, like ungulates and primates, melatonin in this species is regulated by a posttranscriptional mechanism. Compared with AANAT from all rodents examined to date, the predicted amino acid sequence of degu AANAT is phylogenetically more closely related to ungulate and primate AANAT. These data suggest that Octodon degus may provide an ideal model system for laboratory investigation of mechanisms of melatonin synthesis and secretion in diurnal mammals.

Breakpoint graphs and ancestral genome reconstructions

Recently completed whole-genome sequencing projects marked the transition from gene-based phylogenetic studies to phylogenomics analysis of entire genomes. We developed an algorithm MGRA for reconstructing ancestral genomes and used it to study the rearrangement history of seven mammalian genomes: human, chimpanzee, macaque, mouse, rat, dog, and opossum. MGRA relies on the notion of the multiple breakpoint graphs to overcome some limitations of the existing approaches to ancestral genome reconstructions. MGRA also generates the rearrangement-based characters guiding the phylogenetic tree reconstruction when the phylogeny is unknown.

Rodent phylogeny revised: analysis of six nuclear genes from all major rodent clades

BACKGROUND

Rodentia is the most diverse order of placental mammals, with extant rodent species representing about half of all placental diversity. In spite of many morphological and molecular studies, the family-level relationships among rodents and the location of the rodent root are still debated. Although various datasets have already been analyzed to solve rodent phylogeny at the family level, these are difficult to combine because they involve different taxa and genes.

RESULTS

We present here the largest protein-coding dataset used to study rodent relationships. It comprises six nuclear genes, 41 rodent species, and eight outgroups. Our phylogenetic reconstructions strongly support the division of Rodentia into three clades: (1) a "squirrel-related clade", (2) a "mouse-related clade", and (3) Ctenohystrica. Almost all evolutionary relationships within these clades are also highly supported. The primary remaining uncertainty is the position of the root. The application of various models and techniques aimed to remove non-phylogenetic signal was unable to solve the basal rodent trifurcation.

CONCLUSION

Sequencing and analyzing a large sequence dataset enabled us to resolve most of the evolutionary relationships among Rodentia. Our findings suggest that the uncertainty regarding the position of the rodent root reflects the rapid rodent radiation that occurred in the Paleocene rather than the presence of conflicting phylogenetic and non-phylogenetic signals in the dataset.

Phylogenetic inference from binary sequences reduced by primary DNA sequences

Given a bi-classification of nucleotides, we can obtain a reduced binary sequence of a primary DNA sequence. This binary sequence will undoubtedly retain some biological information and lose the rest. Here we want to know what kind of and how much biological information an individual binary sequence carries. Three classifications of nucleotides are explored in the present article. Phylogenetic trees are built from these binary sequences by the Neighbor-Joining (NJ) method, with evolutionary distance evaluated on the basis of a symbolic sequence complexity. We find that, for all data sets studied, binary sequences reduced by the purine/pyrimidine classification give reliable phylogeny (almost the same as that from the primary sequences), while the other two result in discrepancies at different levels. Some possible reasons and a simple model of sequence evolutionary are introduced to interpret this phenomenon.

A Poisson model of sequence comparison and its application to coronavirus phylogeny

In this paper, we propose two metrics to compare DNA and protein sequences based on a Poisson model of word occurrences. Instead of comparing the frequencies of all fixed-length words in two sequences, we consider (1) the probability of ‘generating’ one sequence under the Poisson model estimated from the other; (2) their different expression levels of words. Phylogenetic trees of 25 viruses including SARS-CoVs are constructed to illustrate our approach.

Markov model plus k-word distributions: a synergy that produces novel statistical measures for sequence comparison

MOTIVATION

Many proposed statistical measures can efficiently compare biological sequences to further infer their structures, functions and evolutionary information. They are related in spirit because all the ideas for sequence comparison try to use the information on the k-word distributions, Markov model or both. Motivated by adding k-word distributions to Markov model directly, we investigated two novel statistical measures for sequence comparison, called wre.k.r and S2.k.r.

RESULTS

The proposed measures were tested by similarity search, evaluation on functionally related regulatory sequences and phylogenetic analysis. This offers the systematic and quantitative experimental assessment of our measures. Moreover, we compared our achievements with these based on alignment or alignment-free. We grouped our experiments into two sets. The first one, performed via ROC (receiver operating curve) analysis, aims at assessing the intrinsic ability of our statistical measures to search for similar sequences from a database and discriminate functionally related regulatory sequences from unrelated sequences. The second one aims at assessing how well our statistical measure is used for phylogenetic analysis. The experimental assessment demonstrates that our similarity measures intending to incorporate k-word distributions into Markov model are more efficient.

Analysis of septins across kingdoms reveals orthology and new motifs

Background

Septins are cytoskeletal GTPase proteins first discovered in the fungus Saccharomyces cerevisiae where they organize the septum and link nuclear division with cell division. More recently septins have been found in animals where they are important in processes ranging from actin and microtubule organization to embryonic patterning and where defects in septins have been implicated in human disease. Previous studies suggested that many animal septins fell into independent evolutionary groups, confounding cross-kingdom comparison.

Results

In the current work, we identified 162 septins from fungi, microsporidia and animals and analyzed their phylogenetic relationships. There was support for five groups of septins with orthology between kingdoms. Group 1 (which includes S. cerevisiae Cdc10p and human Sept9) and Group 2 (which includes S. cerevisiae Cdc3p and human Sept7) contain sequences from fungi and animals. Group 3 (which includes S. cerevisiae Cdc11p) and Group 4 (which includes S. cerevisiae Cdc12p) contain sequences from fungi and microsporidia. Group 5 (which includes Aspergillus nidulans AspE) contains sequences from filamentous fungi. We suggest a modified nomenclature based on these phylogenetic relationships. Comparative sequence alignments revealed septin derivatives of already known G1, G3 and G4 GTPase motifs, four new motifs from two to twelve amino acids long and six conserved single amino acid positions. One of these new motifs is septin-specific and several are group specific.

Conclusion

Our studies provide an evolutionary history for this important family of proteins and a framework and consistent nomenclature for comparison of septin orthologs across kingdoms.

Site specific rates of mitochondrial genomes and the phylogeny of eutheria

BACKGROUND

Traditionally, most studies employing data from whole mitochondrial genomes to diagnose relationships among the major lineages of mammals have attempted to exclude regions that potentially complicate phylogenetic analysis. Components generally excluded are 3rd codon positions of protein-encoding genes, the control region, rRNAs, tRNAs, and the ND6 gene (encoded on the opposite strand). We present an approach that includes all the data, with the exception of the control region. This approach is based on a site-specific rate model that accommodates excessive homoplasy and that utilizes secondary structure as a reference for proper alignment of rRNAs and tRNAs.

RESULTS

Mitochondrial genomic data for 78 eutherian mammals, 8 metatherians, and 3 monotremes were analyzed with a Bayesian analysis and our site specific rate model. The resultant phylogeny revealed strong support for most nodes and was highly congruent with more recent phylogenies based on nuclear DNA sequences. In addition, many of the conflicting relationships observed by earlier mitochondrial-based analyses were resolved without need for the exclusion of large subsets of the data.

CONCLUSION

Rather than exclusion of data to minimize presumed noise associated with non-protein encoding genes in the mitochondrial genome, our results indicate that selection of an appropriate model that accommodates rate heterogeneity across data partitions and proper treatment of RNA genes can result in a mitochondrial genome-based phylogeny of eutherian mammals that is reasonably congruent with recent phylogenies derived from nuclear genes.

A phylogenomic study of human, dog, and mouse

In recent years the phylogenetic relationship of mammalian orders has been addressed in a number of molecular studies. These analyses have frequently yielded inconsistent results with respect to some basal ordinal relationships. For example, the relative placement of primates, rodents, and carnivores has differed in various studies. Here, we attempt to resolve this phylogenetic problem by using data from completely sequenced nuclear genomes to base the analyses on the largest possible amount of data. To minimize the risk of reconstruction artifacts, the trees were reconstructed under different criteria-distance, parsimony, and likelihood. For the distance trees, distance metrics that measure independent phenomena (amino acid replacement, synonymous substitution, and gene reordering) were used, as it is highly improbable that all of the trees would be affected the same way by any reconstruction artifact. In contradiction to the currently favored classification, our results based on full-genome analysis of the phylogenetic relationship between human, dog, and mouse yielded overwhelming support for a primate-carnivore clade with the exclusion of rodents.

PTreeRec: Phylogenetic Tree Reconstruction based on genome BLAST distance

Phylogenetic Tree Reconstruction (PTreeRec) is a web-based tool for automatic phylogeny inferences from whole-genome sequences, which accepts files of DNA sequences in the FASTA format and allows users to save the output tree file, and displays the inferred tree through an applet in a web browser. PTreeRec involves three basic steps. First, regions of maximal segment pairs (MSPs) based on an all-against-all pairwise comparison of genomes are located. Second, a distance matrix is calculated from MSP scores or coverage. Finally, a phylogenetic tree is reconstructed by the neighbor-joining method.

LZ complexity distance of DNA sequences and its application in phylogenetic tree reconstruction

DNA sequences can be treated as finite-length symbol strings over a four-letter alphabet (A, C, T, G). As a universal and computable complexity measure, LZ complexity is valid to describe the complexity of DNA sequences. In this study, a concept of conditional LZ complexity between two sequences is proposed according to the principle of LZ complexity measure. An LZ complexity distance metric between two nonnull sequences is defined by utilizing conditional LZ complexity. Based on LZ complexity distance, a phylogenetic tree of 26 species of placental mammals (Eutheria) with three outgroup species was reconstructed from their complete mitochondrial genomes. On the debate that which two of the three main groups of placental mammals, namely Primates, Ferungulates, and Rodents, are more closely related, the phylogenetic tree reconstructed based on LZ complexity distance supports the suggestion that Primates and Ferungulates are more closely related.

Retroposed elements as archives for the evolutionary history of placental mammals

Reconstruction of the placental mammalian (eutherian) evolutionary tree has undergone diverse revisions, and numerous aspects remain hotly debated. Initial hierarchical divisions based on morphology contained many misgroupings due to features that evolved independently by similar selection processes. Molecular analyses corrected many of these misgroupings and the superordinal hierarchy of placental mammals was recently assembled into four clades. However, long or rapid evolutionary periods, as well as directional mutation pressure, can produce molecular homoplasies, similar characteristics lacking common ancestors. Retroposed elements, by contrast, integrate randomly into genomes with negligible probabilities of the same element integrating independently into orthologous positions in different species. Thus, presence/absence analyses of these elements are a superior strategy for molecular systematics. By computationally scanning more than 160,000 chromosomal loci and judiciously selecting from only phylogenetically informative retroposons for experimental high-throughput PCR applications, we recovered 28 clear, independent monophyly markers that conclusively verify the earliest divergences in placental mammalian evolution. Using tests that take into account ancestral polymorphisms, multiple long interspersed elements and long terminal repeat element insertions provide highly significant evidence for the monophyletic clades Boreotheria (synonymous with Boreoeutheria), Supraprimates (synonymous with Euarchontoglires), and Laurasiatheria. More importantly, two retropositions provide new support for a prior scenario of early mammalian evolution that places the basal placental divergence between Xenarthra and Epitheria, the latter comprising all remaining placentals. Due to its virtually homoplasy-free nature, the analysis of retroposon presence/absence patterns avoids the pitfalls of other molecular methodologies and provides a rapid, unequivocal means for revealing the evolutionary history of organisms.

The authors identified and sequenced retroposons in mammalian genomes. The presence and absence of these retroposons provided evolutionary markers from which the authors reconstructed the phylogenetic history of placental mammals.

The Average Common Substring Approach to Phylogenomic Reconstruction

We describe a novel method for efficient reconstruction of phylogenetic trees, based on sequences of whole genomes or proteomes, whose lengths may greatly vary. The core of our method is a new measure of pairwise distances between sequences. This measure is based on computing the average lengths of maximum common substrings, which is intrinsically related to information theoretic tools (Kullback-Leibler relative entropy). We present an algorithm for efficiently computing these distances. In principle, the distance of two l long sequences can be calculated in O(l) time. We implemented the algorithm using suffix arrays our implementation is fast enough to enable the construction of the proteome phylogenomic tree for hundreds of species and the genome phylogenomic forest for almost two thousand viruses. An initial analysis of the results exhibits a remarkable agreement with "acceptable phylogenetic and taxonomic truth." To assess our approach, our results were compared to the traditional (single-gene or protein-based) maximum likelihood method. The obtained trees were compared to implementations of a number of alternative approaches, including two that were previously published in the literature, and to the published results of a third approach. Comparing their outcome and running time to ours, using a "traditional" trees and a standard tree comparison method, our algorithm improved upon the "competition" by a substantial margin. The simplicity and speed of our method allows for a whole genome analysis with the greatest scope attempted so far. We describe here five different applications of the method, which not only show the validity of the method, but also suggest a number of novel phylogenetic insights.

Geographical distance and physical barriers shape the genetic structure of Eurasian red squirrels (Sciurus vulgaris) in the Italian Alps

Red squirrels (Sciurus vulgaris) are widely distributed throughout Eurasia, occurring in many types of coniferous and mixed-deciduous forests. In fragmented landscapes, small and partly isolated populations with low immigration rates show reduced genetic diversity, but reforestation can increase gene flow and restore levels of genetic variation in a few decades. No studies have so far investigated the genetic structure of red squirrel in large, continuous forests. The Italian Alps are presently characterized by almost continuous, recently reconnected forest habitats, that were affected by deep landscape changes during last glaciations but remained mostly unchanged between 10 000 and 200 years bp, when forest cover was heavily reduced. In this study we analyse patterns of genetic variability of red squirrels in and between seven sites distributed over 250 km of Alpine habitat, using mitochondrial DNA (mtDNA) and microsatellites. We use isolation-by-distance (IBD) models to investigate the relative importance that past (Pleistocene glaciations) and recent (fragmentation, bottlenecks) events had on the present genetic situation. Both nuclear and mtDNA data indicate a significant differentiation among study sites and a significant correlation between genetic and geographical distance only over a large scale. No recent bottlenecks are recorded through microsatellites and demographic models strongly support equilibrium between gene flow and drift; however, mtDNA suggests that there may have been local demographic crashes, probably in correspondence with the 19th-century forest fragmentation. These findings indicate that local landscape factors other than geographical distance per se, such as barriers of unsuitable habitat, affect gene flow and determine differentiation.

Phylogenetic position of the small Kashmir flying squirrel, Hylopetes fimbriatus (≡ Eoglaucomys fimbriatus), in the subfamily Pteromyinae

The phylogenetic relationships of flying squirrels (Pteromyinae) were studied by obtaining complete sequence data from the mitochondrial cytochrome b gene of eight Old World and two New World flying squirrel species, with special reference to the systematic and phylogenetic status among Hylopetes fimbriatus (Gray, 1837) (≡ Eoglaucomys fimbriatus (Gray, 1837)) from Pakistan, two Glaucomys Thomas, 1908 species from North America, and two Hylopetes Thomas, 1908 species from Southeast Asia. Phylogenetic trees supported clustering of (i) Belomys pearsonii (Gray, 1842), (ii) H. fimbriatus, the Glaucomys species, Hylopetes lepidus (Horsfield, 1823), and Hylopetes phayrei (Blyth, 1859), (iii) species of Pteromys G. Cuvier, 1800, and (iv) species of Petaurista Link, 1795. Early polytomic divergence among the flying squirrel genera could have taken place in the northern part of the Eurasian continent. The unclear divergence between the Old and New World flying squirrels shows that divergence among flying squirrel genera could have occurred before the formation of the Bering Strait. Hylopetes fimbriatus was more closely related to the two Glaucomys species than to H. lepidus or H. phayrei, supporting placement of the species fimbriatus in the monotypic genus Eoglaucomys Howell, 1915.

Presence of a “CAGA box” in the APP gene unique to amyloid plaque‐forming species and absent in all APLP ‐1/2 genes: implications in Alzheimer's disease

Potentially toxic amyloid beta-peptide (Abeta) in Alzheimer's disease (AD) is generated from a family of Abeta-containing precursor proteins (APP), which is regulated via the 5'-untranslated region (5'-UTR) of its mRNA. We analyzed 5'-UTRs of the APP superfamily, including amyloid plaque-forming and non-amyloid plaque-forming species, and of prions (27 different DNA sequences). A "CAGA" sequence proximal to the "ATG" start codon was present in a location unique to APP genes of amyloid plaque-forming species and absent in all other genes surveyed. This CAGA box is immediately upstream of an interleukin-1-responsive element (acute box). In addition, the proximal CAGA box is predicted to appear on a stem-loop structure in both human and guinea pig APP mRNA. This stem-loop is part of a predicted bulge-loop that encompasses a known iron regulatory element (IRE). Electrophoretic mobility shift with segments of the APP 5'-UTR showed that a region with the proximal CAGA sequence binds nuclear proteins, and this UTR fragment is active in a reporter gene functional assay. Thus, the 5'-UTR in the human APP but not those of APP-like proteins contains a specific region that may participate in APP regulation and may determine a more general model for amyloid generation as seen in AD. The 5'-UTR of human APP contains several interesting control elements, such as an acute box element, a CAGA box, an IRE, and a transforming growth factor-beta-responsive element, that could control APP expression and provide suitable and specific drug targets for AD.

Mapping of spinal cord circuits controlling the bladder and external urethral sphincter functions in the rabbit

AIMS

A primary purpose of this study was to evaluate the rabbit as a model for studying the spinal circuitry controlling the bladder emptying. We aimed to map the locations of the neuronal circuitry controlling the external urethral sphincter (EUS) and the detrusor by stimulating at different spinal cord locations with a microelectrode, while recording the responses from these muscles.

METHODS

Spinal cord microstimulation was performed in the intermediate zone of the gray matter at the L7-S4 spinal cord levels in eight rabbits with empty and full bladders. Bladder activity was measured as intravesical pressure (IVP) changes and EUS activity was measured via electromyographic (EMG) electrodes positioned within the urethra.

RESULTS

Under both bladder conditions, EUS activation was produced from similar locations in the spinal cord comprising a continuous area in the intermediate zone of the S2-S3 spinal cord. This region extended 25 mm in the rostrocaudal dimension, at least 1 mm lateral to the midline, and 0.5-1 mm in the dorsoventral dimension at a depth of 2-3 mm beneath the dorsal surface. No locations in the intermediate zone produced EUS inhibition. The S2-S3 spinal region, stimulation of which produced the strongest EUS activation, also produced modest bladder contractions.

CONCLUSIONS

Overall, the results indicate that spinal cord networks controlling bladder and EUS activation in the rabbit are overlapping and clustered into columns extending rostrocaudally. The lack of spinal locations producing EUS inhibition and large bladder contractions make the rabbit an unattractive model for studies of neuroprosthetic spinal control of micturition.

A memetic-aided approach to hierarchical clustering from distance matrices: application to gene expression clustering and phylogeny

We propose a heuristic approach to hierarchical clustering from distance matrices based on the use of memetic algorithms (MAs). By using MAs to solve some variants of the Minimum Weight Hamiltonian Path Problem on the input matrix, a sequence of the individual elements to be clustered (referred to as patterns) is first obtained. While this problem is also NP-hard, a probably optimal sequence is easy to find with the current advances for this problem and helps to prune the space of possible solutions and/or to guide the search performed by an actual clustering algorithm. This technique has been successfully applied to both a Branch-and-Bound algorithm, and to evolutionary algorithms and MAs. Experimental results are given in the context of phylogenetic inference and in the hierarchical clustering of gene expression data.

Receptive field properties and laminar organization of lateral geniculate nucleus in the gray squirrel (Sciurus carolinensis)

Physiological studies of the lateral geniculate nucleus (LGN) have revealed three classes of relay neurons, called X, Y, and W cells in carnivores and parvocellular (P), magnocellular (M), and koniocellular (K) in primates. The homological relationships among these cell classes and how receptive field (RF) properties of these cells compare with LGN cells in other mammals are poorly understood. To address these questions, we have characterized RF properties and laminar organization in LGN of a highly visual diurnal rodent, the gray squirrel, under isoflurane anesthesia. We identified three classes of LGN cells. One class found in layers 1 and 2 showed sustained, reliable firing, center-surround organization, and was almost exclusively linear in spatial summation. Another class, found in layer 3, showed short response latencies, transient and reliable firing, center-surround organization, and could show either linear (76%) or nonlinear (24%) spatial summation. A third, heterogeneous class found throughout the LGN but primarily in layer 3 showed highly variable responses, a variety of response latencies and could show either center-surround or noncenter-surround receptive field organization and either linear (77%) or nonlinear (23%) spatial summation. RF sizes of all cell classes showed little dependency on eccentricity, and all of these classes showed low contrast gains. When compared with LGN cells in other mammals, our data are consistent with the idea that all mammals contain three basic classes of LGN neurons, one showing reliable, sustained responses, and center-surround organization (X or P); another showing transient but reliable responses, short latencies, and center-surround organization (Y or M); and a third, highly variable and heterogeneous class of cells (W or K). Other properties such as dependency of receptive field size on eccentricity, linearity of spatial summation, and contrast gain appear to vary from species to species.

Mitochondrial phylogeny of hedgehogs and monophyly of Eulipotyphla

We sequenced the complete mitochondrial (mt) genomes of three insectivores: the long-eared hedgehog Hemiechinus auritus, the Japanese mole Mogera wogura, and the greater Japanese shrew-mole Urotrichus talpoides. These mtDNA data together with other previously sequenced mtDNAs were analyzed using a maximum likelihood method to infer their phylogenetic relationships among eutherians. Previous mitochondrial protein analyses used a simple model that did not consider site-heterogeneity, and Erinaceoidea (hedgehogs and moonrats) was placed at the basal eutherian position that is separated from Soricoidea (shrews) and Talpoidea (moles), suggesting the exclusion of the Erinaceoidea-Eulipotyphla tree. By including the new mtDNA sequences and introducing site-heterogeneity into the model, the Erinaceoidea-Eulipotyphla tree emerges as the best tree or as a tree with a log-likelihood score indistinguishable from that of the best tree. However, this conclusion depends on species sampling in Erinaceoidea, demonstrating the importance of both species sampling and use of an appropriate substitution model when inferring phylogenetic relationships.