Next-generation phylogeography of the cockle Cerastoderma glaucum: Highly heterogeneous genetic differentiation in a lagoon species

AIM

Coastal lagoons form an intriguing example of fragmented marine habitats. Restricted gene flow among isolated populations of lagoon species may promote their genetic divergence and may thus provide a first step toward speciation. In the present study, the population genetic structure of the lagoon cockle Cerastoderma glaucum has been investigated to clarify the complex phylogeographic pattern found in previous studies, to localize major genetic breaks, and to discuss their origin and maintenance.

LOCATION

The Atlantic and Mediterranean coasts, including the Baltic, North Sea, and Black Sea.

METHODS

A total of 204 C. glaucum individuals from 14 populations were genotyped using restriction site-associated DNA sequencing (RADseq). The genetic diversity, divergence, and structure were analyzed using genome-wide single nucleotide polymorphisms (SNPs). Phylogenetic relationships were inferred under a coalescent model using svdquartets.

RESULTS

The RADseq approach allowed inferring phylogeographic relationships with an unprecedented resolution. Three deeply divergent lineages were identified within C. glaucum that are separated by many genetic barriers: one lineage in the Aegean-Black Sea region, one in the Ionian Sea, and the last one widely distributed from the Western Mediterranean to the Baltic Sea. The nested branching pattern displayed on the species tree largely agrees with the likely scenario of C. glaucum postglacial expansion from the Mediterranean to the Baltic Sea.

MAIN CONCLUSION

The genetic differentiations between geographically separated lagoons proved to be strong, highlighting the evolutionary influence of these naturally fragmented habitats. The postglacial expansion created complex patterns of spatial segregation of genetic diversity with allele frequency gradients in many outlier loci, but also discrepancies between the nuclear and mitochondrial genetic markers that probably arose from genetic surfing of mitochondrial variation.

Phylogeographic patterns and conservation implications of the endangered Chinese giant salamander

Understanding genetic diversity patterns of endangered species is an important premise for biodiversity conservation. The critically endangered salamander Andrias davidianus, endemic to central and southern mainland in China, has suffered from sharp range and population size declines over the past three decades. However, the levels and patterns of genetic diversity of A. davidianus populations in wild remain poorly understood. Herein, we explore the levels and phylogeographic patterns of genetic diversity of wild-caught A. davidianus using larvae and adult collection with the aid of sequence variation in (a) the mitochondrial DNA (mtDNA) fragments (n = 320 individuals; 33 localities), (b) 19 whole mtDNA genomes, and (c) nuclear recombinase activating gene 2 (RAG2; n = 88 individuals; 19 localities). Phylogenetic analyses based on mtDNA datasets uncovered seven divergent mitochondrial clades (A-G), which likely originated in association with the uplifting of mountains during the Late Miocene, specific habitat requirements, barriers including mountains and drainages and lower dispersal ability. The distributions of clades were geographic partitioned and confined in neighboring regions. Furthermore, we discovered some mountains, rivers, and provinces harbored more than one clades. RAG2 analyses revealed no obvious geographic patterns among the five alleles detected. Our study depicts a relatively intact distribution map of A. davidianus clades in natural species range and provides important knowledge that can be used to improve monitoring programs and develop a conservation strategy for this critically endangered organism.

Genetic diversity and phylogeography of Daphnia similoides sinensis located in the middle and lower reaches of the Yangtze River

Geographical patterns, climate, and environmental change have important influences on the distribution and spread of aquatic organisms. However, the relationships between the geographical pattern and phylogenetics of Daphnia as well as environmental change are not well known. The genetic diversity and phylogeography of seven D. similoides sinensis populations located in the middle and lower reaches of the Yangtze River were investigated based on the combination of mitochondrial (COI gene) and nuclear (14 microsatellite primers) markers. Based on the mitochondrial gene markers, D. similoides sinensis from the middle and lower reaches of the Yangtze River had one ancestral haplotype and two evolutionary clades. In addition, D. similoides sinensis population deviated from neutral evolution, showing signs of a bottleneck effect followed by population expansion. Based on the microsatellite markers, the seven D. similoides sinensis populations formed three main groups. The dendrogram (NJ/ME) showed that D. similoides sinensis based on the mitochondrial genes marker were obviously clustered two main clades, whereas there were three clades based on the microsatellite markers. Our results suggested that the habitat fragmentation due to the barrier of the dams and sluices promoted the genetic differentiation and phylogeography of D. similoides sinensis populations in the middle and lower reaches of the Yangtze River.

