Phylogeography of the Tibetan hamster Cricetulus kamensis in response to uplift and environmental change in the Qinghai-Tibet Plateau

AIM

The evolutionary process of an organism provides valuable data toward an understanding of the Earth evolution history. To investigate the relationship between the uplift of the Qinghai-Tibet Plateau (QTP) and mammalian evolution since the late Cenozoic, the geographic distribution of genetic variations in the Tibetan hamster Cricetulus kamensis was investigated using phylogeographical methods. In particular, population divergence, demographic history, genetic variation, and the prediction of species distribution area were investigated.

LOCATION

The Qinghai-Tibet Plateau.

METHODS

A total of 53 specimens, representing 13 geographic populations, were collected from the QTP. The phylogeographical pattern and demographic history of C. kamensis were analyzed, and the probable factors in the QTP uplift and the Quaternary glacial periods were inferred from one nuclear and four mitochondrial genes. Furthermore, the species distribution model (SDM) was used to predict changes in potentially suitable habitats since the last Interglacial.

RESULTS

Phylogenetic analysis demonstrated that two major genetic differentiations of the C. kamensis population occurred during the Early Pleistocene that were influenced by the Qing-Zang tectonic movement from the Middle Pliocene to the Early Pleistocene. Genetic distance between two major clades indicated low genetic divergence. Demographic history analysis showed that the C. kamensis population was affected by the Quaternary glacial period. SDM analysis indicated that C. kamensis was endemic to the QTP and the suitable habitat was affected by climate change, especially during the Last Glacial Maximum (LGM).

MAIN CONCLUSION

Our results indicated that the QTP uplift led to the population divergence of C. kamensis, and vicariance well accounted for the geographic distribution of genetic variation in C. kamensis as a result of genetic divergence and lack of gene flow. The genetic distance shows that C. alticola may be a subspecies of C. kamensis. Demographic history analysis suggests that the QTP was affected by the last glacial period. SDM analysis supports that almost the entire QTP is covered by a huge ice sheet during the LGM.

Conservation of genetic uniqueness in remaining populations of red squirrels (Sciurus vulgaris L.) in the South of England

The Eurasian red squirrel (Sciurus vulgaris) is an emblematic species for conservation, and its decline in the British Isles exemplifies the impact that alien introductions can have on native ecosystems. Indeed, red squirrels in this region have declined dramatically over the last 60 years due to the spread of squirrelpox virus following the introduction of the gray squirrel (Sciurus carolinensis). Currently, red squirrel populations in Britain are fragmented and need to be closely monitored in order to assess their viability and the effectiveness of conservation efforts. The situation is even more dramatic in the South of England, where S. vulgaris survives only on islands (Brownsea Island, Furzey Island, and the Isle of Wight). Using the D-loop, we investigated the genetic diversity and putative ancestry of the squirrels from Southern England and compared them to a European dataset composed of 1,016 samples from 54 populations. We found that our three populations were more closely related to other squirrels from the British Isles than squirrels from Europe, showed low genetic diversity, and also harbored several private haplotypes. Our study demonstrates how genetically unique the Southern English populations are in comparison with squirrels from the continental European range. We report the presence of four private haplotypes, suggesting that these populations may potentially harbor distinct genetic lineages. Our results emphasize the importance of preserving these isolated red squirrel populations for the conservation of the species.

Ecotypic divergences of the alpine herb Potentilla matsumurae adapted to fellfield-snowbed habitats across a series of mountain sky islands

PREMISE

Divergent selection due to environmental heterogeneity can lead to local adaptation. However, the ecological and evolutionary processes of local adaptation that occurs across multiple regions are often unknown. Our previous studies reported on the ecotypic divergence within a local area of variation of Potentilla matsumurae, an alpine herb adapted to the fellfield-snowbed environment. Here we investigated large-scale geographic patterns of ecotypic differentiation in this species to infer local adaptation and selective forces across multiple regions.

METHODS

We compiled information on the overall distributions of fellfield and snowbed habitats on the mountains in Japan across the distribution of the species. Next, we conducted common garden experiments to test the adaptive divergence of the fellfield-snowbed plants derived from multiple regions. Finally, we evaluated phylogeographic structures based on cpDNA and allozyme variations and inferred the evolutionary history of ecotype differentiation.

RESULTS

