Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Shrinkage-based Random Local Clocks with Scalable Inference

Molecular clock models undergird modern methods of divergence-time estimation. Local clock models propose that the rate of molecular evolution is constant within phylogenetic subtrees. Current local clock inference procedures exhibit one or more weaknesses, namely they achieve limited scalability to trees with large numbers of taxa, impose model misspecification, or require a priori knowledge of the existence and location of clocks. To overcome these challenges, we present an autocorrelated, Bayesian model of heritable clock rate evolution that leverages heavy-tailed priors with mean zero to shrink increments of change between branch-specific clocks. We further develop an efficient Hamiltonian Monte Carlo sampler that exploits closed form gradient computations to scale our model to large trees. Inference under our shrinkage clock exhibits a speed-up compared to the popular random local clock when estimating branch-specific clock rates on a variety of simulated datasets. This speed-up increases with the size of the problem. We further show our shrinkage clock recovers known local clocks within a rodent and mammalian phylogeny. Finally, in a problem that once appeared computationally impractical, we investigate the heritable clock structure of various surface glycoproteins of influenza A virus in the absence of prior knowledge about clock placement. We implement our shrinkage clock and make it publicly available in the BEAST software package.

Adaptive dating and fast proposals: Revisiting the phylogenetic relaxed clock model

Relaxed clock models enable estimation of molecular substitution rates across lineages and are widely used in phylogenetics for dating evolutionary divergence times. Under the (uncorrelated) relaxed clock model, tree branches are associated with molecular substitution rates which are independently and identically distributed. In this article we delved into the internal complexities of the relaxed clock model in order to develop efficient MCMC operators for Bayesian phylogenetic inference. We compared three substitution rate parameterisations, introduced an adaptive operator which learns the weights of other operators during MCMC, and we explored how relaxed clock model estimation can benefit from two cutting-edge proposal kernels: the AVMVN and Bactrian kernels. This work has produced an operator scheme that is up to 65 times more efficient at exploring continuous relaxed clock parameters compared with previous setups, depending on the dataset. Finally, we explored variants of the standard narrow exchange operator which are specifically designed for the relaxed clock model. In the most extreme case, this new operator traversed tree space 40% more efficiently than narrow exchange. The methodologies introduced are adaptive and highly effective on short as well as long alignments. The results are available via the open source optimised relaxed clock (ORC) package for BEAST 2 under a GNU licence (https://github.com/jordandouglas/ORC).

Effect of Dinotefuran, Permethrin, and Pyriproxyfen (Vectra® 3D) on the Foraging and Blood-Feeding Behaviors of Aedes albopictus Using Laboratory Rodent Model

Simple Summary

Tiger mosquito (Aedes albopictus) is a harmful vector involved in the transmission of several diseases to humans and their pets. Currently, several veterinary products are used to prevent pets against bites of arthropod vectors. However, there is no available information on the effect of these products on feeding and host choice behaviors of Aedes albopictus in the presence of treated and untreated hosts, as is the case of treated dogs present in close physical contact with their owners. The present study investigated the effect of a spot-on product (Vectra^® 3D) on the feeding and host choice behaviours of Aedes albopictus when treated and untreated hosts are presents. Laboratory rodent model was performed to simulate the natural conditions. Rat and mouse hosts were alternately treated with Vectra^® 3D and exposed simultaneously to starved mosquitoes. Results showed that Vectra^® 3D-treated hosts are perfectly protected against up to 82% of mosquitoes. While up to 21% of mosquitoes were repelled from untreated hosts when these latter are present in close physical contact (30 cm) with treated ones suggesting an indirect protection that can allowed the protection of owners who treat their pets with Vectra^® 3D.

Abstract

Dinotefuran-Permethrin-Pyriproxyfen (DPP) is used to kill and repel mosquitoes from dogs. However, the influence of the product on the host-seeking behavior of mosquitoes remains unknown. The interference of DPP with the host selection of unfed female Aedes albopictus was investigated. A total of 18 animals (9 mice and 9 rats) were divided into three groups of six animals each. DU: DPP treated rats (n = 3) with untreated mice (n = 3), UD: DPP treated mice (n = 3) with untreated rats (n = 3) and control UU: untreated mice (n = 3) and untreated rats (n = 3). In each group, the rats and mice were placed 30 cm apart. After sedation, the animals in each group were exposed twice (Day 1 and Day 7 post-treatment) for one hour to 71 ± 3 female mosquitoes. Mosquitoes were categorized after the 2-h post-exposure period as dead or alive. Blood-meal origin was determined from mosquitoes using a newly customized duplex qPCR. The highest values of forage ratio (1.36 ≥ wi ≤ 1.88) and selection index (0.63 ≥ Bi ≤ 0.94) for rat hosts indicates a preference of mosquitoes for this species as compared to mice when co-housed during the exposure. The mosquitoes only seldom fed on mice, even in the untreated group. The anti-feeding effect of DPP was therefore only assessed on rat’s hosts. The results showed that DPP, when directly applied on rats, provided a direct protection of 82% and 61% on Day 1 and Day 7, respectively, while when applied on mice hosts (UD), the DPP provided an indirect protection of 21% and 10% on Day 1 and Day 7, respectively. The results showed also that DPP, when applied on rats, provided a direct protection against Ae. albopictus bites. This effect did not result in increased exposure of the untreated host placed in the same cage at a distance of 30 cm.

Rapid host expansion of an introduced parasite, the spiny rat louse Polyplax spinulosa (Psocodea: Phthiraptera: Polyplacidae), among endemic rodents in Australia

BACKGROUND

Historical European exploration and colonization resulted in the introduction of four species of rodents to the Australian continent from Eurasia: the brown rat, Rattus norvegicus, the black rat, R. rattus, the Pacific rat, R. exulans, and the house mouse, Mus musculus. The spread of these rodents created opportunities for their co-introduced sucking lice to parasitize and adapt to endemic rodents in Australia.

METHODS

We collected sucking lice from rodent specimens in seven museums across Australia. We identified the spiny rat louse, Polyplax spinulosa, based on morphology. We sequenced the mitochondrial cox1 and rrnL genes of P. spinulosa specimens and constructed a phylogenetic tree with rrnL sequences.

RESULTS

We examined 989 rodent specimens of 54 species and collected 2111 adult sucking lice and 1064 nymphal sucking lice. We found that P. spinulosa had nearly doubled its host range by parasitizing at least six endemic rodent species in Australia. The other two introduced lice, P. serrata and Hoplopleura pacifica, however, have apparently failed to expand to any endemic rodents in Australia. Our analysis of mitochondrial rrnL gene sequences divided P. spinulosa into two genotypes (European vs Southeast Asian), which differ by 7.5%; both genotypes were introduced into Australia and then expanded their host ranges to include endemic rodents.

CONCLUSIONS

The earliest record of a European ship landing in Australia was in 1606, followed by British settlement in 1788. The expansion of P. spinulosa to at least six endemic rodent species in Australia has therefore occurred in the time frame of 200 to 400 years, which is extremely rapid relative to its host expansion to eight native rat species in Eurasia in ~ 16 million years since it diverged from P. serrata. The host expansion of P. spinulosa is remarkable for a blood-sucking louse and is in stark contrast to the absence of host expansion by P. serrata and H. pacifica. Comparison among these three introduced sucking lice indicated that both louse-specific factors and host-specific factors can contribute to the success or failure of host expansion.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

Fast Dating Using Least-Squares Criteria and Algorithms

Phylogenies provide a useful way to understand the evolutionary history of genetic samples, and data sets with more than a thousand taxa are becoming increasingly common, notably with viruses (e.g., human immunodeficiency virus (HIV)). Dating ancestral events is one of the first, essential goals with such data. However, current sophisticated probabilistic approaches struggle to handle data sets of this size. Here, we present very fast dating algorithms, based on a Gaussian model closely related to the Langley–Fitch molecular-clock model. We show that this model is robust to uncorrelated violations of the molecular clock. Our algorithms apply to serial data, where the tips of the tree have been sampled through times. They estimate the substitution rate and the dates of all ancestral nodes. When the input tree is unrooted, they can provide an estimate for the root position, thus representing a new, practical alternative to the standard rooting methods (e.g., midpoint). Our algorithms exploit the tree (recursive) structure of the problem at hand, and the close relationships between least-squares and linear algebra. We distinguish between an unconstrained setting and the case where the temporal precedence constraint (i.e., an ancestral node must be older that its daughter nodes) is accounted for. With rooted trees, the former is solved using linear algebra in linear computing time (i.e., proportional to the number of taxa), while the resolution of the latter, constrained setting, is based on an active-set method that runs in nearly linear time. With unrooted trees the computing time becomes (nearly) quadratic (i.e., proportional to the square of the number of taxa). In all cases, very large input trees (>10,000 taxa) can easily be processed and transformed into time-scaled trees. We compare these algorithms to standard methods (root-to-tip, r8s version of Langley–Fitch method, and BEAST). Using simulated data, we show that their estimation accuracy is similar to that of the most sophisticated methods, while their computing time is much faster. We apply these algorithms on a large data set comprising 1194 strains of Influenza virus from the pdm09 H1N1 Human pandemic. Again the results show that these algorithms provide a very fast alternative with results similar to those of other computer programs. These algorithms are implemented in the LSD software (least-squares dating), which can be downloaded from http://www.atgc-montpellier.fr/LSD/, along with all our data sets and detailed results. An Online Appendix, providing additional algorithm descriptions, tables, and figures can be found in the Supplementary Material available on Dryad at http://dx.doi.org/10.5061/dryad.968t3.

Global pattern of phylogenetic species composition of shark and its conservation priority

The diversity of marine communities is in striking contrast with the diversity of terrestrial communities. In all oceans, species richness is low in tropical areas and high at latitudes between 20 and 40°. While species richness is a primary metric used in conservation and management strategies, it is important to take into account the complex phylogenetic patterns of species compositions within communities. We measured the phylogenetic skew and diversity of shark communities throughout the world. We found that shark communities in tropical seas were highly phylogenetically skewed, whereas temperate sea communities had phylogenetically diversified species compositions. Interestingly, although geographically distant from one another, tropical sea communities were all highly skewed toward requiem sharks (Carcharhinidae), hammerhead sharks (Sphyrnidae), and whale sharks (Rhincodon typus). Worldwide, the greatest phylogenetic evenness in terms of clades was found in the North Sea and coastal regions of countries in temperate zones, such as the United Kingdom, Ireland, southern Australia, and Chile. This study is the first to examine patterns of phylogenetic diversity of shark communities on a global scale. Our findings suggest that when establishing conservation activities, it is important to take full account of phylogenetic patterns of species composition and not solely use species richness as a target. Protecting areas of high phylogenetic diversity in sharks, which were identified in this study, could form a broader strategy for protecting other threatened marine species.

Flying rocks and flying clocks: disparity in fossil and molecular dates for birds

Major disparities are recognized between molecular divergence dates and fossil ages for critical nodes in the Tree of Life, but broad patterns and underlying drivers remain elusive. We harvested 458 molecular age estimates for the stem and crown divergences of 67 avian clades to explore empirical patterns between these alternate sources of temporal information. These divergence estimates were, on average, over twice the age of the oldest fossil in these clades. Mitochondrial studies yielded older ages than nuclear studies for the vast majority of clades. Unexpectedly, disparity between molecular estimates and the fossil record was higher for divergences within major clades (crown divergences) than divergences between major clades (stem divergences). Comparisons of dates from studies classed by analytical methods revealed few significant differences. Because true divergence ages can never be known with certainty, our study does not answer the question of whether fossil gaps or molecular dating error account for a greater proportion of observed disparity. However, empirical patterns observed here suggest systemic overestimates for shallow nodes in existing molecular divergence dates for birds. We discuss underlying biases that may drive these patterns.

Phylogenetic skew: an index of community diversity

The distribution of divergence times between member species of a community reflects the pattern of species composition. In this study, we contrast the species composition of a community against the meta-community, which we define as the species composition of a set of target communities. We regard the collection of species that comprise a community as a sample from the set of member species of the meta-community, and interpret the pattern of the community species composition in terms of the type of species sampled from the meta-community. A newly defined effective species sampling proportion explains the amount of the difference between the divergence time distributions of the community and that of the meta-community, assuming random sampling. We propose a new index of phylogenetic skew (PS), as the ratio of the maximum-likelihood estimate of the effective species sampling proportion to the observed sampling proportion. A PS value of 1 is interpreted as random sampling. If the value is >1, the sampling is suspected to be phylogenetically skewed. If it is <1, systematic thinning of species is likely. Unlike other indices, the PS does not depend on species richness as long as the community has more than a few members of a species. Because it is possible to compare partially observed communities, the index may be effectively used in exploratory analysis to detect candidate communities with unique species compositions from a large number of communities.

Ocular comparative anatomy of the family Rodentia

There is little information regarding ocular anatomy and histology in many of the rodent species. Histological analyses for morphologic features were performed in 31 globes from 18 rodent species submitted to and archived at the Comparative Ocular Pathology Laboratory of Wisconsin. The following measurements were taken: thickness of the cornea, corneal epithelium, corneal stroma, Descemet's membrane, and retina. H&E sections were evaluated for the following anatomical features: presence of pigmented epithelial cells in the peripheral cornea, presence and location of Schlemm's canal, presence of iridal sphincter and dilator and ciliary body muscles, presence of pars plicata and plana, presence of retinal vessels, presence of lamina cribrosa, and presence of tapetum lucidum. The springhaas was the only rodent in our collection that presented a well-developed tapetum lucidum fibrosum. The presence of retinal vessels was variable: vessels were observed in all of the members of the mouse-related clade, except the springhaas and the beaver, in all of the squirrel-related clade members, and in none of the Ctenohystrica. In the flying squirrels, blood vessels extended to the outer limiting membrane in the photoreceptor layer. Beavers, chinchillas, capybara, and guinea pigs lacked vessels within the retina; however, they had vessels within the optic nerve head. Ground squirrels have an optic nerve head, which is linear in the horizontal plane and an asymmetric retina. The tree-dwelling squirrels have a rounded but still elongated optic nerve, and the flying squirrel has a round optic nerve head like all the other rodents.

Accelerated evolution after gene duplication: a time-dependent process affecting just one copy

Gene duplication is widely regarded as a major mechanism modeling genome evolution and function. However, the mechanisms that drive the evolution of the two, initially redundant, gene copies are still ill defined. Many gene duplicates experience evolutionary rate acceleration, but the relative contribution of positive selection and random drift to the retention and subsequent evolution of gene duplicates, and for how long the molecular clock may be distorted by these processes, remains unclear. Focusing on rodent genes that duplicated before and after the mouse and rat split, we find significantly increased sequence divergence after duplication in only one of the copies, which in nearly all cases corresponds to the novel daughter copy, independent of the mechanism of duplication. We observe that the evolutionary rate of the accelerated copy, measured as the ratio of nonsynonymous to synonymous substitutions, is on average 5-fold higher in the period spanning 4-12 My after the duplication than it was before the duplication. This increase can be explained, at least in part, by the action of positive selection according to the results of the maximum likelihood-based branch-site test. Subsequently, the rate decelerates until purifying selection completely returns to preduplication levels. Reversion to the original rates has already been accomplished 40.5 My after the duplication event, corresponding to a genetic distance of about 0.28 synonymous substitutions per site. Differences in tissue gene expression patterns parallel those of substitution rates, reinforcing the role of neofunctionalization in explaining the evolution of young gene duplicates.

Characterization of thymosin β4 in mammals' saliva

Thymosin β4 (Tβ4) is a low molecular weight peptide found in several mammalian tissues and is known mainly by its ability to bind cytoskeletal actin, influencing cell migration and differentiation, and promoting tissue repair. Considering the functional role of this peptide, the main goal of this work was to characterize Tβ4 in mammals' saliva by using evolutionary and proteomic tools. For this, mammalian Tβ4 sequences were retrieved from NCBI, SwissProt and Ensembl databases. The alignment of Tβ4 amino acid sequences showed a high degree of conservation between species. The gene seems to be evolving under negative selection as indicated by a dN/dS ratio of 0.05. Whole saliva was collected from dog, human, rabbit, cow, horse and sheep and the salivary peptides were isolated through filtration and analyzed by LC-MS/MS. Spectra was processed against the database constructed with the retrieved Tβ4 sequences. For the first time, the identification of this peptide was achieved in rat, dog, horse and bovine saliva. Detection in these mammal species and its amino acid conservation suggest an important role of Tβ4 in the homeostasis of the mammalian oral cavity.

Molecular and paleontological evidence for a post-Cretaceous origin of rodents

The timing of the origin and diversification of rodents remains controversial, due to conflicting results from molecular clocks and paleontological data. The fossil record tends to support an early Cenozoic origin of crown-group rodents. In contrast, most molecular studies place the origin and initial diversification of crown-Rodentia deep in the Cretaceous, although some molecular analyses have recovered estimated divergence times that are more compatible with the fossil record. Here we attempt to resolve this conflict by carrying out a molecular clock investigation based on a nine-gene sequence dataset and a novel set of seven fossil constraints, including two new rodent records (the earliest known representatives of Cardiocraniinae and Dipodinae). Our results indicate that rodents originated around 61.7–62.4 Ma, shortly after the Cretaceous/Paleogene (K/Pg) boundary, and diversified at the intraordinal level around 57.7–58.9 Ma. These estimates are broadly consistent with the paleontological record, but challenge previous molecular studies that place the origin and early diversification of rodents in the Cretaceous. This study demonstrates that, with reliable fossil constraints, the incompatibility between paleontological and molecular estimates of rodent divergence times can be eliminated using currently available tools and genetic markers. Similar conflicts between molecular and paleontological evidence bedevil attempts to establish the origination times of other placental groups. The example of the present study suggests that more reliable fossil calibration points may represent the key to resolving these controversies.

The diversity of class II transposable elements in mammalian genomes has arisen from ancestral phylogenetic splits during ancient waves of proliferation through the genome

DNA transposons make up 3% of the human genome, approximately the same percentage as genes. However, because of their inactivity, they are often ignored in favor of the more abundant, active, retroelements. Despite this relative ignominy, there are a number of interesting questions to be asked of these transposon families. One particular question relates to the timing of proliferation and inactivation of elements in a family. Does an ongoing process of turnover occur, or is the process more akin to a life cycle for the family, with elements proliferating rapidly before deactivation at a later date? We answer this question by tracing back to the most recent common ancestor (MRCA) of each modern transposon family, using two different methods. The first method identifies the MRCA of the species in which a family of transposon fossils can still be found, which we assume will have existed soon after the true origin date of the transposon family. The second method uses molecular dating techniques to predict the age of the MRCA element from which all elements found in a modern genome are descended. Independent data from five pairs of species are used in the molecular dating analysis: human–chimpanzee, human–orangutan, dog–panda, dog–cat, and cow–pig. Orthologous pairs of elements from host species pairs are included, and the divergence dates of these species are used to constrain the analysis. We discover that, in general, the times to element common ancestry for a given family are the same for the different species pairs, suggesting that there has been no order-specific process of turnover. Furthermore, for most families, the ages of the common ancestor of the host species and of that of the elements are similar, suggesting a life cycle model for the proliferation of transposons. Where these two ages differ, in families found only in Primates and Rodentia, for example, we find that the host species date is later than that of the common ancestor of the elements, implying that there may be large deletions of elements from host species, examples of which were found in their ancestors.

Rapid turnover of long noncoding RNAs and the evolution of gene expression

A large proportion of functional sequence within mammalian genomes falls outside protein-coding exons and can be transcribed into long RNAs. However, the roles in mammalian biology of long noncoding RNA (lncRNA) are not well understood. Few lncRNAs have experimentally determined roles, with some of these being lineage-specific. Determining the extent by which transcription of lncRNA loci is retained or lost across multiple evolutionary lineages is essential if we are to understand their contribution to mammalian biology and to lineage-specific traits. Here, we experimentally investigated the conservation of lncRNA expression among closely related rodent species, allowing the evolution of DNA sequence to be uncoupled from evolution of transcript expression. We generated total RNA (RNAseq) and H3K4me3-bound (ChIPseq) DNA data, and combined both to construct catalogues of transcripts expressed in the adult liver of Mus musculus domesticus (C57BL/6J), Mus musculus castaneus, and Rattus norvegicus. We estimated the rate of transcriptional turnover of lncRNAs and investigated the effects of their lineage-specific birth or death. LncRNA transcription showed considerably greater gain and loss during rodent evolution, compared with protein-coding genes. Nucleotide substitution rates were found to mirror the in vivo transcriptional conservation of intergenic lncRNAs between rodents: only the sequences of noncoding loci with conserved transcription were constrained. Finally, we found that lineage-specific intergenic lncRNAs appear to be associated with modestly elevated expression of genomically neighbouring protein-coding genes. Our findings show that nearly half of intergenic lncRNA loci have been gained or lost since the last common ancestor of mouse and rat, and they predict that such rapid transcriptional turnover contributes to the evolution of tissue- and lineage-specific gene expression.

The best-understood portion of mammalian genomes contains genes transcribed into RNAs, which are subsequently translated into proteins. These genes are generally under high selective pressure and deeply conserved between species. Recent publications have revealed novel classes of genes, which are also transcribed into RNA but are not subsequently translated into proteins. One such novel class are long noncoding RNA (lncRNA). LncRNA loci are controlled in a similar manner to protein-coding genes, yet are more often expressed tissue-specifically, and their conservation and function(s) are mostly unknown. Previous reports suggest that lncRNAs can affect the expression of nearby protein-coding genes or act at a distance to control broader biological processes. Also, lncRNA sequence is poorly conserved between mammals compared with protein-coding genes, but how rapidly their transcription evolves, particularly between closely related species, remains unknown. By comparing lncRNA expression between homologous tissues in two species of mouse and in rat, we discovered that lncRNA genes are “born” or “die” more rapidly than protein-coding genes and that this rapid evolution impacts the expression levels of nearby coding genes. This local regulation of gene expression reveals a functional role for the rapid evolution of lncRNAs, which may contribute to biological differences between species.

A glimpse on the pattern of rodent diversification: a phylogenetic approach

Background

Development of phylogenetic methods that do not rely on fossils for the study of evolutionary processes through time have revolutionized the field of evolutionary biology and resulted in an unprecedented expansion of our knowledge about the tree of life. These methods have helped to shed light on the macroevolution of many taxonomic groups such as the placentals (Mammalia). However, despite the increase of studies addressing the diversification patterns of organisms, no synthesis has addressed the case of the most diversified mammalian clade: the Rodentia.

Results

Here we present a rodent maximum likelihood phylogeny inferred from a molecular supermatrix. It is based on 11 mitochondrial and nuclear genes that covers 1,265 species, i.e., respectively 56% and 81% of the known specific and generic rodent diversity. The inferred topology recovered all Rodentia clades proposed by recent molecular works. A relaxed molecular clock dating approach provided a time framework for speciation events. We found that the Myomorpha clade shows a greater degree of variation in diversification rates than Sciuroidea, Caviomorpha, Castorimorpha and Anomaluromorpha. We identified a number of shifts in diversification rates within the major clades: two in Castorimorpha, three in Ctenohystrica, 6 within the squirrel-related clade and 24 in the Myomorpha clade. The majority of these shifts occurred within the most recent familial rodent radiations: the Cricetidae and Muridae clades. Using the topological imbalances and the time line we discuss the potential role of different diversification factors that might have shaped the rodents radiation.

Conclusions

The present glimpse on the diversification pattern of rodents can be used for further comparative meta-analyses. Muroid lineages have a greater degree of variation in their diversification rates than any other rodent group. Different topological signatures suggest distinct diversification processes among rodent lineages. In particular, Muroidea and Sciuroidea display widespread distribution and have undergone evolutionary and adaptive radiation on most of the continents. Our results show that rodents experienced shifts in diversification rate regularly through the Tertiary, but at different periods for each clade. A comparison between the rodent fossil record and our results suggest that extinction led to the loss of diversification signal for most of the Paleogene nodes.

Joint reconstruction of divergence times and life-history evolution in placental mammals using a phylogenetic covariance model

Violation of the molecular clock has been amply documented, and is now routinely taken into account by molecular dating methods. Comparative analyses have revealed a systematic component in rate variation, relating it to the evolution of life-history traits, such as body size or generation time. Life-history evolution can be reconstructed using Brownian models. However, the resulting estimates are typically uncertain, and potentially sensitive to the underlying assumptions. As a way of obtaining more accurate ancestral trait and divergence time reconstructions, correlations between life-history traits and substitution rates could be used as an additional source of information. In this direction, a Bayesian framework for jointly reconstructing rates, traits, and dates was previously introduced. Here, we apply this model to a 17 protein-coding gene alignment for 73 placental taxa. Our analysis indicates that the coupling between molecules and life history can lead to a reevaluation of ancestral life-history profiles, in particular for groups displaying convergent evolution in body size. However, reconstructions are sensitive to fossil calibrations and to the Brownian assumption. Altogether, our analysis suggests that further integrating inference of rates and traits might be particularly useful for neontological macroevolutionary comparative studies.

Investigation of the origin and spread of a Mammalian transposable element based on current sequence diversity

Almost half the human genome consists of mobile DNA elements, and their analysis is a vital part of understanding the human genome as a whole. Many of these elements are ancient and have persisted in the genome for tens or hundreds of millions of years, providing a window into the evolution of modern mammals. The Golem family have been used as model transposons to highlight computational analyses which can be used to investigate these elements, particularly the use of molecular dating with large transposon families. Whole-genome searches found Golem sequences in 20 mammalian species. Golem A and B subsequences were only found in primates and squirrel. Interestingly, the full-length Golem, found as a few copies in many mammalian genomes, was found abundantly in horse. A phylogenetic profile suggested that Golem originated after the eutherian–metatherian divergence and that the A and B subfamilies originated at a much later date. Molecular dating based on sequence diversity suggests an early age, of 175 Mya, for the origin of the family and that the A and B lineages originated much earlier than expected from their current taxonomic distribution and have subsequently been lost in some lineages. Using publically available data, it is possible to investigate the evolutionary history of transposon families. Determining in which organisms a transposon can be found is often used to date the origin and expansion of the families. However, in this analysis, molecular dating, commonly used for determining the age of gene sequences, has been used, reducing the likelihood of errors from deleted lineages.

Updated neuronal scaling rules for the brains of Glires (rodents/lagomorphs)

Brain size scales as different functions of its number of neurons across mammalian orders such as rodents, primates, and insectivores. In rodents, we have previously shown that, across a sample of 6 species, from mouse to capybara, the cerebral cortex, cerebellum and the remaining brain structures increase in size faster than they gain neurons, with an accompanying decrease in neuronal density in these structures [Herculano-Houzel et al.: Proc Natl Acad Sci USA 2006;103:12138-12143]. Important remaining questions are whether such neuronal scaling rules within an order apply equally to all pertaining species, and whether they extend to closely related taxa. Here, we examine whether 4 other species of Rodentia, as well as the closely related rabbit (Lagomorpha), conform to the scaling rules identified previously for rodents. We report the updated neuronal scaling rules obtained for the average values of each species in a way that is directly comparable to the scaling rules that apply to primates [Gabi et al.: Brain Behav Evol 2010;76:32-44], and examine whether the scaling relationships are affected when phylogenetic relatedness in the dataset is accounted for. We have found that the brains of the spiny rat, squirrel, prairie dog and rabbit conform to the neuronal scaling rules that apply to the previous sample of rodents. The conformity to the previous rules of the new set of species, which includes the rabbit, suggests that the cellular scaling rules we have identified apply to rodents in general, and probably to Glires as a whole (rodents/lagomorphs), with one notable exception: the naked mole-rat brain is apparently an outlier, with only about half of the neurons expected from its brain size in its cerebral cortex and cerebellum.

The cecal appendix: one more immune component with a function disturbed by post-industrial culture

This review assesses the current state of knowledge regarding the cecal appendix, its apparent function, and its evolution. The association of the cecal appendix with substantial amounts of immune tissue has long been taken as an indicator that the appendix may have some immune function. Recently, an improved understanding of the interactions between the normal gut flora and the immune system has led to the identification of the appendix as an apparent safe-house for normal gut bacteria. Further, a variety of observations related to the evolution and morphology of the appendix, including the identification of the structure as a "recurrent trait" in some clades, the presence of appendix-like structures in monotremes and some non-mammalian species, and consistent features of the cecal appendix such as its narrow diameter, provide direct support for an important function of the appendix. This bacterial safe-house, which is likely important in the event of diarrheal illness, is presumably of minimal importance to humans living with abundant nutritional resources, modern medicine and modern hygiene practices that include clean drinking water. Consistent with this idea, epidemiologic studies demonstrate that diarrheal illness is indeed a major source of selection pressure in developing countries but not in developed countries, whereas appendicitis shows the opposite trend, being associated with modern hygiene and medicine. The cecal appendix may thus be viewed as a part of the immune system that, like those immune compartments that cause allergy, is vital to life in a "natural" environment, but which is poorly suited to post-industrialized societies.

Bayesian random local clocks, or one rate to rule them all

Background

Relaxed molecular clock models allow divergence time dating and "relaxed phylogenetic" inference, in which a time tree is estimated in the face of unequal rates across lineages. We present a new method for relaxing the assumption of a strict molecular clock using Markov chain Monte Carlo to implement Bayesian modeling averaging over random local molecular clocks. The new method approaches the problem of rate variation among lineages by proposing a series of local molecular clocks, each extending over a subregion of the full phylogeny. Each branch in a phylogeny (subtending a clade) is a possible location for a change of rate from one local clock to a new one. Thus, including both the global molecular clock and the unconstrained model results, there are a total of 2^2n-2 possible rate models available for averaging with 1, 2, ..., 2n - 2 different rate categories.

Results

We propose an efficient method to sample this model space while simultaneously estimating the phylogeny. The new method conveniently allows a direct test of the strict molecular clock, in which one rate rules them all, against a large array of alternative local molecular clock models. We illustrate the method's utility on three example data sets involving mammal, primate and influenza evolution. Finally, we explore methods to visualize the complex posterior distribution that results from inference under such models.

Conclusions

The examples suggest that large sequence datasets may only require a small number of local molecular clocks to reconcile their branch lengths with a time scale. All of the analyses described here are implemented in the open access software package BEAST 1.5.4 (http://beast-mcmc.googlecode.com/).

Massive turnover of functional sequence in human and other mammalian genomes

Despite the availability of dozens of animal genome sequences, two key questions remain unanswered: First, what fraction of any species' genome confers biological function, and second, are apparent differences in organismal complexity reflected in an objective measure of genomic complexity? Here, we address both questions by applying, across the mammalian phylogeny, an evolutionary model that estimates the amount of functional DNA that is shared between two species' genomes. Our main findings are, first, that as the divergence between mammalian species increases, the predicted amount of pairwise shared functional sequence drops off dramatically. We show by simulations that this is not an artifact of the method, but rather indicates that functional (and mostly noncoding) sequence is turning over at a very high rate. We estimate that between 200 and 300 Mb (∼6.5%-10%) of the human genome is under functional constraint, which includes five to eight times as many constrained noncoding bases than bases that code for protein. In contrast, in D. melanogaster we estimate only 56-66 Mb to be constrained, implying a ratio of noncoding to coding constrained bases of about 2. This suggests that, rather than genome size or protein-coding gene complement, it is the number of functional bases that might best mirror our naïve preconceptions of organismal complexity.

An application of supertree methods to Mammalian mitogenomic sequences

TWO DIFFERENT APPROACHES CAN BE USED IN PHYLOGENOMICS: combined or separate analysis. In the first approach, different datasets are combined in a concatenated supermatrix. In the second, datasets are analyzed separately and the phylogenetic trees are then combined in a supertree. The supertree method is an interesting alternative to avoid missing data, since datasets that are analyzed separately do not need to represent identical taxa. However, the supertree approach and the corresponding consensus methods have been highly criticized for not providing valid phylogenetic hypotheses. In this study, congruence of trees estimated by consensus and supertree approaches were compared to model trees obtained from a combined analysis of complete mitochondrial sequences of 102 species representing 93 mammal families. The consensus methods produced poorly resolved consensus trees and did not perform well, except for the majority rule consensus with compatible groupings. The weighted supertree and matrix representation with parsimony methods performed equally well and were highly congruent with the model trees. The most similar supertree method was the least congruent with the model trees. We conclude that some of the methods tested are worth considering in a phylogenomic context.

Divergence time estimates of mammals from molecular clocks and fossils: relevance of new fossil finds from India

This paper presents a brief review of recent advances in the classification of mammals at higher levels using fossils and molecular clocks. It also discusses latest fossil discoveries from the Cretaceous - Eocene (66-55 m.y.) rocks of India and their relevance to our current understanding of placental mammal origins and diversifications.

Co-evolutionary networks of genes and cellular processes across fungal species

BACKGROUND

The introduction of measures such as evolutionary rate and propensity for gene loss have significantly advanced our knowledge of the evolutionary history and selection forces acting upon individual genes and cellular processes.

RESULTS

We present two new measures, the 'relative evolutionary rate pattern' (rERP), which records the relative evolutionary rates of conserved genes across the different branches of a species' phylogenetic tree, and the 'copy number pattern' (CNP), which quantifies the rate of gene loss of less conserved genes. Together, these measures yield a high-resolution study of the co-evolution of genes in 9 fungal species, spanning 3,540 sets of orthologs. We find that the evolutionary tempo of conserved genes varies in different evolutionary periods. The co-evolution of genes' Gene Ontology categories exhibits a significant correlation with their functional distance in the Gene Ontology hierarchy, but not with their location on chromosomes, showing that cellular functions are a more important driving force in gene co-evolution than their chromosomal proximity. Two fundamental patterns of co-evolution of conserved genes, cooperative and reciprocal, are identified; only genes co-evolving cooperatively functionally back each other up. The co-evolution of conserved and less conserved genes exhibits both commonalities and differences; DNA metabolism is positively correlated with nuclear traffic, transcription processes and vacuolar biology in both analyses.

CONCLUSIONS

Overall, this study charts the first global network view of gene co-evolution in fungi. The future application of the approach presented here to other phylogenetic trees holds much promise in characterizing the forces that shape cellular co-evolution.

Identification of myo-inositol-3-phosphate synthase isoforms: characterization, expression, and putative role of a 16-kDa gamma(c) isoform

Myo-inositol is an important constituent of membrane phospholipids and is a precursor for the phosphoinositide signaling pathway. It is synthesized from glucose 6-phosphate by myo-inositol-3-phosphate synthase (IP synthase), a homotrimer composed of a 68-kDa polypeptide in most mammalian tissues. It is a putative target for mood-stabilizing drugs such as lithium and valproate. Here, we show that the rat gene (Isyna1) encoding this enzyme generates a number of alternatively spliced transcripts in addition to the fully spliced form that encodes the 68-kDa subunit (the alpha isoform). Specifically, we identify a small 16-kDa subunit (the gamma(c) isoform) derived by an intron retention mechanism and provide evidence for its existence in rat tissues. The gamma(c) isoform is highly conserved in mammals, but it lacks the catalytic domain while retaining the NAD(+) binding domain. Both alpha and gamma(c) isoforms are predominantly expressed in many rat tissues and display apparent stoichiometry in purified enzyme preparations. An IP synthase polyclonal antibody not only detects the alpha and gamma(c) isoforms but also several other isoforms in pancreas, intestine, and testis suggesting that the holoenzyme is composed of unique subunits in various tissues. Interestingly, the alpha isoform is not expressed in the intestine. IP synthase activity assays using purified alpha and gamma(c) isoforms indicate that the latter negatively modulates alpha isoform activity, possibly by competing for NAD(+) molecules. Our findings have important ramifications for understanding the mood stabilization process and suggest that inositol biosynthesis is a highly regulated and dynamic process.

Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura)

Cospeciation between hosts and parasites offers a unique opportunity to use information from parasites to infer events in host evolutionary history. Although lice (Insecta: Phthiraptera) are known to cospeciate with their hosts and have frequently served as important markers to infer host evolutionary history, most molecular studies are based on only one or two markers. Resulting phylogenies may, therefore, represent gene histories (rather than species histories), and analyses of multiple molecular markers are needed to increase confidence in the results of phylogenetic analyses. Herein, we phylogenetically examine nine molecular markers in primate sucking lice (Phthiraptera: Anoplura) and we use these markers to estimate divergence times among louse lineages. Individual and combined analyses of these nine markers are, for the most part, congruent, supporting relationships hypothesized in previous studies. Only one marker, the nuclear protein-coding gene Histone 3, has a significantly different tree topology compared to the other markers. The disparate evolutionary history of this marker, however, has no significant effect on topology or nodal support in the combined phylogenetic analyses. Therefore, phylogenetic results from the combined data set likely represent a solid hypothesis of species relationships. Additionally, we find that simultaneous use of multiple markers and calibration points provides the most reliable estimates of louse divergence times, in agreement with previous studies estimating divergences among species. Estimates of phylogenies and divergence times also allow us to verify the results of [Reed, D.L., Light, J.E., Allen, J.M., Kirchman, J.J., 2007. Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 5, 7.]; there was probable contact between gorilla and archaic hominids roughly 3 Ma resulting in a host switch of Pthirus lice from gorillas to archaic hominids. Thus, these results provide further evidence that data from cospeciating organisms can yield important information about the evolutionary history of their hosts.

Early Eocene lagomorph (Mammalia) from Western India and the early diversification of Lagomorpha

We report the oldest known record of Lagomorpha, based on distinctive, small ankle bones (calcaneus and talus) from Early Eocene deposits (Middle Ypresian equivalent, ca 53 Myr ago) of Gujarat, west-central India. The fossils predate the oldest previously known crown lagomorphs by several million years and extend the record of lagomorphs on the Indian subcontinent by 35 Myr. The bones show a mosaic of derived cursorial adaptations found in gracile Leporidae (rabbits and hares) and primitive traits characteristic of extant Ochotonidae (pikas) and more robust leporids. Together with gracile and robust calcanei from the Middle Eocene of Shanghuang, China, also reported here, the Indian fossils suggest that diversification within crown Lagomorpha and possibly divergence of the family Leporidae were already underway in the Early Eocene.

Correlates of substitution rate variation in mammalian protein-coding sequences

BACKGROUND

Rates of molecular evolution in different lineages can vary widely, and some of this variation might be predictable from aspects of species' biology. Investigating such predictable rate variation can help us to understand the causes of molecular evolution, and could also help to improve molecular dating methods. Here we present a comprehensive study of the life history correlates of substitution rate variation across the mammals, comparing results for mitochondrial and nuclear loci, and for synonymous and non-synonymous sites. We use phylogenetic comparative methods, refined to take into account the special nature of substitution rate data. Particular attention is paid to the widespread correlations between the components of mammalian life history, which can complicate the interpretation of results.

RESULTS

We find that mitochondrial synonymous substitution rates, estimated from the 9 longest mitochondrial genes, show strong negative correlations with body mass and with maximum recorded lifespan. But lifespan is the sole variable to remain after multiple regression and model simplification. Nuclear synonymous substitution rates, estimated from 6 genes, show strong negative correlations with body mass and generation time, and a strong positive correlation with fecundity. In contrast to the mitochondrial results, the same trends are evident in rates of nonsynonymous substitution.

CONCLUSION

A substantial proportion of variation in mammalian substitution rates can be explained by aspects of their life history, implying that molecular and life history evolution are closely interlinked in this group. The strength and consistency of the nuclear body mass effect suggests that molecular dating studies may have been systematically misled, but also that methods could be improved by incorporating the finding as a priori information. Mitochondrial synonymous rates also show the body mass effect, but for apparently quite different reasons, and the strength of the relationship with maximum lifespan provides support for the hypothesis that mtDNA damage is causally linked to aging.

Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds

A limiting factor in many molecular dating studies is shortage of reliable calibrations. Current methods for choosing calibrations (e.g. cross-validation) treat them as either correct or incorrect, whereas calibrations probably lie on a continuum from highly accurate to very poor. Bayesian relaxed clock analysis permits inclusion of numerous candidate calibrations as priors: provided most calibrations are reliable, the model appropriate and the data informative, the accuracy of each calibration prior can be evaluated. If a calibration is accurate, then the analysis will support the prior so that the posterior estimate reflects the prior; if a calibration is poor, the posterior will be forced away from the prior. We use this approach to test two fossil dates recently proposed as standard calibrations within vertebrates. The proposed bird-crocodile calibration (approx. 247Myr ago) appears to be accurate, but the proposed bird-lizard calibration (approx. 255Myr ago) is substantially too recent.

Accelerated rate of gene gain and loss in primates

The molecular changes responsible for the evolution of modern humans have primarily been discussed in terms of individual nucleotide substitutions in regulatory or protein coding sequences. However, rates of nucleotide substitution are slowed in primates, and thus humans and chimpanzees are highly similar at the nucleotide level. We find that a third source of molecular evolution, gene gain and loss, is accelerated in primates relative to other mammals. Using a novel method that allows estimation of rate heterogeneity among lineages, we find that the rate of gene turnover in humans is more than 2.5 times faster than in other mammals and may be due to both mutational and selective forces. By reconciling the gene trees for all of the gene families included in the analysis, we are able to independently verify the numbers of inferred duplications. We also use two methods based on the genome assembly of rhesus macaque to further verify our results. Our analyses identify several gene families that have expanded or contracted more rapidly than is expected even after accounting for an overall rate acceleration in primates, including brain-related families that have more than doubled in size in humans. Many of the families showing large expansions also show evidence for positive selection on their nucleotide sequences, suggesting that selection has been important in shaping copy-number differences among mammals. These findings may help explain why humans and chimpanzees show high similarity between orthologous nucleotides yet great morphological and behavioral differences.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Evolution of the CYP2ABFGST gene cluster in rat, and a fine-scale comparison among rodent and primate species

The evolution of gene families can be best understood by studying the modern organization and functions of family members, and by comparing parallel families in different species. In this study, the CYP2ABFGST gene cluster has been characterized in rat and compared to the syntenic clusters in mouse and human, providing an interesting example of gene family evolution. In the rat, 18 loci from six subfamilies have been identified by specifically amplifying and sequencing gene fragments from cloned DNA, and have been exactly placed on chromosome 1. The overall organization of the gene cluster in rat is relatively simple, with genes from each subfamily in tandem, and is more similar to the mouse than to the human cluster. We have reconstructed the probable structure of the CYP2ABFGST cluster in the common ancestor of primates and rodents, and inferred a model of the evolution of this gene cluster in the three species. Numerous nontandem and block duplications, inversions, and translocations have occurred entirely inside the cluster, indicating that pairing between duplicate genes is keeping the rearrangements within the cluster region. The initial tandem duplication of a CYP2 gene in an early mammalian ancestor has made this region particularly subject to such localized rearrangements. Even if duplicated genes do not have a large-scale effect on chromosomal rearrangements, on a local level clustered gene families may have contributed significantly to the genomic complexity of modern mammals.

Cretaceous eutherians and Laurasian origin for placental mammals near the K/T boundary

Estimates of the time of origin for placental mammals from DNA studies span nearly the duration of the Cretaceous period (145 to 65 million years ago), with a maximum of 129 million years ago and a minimum of 78 million years ago. Palaeontologists too are divided on the timing. Some support a deep Cretaceous origin by allying certain middle Cretaceous fossils (97-90 million years old) from Uzbekistan with modern placental lineages, whereas others support the origin of crown group Placentalia near the close of the Cretaceous. This controversy has yet to be addressed by a comprehensive phylogenetic analysis that includes all well-known Cretaceous fossils and a wide sample of morphology among Tertiary and recent placentals. Here we report the discovery of a new well-preserved mammal from the Late Cretaceous of Mongolia and a broad-scale phylogenetic analysis. Our results exclude Cretaceous fossils from Placentalia, place the origin of Placentalia near the Cretaceous/Tertiary (K/T) boundary in Laurasia rather than much earlier within the Cretaceous in the Southern Hemisphere, and place afrotherians and xenarthrans in a nested rather than a basal position within Placentalia.

Evolution and pathology in chagas disease--a review

Trypanosoma cruzi acute infections often go unperceived, but one third of chronically infected individuals die of Chagas disease, showing diverse manifestations affecting the heart, intestines, and nervous systems. A common denominator of pathology in Chagas disease is the minimal rejection unit, whereby parasite-free target host cells are destroyed by immune system mononuclear effectors cells infiltrates. Another key feature stemming from T. cruzi infection is the integration of kDNA minicircles into the vertebrate host genome; horizontal transfer of the parasite DNA can undergo vertical transmission to the progeny of mammals and birds. kDNA integration-induced mutations can enter multiple loci in diverse chromosomes, generating new genes, pseudo genes and knock-outs, and resulting in genomic shuffling and remodeling over time. As a result of the juxtaposition of kDNA insertions with host open reading frames, novel chimeric products may be generated. Germ line transmission of kDNA-mutations determined the appearance of lesions in birds that are indistinguishable from those seen in Chagas disease patients. The production of tissue lesions showing typical minimal rejection units in birds' refractory to T. cruzi infection is consistent with the hypothesis that autoimmunity, likely triggered by integration-induced phenotypic alterations, plays a major role in the pathogenesis of Chagas disease.

Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent Anomalurus

Background

Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious.

Results

We have sequenced and performed extensive evolutionary analyses on the mitochondrial genome of the scaly-tailed flying squirrel Anomalurus sp., an enigmatic rodent whose phylogenetic affinities have been obscure and extensively debated. Our phylogenetic analyses of the coding regions of available complete mitochondrial genome sequences from Euarchontoglires suggest that Anomalurus is a sister taxon to the Hystricognathi, and that this clade represents the most basal divergence among sampled Rodentia. Bayesian dating methods incorporating a relaxed molecular clock provide divergence-time estimates which are consistently in agreement with the fossil record and which indicate a rapid radiation within Glires around 60 million years ago.

Conclusion

Taken together, the data presented provide a working hypothesis as to the phylogenetic placement of Anomalurus, underline the utility of mitochondrial sequences in the resolution of even relatively deep divergences and go some way to explaining the difficulty of conclusively resolving higher-level relationships within Glires with available data and methodologies.

Bonnet monkey (Macaca radiata) ovaries, like human oocytes, express four zona pellucida glycoproteins

Recent studies document that the zona pellucida matrix of human oocytes is composed of four glycoproteins designated as ZP1, ZP2, ZP3, and ZP4 instead of three as proposed in mouse model. In the present study, investigations were carried out to find the presence of the fourth ZP glycoprotein in the ovaries of non-human primates namely bonnet monkey (Macaca radiata). Employing total RNA isolated from bonnet monkey ovaries, the complementary deoxyribonucleic acid (cDNA) encoding bonnet monkey ZP1 (up to furin cleavage site) was successfully amplified by reverse transcribed polymerase chain reaction (RT-PCR). The deduced amino acid (aa) sequence of bonnet monkey ZP1 revealed 96.0% identity with human ZP1. The 21 cysteine residues present in bonnet monkey ZP1 were conserved in human, mouse, rat, quail, and chicken. Simultaneously, polyclonal antibodies were generated in mice against synthetic peptides corresponding to human ZP1 (P1, 137-150 aa; P2, 223-235 aa; P3, 237-251 aa; P4, 413-432 aa; and P5, 433-451 aa). Employing anti-peptide antibodies that were devoid of cross-reactivity as determined by ELISA with human/bonnet monkey recombinant ZP2, ZP3, and ZP4, the presence of ZP1 in the ovaries of bonnet monkey and human oocytes was demonstrated by indirect immunofluorescence. The antibodies against peptides P3 and P4 reacted only with the ZP of bonnet monkey ovaries and not other ovarian associated cells. The data presented in this manuscript provide evidence, for the first time, that the bonnet monkey ZP matrix is composed of four glycoproteins.

Bovine and human cathelicidin cationic host defense peptides similarly suppress transcriptional responses to bacterial lipopolysaccharide

Genomic approaches can be exploited to expose the complexities and conservation of biological systems such as the immune network across various mammalian species. In this study, temporal transcriptional expression profiles were analyzed in human and bovine monocytic cells in response to the TLR-4 agonist, LPS, in the presence or absence of their respective host defense peptides. The cathelicidin peptides, human LL-37 and bovine myeloid antimicrobial peptide-27 (BMAP-27), are homologs, yet they have diverged notably in terms of sequence similarity. In spite of their low sequence similarities, both of these cathelicidin peptides demonstrated potent, antiendotoxin activity in monocytic cells at low, physiologically relevant concentrations. Microarray studies indicated that 10 ng/ml LPS led to the up-regulation of 125 genes in human monocytes, 106 of which were suppressed in the presence of 5 mug/ml of the human peptide LL-37. To confirm and extend these data, temporal transcriptional responses to LPS were assessed in the presence or absence of the species-specific host defense peptides by quantitative real-time PCR. The transcriptional trends of 20 LPS-induced genes were analyzed in bovine and human monocytic cells. These studies demonstrated conserved trends of gene responses in that both peptides were able to profoundly suppress many LPS-induced genes. Consistent with this, the human and bovine peptides suppressed LPS-induced translocation of NF-kappaB subunits p50 and p65 into the nucleus of monocytic cells. However, there were also distinct differences in responses to LPS and the peptides; for example, treatment with 5 mug/ml BMAP-27 alone tended to influence gene expression (RELA, TNF-alpha-induced protein 2, MAPK phosphatase 1/dual specificity phosphatase 1, IkappaBkappaB, NFkappaBIL1, TNF receptor-associated factor 2) to a greater extent than did the same amount of human LL-37. We hypothesize that the immunomodulatory effects of the species-specific host defense peptides play a critical role in regulating inflammation and represent an evolutionarily conserved mechanism for maintaining homeostasis, although the sequence divergence of these peptides is substantial.

High amino acid variation in the intracellular domain of the pig prolactin receptor (PRLR) and its relation to ovulation rate and piglet survival traits1

Two polymorphisms of the porcine prolactin receptor (PRLR) gene were previously related to litter size by several authors; however, the magnitude and direction of such effects varied depending on the population analyzed. We have sequenced the complete coding region of the porcine PRLR gene and found 6 nonconservative SNP: C1217T (Leu/Pro406), C1283A (Asp/Ala428), G1439A (Lys/Arg480), T1528A (Met/Leu510), G1600A (Gly/Ser534), and G1789A (Gly/Ser597), within exon 10 of the gene, which encodes the entire intracytoplasmic domain of the protein. Eight haplotypes were found and were segregating at different frequencies in 6 porcine breeds. The effects of each individual SNP and haplotype were evaluated in an Iberian x Meishan F2 population using a univariate mixed-inheritance animal model. Significant effects on the number of corpora lutea were found for PRLR haplotypes (P < 0.012), confirming the previously reported associations of PRLR in this process and highlighting the importance of performing analysis of haplotypes rather than of individual SNP. Suggestive effects or tendencies were found for heart rate at birth (P < 0.10), rectal temperature (P < 0.05), and oxygen saturation (P < 0.10) 1 h after birth, and time to first suckle (P < 0.10). We found greater than expected levels of amino acid variability within the intracellular domain of the porcine PRLR, which have been associated with differences in the number of corpus lutea of sows and the preweaning survivability of piglets.