Possible female philopatry of the smooth hammerhead shark Sphyrna zygaena revealed by genetic structure patterns

We assessed the spatial pattern of genetic structure of smooth hammerhead shark Sphyrna zygaena in 10 localities from the Northern Mexican Pacific. A total of 35 haplotypes were identified in 129 sequences of the mtDNA control region. The results showed slight but significant genetic structure among localities (Φ_ST = 0.044, P < 0.001). In addition, the localities with highest number of juveniles were genetically different (Φ_ST = 0.058, P < 0.024), which may be representative of nursery areas. The genetic differentiation pattern can be associated to female philopatry and preference for particular birthing sites. Finally, historical demography shows that S. zygaena populations present a recent demographic expansion that occurred during glacial events in the late Pleistocene to early Holocene.

Phylogeography of the red-bellied lizard, Darevskia parvula in Turkey

The groups of red-bellied lizards had a small distribution area in the Pontic zone. The several studies performed on these lizard groups are based on taxonomy and systematics. Although there were several taxonomic or systematic researches on some species of this group, the phylogeographical pattern and species disturbing boundaries of this group is still not clear. In the present study, we aimed to resolve the taxonomic and phylogenetic relationships of the red-bellied lizards in Turkey, based on two combined mitochondrial gene fragments and one protein-coding nuclear gene (rag1). Also, we evaluated ecological niches differentiations among red-bellied lizard groups. The mitochondrial DNA genes were found to be highly polymorphic in this group. One hundred and one variable nucleotide sites were detected on the combined gene sequences. According to phylogenetic trees based on the maximum likelihood (ML) and Bayesian inference (BI), the red-bellied lizards group have three species groups; Darevskia parvula, D. adjarica and unnamed Darevskia sp. (candidate species for Darevskia genus). This situation was supported by high bootstrap and posterior probability values in the trees of mitochondrial DNA gene fragments. However, no genetic variation was detected according to nuclear DNA (rag1) sequence. Because the species groups have no overlaps in terms of their ecological niches, ecological niche modelling (ENM) results revealed differences among the groups of D. parvula, D. adjarica, and unnamed Darevskia sp. Besides, we detected no geographical overlaps among three species groups, since there were geographical isolation zones among the species groups of red-bellied lizard.

Mountains as Evolutionary Arenas: Patterns, Emerging Approaches, Paradigm Shifts, and Their Implications for Plant Phylogeographic Research in the Tibeto-Himalayan Region

Recently, the "mountain-geobiodiversity hypothesis" (MGH) was proposed as a key concept for explaining the high levels of biodiversity found in mountain systems of the Tibeto-Himalayan region (THR), which comprises the Qinghai-Tibetan Plateau, the Himalayas, and the biodiversity hotspot known as the "Mountains of Southwest China" (Hengduan Mountains region). In addition to the MGH, which covers the entire life span of a mountain system, a complementary concept, the so-called "flickering connectivity system" (FCS), was recently proposed for the period of the Quaternary. The FCS focuses on connectivity dynamics in alpine ecosystems caused by the drastic climatic changes during the past ca. 2.6 million years, emphasizing that range fragmentation and allopatric speciation are not the sole factors for accelerated evolution of species richness and endemism in mountains. I here provide a review of the current state of knowledge concerning geological uplift, Quaternary glaciation, and the main phylogeographic patterns ("contraction/recolonization," "platform refugia/local expansion," and "microrefugia") of seed plant species in the THR. In addition, I make specific suggestions as to which factors future avenues of phylogeographic research should take into account based on the fundamentals presented by the MGH and FCS, and associated complementary paradigm shifts.

Molecular Phylogeography and Ecological Niche Modeling of Sibbaldia procumbens s.l. (Rosaceae)

The phylogeographical analysis and ecological niche modeling (ENM) of the widely distributed Northern Hemisphere Sibbaldia procumbens s.l. can help evaluate how tectonic motion and climate change helped shape the current distribution patterns of this species. Three chloroplast regions (the atpI-atpH and trnL-trnF intergenic spacers and the trnL intron) were obtained from 332 (156 from present study and 176 from the previous study) individuals of S. procumbens s.l. An unrooted haplotype network was constructed using the software NETWORK, while BEAST was used to estimate the divergence times among haplotypes. ENM was performed by MAXENT to explore the historical dynamic distribution of S. procumbens s.l. The haplotype distribution demonstrates significant phylogeographical structure (N_ST > G_ST; P < 0.01). The best partitioning of genetic diversity by SAMOVA produced three groups, while the time to the most recent common ancestor of all haplotypes was estimated to originate during the Miocene, with most of the haplotype diversity having occurred during the Quaternary. The MAXENT analysis showed S. procumbens s.l. had a wider distribution range during the last glacial maximum and a narrower distribution range during the last interglacial, with predictions into the future showing the distribution range of S. procumbens s.l. shrinking.