The mosaic distribution of the fellfield-snowbed ecotypes across isolated mountaintops constitutes indirect evidence for habitat-specific natural selection. The significant difference in survivorship between the ecotypes observed in a controlled snow environment provides more substantial evidence of local selection. Phylogeographic structures support the hypothesis that ecotypic divergence events from fellfield to snowbed populations occurred independently in at least two distinct regions.

CONCLUSIONS

Ecotypic divergence of P. matsumurae has occurred across a series of mountain sky islands. Local selection in snowy environments is a driving force that maintains the divergent ecotypes across multiple mountain regions and can contribute to the diversification of plants in heavy-snow regions.

The comparative biogeography of Philippine geckos challenges predictions from a paradigm of climate-driven vicariant diversification across an island archipelago

A primary goal of biogeography is to understand how large-scale environmental processes, like climate change, affect diversification. One often-invoked but seldom tested process is the "species-pump" model, in which repeated bouts of cospeciation are driven by oscillating climate-induced habitat connectivity cycles. For example, over the past three million years, the landscape of the Philippine Islands has repeatedly coalesced and fragmented due to sea-level changes associated with glacial cycles. This repeated climate-driven vicariance has been proposed as a model of speciation across evolutionary lineages codistributed throughout the islands. This model predicts speciation times that are temporally clustered around the times when interglacial rises in sea level fragmented the islands. To test this prediction, we collected comparative genomic data from 16 pairs of insular gecko populations. We analyze these data in a full-likelihood, Bayesian model-choice framework to test for shared divergence times among the pairs. Our results provide support against the species-pump model prediction in favor of an alternative interpretation, namely that each pair of gecko populations diverged independently. These results suggest the repeated bouts of climate-driven landscape fragmentation have not been an important mechanism of speciation for gekkonid lizards across the Philippine Archipelago.

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

Phylogeography and species limits in the red-shouldered hawk (Buteo lineatus): Characterization of the Northern Florida Suture Zone in birds

The North American red-shouldered hawk, Buteo lineatus, is comprised of two widely allopatric eastern and western populations with an additional well-marked subspecies in the Florida peninsula. The two eastern populations meet in northern Florida, the location of a well-known suture zone in many nonavian organisms. We sequenced the complete mitochondrial ND2 gene and two nuclear introns to investigate its genetic population structure and species status. No mitochondrial haplotypes were shared between the eastern and western populations, and genetic variance among 14 populations was 0.42; almost all of this (0.40) was distributed among the three regions. A clade of haplotypes very common in the Florida peninsula decreased in frequency elsewhere and, when modeled as a hybrid zone, had an estimated width of 1,158 km with a center near Ocala, FL. Ecological niche modeling suggests the western, eastern, and Florida peninsula populations were geographically isolated during the last glacial maximum. We consider these to represent three phylogenetic species. A coalescent analysis incorporating incomplete lineage sorting and gene tree uncertainty also suggested the divergence between the western and eastern populations is consistent with species-level divergence. With the addition of this hawk, four avian species are now known to hybridize along the Gulf Coast of the United States in or near the Northern Florida Suture Zone. The widths of these avian zones vary substantially (176-1,158 km) and appear to reflect magnitude of gene flow, rather than extent of genetic differentiation. None of these birds was suggested as possible exemplars in the original description of the suture zone. Of the six species that were so identified, three have been surveyed to date, but none of those was found to be genetically differentiated.

Phylogeography of Pulsatilla cernua (Ranunculaceae), a grassland species, in Japan

The genetic diversity and structure of Pulsatilla cernua, a continental-grassland relict, were investigated using variations in chloroplast DNA (cpDNA) and microsatellites of nuclear DNA. In the analyses of three cpDNA regions, 17 haplotypes were found in 24 populations of P. cernua from Japan, Korea, and Russia. Although the route and time of migration between the continent of Asia and Japan could not be well resolved, the cpDNA haplotype network suggests the existence of several ancient lineages in Japan and a recent secondary migration from Japan to the continent. Microsatellite analyses did not indicate genetic structure among the Japanese populations, indicating the existence of gene flow across the distribution area until recently. These results indicate that the present fragmentation of P. cernua in Japan may reflect a rapid, recent reduction from a previously large, continuous distribution.

