Chromosomal Fusions Facilitate Adaptation to Divergent Environments in Threespine Stickleback.

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.

Molecular mechanisms of Eda-mediated adaptation to freshwater in threespine stickleback.

A main goal of evolutionary biology is to understand the genetic basis of adaptive evolution. Although the genes that underlie some adaptive phenotypes are now known, the molecular pathways and regulatory mechanisms mediating the phenotypic effects of those genes often remain a black box. Unveiling this black box is necessary to fully understand the genetic basis of adaptive phenotypes, and to understand why particular genes might be used during phenotypic evolution. Here, we investigated which genes and regulatory mechanisms are mediating the phenotypic effects of the Eda haplotype, a locus responsible for the loss of lateral plates and changes in the sensory lateral line of freshwater threespine stickleback (Gasterosteus aculeatus) populations. Using a combination of RNAseq and a cross design that isolated the Eda haplotype on a fixed genomic background, we found that the Eda haplotype affects both gene expression and alternative splicing of genes related to bone development, neuronal development and immunity. These include genes in conserved pathways, like the BMP, netrin and bradykinin signalling pathways, known to play a role in these biological processes. Furthermore, we found that differentially expressed and differentially spliced genes had different levels of connectivity and expression, suggesting that these factors might influence which regulatory mechanisms are used during phenotypic evolution. Taken together, these results provide a better understanding of the mechanisms mediating the effects of an important adaptive locus in stickleback and suggest that alternative splicing could be an important regulatory mechanism mediating adaptive phenotypes.

Field vole-associated Traemmersee hantavirus from Germany represents a novel hantavirus species.

Vole-associated hantaviruses occur in the Old and New World. Tula orthohantavirus (TULV) is widely distributed throughout the European continent in its reservoir, the common vole (Microtus arvalis), but the virus was also frequently detected in field voles (Microtus agrestis) and other vole species. TULV and common voles are absent from Great Britain. However, field voles there harbor Tatenale and Kielder hantaviruses. Here we screened 126 field voles and 13 common voles from Brandenburg, Germany, for hantavirus infections. One common vole and four field voles were anti-TULV antibody and/or TULV RNA positive. In one additional, seropositive field vole a novel hantavirus sequence was detected. The partial S and L segment nucleotide sequences were only 61.1% and 75.6% identical to sympatrically occurring TULV sequences, but showed highest similarity of approximately 80% to British Tatenale and Kielder hantaviruses. Subsequent determination of the entire nucleocapsid (N), glycoprotein (GPC), and RNA-dependent RNA polymerase encoding sequences and determination of the pairwise evolutionary distance (PED) value for the concatenated N and GPC amino acid sequences confirmed a novel orthohantavirus species, tentatively named Traemmersee orthohantavirus. The identification of this novel hantavirus in a field vole from eastern Germany underlines the necessity of a large-scale, broad geographical hantavirus screening of voles to understand evolutionary processes of virus-host associations and host switches.

Cocirculation of Leptospira spp. and multiple orthohantaviruses in rodents, Lithuania, Northern Europe.

In Europe, zoonotic Leptospira spp. and orthohantaviruses are mainly associated with specific rodent hosts. These pathogens cause febrile human diseases with similar symptoms and disease progression. In Lithuania, the presence of Dobrava-Belgrade orthohantavirus (DOBV), Tula orthohantavirus (TULV) and Leptospira spp. in rodent reservoirs is still unknown, and Puumala orthohantavirus (PUUV) was detected in bank voles (Clethrionomys glareolus) at only one site. Therefore, we collected and screened 1617 rodents and insectivores from Lithuania for zoonotic (re-)emerging Leptospira and orthohantaviruses. We detected Leptospira DNA in six rodent species, namely striped field mouse (Apodemus agrarius), yellow-necked mouse (Apodemus flavicollis), bank vole, common vole (Microtus arvalis), field vole (Microtus agrestis) and root vole (Microtus oeconomus). Leptospira DNA was detected with an overall mean prevalence of 4.4% (range 3.7%-7.9% per rodent species). We detected DOBV RNA in 5.6% of the striped field mice, PUUV RNA in 1% of bank voles and TULV RNA in 4.6% of common voles, but no Leptospira DNA in shrews and no hantavirus-Leptospira coinfections in rodents. Based on the complete coding sequences of the three genome segments, two distant DOBV phylogenetic lineages in striped field mice, one PUUV strain in bank voles and two TULV strains in common voles were identified. The Leptospira prevalence for striped field mice and yellow-necked mice indicated a significant negative effect of the distance to water points. The detection of (re-)emerging human pathogenic Leptospira and three orthohantaviruses in rodent reservoirs in Lithuania calls for increased awareness of public health institutions and allows the improvement of molecular diagnostics for pathogen identification.

Identification of a novel hantavirus strain in the root vole (Microtus oeconomus) in Lithuania, Eastern Europe.