Phylogeography of the Assassin Bug Sphedanolestes impressicollis in East Asia Inferred From Mitochondrial and Nuclear Gene Sequences

The assassin bug, Sphedanolestes impressicollis (Hemiptera: Reduviidae), is widely distributed in East Asia. It is an ideal model for evaluating the effects of climatic fluctuation and geographical events on the distribution patterns of East Asian reduviids. Here, we used two mitochondrial genes and one nuclear gene to investigate the phylogeographic pattern of the assassin bug based on comprehensive sampling in China, Japan, South Korea, Vietnam, and Laos. High levels of genetic differentiation were detected among the geographic populations classified into the northern and southern groups. A significant correlation was detected between genetic and geographical distances. The East China Sea land bridge served as a “dispersal corridor” during Pleistocene glaciation. The estimated divergence time indicated that the northern group may have separated from the eastern Chinese populations when the sea level rapidly rose during the “Ryukyu Coral Sea Stage” and the East China Sea land bridge was completely submerged. Demographic history and ecological niche modeling suggested that appropriate climatic conditions may have accounted for the rapid spread across the Korean Peninsula and Japan during the late Pleistocene. Our study underscores the pivotal roles of the Pleistocene sea level changes and climatic fluctuations in determining the distribution patterns of East Asian reduviids.

Genetic diversity and maternal origin of Northeast African and South American donkey populations

To investigate the mtDNA variation and origin of maternal lineages in South American donkeys and to reassess the domestication of donkeys in northeast Africa, we analyzed sequences (489 bp of the D-loop) from 323 domestic donkeys sampled from Peru, Brazil, Ethiopia and Egypt. Altogether, the 323 sequences displayed 53 different haplotypes (45 in Ethiopia, 14 in Egypt, eight in Peru and six in Brazil). Among the four populations, Egyptian donkeys possessed the highest haplotype diversity (0.910 ± 0.032), followed by Brazilian donkeys (0.879 ± 0.060). The Clade I haplotypes dominated in Peruvian donkeys (65%), whereas Clade II haplotypes dominated in Brazilian donkeys (67%). Estimates of F_ST values showed a high genetic differentiation between Peruvian and Brazilian donkey populations (F_ST  = 0.4066), which could be explained by the complex introduction history of South American donkeys. Phylogeographic analysis indicates that northeast Africa could be the most probable domestication center for Clade I donkeys. Analysis of molecular variance confirmed a weak genetic structure in domestic donkey populations among four continents (Europe, Asia, Africa and South America).

Calochortus gunnisonii furthers evidence for the complex genetic legacy of historical climate change in the southern Rocky Mountains

PREMISE OF THE STUDY

Climate cycles of the Quaternary have impacted plants at a global scale, leaving behind a complex genetic legacy. Species of the northern Rocky Mountains of North America were exposed to more uniform glacial patterns than the central and southern ranges, where synergistic relationships between temperature and precipitation caused differences in the timing and extent of glacier onset. We examined the genetic impacts of climate oscillations on Calochortus gunnisonii (Liliaceae) in the central and southern Rocky Mountains.

METHODS

Populations were sampled from disjunct mountain ranges across the basins of Wyoming and northern and central Colorado. Allelic data from nuclear microsatellites and plastid sequences (trnV-ndhC, petA-psbJ, and rpl16) were used to examine patterns of genetic structure between and among populations along the southern Rocky Mountain corridor.

KEY RESULTS

We infer considerable population structure concordant with mountain range of origin. Clustering analysis supports separate north and south genetic clusters on either side of major basins in Wyoming, suggesting that populations were maintained in two distinct refugia. Additionally, populations within the Sierra Madre Range of southern Wyoming show localized, divergent genetic signal indicative of a third potential glacial refugium. By contrast, recent genetic admixture is observed in the Laramie, Medicine Bow, and Front ranges, where population expansion from glacial refugia has likely occurred.

CONCLUSIONS

We conclude that during climate cycles of the Quaternary, C. gunnisonii experienced periods of population expansion and reduction, habitat fragmentation, isolation in three or more refugia, and admixture mirroring genetic impacts of other southern Rocky Mountains organisms.