The Apennines as a cryptic Pleistocene refugium of the bark beetle Pityogenes chalcographus (Coleoptera: Curculionidae)

The Apennine Mountains in Italy are an important biogeographical region and of particular interest in phylogeographical research, because they have been a refugium during Pleistocene glaciation events for numerous European species. We performed a genetic study on the Eurasian bark beetle Pityogenes chalcographus (Linnaeus, 1760), focusing on two Apennine (Italian) and two Central European (Austrian) locations to assess the influence of the Apennines in the evolutionary history of the beetle, particularly during the Pleistocene. We analysed a part of the mitochondrial COI gene and a set of 5470 informative genome-wide markers to understand its biogeography. We found 75 distinct mitochondrial haplotypes, which are structured in three main clades. In general, the Apennine locations harbour a higher number of mitochondrial clades than Central European sites, with one specific clade exclusively detected in the Apennines. Analysis of our genome-wide, multi-locus dataset reveals a clustering of P. chalcographus by geography, with Italian individuals clearly separated from Austrian samples. Our data highlight the significance of the Apennines for the genetic diversity of P. chalcographus and support the hypothesis that this area was an important refugium during unfavourable conditions in the Pleistocene. We discuss additional life-history traits and processes that shaped the evolution of this widespread beetle.

Reconstruction of the Genetic History and the Current Spread of HIV-1 Subtype A in Germany

HIV-1 non-B infections have been increasing in Europe for several years. In Germany, subtype A belongs to the most abundant non-B subtypes showing an increasing prevalence of 8.3% among new infections in 2016. Here we trace the origin and examine the current spread of the German HIV-1 subtype A epidemic. Bayesian coalescence and birth-death analyses were performed with 180 German HIV-1 pol sequences and 528 related and publicly available sequences to reconstruct the population dynamics and fluctuations for each of the transmission groups. Our reconstructions indicate two distinct sources of the German subtype A epidemic, with an Eastern European and an Eastern African lineage both cocirculating in the country. A total of 13 German-origin clusters were identified; among these, 6 clusters showed recent activity. Introductions leading to further countrywide spread originated predominantly from Eastern Africa when introduced before 2005. Since 2005, however, spreading introductions have occurred exclusively within the Eastern European clade. Moreover, we observed changes in the main route of subtype A transmission. The beginning of the German epidemic (1985 to 1995) was dominated by heterosexual transmission of the Eastern African lineage. Since 2005, transmissions among German men who have sex with men (MSM) have been increasing and have been associated with the Eastern European lineage. Infections among people who inject drugs dominated between 1998 and 2005. Our findings on HIV-1 subtype A infections provide new insights into the spread of this virus and extend the understanding of the HIV epidemic in Germany.IMPORTANCE HIV-1 subtype A is the second most prevalent subtype worldwide, with a high prevalence in Eastern Africa and Eastern Europe. However, an increase of non-B infections, including subtype A infections, has been observed in Germany and other European countries. There has simultaneously been an increased flow of refugees into Europe and especially into Germany, raising the question of whether the surge in non-B infections resulted from this increased immigration or whether German transmission chains are mainly involved. This study is the first comprehensive subtype A study from a western European country analyzing in detail its phylogenetic origin, the impact of various transmission routes, and its current spread. The results and conclusions presented provide new and substantial insights for virologists, epidemiologists, and the general public health sector. In this regard, they should be useful to those authorities responsible for developing public health intervention strategies to combat the further spread of HIV/AIDS.

A test of trophic and functional island biogeography theory with the avifauna of a continental archipelago

