Museomics Dissects the Genetic Basis for Adaptive Seasonal Coloration in the Least Weasel.

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat coloration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage color moults across regions with varying winter snow. Whole-genome sequence data were obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris coloration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200-kb genomic region to coloration morph, which was validated by genotyping museum specimens from intermorph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change cosegregating with coloration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favored winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter color variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages.

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

The evolutionary pathways for local adaptation in mountain hares.

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole-genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome-wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area-specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression-driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares' adaptive variants appear predominantly species-specific. However, using coalescent simulations, we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our study shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post-glacial dynamics of species.

A low-latitude species pump: Peripheral isolation, parapatric speciation and mating-system evolution converge in a marine radiation.

Geologically recent radiations can shed light on speciation processes, but incomplete lineage sorting and introgressive gene flow render accurate evolutionary reconstruction and interpretation challenging. Independently evolving metapopulations of low dispersal taxa may provide an additional level of phylogeographic information, given sufficiently broad sampling and genome-wide sequencing. Evolution in the marine brown algal genus Fucus in the south-eastern North Atlantic was shaped by Quaternary climate-driven range shifts. Over this timescale, divergence and speciation occurred against a background of expansion-contraction cycles from multiple refugia, together with mating-system shifts from outcrossing (dioecy) to selfing hermaphroditism. We tested the hypothesis that peripheral isolation of range edge (dioecious) F. vesiculosus led to parapatric speciation and radiation of hermaphrodite lineages. Species tree methods using 876 single-copy nuclear genes and extensive geographic coverage produced conflicting topologies with respect to geographic clades of F. vesiculosus. All methods, however, revealed a new and early diverging hermaphrodite species, Fucus macroguiryi sp. nov. Both the multispecies coalescent and polymorphism-aware models (in contrast to concatenation) support sequential paraphyly in F. vesiculosus resulting from distinct evolutionary processes. Our results support (1) peripheral isolation of the southern F. vesiculosus clade prior to parapatric speciation and radiation of hermaphrodite lineages-a "low-latitude species pump". (2) Directional introgressive gene flow into F. vesiculosus around the present-day secondary contact zone (sympatric-allopatric boundary) between dioecious/hermaphrodite lineages as hermaphrodites expanded northwards, supported by concordance analysis and statistical tests of introgression. (3) Species boundaries in the extensive sympatric range are probably maintained by reproductive system (selfing in hermaphrodites) and reinforcement.

The Legacy of Recurrent Introgression during the Radiation of Hares.

Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1-4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation. [Adaptation; ancient introgression; hybridization; Lepus; phylogenomics.].

Introgression drives repeated evolution of winter coat color polymorphism in hares.

Changing from summer-brown to winter-white pelage or plumage is a crucial adaptation to seasonal snow in more than 20 mammal and bird species. Many of these species maintain nonwhite winter morphs, locally adapted to less snowy conditions, which may have evolved independently. Mountain hares (Lepus timidus) from Fennoscandia were introduced into the Faroe Islands in 1855. While they were initially winter-white, within ∼65 y all Faroese hares became winter-gray, a morph that occurs in the source population at low frequency. The documented population history makes this a valuable model for understanding the genetic basis and evolution of the seasonal trait polymorphism. Through whole-genome scans of differentiation and single-nucleotide polymorphism (SNP) genotyping, we associated winter coat color polymorphism to the genomic region of the pigmentation gene Agouti, previously linked to introgression-driven winter coat color variation in the snowshoe hare (Lepus americanus). Lower Agouti expression in the skin of winter-gray individuals during the autumn molt suggests that regulatory changes may underlie the color polymorphism. Variation in the associated genomic region shows signatures of a selective sweep in the Faroese population, suggesting that positive selection drove the fixation of the variant after the introduction. Whole-genome analyses of several hare species revealed that the winter-gray variant originated through introgression from a noncolor changing species, in keeping with the history of ancient hybridization between the species. Our findings show the recurrent role of introgression in generating winter coat color variation by repeatedly recruiting the regulatory region of Agouti to modulate seasonal coat color change.

Nucleotide diversity of functionally different groups of immune response genes in Old World camels based on newly annotated and reference-guided assemblies.