The speciation history of northern- and southern-sourced Eranthis (Ranunculaceae) species on the Korean peninsula and surrounding areas

The temporal and spatial origins and evolution of the genus Eranthis have not been previously studied. We investigated the speciation and establishment histories of four Eranthis species: Eranthis byunsanensis, E. pungdoensis, E. stellata, and E. pinnatifida. The sampling localities were Korea, Japan, Jilin in China, and the area near Vladivostok in Primorskiy, Russia. We used 12 chloroplast microsatellite loci (n = 935 individuals) and two chloroplast noncoding regions (rpl16 intron, petL-psbE intergenic spacer; n = 33 individuals). The genetic diversity, genetic structure, phylogenetic relationships of the four species were analyzed, and their ancestral areas were reconstructed. The high genetic diversity of the Jeju island population of E. byunsanensis and Russian populations of E. stellata indicated these species' northward and southward dispersal, respectively. The genetic structure analyses suggest that the populations in these four species have limited geographical structure, except for the Chinese E. stellata population (SCP). The phylogenetic analyses suggest that E. byunsanensis and E. pinnatifida are sister species and that Chinese SCP may not belong to E. stellata. The ancestral area reconstruction revealed that the most recent common ancestor of the four species existed in the current Chinese habitat of E. stellata. This study shows that E. byunsanensis and E. pinnatifida originated from a southern Eranthis species and speciated into their current forms near Jeju island and near western regions of Japan, respectively, during the Miocene. E. stellata may have dispersed southward on and near the Korean peninsula, though its specific origin remains unclear. Interestingly, the Chinese E. stellata population SCP suggests that the Chinese population might be most ancient among all the four Eranthis species. E. pungdoensis may have allopatrically speciated from E. byunsanensis during the Holocene. The Korean peninsula and the surrounding areas can be considered interesting regions which provide the opportunity to observe both northern- and southern-sourced Eranthis species.

Species distribution modelling and cpSSR reveal population history of the Neotropical annual herb Tithonia rotundifolia (Asteraceae)

The impacts of the historical geologic and climatic events on the diversity and genetic structure of Neotropical taxa have recently become a subject of study. However, annual plants associated with tropical dry forests remain under-studied. The exploration of additional taxa in contrasting environments will improve the current understanding of responses of the Neotropical biota to these events. Here, we explore the species distribution and geographic structure of the annual herb Tithonia rotundifolia. We sampled 175 individuals from 19 populations of T. rotundifolia. Species distribution modelling and six microsatellite chloroplast loci were used to infer its population history. We identified areas of historical climate suitability and then tested if there is genetic structuring among these areas. Haplotypes showed strong phylogeographic structure. Historical climatic suitability areas were found along the Pacific coast; however, a gap was found at the Isthmus of Tehuantepec (IT). Although Bayesian analysis showed population structuring, amova revealed that the IT is not its main driver. Instead, a subdivision into a higher number of regions had higher F_CT values. Also, populations to the east of the IT showed evidence of recent population expansion and migration in a south-north direction. Pleistocene climate fluctuations partially explain the geographic structure of T. rotundifolia. However, life-history characteristics such as limited seed dispersal and the patchy distribution of suitable habitats explain the high haplotype diversity and population sub-structuring and diversity. Lastly, the absence of geographic structure of some haplotypes may indicate long-distance dispersal, or hybridisation with the closely related T. tubaeformis.

Increasing importance of European lineages in seeding the hepatitis C virus subtype 1a epidemic in Spain

BackgroundReducing the burden of the hepatitis C virus (HCV) requires large-scale deployment of intervention programmes, which can be informed by the dynamic pattern of HCV spread. In Spain, ongoing transmission of HCV is mostly fuelled by people who inject drugs (PWID) infected with subtype 1a (HCV1a).AimOur aim was to map how infections spread within and between populations, which could help formulate more effective intervention programmes to halt the HCV1a epidemic in Spain.MethodsEpidemiological links between HCV1a viruses from a convenience sample of 283 patients in Spain, mostly PWID, collected between 2014 and 2016, and 1,317, 1,291 and 1,009 samples collected abroad between 1989 and 2016 were reconstructed using sequences covering the NS3, NS5A and NS5B genes. To efficiently do so, fast maximum likelihood-based tree estimation was coupled to a flexible Bayesian discrete phylogeographic inference method.ResultsThe transmission network structure of the Spanish HCV1a epidemic was shaped by continuous seeding of HCV1a into Spain, almost exclusively from North America and European countries. The latter became increasingly relevant and have dominated in recent times. Export from Spain to other countries in Europe was also strongly supported, although Spain was a net sink for European HCV1a lineages. Spatial reconstructions showed that the epidemic in Spain is diffuse, without large, dominant within-country networks.ConclusionTo boost the effectiveness of local intervention efforts, concerted supra-national strategies to control HCV1a transmission are needed, with a strong focus on the most important drivers of ongoing transmission, i.e. PWID and other high-risk populations.