The classical MacArthur-Wilson theory of island biogeography (TIB) emphasizes the role of island area and isolation in determining island biotas, but is neutral with respect to species differences that could affect community assembly and persistence. Recent extensions of island biogeography theory address how functional differences among species may lead to non-random community assembly processes and different diversity-area scaling patterns. First, the trophic TIB considers how diversity scaling varies across trophic position in a community, with species at higher trophic levels being most strongly influenced by island area. Second, further extensions have predicted how trait distributions, and hence functional diversity, should scale with area. Trait-based theory predicts richness-corrected functional diversity should be low on small islands but converge to null on larger islands. Conversely, competitive assembly predicts high diversity on small islands converging to null with increasing size. However, despite mounting interest in diversity-area relationships across different dimensions of diversity, these predictions derived from theory have not been extensively tested across taxa and island systems. Here, we develop and test predictions of the trophic TIB and extensions to functional traits, by examining the diversity-area relationship across multiple trophic ranks and dimensions of avian biodiversity in the Ryūkyū archipelago of Japan. We find evidence for a positive species- and phylogenetic diversity-area relationship, but functional diversity was not strongly affected by island area. Counter to the trophic TIB, we found no differences in the slopes of species-area relationships among trophic ranks, although slopes varied among trophic guilds at the same rank. We revealed differential assembly of trophic ranks, with evidence of trait-based assembly of intermediate predators but otherwise neutral community assembly. Our results suggest that niche space differs among trophic guilds of birds, but that differences are mostly not predicted by current extensions of island biogeography theory. While predicted patterns do not fit the empirical data well in this case, the development of such theory provides a useful framework to analyse island patterns from new perspectives. The application of empirical datasets such as ours should help provide a basis for developing further iterations of island biogeography theory.

Distinct YFV Lineages Co-circulated in the Central-Western and Southeastern Brazilian Regions From 2015 to 2018

The current outbreak of yellow fever virus (YFV) that is afflicting Brazil since the end of 2016 probably originated from a re-introduction of YFV from endemic areas into the non-endemic Southeastern Brazil. However, the lack of genomic sequences from endemic regions hinders the tracking of YFV's dissemination routes. We assessed the origin and spread of the ongoing YFV Brazilian outbreak analyzing a new set of YFV strains infecting humans, non-human primates (NHPs) and mosquitoes sampled across five Brazilian states from endemic and non-endemic regions between 2015 and 2018. We found two YFV sub-clade 1E lineages circulating in NHP from Goiás state (GO), resulting from independent viral introductions into the Araguaia tributary river basin: while one strain from 2017 clustered intermingled with Venezuelan YFV strains from 2000, the other YFV strains sampled in 2015 and 2017 clustered with sequences of the current YFV outbreak in the Brazilian Southeastern region (named YFV2015-2018 lineage), displaying the same molecular signature associated to the current YFV outbreak. After its introduction in GO at around mid-2014, the YFV2015-2018 lineage followed two paths of dissemination outside GO, originating two major YFV sub-lineages: (1) the YFVMG/ES/RJ sub-lineage spread sequentially from the eastern area of Minas Gerais state to Espírito Santo and then to Rio de Janeiro states, following the Southeast Atlantic basin; (2) the YFVMG/SP sub-lineage spread from the southwestern area of Minas Gerais to the metropolitan region of São Paulo state, following the Paraná basin. These results indicate the ongoing YFV outbreak in Southeastern Brazil originated from a dissemination event from GO almost 2 years before its recognition at the end of 2016. From GO this lineage was introduced in Minas Gerais state at least two times, originating two sub-lineages that followed different routes toward densely populated areas. The spread of YFV outside endemic regions for at least 4 years stresses the imperative importance of the continuous monitoring of YFV to aid decision-making for effective control policies aiming the increase of vaccination coverage to avoid the YFV transmission in densely populated urban centers.

Climate change, extinction, and Sky Island biogeography in a montane lizard

Around the world, many species are confined to "Sky Islands," with different populations in isolated patches of montane habitat. How does this pattern arise? One scenario is that montane species were widespread in lowlands when climates were cooler, and were isolated by local extinction caused by warming conditions. This scenario implies that many montane species may be highly susceptible to anthropogenic warming. Here, we test this scenario in a montane lizard (Sceloporus jarrovii) from the Madrean Sky Islands of southeastern Arizona. We combined data from field surveys, climate, population genomics, and physiology. Overall, our results support the hypothesis that this species' current distribution is explained by local extinction caused by past climate change. However, our results for this species differ from simple expectations in several ways: (a) their absence at lower elevations is related to warm winter temperatures, not hot summer temperatures; (b) they appear to exclude a low-elevation congener from higher elevations, not the converse; (c) they are apparently absent from many climatically suitable but low mountain ranges, seemingly "pushed off the top" by climates even warmer than those today; (d) despite the potential for dispersal among ranges during recent glacial periods (~18,000 years ago), populations in different ranges diverged ~4.5-0.5 million years ago and remained largely distinct; and (e) body temperatures are inversely related to climatic temperatures among sites. These results may have implications for many other Sky Island systems. More broadly, we suggest that Sky Island species may be relevant for predicting responses to future warming.