Hantaviruses are zoonotic pathogens that can cause subclinical to lethal infections in humans. In Europe, five orthohantaviruses are present in rodents: Myodes-associated Puumala orthohantavirus (PUUV), Microtus-associated Tula orthohantavirus, Traemmersee hantavirus (TRAV)/ Tatenale hantavirus (TATV)/ Kielder hantavirus, rat-borne Seoul orthohantavirus, and Apodemus-associated Dobrava-Belgrade orthohantavirus (DOBV). Human PUUV and DOBV infections were detected previously in Lithuania, but the presence of Microtus-associated hantaviruses is not known. For this study we screened 234 Microtus voles, including root voles (Microtus oeconomus), field voles (Microtus agrestis) and common voles (Microtus arvalis) from Lithuania for hantavirus infections. This initial screening was based on reverse transcription-polymerase chain reaction (RT-PCR) targeting the S segment and serological analysis. A novel hantavirus was detected in eight of 79 root voles tentatively named "Rusne virus" according to the capture location and complete genome sequences were determined. In the coding regions of all three genome segments, Rusne virus showed high sequence similarity to TRAV and TATV and clustered with Kielder hantavirus in phylogenetic analyses of partial S and L segment sequences. Pairwise evolutionary distance analysis confirmed Rusne virus as a strain of the species TRAV/TATV. Moreover, we synthesized the entire nucleocapsid (N) protein of Rusne virus in Saccharomyces cerevisiae. We observed cross-reactivity of antibodies raised against other hantaviruses, including PUUV, with this new N protein. ELISA investigation of all 234 voles detected Rusne virus-reactive antibodies exclusively in four of 79 root voles, all being also RNA positive, but not in any other vole species. In conclusion, the detection of Rusne virus RNA in multiple root voles at the same trapping site during three years and its absence in sympatric field voles suggests root voles as the reservoir host of this novel virus. Future investigations should evaluate host association of TRAV, TATV, Kielder virus and the novel Rusne virus and their evolutionary relationships.

Assessing Genome-Wide Diversity in European Hantaviruses through Sequence Capture from Natural Host Samples.

Research on the ecology and evolution of viruses is often hampered by the limitation of sequence information to short parts of the genomes or single genomes derived from cultures. In this study, we use hybrid sequence capture enrichment in combination with high-throughput sequencing to provide efficient access to full genomes of European hantaviruses from rodent samples obtained in the field. We applied this methodology to Tula (TULV) and Puumala (PUUV) orthohantaviruses for which analyses from natural host samples are typically restricted to partial sequences of their tri-segmented RNA genome. We assembled a total of ten novel hantavirus genomes de novo with very high coverage (on average >99%) and sequencing depth (average >247×). A comparison with partial Sanger sequences indicated an accuracy of >99.9% for the assemblies. An analysis of two common vole (Microtus arvalis) samples infected with two TULV strains each allowed for the de novo assembly of all four TULV genomes. Combining the novel sequences with all available TULV and PUUV genomes revealed very similar patterns of sequence diversity along the genomes, except for remarkably higher diversity in the non-coding region of the S-segment in PUUV. The genomic distribution of polymorphisms in the coding sequence was similar between the species, but differed between the segments with the highest sequence divergence of 0.274 for the M-segment, 0.265 for the S-segment, and 0.248 for the L-segment (overall 0.258). Phylogenetic analyses showed the clustering of genome sequences consistent with their geographic distribution within each species. Genome-wide data yielded extremely high node support values, despite the impact of strong mutational saturation that is expected for hantavirus sequences obtained over large spatial distances. We conclude that genome sequencing based on capture enrichment protocols provides an efficient means for ecological and evolutionary investigations of hantaviruses at an unprecedented completeness and depth.

Spatial and Temporal Evolutionary Patterns in Puumala Orthohantavirus (PUUV) S Segment.

The S segment of bank vole (Clethrionomys glareolus)-associated Puumala orthohantavirus (PUUV) contains two overlapping open reading frames coding for the nucleocapsid (N) and a non-structural (NSs) protein. To identify the influence of bank vole population dynamics on PUUV S segment sequence evolution and test for spillover infections in sympatric rodent species, during 2010-2014, 883 bank voles, 357 yellow-necked mice (Apodemus flavicollis), 62 wood mice (A. sylvaticus), 149 common voles (Microtus arvalis) and 8 field voles (M. agrestis) were collected in Baden-Wuerttemberg and North Rhine-Westphalia, Germany. In total, 27.9% and 22.3% of bank voles were positive for PUUV-reactive antibodies and PUUV-specific RNA, respectively. One of eight field voles was PUUV RNA-positive, indicating a spillover infection, but none of the other species showed evidence of PUUV infection. Phylogenetic and isolation-by-distance analyses demonstrated a spatial clustering of PUUV S segment sequences. In the hantavirus outbreak years 2010 and 2012, PUUV RNA prevalence was higher in our study regions compared to non-outbreak years 2011, 2013 and 2014. NSs amino acid and nucleotide sequence types showed temporal and/or local variation, whereas the N protein was highly conserved in the NSs overlapping region and, to a lower rate, in the N alone coding part.