Genetic diversity and drivers of dwarfism in extinct island emu populations

Australia's iconic emu (Dromaius novaehollandiae novaehollandiae) is the only living representative of its genus, but fossil evidence and reports from early European explorers suggest that three island forms (at least two of which were dwarfs) became extinct during the nineteenth century. While one of these-the King Island emu-has been found to be conspecific with Australian mainland emus, little is known about how the other two forms-Kangaroo Island and Tasmanian emus-relate to the others, or even the size of Tasmanian emus. We present a comprehensive genetic and morphological analysis of Dromaius diversity, including data from one of the few definitively genuine Tasmanian emu specimens known. Our genetic analyses suggest that all the island populations represent sub-populations of mainland Dnovaehollandiae Further, the size of island emus and those on the mainland appears to scale linearly with island size but not time since isolation, suggesting that island size-and presumably concomitant limitations on resource availability-may be a more important driver of dwarfism in island emus, though its precise contribution to emu dwarfism remains to be confirmed.

Incorporating sampling uncertainty in the geospatial assignment of taxa for virus phylogeography

Discrete phylogeography using software such as BEAST considers the sampling location of each taxon as fixed; often to a single location without uncertainty. When studying viruses, this implies that there is no possibility that the location of the infected host for that taxa is somewhere else. Here, we relaxed this strong assumption and allowed for analytic integration of uncertainty for discrete virus phylogeography. We used automatic language processing methods to find and assign uncertainty to alternative potential locations. We considered two influenza case studies: H5N1 in Egypt; H1N1 pdm09 in North America. For each, we implemented scenarios in which 25 per cent of the taxa had different amounts of sampling uncertainty including 10, 30, and 50 per cent uncertainty and varied how it was distributed for each taxon. This includes scenarios that: (i) placed a specific amount of uncertainty on one location while uniformly distributing the remaining amount across all other candidate locations (correspondingly labeled 10, 30, and 50); (ii) assigned the remaining uncertainty to just one other location; thus ‘splitting’ the uncertainty among two locations (i.e. 10/90, 30/70, and 50/50); and (iii) eliminated uncertainty via two predefined heuristic approaches: assignment to a centroid location (CNTR) or the largest population in the country (POP). We compared all scenarios to a reference standard (RS) in which all taxa had known (absolutely certain) locations. From this, we implemented five random selections of 25 per cent of the taxa and used these for specifying uncertainty. We performed posterior analyses for each scenario, including: (a) virus persistence, (b) migration rates, (c) trunk rewards, and (d) the posterior probability of the root state. The scenarios with sampling uncertainty were closer to the RS than CNTR and POP. For H5N1, the absolute error of virus persistence had a median range of 0.005–0.047 for scenarios with sampling uncertainty—(i) and (ii) above—versus a range of 0.063–0.075 for CNTR and POP. Persistence for the pdm09 case study followed a similar trend as did our analyses of migration rates across scenarios (i) and (ii). When considering the posterior probability of the root state, we found all but one of the H5N1 scenarios with sampling uncertainty had agreement with the RS on the origin of the outbreak whereas both CNTR and POP disagreed. Our results suggest that assigning geospatial uncertainty to taxa benefits estimation of virus phylogeography as compared to ad-hoc heuristics. We also found that, in general, there was limited difference in results regardless of how the sampling uncertainty was assigned; uniform distribution or split between two locations did not greatly impact posterior results. This framework is available in BEAST v.1.10. In future work, we will explore viruses beyond influenza. We will also develop a web interface for researchers to use our language processing methods to find and assign uncertainty to alternative potential locations for virus phylogeography.

Species delimitation in endangered groundwater salamanders: Implications for aquifer management and biodiversity conservation