Admixture analyses and phylogeographic relationships reveal complete genetic distinctiveness of Polish farm and wild red foxes (Vulpes vulpes) and the North American origin of farm-bred individuals

A number of studies showed that many mtDNA haplotypes were shared among contemporary farm red foxes bred on different continents and the historical wild red foxes of North American origin. Therefore, in this study, the population genetic structure and phylogeographic relationships of Polish red foxes kept on fur farms and their wild conspecifics were investigated to assess the ancestry of the farm red foxes in Poland. A total of 330 tissue samples (200 from farm foxes and 130 from wild foxes) were used for the genetic analyses. Thirty microsatellite loci and two regions of mtDNA were used to assess the level of admixture between farm- and wild red foxes, to construct haplotype networks and create a phylogenetic tree. The genetic structure analysis clearly indicated two genetic clusters as being the most probable number of genetically distinct populations. The fixation index revealed a significant genetic distance between the farm- and wild red fox populations (F_ST  = 0.27, p < 0.05). Haplotype networks based on frequencies showing relationships between concatenated haplotypes of Polish farm- and wild red foxes and the constructed phylogenetic tree clearly indicated two genetically distinct groups. The results of this study provide strong evidence confirming the North American origin of red foxes bred on Polish farms and the genetic distinctiveness of both studied populations.

Cryptic diversity in Lithobateswarszewitschii (Amphibia, Anura, Ranidae)

Lithobateswarszewitschii is a species of ranid frog distributed from southern Honduras to Panama. This species suffered severe population declines at higher elevations (above 500 m a.s.l.) from the 1980s to early 1990s, but there is more recent evidence of recovery in parts of its range. Here we advocate for the status of Lithobateswarszewitschii as a candidate cryptic species complex based on sequence data from mitochondrial genes CO1 and 16S. Using concatenated phylogenies, nucleotide diversity (K2P-π), net between group mean distance (NBGMD) (π_net) and species delimitation methods, we further elucidate cryptic diversity within this species. All phylogenies display polyphyletic lineages within Costa Rica and Panama. At both loci, observed genetic polymorphism (K2P-π) is also high within and between geographic populations, surpassing proposed species threshold values for amphibians. Additionally, patterns of phylogeographic structure are complicated for this species, and do not appear to be explained by geographic barriers or isolation by distance. These preliminary findings suggest L.warszewitschii is a wide-ranging species complex. Therefore, we propose further research within its wider range, and recommend integrative taxonomic assessment is merited to assess species status.

Genetic turnovers and northern survival during the last glacial maximum in European brown bears

The current phylogeographic pattern of European brown bears (Ursus arctos) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid-Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter-specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought.

Genetic diversity and population structure of the blue jack mackerel Trachurus picturatus across its western distribution

Blue jack mackerel Trachurus picturatus collected at six sampling locations of the north-east Atlantic Ocean (Azores, Madeira, Canaries, and Matosinhos, Peniche and Portimão, mainland Portugal) and one location in the Mediterranean (Sicily), were used to examine the genetic structure of this species. Three mitochondrial gene regions (cytochrome c oxidase subunit I, cytochrome b and control region) were used to study the genetic structure of the species in Macaronesia, as well as to compare the genetic diversity of this region with published results from its eastern distribution. All markers indicated the absence of genetic structure among populations, with high indices of genetic diversity. These results suggest that the species went through a bottleneck event, followed by a recent population expansion. Moreover, the comparison with previously published results from the T. picturatus Mediterranean distribution suggests the existence of a single panmictic population throughout the species' full range. This was, however, an unexpected result since other methodologies have shown the presence of, at least, three different population-units in the NE Atlantic Ocean.

Full Bayesian Comparative Phylogeography from Genomic Data

