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.

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.

Among the protein translocation pathways of the thylakoid membrane in chloroplasts, the DeltapH/TAT pathway is unique in several aspects. In vitro transport assays with isolated chloroplasts or thylakoids have defined the trans-thylakoidal proton gradient as the sole requirement for effecting transport. From these studies, evidence has also accumulated indicating that, in contrast to the remaining protein transport pathways present in the thylakoid membrane, the DeltapH/TAT pathway is able to mediate the transport of folded proteins. The present work has established a novel approach to demonstrate the transport of folded proteins by this pathway in vivo. For this purpose, Arabidopsis thaliana plants were stably transformed with gene constructs expressing enhanced green fluorescent protein (EGFP) alone or fused to the transit peptides of different chloroplast proteins under the control of the 35S CAMV promoter. The intracellular and intraorganellar distribution of EGFP in the resulting transformants showed that while all the chloroplast transit peptides efficiently mediated the transport of EGFP into plastids, only those specific for the DeltapH/TAT pathway were able to direct the protein into the thylakoid lumen as well. This could be demonstrated both by fluorescence and immunoelectron microscopy. Analysis of isolated and fractionated chloroplasts using western blot and spectrofluorometric assays confirmed the presence of folded EGFP solely within the thylakoid lumen of these lines. These results strongly suggest that the protein adopts a folded state in the chloroplast stroma and thus, can only be translocated further into the chloroplast lumen by the DeltapH/TAT pathway.

We have tested the potential of EGFP, a derivative of the green fluorescent protein (GFP), as a passenger protein for the analysis of protein transport processes across the thylakoid membranes in chloroplasts. In contrast to the majority of fusion proteins commonly used in such studies, EGFP is not of plant origin and can therefore be assumed to behave like a "neutral" passenger protein that is unaffected by any internal plant regulatory circuits. Our in vitro transport experiments clearly demonstrate that EGFP is a suitable passenger protein that can be correctly targeted either to the stroma or to the thylakoid lumen if fused to the appropriate transit peptide. The transport of EGFP across the thylakoid membrane shows, however, a clear pathway preference. While the protein is efficiently targeted by the deltapH/TAT pathway, transport by the Sec pathway is barely detectable, either with isolated thylakoids or with intact chloroplasts. This pathway specificity suggests that EGFP is folded immediately after import into the chloroplast stroma, thus preventing further translocation across the thylakoid membrane by the Sec translocase. The data obtained provide a good basis for the development of molecular tools for transport studies using EGFP as a passenger protein. Furthermore, plant lines expressing corresponding EGFP chimeras are expected to allow in vivo studies on the transport and sorting mechanisms involved in the biogenesis of the chloroplast.

Infection of Balb/c mice with Echinococcus granulosus protoscoleces constitutes the model for secondary hydatid infection. The immune response of Balb/c mice infected with E. granulosus is characterized by secretion of antibodies specific for carbohydrate epitopes and production of type-2 cytokines. A role for glycoconjugates in the induction of type-2 responses has been suggested in other host--parasite systems. Although glycoconjugates are immunogenic in E. granulosus infection, the role of these molecules in the establishment of the type-2 response has never been analysed. In this study, a carbohydrate rich fraction (E4+) from E. granulosus protoscoleces was obtained using the monoclonal antibody E492/G1 specific for the moiety Galalpha(1,4)Gal which is widely represented in protoscoleces and other E. granulosus antigenic preparations. The results showed that E4+ was immunogenic in Balb/c mice evoking an antibody response mainly directed against carbohydrate epitopes. In addition, splenocytes from E4+-immunized mice showed suppressed proliferative responses to Con A and E4+ induced IL-10 secretion by E4+-primed and naive splenocytes. The fraction E4+ also was immunogenic in infected mice during early infection. In this case also, splenocytes from infected mice as well as peritoneal cells from infected or naive mice, when stimulated in vitro with E4+, secreted IL-10. Collectively, these results suggest that E4+ may be involved in immunosuppression phenomena and, by stimulating IL-10 secretion, may contribute to the induction and sustaining of the type-2 cytokine response established in early experimental infection.

Introgressive hybridization is an important and widespread evolutionary process, but the relative roles of neutral demography and natural selection in promoting massive introgression are difficult to assess and an important matter of debate. Hares from the Iberian Peninsula provide an appropriate system to study this question. In its northern range, the Iberian hare, Lepus granatensis, shows a northwards gradient of increasing mitochondrial DNA (mtDNA) introgression from the arctic/boreal L. timidus, which it presumably replaced after the last glacial maximum. Here, we asked whether a south-north expansion wave of L. granatensis into L. timidus territory could underlie mtDNA introgression, and whether nuclear genes interacting with mitochondria (“mitonuc” genes) were affected. We extended previous RNA-sequencing and produced a comprehensive annotated transcriptome assembly for L. granatensis. We then genotyped 100 discovered nuclear SNPs in 317 specimens spanning the species range. The distribution of allele frequencies across populations suggests a northwards range expansion, particularly in the region of mtDNA introgression. We found no correlation between variants at 39 mitonuc genes and mtDNA introgression frequency. Whether the nuclear and mitochondrial genomes coevolved will need a thorough investigation of the hundreds of mitonuc genes, but range expansion and species replacement likely promoted massive mtDNA introgression.

Vertebrate genomes encode a diversity of G protein-coupled receptor (GPCR) that belong to large gene families and are used by olfactory systems to detect chemical cues found in the environment. It is not clear however, if individual receptors from these large gene families have evolved roles that are specific to certain life stages. Here, we used deep sequencing to identify differentially expressed receptor transcripts in the olfactory epithelia (OE) of freshwater, seawater and sexually mature male eels (Anguilla anguilla). This species is particularly intriguing because of its complex life cycle, extreme long-distance migrations and early-branching position within the teleost phylogeny. In the A. anguillaOE, we identified full-length transcripts for 13, 112, 6 and 38 trace amine-associated receptors, odorant receptors (OR) and type I and type II vomeronasal receptors (V1R and V2R). Most of these receptors were expressed at similar levels at different life stages and a subset of OR and V2R-like transcripts was more abundant in sexually mature males suggesting that ORs and V2R-like genes are important for reproduction. We also identified a set of GPCR signal transduction genes that were differentially expressed indicating that eels make use of different GPCR signal transduction genes at different life stages. The finding that a diversity of chemosensory receptors is expressed in the olfactory epithelium and that a subset is differentially expressed suggests that most receptors belonging to large chemosensory gene families have functions that are important at multiple life stages, while a subset has evolved specific functions at different life stages.

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.