Groundwater-dependent species are among the least-known components of global biodiversity, as well as some of the most vulnerable because of rapid groundwater depletion at regional and global scales. The karstic Edwards-Trinity aquifer system of west-central Texas is one of the most species-rich groundwater systems in the world, represented by dozens of endemic groundwater-obligate species with narrow, naturally fragmented distributions. Here, we examine how geomorphological and hydrogeological processes have driven population divergence and speciation in a radiation of salamanders (Eurycea) endemic to the Edwards-Trinity system using phylogenetic and population genetic analysis of genome-wide DNA sequence data. Results revealed complex patterns of isolation and reconnection driven by surface and subsurface hydrology, resulting in both adaptive and nonadaptive population divergence and speciation. Our results uncover cryptic species diversity and refine the borders of several threatened and endangered species. The US Endangered Species Act has been used to bring state regulation to unrestricted groundwater withdrawals in the Edwards (Balcones Fault Zone) Aquifer, where listed species are found. However, the Trinity and Edwards-Trinity (Plateau) aquifers harbor additional species with similarly small ranges that currently receive no protection from regulatory programs designed to prevent groundwater depletion. Based on regional climate models that predict increased air temperature, together with hydrologic models that project decreased springflow, we conclude that Edwards-Trinity salamanders and other codistributed groundwater-dependent organisms are highly vulnerable to extinction within the next century.

Candida auris Infection Leading to Nosocomial Transmission, Israel, 2017

A patient transferred from South Africa to Israel acquired a Candida auris infection. Phylogenetic analysis showed resemblance of C. auris to isolates from South Africa but not Israel, suggesting travel-associated infection. C. auris infection occurred weeks later in another patient at the same hospital, suggesting prolonged environmental persistence.

Mycobacterium ulcerans Population Genomics To Inform on the Spread of Buruli Ulcer across Central Africa

Buruli ulcer is a neglected tropical disease of skin and subcutaneous tissue caused by infection with the pathogen Mycobacterium ulcerans Many critical issues for disease control, such as understanding the mode of transmission and identifying source reservoirs of M. ulcerans, are still largely unknown. Here, we used genomics to reconstruct in detail the evolutionary trajectory and dynamics of M. ulcerans populations at a central African scale and at smaller geographical village scales. Whole-genome sequencing (WGS) data were analyzed from 179 M. ulcerans strains isolated from all Buruli ulcer foci in the Democratic Republic of the Congo, The Republic of Congo, and Angola that have ever yielded positive M. ulcerans cultures. We used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our phylogeographic analysis revealed one almost exclusively predominant sublineage of M. ulcerans that arose in Central Africa and proliferated in its different regions of endemicity during the Age of Discovery. We observed how the best sampled endemic hot spot, the Songololo territory, became an area of endemicity while the region was being colonized by Belgium (1880s). We furthermore identified temporal parallels between the observed past population fluxes of M. ulcerans from the Songololo territory and the timing of health policy changes toward control of the Buruli ulcer epidemic in that region. These findings suggest that an intervention based on detecting and treating human cases in an area of endemicity might be sufficient to break disease transmission chains, irrespective of other reservoirs of the bacterium.IMPORTANCE Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans The disease is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Currently, the major hurdles facing disease control are incomplete understandings of both the mode of transmission and environmental reservoirs of M. ulcerans As decades of spasmodic environmental sampling surveys have not brought us much closer to overcoming these hurdles, the Buruli ulcer research community has recently switched to using comparative genomics. The significance of our research is in how we used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our approach shows that it might be possible to use bacterial population genomics to assess the impact of health interventions, providing valuable feedback for managers of disease control programs in areas where health surveillance infrastructure is poor.

Analysis of the Genetic Diversity and Structure of the Grey-Faced Buzzard (Butastur indicus) in Japan, Based on mtDNA

The grey-faced buzzard (Butastur indicus) is a raptor that inhabits East Asia, including Japan. Because the number of individuals has decreased by 75% over the last 40 years, this species is classified as vulnerable (VU) in Japan. In the present study, wesought to reveal the genetic structure of the Japanese grey-faced buzzard population at several breeding sites, and to assess the levels of genetic diversity within the Japanese population. We sequenced 555 bp of the mitochondrial DNA of 96 individuals sampled during the breeding season at 18 sites, and 11 individuals sampled during the winter season at one site. In total, 21 variable sites were found in the control region, and we detected 26 haplotypes among the 107 individuals. Fukuoka represented the core breeding area for grey-faced buzzards, as half of all haplotypes were detected there. Four unique haplotypes were detected in the overwintering area. The results of the network and mismatch distribution analyses indicated that the grey-faced buzzard has not experienced a genetic bottleneck in the past, but did experience recent population expansion. In addition, comparisons with other raptors revealed rich genetic diversity in the grey-faced buzzard population. Our results indicate that conservation of both breeding and wintering areas is important for the protection of the grey-faced buzzard.