A challenge to understanding biological diversification is accounting for community-scale processes that cause multiple, co-distributed lineages to co-speciate. Such processes predict non-independent, temporally clustered divergences across taxa. Approximate-likelihood Bayesian computation (ABC) approaches to inferring such patterns from comparative genetic data are very sensitive to prior assumptions and often biased toward estimating shared divergences. We introduce a full-likelihood Bayesian approach, ecoevolity, which takes full advantage of information in genomic data. By analytically integrating over gene trees, we are able to directly calculate the likelihood of the population history from genomic data, and efficiently sample the model-averaged posterior via Markov chain Monte Carlo algorithms. Using simulations, we find that the new method is much more accurate and precise at estimating the number and timing of divergence events across pairs of populations than existing approximate-likelihood methods. Our full Bayesian approach also requires several orders of magnitude less computational time than existing ABC approaches. We find that despite assuming unlinked characters (e.g., unlinked single-nucleotide polymorphisms), the new method performs better if this assumption is violated in order to retain the constant characters of whole linked loci. In fact, retaining constant characters allows the new method to robustly estimate the correct number of divergence events with high posterior probability in the face of character-acquisition biases, which commonly plague loci assembled from reduced-representation genomic libraries. We apply our method to genomic data from four pairs of insular populations of Gekko lizards from the Philippines that are not expected to have co-diverged. Despite all four pairs diverging very recently, our method strongly supports that they diverged independently, and these results are robust to very disparate prior assumptions.

Evolution and phylogeography analysis of diploid and polyploid Misgurnus anguillicaudatus populations across China

The origin and evolution of polyploid organisms have been extensively studied in plants, but this topic remains only partially understood in vertebrates, where polyploidy is relatively rare. In this study, we used Misgurnus anguillicaudatus, a fish that comprises five ploidy levels in nature, as a model animal to improve our understanding of biogeographic history and evolution of polyploid vertebrates. After collecting samples from different geographical populations in China, their ploidy levels were determined using flow cytometry. Two mitochondrial markers ( cytochrome b and control region) were then used for phylogeographic analyses to unravel the possible origins of diploids and tetraploids in China. The results showed that diploids have wider geographical distribution than tetraploids and triploids. There was no clear allopatric geographical range or boundary to divide diploid and polyploid populations. Rather, the analysis of mitochondrial DNA sequences indicated that tetraploids were autopolyploids, with lower genetic diversity than diploids. This suggests that tetraploids originated from sympatric diploids via multiple independent polyploidization events. Genetic structure patterns were similar between diploids and tetraploids, whereas complex genetic differentiation was found among different regions. The potential origin of M. anguillicaudatus was deduced to be in the Pearl River basin, which exhibited the highest nucleotide diversity and genetic differentiation. These findings provide insights into the evolution of polyploidy in vertebrates.

Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake

Comparing mitochondrial and genomic phylogenies is an essential tool for investigating speciation processes, because each genome carries different inheritance properties and evolutionary characteristics. Furthermore, mitonuclear discordance may arise from ecological adaptation, historic isolation, population size changes, and sex-biased dispersal. Closely related taxa are expected to experience gene flow; however, this may not be true for insular populations or populations isolated in refugia. The four-lined snake Elaphe quatuorlineata has a fragmented distribution, separating populations of the Italian and Balkan Peninsulas, whereas several insular Aegean populations of significantly smaller body size (Cyclades island group and Skyros Island, Greece) are currently considered distinct subspecies. We constructed the species-tree phylogeny of this species utilizing genome-wide single nucleotide polymorphisms and a gene-tree based on complete cytochrome b sequences, aiming to detect convergence and discrepancies between biparentally and maternally inherited genomes. Population structuring, phylogenetic patterns and migration events among geographically defined lineages supported our hypothesis of isolation in multiple sub-refugia. Where biogeographical barriers did not restrict subsequent dispersal, extensive genetic exchange occurred between mainland Balkan populations. This process has led to the mitochondrial sweep of an ancestral mitolineage that survived only in peripheral (East Greece) and insular populations (North Cyclades and Skyros). The Central Cyclades represent an ancient lineage for both molecular markers that emerged almost 3.3 Mya. Considering their distinct morphology, insular E. quatuorlineata populations should be the future focus of an extensive sampling, especially since the mitonuclear discordance observed in this species could be related to ecological adaptations, such as the island-dwarfism phenomenon.

Shedding Light on a Secretive Tertiary urodelean Relict: Hynobiid salamanders (Paradactylodon persicus s.l.) from Iran, Illuminated by Phylogeographic, Developmental and Transcriptomic Data