BACKGROUND: Immune-response (IR) genes have an important role in the defense against highly variable pathogens, and therefore, diversity in these genomic regions is essential for species' survival and adaptation. Although current genome assemblies from Old World camelids are very useful for investigating genome-wide diversity, demography and population structure, they have inconsistencies and gaps that limit analyses at local genomic scales. Improved and more accurate genome assemblies and annotations are needed to study complex genomic regions like adaptive and innate IR genes. RESULTS: In this work, we improved the genome assemblies of the three Old World camel species - domestic dromedary and Bactrian camel, and the two-humped wild camel - via different computational methods. The newly annotated dromedary genome assembly CamDro3 served as reference to scaffold the NCBI RefSeq genomes of domestic Bactrian and wild camels. These upgraded assemblies were then used to assess nucleotide diversity of IR genes within and between species, and to compare the diversity found in immune genes and the rest of the genes in the genome. We detected differences in the nucleotide diversity among the three Old World camelid species and between IR gene groups, i.e., innate versus adaptive. Among the three species, domestic Bactrian camels showed the highest mean nucleotide diversity. Among the functionally different IR gene groups, the highest mean nucleotide diversity was observed in the major histocompatibility complex. CONCLUSIONS: The new camel genome assemblies were greatly improved in terms of contiguity and increased size with fewer scaffolds, which is of general value for the scientific community. This allowed us to perform in-depth studies on genetic diversity in immunity-related regions of the genome. Our results suggest that differences of diversity across classes of genes appear compatible with a combined role of population history and differential exposures to pathogens, and consequent different selective pressures.

Analysis of substitution rates showed that TLR5 is evolving at different rates among mammalian groups.

BACKGROUND: Toll-like receptors (TLRs) are the most widely studied innate immunity receptors responsible for recognition of invading pathogens. Among the TLR family, TLR5 is the only that senses and recognizes flagellin, the major protein of bacterial flagella. TLR5 has been reported to be under overall purifying selection in mammals, with a small proportion of codons under positive selection. However, the variation of substitution rates among major mammalian groups has been neglected. Here, we studied the evolution of TLR5 in mammals, comparing the substitution rates among groups. RESULTS: In this study we analysed the TLR5 substitution rates in Euungulata, Carnivora, Chiroptera, Primata, Rodentia and Lagomorpha, groups. For that, Tajima's relative rate test, Bayesian inference of evolutionary rates and genetic distances were estimated with CODEML's branch model and RELAX. The combined results showed that in the Lagomorpha, Rodentia, Carnivora and Chiroptera lineages TLR5 is evolving at a higher substitution rate. The RELAX analysis further suggested a significant relaxation of selective pressures for the Lagomorpha (K = 0.22, p < 0.01), Rodentia (K = 0.58, p < 0.01) and Chiroptera (K = 0.65, p < 0.01) lineages and for the Carnivora ancestral branches (K = 0.13, p < 0.01). CONCLUSIONS: Our results show that the TLR5 substitution rate is not uniform among mammals. In fact, among the different mammal groups studied, the Lagomorpha, Rodentia, Carnivora and Chiroptera are evolving faster. This evolutionary pattern could be explained by 1) the acquisition of new functions of TLR5 in the groups with higher substitution rate, i.e. TLR5 neofunctionalization, 2) by the beginning of a TLR5 pseudogenization in these groups due to some redundancy between the TLRs genes, or 3) an arms race between TLR5 and species-specific parasites.

The evolutionary history of the Cape hare (Lepus capensis sensu lato): insights for systematics and biogeography.

Inferring the phylogeography of species with large distributions helps deciphering major diversification patterns that may occur in parallel across taxa. Here, we infer the evolutionary history of the Cape hare, Lepus capensis sensu lato, a species distributed from southern Africa to Asia, by analyzing variation at 18 microsatellites and 9 DNA (1 mitochondrial and 8 nuclear) sequenced loci, from field and museum-collected samples. Using a combination of assignment and coalescent-based methods, we show that the Cape hare is composed of five evolutionary lineages, distributed in distinct biogeographic regions-north-western Africa, eastern Africa, southern Africa, the Near East and the Arabian Peninsula. A deep phylogenetic break possibly dating to the Early Pleistocene was inferred between the African and Asian L. capensis groups, and the latter appear more closely related to other Eurasian hare species than to African Cape hares. The inferred phylogeographic structure is shared by numerous taxa distributed across the studied range, suggesting that environmental changes, such as the progressive aridification of the Saharo-Arabian desert and the fluctuations of savannah habitats in Sub-Saharan Africa, had comparable impacts across species. Fine-scale analyses of the western Sahara-Sahel populations showed rich fragmentation patterns for mitochondrial DNA but not for microsatellites, compatible with the environmental heterogeneity of the region and female philopatry. The complex evolutionary history of L. capensis sensu lato, which possibly includes interspecific gene flow, is not reflected by taxonomy. Integrating evolutionary inference contributes to an improved characterization of biodiversity, which is fundamental to foster the conservation of relevant evolutionary units.

The transcriptional landscape of seasonal coat colour moult in the snowshoe hare.