Extensive secondary contact among three glacial lineages of Arctic Char (Salvelinus alpinus) in Labrador and Newfoundland

AIM

The Pleistocene glaciation event prompted the allopatric divergence of multiple glacial lineages of Arctic char (Salvelinus alpinus), some of which have come into secondary contact upon their recolonization of the Holarctic. While three glacial lineages (Arctic, Atlantic, and Acadian) are known to have recolonized the western Atlantic, the degree of overlap of these three lineages is largely unknown. We sought to determine the distribution of these three glacial lineages in Labrador and Newfoundland at a fine spatial scale to assess their potential for introgression and their relative contribution to local fisheries.

LOCATION

Labrador and Newfoundland, Canada.

METHODS

We sequenced a portion of the D-loop region in over 1,000 Arctic char (S. alpinus) samples from 67 locations across Labrador and Newfoundland.

RESULTS

Within Labrador, the Arctic and Atlantic lineages were widespread. Two locations (one landlocked and one with access to the sea) also contained individuals of the Acadian lineage, constituting the first record of this lineage in Labrador. Atlantic and Acadian lineage individuals were found in both eastern and western Newfoundland. Multiple sampling locations in Labrador and Newfoundland contained fish of two or more different glacial lineages, implying their introgression. Glacial lineage did not appear to dictate contemporary genetic divergence between the pale and dark morph of char present in Gander Lake, Newfoundland. Both were predominately of the Atlantic lineage, suggesting the potential for their divergence in sympatry.

MAIN CONCLUSIONS

Our study reveals Labrador and Newfoundland to be a unique junction of three glacial lineages which have likely hybridized extensively in this region.

Phylogeography of lethal male fighting in a social spider mite

When males fight for access to females, such conflict rarely escalates into lethal fight because the risks and costs involved, that is, severe injury or death, are too high. The social spider mite, Stigmaeopsis miscanthi, does exhibit lethal male fights, and this male-male aggressiveness varies among populations. To understand the evolution of lethal fighting, we investigated aggressiveness in 42 populations and phylogenetic relationships in 47 populations along the Japanese archipelago. By analysis of the male weapon morph, a proxy for aggressiveness, we confirmed the existence of a mildly aggressive (ML) form, besides the low aggression (LW) and high aggression (HG) forms reported earlier. To evaluate demographic history of these three forms, we employed the approximate Bayesian computation approach using mtCOI sequences and taking into consideration the postlast glacial expansion history of the host plant, Miscanthus sinensis. As results, hierarchical split models are more likely to explain the observed genetic pattern than admixture models, and the ML form in the subtropical region was considered the ancestral group. The inferred demographic history was consistent with the one reconstructed for the host plant in a previous study. The LW form was split from the ML form during the last glacial period (20,000-40,000 years BP), and subsequently, the HG form was split from the ML form at the end of or after the last glacial period (5,494-10,988 years BP). The results also suggest that the mite invaded Japan more than once, resulting in the present parapatric distribution of LW and HG forms in eastern Japan.

Hitchhiking the high seas: Global genomics of rafting crabs

Population differentiation and diversification depend in large part on the ability and propensity of organisms to successfully disperse. However, our understanding of these processes in organisms with high dispersal ability is biased by the limited genetic resolution offered by traditional genotypic markers. Many neustonic animals disperse not only as pelagic larvae, but also as juveniles and adults while drifting or rafting at the surface of the open ocean. In theory, the heightened dispersal ability of these animals should limit opportunities for species diversification and population differentiation. To test these predictions, we used next-generation sequencing of genomewide restriction-site-associated DNA tags (RADseq) and traditional mitochondrial DNA sequencing, to investigate the species-level relationships and global population structure of Planes crabs collected from oceanic flotsam and sea turtles. Our results indicate that species diversity in this clade is low-likely three closely related species-with no evidence of cryptic or undescribed species. Moreover, our results indicate weak population differentiation among widely separated aggregations with genetic indices showing only subtle genetic discontinuities across all oceans of the world (RADseq F ST = 0.08-0.16). The results of this study provide unprecedented resolution of the systematics and global biogeography of this group and contribute valuable information to our understanding of how theoretical dispersal potential relates to actual population differentiation and diversification among marine organisms. Moreover, these results demonstrate the limitations of single gene analyses and the value of genomic-level resolution for estimating contemporary population structure in organisms with large, highly connected populations.

Genetic differentiation and asymmetric gene flow among Carpathian brown bear (Ursus arctos) populations-Implications for conservation of transboundary populations