The Hyrcanian Forests present a unique Tertiary relict ecosystem, covering the northern Elburz and Talysh Ranges (Iran, Azerbaijan), a poorly investigated, unique biodiversity hotspot with many cryptic species. Since the 1970s, two nominal species of Urodela, Hynobiidae, Batrachuperus (later: Paradactylodon) have been described: Paradactylodon persicus from northwestern and P. gorganensis from northeastern Iran. Although P. gorganensis has been involved in studies on phylogeny and development, there is little data on the phylogeography, systematics, and development of the genus throughout the Hyrcanian Forests; genome-wide resources have been entirely missing. Given the huge genome size of hynobiids, making whole genome sequencing hardly affordable, we aimed to publish the first transcriptomic resources for Paradactylodon from an embryo and a larva (9.17 Gb RNA sequences; assembled to 78,918 unigenes). We also listed 32 genes involved in vertebrate sexual development and sex determination. Photographic documentation of the development from egg sacs across several embryonal and larval stages until metamorphosis enabled, for the first time, comparison of the ontogeny with that of other hynobiids and new histological and transcriptomic insights into early gonads and timing of their differentiation. Transcriptomes from central Elburz, next-generation sequencing (NGS) libraries of archival DNA of topotypic P. persicus, and GenBank-sequences of eastern P. gorganensis allowed phylogenetic analysis with three mitochondrial genomes, supplemented by PCR-amplified mtDNA-fragments from 17 museum specimens, documenting <2% uncorrected intraspecific genetic distance. Our data suggest that these rare salamanders belong to a single species P. persicus s.l. Humankind has a great responsibility to protect this species and the unique biodiversity of the Hyrcanian Forest ecosystems.

Deep Phylogeographic Breaks in Magpie Pica pica Across the Holarctic: Concordance with Bioacoustics and Phenotypes

We examined sequences of the mitochondrial control region in magpies (Pica pica) from the entire distribution range and found deep genetic splits into four major lineages: (1) group West (Europe-Siberia), (2) group East (southern Far East), (3) P. p. mauritanica (North Africa), and (4) P. p. hudsonia (North America). These lineages show a geographic pattern corresponding to known subspecies or subspecies groups. Genetic variation within the widely-distributed group West is low and neutrality tests supported a recent expansion scenario. The haplotypes from Kamchatka, representing a separated sublineage with clear affinity to the European-Siberian group, are almost identical, implying a recent bottleneck. Group East contained two subclades without clear geographic pattern, presumably due to admixing of populations that had diverged in Pleistocene refuges. The homogeneity of the Kyushu population supports historical reports of introduction of the species from Korea. In contrast, the high variation in the recently established Hokkaido population may reflect an ongoing invasion from several populations of the Far Eastern mainland. Bioacoustic data based on chatter call differentiate groups of subspecies and reflect phylogeographic patterns, i.e., mitochondrial lineages. Furthermore, we report the fast spreading of P. p. jankowskii towards the west along the upper Amur River, and a slower shifting of P. p. leucoptera in the opposite direction thus yielding a new contact zone. Overall, our data support a scenario of divergence in geographic isolation, but the ongoing expansion of distribution ranges may lead to major changes in phylogeographic patterns.

Human-Origin Influenza A(H3N2) Reassortant Viruses in Swine, Southeast Mexico

The genetic diversity of influenza A viruses circulating in swine in Mexico complicates control efforts in animals and presents a threat to humans, as shown by influenza A(H1N1)pdm09 virus. To describe evolution of swine influenza A viruses in Mexico and evaluate strains for vaccine development, we sequenced the genomes of 59 viruses and performed antigenic cartography on strains from 5 regions. We found that genetic and antigenic diversity were particularly high in southeast Mexico because of repeated introductions of viruses from humans and swine in other regions in Mexico. We identified novel reassortant H3N2 viruses with genome segments derived from 2 different viruses that were independently introduced from humans into swine: pandemic H1N1 viruses and seasonal H3N2 viruses. The Mexico swine viruses are antigenically distinct from US swine lineages. Protection against these viruses is unlikely to be afforded by US virus vaccines and would require development of new vaccines specifically targeting these diverse strains.

Divergent mitochondrial lineages arose within a large, panmictic population of the Savannah sparrow (Passerculus sandwichensis)