Seasonal coat colour change is an important adaptation to seasonally changing environments but the evolution of this and other circannual traits remains poorly understood. In this study, we use gene expression to understand seasonal coat colour moulting in wild snowshoe hares (Lepus americanus). We used hair colour to follow the progression of the moult, simultaneously sampling skin from three moulting stages in hares collected during the peak of the spring moult from white winter to brown summer pelage. Using RNA sequencing, we tested whether patterns of expression were consistent with predictions based on the established phases of the hair growth cycle. We found functionally consistent clustering across skin types, with 766 genes differentially expressed between moult stages. "White" pelage showed more differentially expressed genes that were upregulated relative to other skin types, involved in the transition between late telogen (quiescent stage) and the onset of anagen (proliferative stage). Skin samples from transitional "intermediate" and "brown" pelage were transcriptionally similar and resembled the regressive transition to catagen (regressive stage). We also detected differential expression of several key circadian clock and pigmentation genes, providing important means to dissect the bases of alternate seasonal colour morphs. Our results reveal that pelage colour is a useful biomarker for seasonal change but that there is a consistent lag between the main gene expression waves and change in visible coat colour. These experiments establish that developmental sampling from natural populations of nonmodel organisms can provide a crucial resource to dissect the genetic basis and evolution of complex seasonally changing traits.

The genomic architecture of population divergence between subspecies of the European rabbit.

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

Reference-free population genomics from next-generation transcriptome data and the vertebrate-invertebrate gap.

In animals, the population genomic literature is dominated by two taxa, namely mammals and drosophilids, in which fully sequenced, well-annotated genomes have been available for years. Data from other metazoan phyla are scarce, probably because the vast majority of living species still lack a closely related reference genome. Here we achieve de novo, reference-free population genomic analysis from wild samples in five non-model animal species, based on next-generation sequencing transcriptome data. We introduce a pipe-line for cDNA assembly, read mapping, SNP/genotype calling, and data cleaning, with specific focus on the issue of hidden paralogy detection. In two species for which a reference genome is available, similar results were obtained whether the reference was used or not, demonstrating the robustness of our de novo inferences. The population genomic profile of a hare, a turtle, an oyster, a tunicate, and a termite were found to be intermediate between those of human and Drosophila, indicating that the discordant genomic diversity patterns that have been reported between these two species do not reflect a generalized vertebrate versus invertebrate gap. The genomic average diversity was generally higher in invertebrates than in vertebrates (with the notable exception of termite), in agreement with the notion that population size tends to be larger in the former than in the latter. The non-synonymous to synonymous ratio, however, did not differ significantly between vertebrates and invertebrates, even though it was negatively correlated with genetic diversity within each of the two groups. This study opens promising perspective regarding genome-wide population analyses of non-model organisms and the influence of population size on non-synonymous versus synonymous diversity.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach.

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

Sequencing of Sylvilagus VDJ genes reveals a new VHa allelic lineage and shows that ancient VH lineages were retained differently in leporids.

Antigen recognition by immunoglobulins depends upon initial rearrangements of heavy chain V, D, and J genes. In leporids, a unique system exists for the VH genes usage that exhibit highly divergent lineages: the VHa allotypes, the Lepus sL lineage and the VHn genes. For the European rabbit (Oryctolagus cuniculus), four VHa lineages have been described, the a1, a2, a3 and a4. For hares (Lepus sp.), one VHa lineage was described, the a2L, as well as a more ancient sL lineage. Both genera use the VHn genes in a low frequency of their VDJ rearrangements. To address the hypothesis that the VH specificities could be associated with different environments, we sequenced VDJ genes from a third leporid genus, Sylvilagus. We found a fifth and equally divergent VHa lineage, the a5, and an ancient lineage, the sS, related to the hares' sL, but failed to obtain VHn genes. These results show that the studied leporids employ different VH lineages in the generation of the antibody repertoire, suggesting that the leporid VH genes are subject to strong selective pressure likely imposed by specific pathogens.

The hidden history of the snowshoe hare, Lepus americanus: extensive mitochondrial DNA introgression inferred from multilocus genetic variation.

Hybridization drives the evolutionary trajectory of many species or local populations, and assessing the geographic extent and genetic impact of interspecific gene flow may provide invaluable clues to understand population divergence or the adaptive relevance of admixture. In North America, hares (Lepus spp.) are key species for ecosystem dynamics and their evolutionary history may have been affected by hybridization. Here we reconstructed the speciation history of the three most widespread hares in North America - the snowshoe hare (Lepus americanus), the white-tailed jackrabbit (L. townsendii) and the black-tailed jackrabbit (L. californicus) - by analysing sequence variation at eight nuclear markers and one mitochondrial DNA (mtDNA) locus (6240 bp; 94 specimens). A multilocus-multispecies coalescent-based phylogeny suggests that L. americanus diverged ~2.7 Ma and that L. californicus and L. townsendii split more recently (~1.2 Ma). Within L. americanus, a deep history of cryptic divergence (~2.0 Ma) was inferred, which coincides with major speciation events in other North American species. While the isolation-with-migration model suggested that nuclear gene flow was generally rare or absent among species or major genetic groups, coalescent simulations of mtDNA divergence revealed historical mtDNA introgression from L. californicus into the Pacific Northwest populations of L. americanus. This finding marks a history of past reticulation between these species, which may have affected other parts of the genome and influence the adaptive potential of hares during climate change.

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare.

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.