The abundance and distribution of large carnivores in Europe have been historically reduced. Their recovery requires multilevel coordination, especially regarding transboundary populations. Here, we apply nuclear and mitochondrial genetic markers to test for admixture level and its impact on population genetic structure of contemporary brown bears (Ursus arctos) from the Eastern, Southern, and Western Carpathians. Carpathian Mountains (Europe). Nearly 400 noninvasive brown bear DNA samples from the Western (Poland) and Eastern Carpathians (Bieszczady Mountains in Poland, Slovakia, Ukraine) were collected. Together with DNA isolates from Slovakia and Romania, they were analyzed using the set of eight microsatellite loci and two mtDNA regions (control region and cytochrome b). A set of 113 individuals with complete genotypes was used to investigate genetic differentiation across national boundaries, genetic structuring within and between populations, and movement between populations. Transboundary brown bear subpopulations (Slovakia and Poland) did not show significant internal genetic structure, and thus were treated as cohesive units. All brown bears from the Western Carpathians carried mitochondrial haplotypes from the Eastern lineage, while the Western lineage prevailed in the brown bears from the Bieszczady Mountains. Despite similar levels of microsatellite variability, we documented significant differentiation among the studied populations for nuclear markers and mtDNA. We also detected male-biased and asymmetrical movement into the Bieszczady Mountains population from the Western Carpathians. Our findings suggest initial colonization of the Western Carpathians by brown bears possessing mtDNA from the Eastern lineage. Genetic structuring among populations at microsatellite loci could be a result of human-mediated alterations. Detected asymmetric gene flow suggests ongoing expansion from more abundant populations into the Bieszczady Mountains and thus supports a metapopulation model. The knowledge concerning this complex pattern can be implemented in a joint Carpathian brown bear management plan that should allow population mixing by dispersing males.

Ribosomal DNA polymorphisms reveal genetic structure and a phylogeographic pattern in the Burgundy truffle Tuber aestivum Vittad

Ectomycorrhizal ascomycetes belonging to the genus Tuber produce edible fruiting bodies known as truffles. Tuber aestivum, in particular, is a fungus appreciated worldwide and has a natural distribution throughout Europe. Most of the molecular studies conducted on this species have been focused on the question as to whether or not T. aestivum and the morphologically similar T. uncinatum are conspecific. Conversely, only a handful of studies have assessed the level and distribution of genetic diversity and occurrence of phylogeographic patterns in this species. Here, we analyzed the genetic diversity of T. aestivum over a wide geographic range, performing an extensive sampling of specimens from Turkey, which is novel, to the best of our knowledge. We compared the internal transcribed spacer (ITS) profiles of 45 samples collected in different Turkish areas with those of 144 samples from all over Europe. We identified 63 haplotypes, 32 of which were exclusively present in Turkey. The majority of these haplotyes were also population specific. Haplotype network analysis and statistical tests highlighted the presence of a genetic structure and phylogeographic pattern, with three spatially distinct genetic clusters (northeastern Europe, southern Europe, and Turkey), with Turkey representing a diversity hotspot. Based on these results, we hypothesize the presence of glacial refugia for T. aestivum in Turkey, whereas European populations likely experienced a population bottleneck. The possible occurrence of cryptic species among Turkish T. aestivum samples also emerged. Our results are of practical relevance for the marketing of T. aestivum truffles and mycorrhizal seedlings and the preservation of the biodiversity of this species.

Invasion Dynamics of A Termite, Reticulitermes flavipes, at Different Spatial Scales in France

Termites are social insects that can also be major pests. A well-known problem species is the subterranean termite, Reticulitermes flavipes. It is invasive in France and is thought to have arrived from Louisiana during the 18th century. While the putative source of French populations has been identified, little is known about how the termite spread following its establishment. Here, we examined expansion patterns at different spatial scales in urban areas to clarify how R. flavipes spread in France. Based on our analyses of phylogeography and population genetics, results suggest a scenario of successive introductions into the Charente-Maritime region, on the Atlantic Coast. Two major expansion fronts formed: one that spread toward the northeast and the other toward the southeast. At the regional scale, different spatial and genetic distribution patterns were observed: there was heterogeneity in Île-de-France and aggregation in Centre-Val de Loire. At the local scale, we found that our three focal urban sites each formed a single large colony that contained several secondary reproductives. Our findings represent a second step in efforts to reconstruct termite’s invasion dynamics. They also highlight the role that may have been played by the French railway network in transporting termites over long distances.