Unusual patterns of mtDNA diversity can reveal interesting aspects of a species' biology. However, making such inferences requires discerning among the many alternative scenarios that could underlie any given mtDNA pattern. Next-generation sequencing methods provide large, multilocus data sets with increased power to resolve unusual mtDNA patterns. A mtDNA-based phylogeography of the Savannah sparrow (Passerculus sandwichensis) previously identified two sympatric, but divergent (~2%) clades within the nominate subspecies group and a third clade that consisted of birds sampled from northwest Mexico. We revisited the phylogeography of this species using a population genomic data set to resolve the processes leading to the evolution of sympatric and divergent mtDNA lineages. We identified two genetic clusters in the genomic data set corresponding to (a) the nominate subspecies group and (b) northwestern Mexico birds. Following divergence, the nominate clade maintained a large, stable population, indicating that divergent mitochondrial lineages arose within a panmictic population. Simulations based on parameter estimates from this model further confirmed that this demographic history could produce observed levels of mtDNA diversity. Patterns of divergent, sympatric mtDNA lineages are frequently interpreted as admixture of historically isolated lineages. Our analyses reject this interpretation for Savannah sparrows and underscore the need for genomic data sets to resolve the evolutionary mechanisms behind anomalous, locus-specific patterns.

Patterns of mtDNA variation reveal complex evolutionary history of relict and endangered peat bog pine (Pinus uliginosa)

Estimates of genetic differentiation at intra- and interspecific level are often hindered by the lack of suitable molecular markers. Low phylogeographic resolution limits development of appropriate conservation strategies especially in case of endangered forest tree species with small and disjunct distribution. In this study, we assessed fine-scale genetic structure of relict and endangered peat bog pine (Pinus uliginosa) and two other closely related European pine species (Pinus mugo and Pinus uncinata) using a set of 15 newly developed maternally inherited and seed-mediated mitochondrial DNA (mtDNA) markers and two previously known polymorphic mtDNA regions (nad1, nad7). Three main groups, corresponding in general to three investigated species were revealed in the haplotype network analysis. However, only P. uncinata was clearly distinct at all levels of analysis, whereas great genetic similarity and haplotype sharing was observed between P. uliginosa and P. mugo. Strong phylogeographic structure was found in P. uliginosa that showed high differentiation at relatively short geographical distance among populations and the existence of mitochondrial lineages of different evolutionary history. Hybridization with other pine species has likely contributed to genetic differentiation of P. uliginosa as indicated by contemporary distribution of mtDNA haplotypes. The research emphasizes the importance of accurate assessments of genetic structure of endangered species with complex evolutionary history for development of efficient conservation strategies.

The role of glacial-interglacial climate change in shaping the genetic structure of eastern subterranean termites in the southern Appalachian Mountains, USA

The eastern subterranean termite, Reticulitermes flavipes, currently inhabits previously glaciated regions of the northeastern U.S., as well as the unglaciated southern Appalachian Mountains and surrounding areas. We hypothesized that Pleistocene climatic fluctuations have influenced the distribution of R. flavipes, and thus the evolutionary history of the species. We estimated contemporary and historical geographic distributions of R. flavipes by constructing Species Distribution Models (SDM). We also inferred the evolutionary and demographic history of the species using mitochondrial (cytochrome oxidase I and II) and nuclear (endo-beta-1,4-glucanase) DNA sequence data. To do this, genetic populations were delineated using Bayesian spatial-genetic clustering, competing hypotheses about population divergence were assessed using approximate Bayesian computation (ABC), and changes in population size were estimated using Bayesian skyline plots. SDMs identified areas in the north with suitable habitat during the transition from the Last Interglacial to the Last Glacial Maximum, as well as an expanding distribution from the mid-Holocene to the present. Genetic analyses identified three geographically cohesive populations, corresponding with northern, central, and southern portions of the study region. Based on ABC analyses, divergence between the Northern and Southern populations was the oldest, estimated to have occurred 64.80 thousand years ago (kya), which corresponds with the timing of available habitat in the north. The Central and Northern populations diverged in the mid-Holocene, 8.63 kya, after which the Central population continued to expand. Accordingly, phylogeographic patterns of R. flavipes in the southern Appalachians appear to have been strongly influenced by glacial-interglacial climate change.

OPEN RESEARCH BADGES

This article has been awarded Open Materials, Open Data Badges. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.5061/dryad.5hr7f31.