Advances in Ecological Speciation: an integrative approach.

The role of natural selection in promoting reproductive isolation has received substantial renewed interest within the last two decades. As a consequence, the study of ecological speciation has become an extremely productive research area in modern evolutionary biology. Recent innovations in sequencing technologies offer an unprecedented opportunity to study the mechanisms involved in ecological speciation. Genome scans provide significant insights but have some important limitations; efforts are needed to integrate them with other approaches to make full use of the sequencing data deluge. An international conference 'Advances in Ecological Speciation' organized by the University of Porto (Portugal) aimed to review current progress in ecological speciation. Using some of the examples presented at the conference, we highlight the benefits of integrating ecological and genomic data and discuss different mechanisms of parallel evolution. Finally, future avenues of research are suggested to advance our knowledge concerning the role of natural selection in the establishment of reproductive isolation during ecological speciation.

Evolutionary and functional characterization of lagomorph guanylate-binding proteins: a story of gain and loss and shedding light on expression, localization and innate immunity-related functions.

Guanylate binding proteins (GBPs) are an evolutionarily ancient family of proteins that are widely distributed among eukaryotes. They belong to the dynamin superfamily of GTPases, and their expression can be partially induced by interferons (IFNs). GBPs are involved in the cell-autonomous innate immune response against bacterial, parasitic and viral infections. Evolutionary studies have shown that GBPs exhibit a pattern of gene gain and loss events, indicative for the birth-and-death model of evolution. Most species harbor large GBP gene clusters that encode multiple paralogs. Previous functional and in-depth evolutionary studies have mainly focused on murine and human GBPs. Since rabbits are another important model system for studying human diseases, we focus here on lagomorphs to broaden our understanding of the multifunctional GBP protein family by conducting evolutionary analyses and performing a molecular and functional characterization of rabbit GBPs. We observed that lagomorphs lack GBP3, 6 and 7. Furthermore, Leporidae experienced a loss of GBP2, a unique duplication of GBP5 and a massive expansion of GBP4. Gene expression analysis by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and transcriptome data revealed that leporid GBP expression varied across tissues. Overexpressed rabbit GBPs localized either uniformly and/or discretely to the cytoplasm and/or to the nucleus. Oryctolagus cuniculus (oc)GBP5L1 and rarely ocGBP5L2 were an exception, colocalizing with the trans-Golgi network (TGN). In addition, four ocGBPs were IFN-inducible and only ocGBP5L2 inhibited furin activity. In conclusion, from an evolutionary perspective, lagomorph GBPs experienced multiple gain and loss events, and the molecular and functional characteristics of ocGBP suggest a role in innate immunity.

Not so pseudo: the evolutionary history of protein phosphatase 1 regulatory subunit 2 and related pseudogenes.

BACKGROUND: Pseudogenes are traditionally considered "dead" genes, therefore lacking biological functions. This view has however been challenged during the last decade. This is the case of the Protein phosphatase 1 regulatory subunit 2 (PPP1R2) or inhibitor-2 gene family, for which several incomplete copies exist scattered throughout the genome. RESULTS: In this study, the pseudogenization process of PPP1R2 was analyzed. Ten PPP1R2-related pseudogenes (PPP1R2P1-P10), highly similar to PPP1R2, were retrieved from the human genome assembly present in the databases. The phylogenetic analysis of mammalian PPP1R2 and related pseudogenes suggested that PPP1R2P7 and PPP1R2P9 retroposons appeared before the great mammalian radiation, while the remaining pseudogenes are primate-specific and retroposed at different times during Primate evolution. Although considered inactive, four of these pseudogenes seem to be transcribed and possibly possess biological functions. Given the role of PPP1R2 in sperm motility, the presence of these proteins was assessed in human sperm, and two PPP1R2-related proteins were detected, PPP1R2P3 and PPP1R2P9. Signatures of negative and positive selection were also detected in PPP1R2P9, further suggesting a role as a functional protein. CONCLUSIONS: The results show that contrary to initial observations PPP1R2-related pseudogenes are not simple bystanders of the evolutionary process but may rather be at the origin of genes with novel functions.

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome.

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.