Assessing pyrethroid resistance in Aedes aegypti from Cordoba Colombia: Implications of kdr mutations.

Resistance to insecticides is one of the great challenges that vector control programs must face. The constant use of pyrethroid-type insecticides worldwide has caused selection pressure in populations of the Aedes aegypti vector, which has promoted the emergence of resistant populations. The resistance mechanism to pyrethroid insecticides most studied to date is target-site mutations that desensitize the voltage-gated sodium channel (VGSC) of the insect to the action of pyrethroids. In the present study, susceptibility to the pyrethroid insecticides permethrin, lambda-cyhalothrin, and deltamethrin was evaluated in fourteen populations from the department of Córdoba, Colombia. The CDC bottle bioassay and WHO tube methods were used. Additionally, the frequencies of the F1534C, V1016I, and V410L mutations were determined, and the association of resistance with the tri-locus haplotypes was examined. The results varied between the two techniques used, with resistance to permethrin observed in thirteen of the fourteen populations, resistance to lambda-cyhalothrin in two populations, and susceptibility to deltamethrin in all the populations under study with the CDC method. In contrast, the WHO method showed resistance to the three insecticides evaluated in all populations. The frequencies of the mutated alleles ranged from 0.05-0.43 for 1016I, 0.94-1.0 for 1534C, and 0.01-0.59 for 410L. The triple homozygous mutant CIL haplotype was associated with resistance to all three pyrethroids evaluated with the WHO bioassay, while with the CDC bioassay, it was only associated with resistance to permethrin. This study highlights the importance of implementing systematic monitoring of kdr mutations, allowing resistance management strategies to be dynamically adjusted to achieve effective control of Aedes aegypti.

VolcanoSV enables accurate and robust structural variant calling in diploid genomes from single-molecule long read sequencing.

Structural variants (SVs) significantly contribute to human genome diversity and play a crucial role in precision medicine. Although advancements in single-molecule long-read sequencing offer a groundbreaking resource for SV detection, identifying SV breakpoints and sequences accurately and robustly remains challenging. We introduce VolcanoSV, an innovative hybrid SV detection pipeline that utilizes both a reference genome and local de novo assembly to generate a phased diploid assembly. VolcanoSV uses phased SNPs and unique k-mer similarity analysis, enabling precise haplotype-resolved SV discovery. VolcanoSV is adept at constructing comprehensive genetic maps encompassing SNPs, small indels, and all types of SVs, making it well-suited for human genomics studies. Our extensive experiments demonstrate that VolcanoSV surpasses state-of-the-art assembly-based tools in the detection of insertion and deletion SVs, exhibiting superior recall, precision, F1 scores, and genotype accuracy across a diverse range of datasets, including low-coverage (10x) datasets. VolcanoSV outperforms assembly-based tools in the identification of complex SVs, including translocations, duplications, and inversions, in both simulated and real cancer data. Moreover, VolcanoSV is robust to various evaluation parameters and accurately identifies breakpoints and SV sequences.

Chromosome-level genome assembly and annotation of Flueggea virosa (Phyllanthaceae).

Flueggea virosa (Roxb. ex Willd.) Royle, an evergreen shrub and small tree in the Phyllanthaceae family, holds significant potential in garden landscaping and pharmacological applications. However, the lack of genomic data has hindered further scientific understanding of its horticultural and medicinal values. In this study, we have assembled a haplotype-resolved genome of F. virosa for the first time. The two haploid genomes, named haplotype A genome and haplotype B genome, are 487.33 Mb and 477.53 Mb in size, respectively, with contig N50 lengths of 31.45 Mb and 32.81 Mb. More than 99% of the assembled sequences were anchored to 13 pairs of pseudo-chromosomes. Furthermore, 21,587 and 21,533 protein-coding genes were predicted in haplotype A and haplotype B genomes, respectively. The availability of this chromosome-level genome fills the gap in genomic data for F. virosa and provides valuable resources for molecular studies of this species, supporting future research on speciation, functional genomics, and comparative genomics within the Phyllanthaceae family.

Chromosome-level and haplotype-resolved genome assembly of Dracaena cambodiana (Asparagaceae).

Dracaena cambodiana Pierre ex Gagn. (Asparagaceae) is the source plant of Dragon's blood and has high ornamental values in gardening. Currently, this species is classified as the second-class state-protected species in the National Key Protected Wild Plants (NKPWP) of China. However, limited genomic data has hindered a more comprehensive scientific understanding of the processes involved in the production of Dragon's blood and the related conservation genomics research. In this study, we assembled a haplotype-resolved genome of D. cambodiana. The haploid genomes, haplotype A and haplotype B, are 1,015.22 Mb and 1,003.13 Mb in size, respectively. The completeness of haplotype A and haplotype B genomes was 98.60% and 98.20%, respectively, using the "embryophyta_10" dataset. Haplotype A and haplotype B genomes contained 27,361 and 27,066 protein-coding genes, respectively, with nearly all being functionally annotated. These findings provide new insights into the genomic characteristics of D. cambodiana and will offer additional genomic resources for studying the biosynthesis mechanism of Dragon's blood and the horticultural application of Dragon trees.

MAPT haplotype-associated transcriptomic changes in progressive supranuclear palsy.

Progressive supranuclear palsy (PSP) is a neurodegenerative movement and cognitive disorder characterized by abnormal accumulation of the microtubule-associated protein tau in the brain. Biochemically, inclusions in PSP are enriched for tau proteoforms with four microtubule-binding domain repeats (4R), an isoform that arises from alternative tau pre-mRNA splicing. While preferential aggregation and reduced degradation of 4R tau protein is thought to play a role in inclusion formation and toxicity, an alternative hypothesis is that altered expression of tau mRNA isoforms plays a causal role. This stems from the observation that PSP is associated with common variation in the tau gene (MAPT) at the 17q21.31 locus which contains low copy number repeats flanking a large recurrent genomic inversion. The complex genomic structural changes at the locus give rise to two dominant haplotypes, termed H1 and H2, that have the potential to markedly influence gene expression. Here, we explored haplotype-dependent differences in gene expression using a bulk RNA-seq dataset derived from human post-mortem brain tissue from PSP (n = 84) and controls (n = 77) using a rigorous computational pipeline, including alternative pre-mRNA splicing. We found 3579 differentially expressed genes in the temporal cortex and 10,011 in the cerebellum. We also found 7214 differential splicing events in the temporal cortex and 18,802 in the cerebellum. In the cerebellum, total tau mRNA levels and the proportion of transcripts encoding 4R tau were significantly increased in PSP compared to controls. In the temporal cortex, the proportion of reads that expressed 4R tau was increased in cases compared to controls. 4R tau mRNA levels were significantly associated with the H1 haplotype in the temporal cortex. Further, we observed a marked haplotype-dependent difference in KANSL1 expression that was strongly associated with H1 in both brain regions. These findings support the hypothesis that sporadic PSP is associated with haplotype-dependent increases in 4R tau mRNA that might play a causal role in this disorder.

Haplotype analysis identifies functional elements in monoclonal gammopathy of unknown significance.

Genome-wide association studies (GWASs) based on common single nucleotide polymorphisms (SNPs) have identified several loci associated with the risk of monoclonal gammopathy of unknown significance (MGUS), a precursor condition for multiple myeloma (MM). We hypothesized that analyzing haplotypes might be more useful than analyzing individual SNPs, as it could identify functional chromosomal units that collectively contribute to MGUS risk. To test this hypothesis, we used data from our previous GWAS on 992 MGUS cases and 2910 controls from three European populations. We identified 23 haplotypes that were associated with the risk of MGUS at the genome-wide significance level (p < 5 × 10-8) and showed consistent results among all three populations. In 10 genomic regions, strong promoter, enhancer and regulatory element-related histone marks and their connections to target genes as well as genome segmentation data supported the importance of these regions in MGUS susceptibility. Several associated haplotypes affected pathways important for MM cell survival such as ubiquitin-proteasome system (RNF186, OTUD3), PI3K/AKT/mTOR (HINT3), innate immunity (SEC14L1, ZBP1), cell death regulation (BID) and NOTCH signaling (RBPJ). These pathways are important current therapeutic targets for MM, which may highlight the advantage of the haplotype approach homing to functional units.

Indian ASD probands with 25(OH)D and vitamin D binding protein deficiency exhibited higher severity.

The severity of autism spectrum disorder (ASD) shows wide variations, though the reason remains unclear. Vitamin D (VitD) deficiency is considered a risk factor for ASD and its supplementation was reported to reduce symptom severity. Since VitD, either synthesized in the skin or absorbed from the food, is transported to the liver by the vitamin D binding protein (DBP), we have analyzed DBP genetic polymorphisms [rs7041 (A/C), rs4588 (G/T), and rs3755967 (C/T)] affecting DBP function [Case = 411; Control = 397], levels of plasma 25(OH)D and DBP [Case = 25; Control = 26], and DBP mRNA expression [Case = 74; Control = 44] in a group of Indo-Caucasoid ASD probands and neurotypical subjects. ASD probands with rs7041'CC', rs4588 'TT', and rs3755967 'TT' genotypes exhibited higher scores for a few traits. Scores for Imitation and Listening response were also higher in the presence of the "A-T" haplotype (rs7041-rs4588). Plasma 25(OH)D and DBP levels as well as DBP mRNA expressions were significantly lower in the ASD probands as compared to the neurotypical subjects. We infer that DBP deficiency, in the presence of risk genetic variants, could be one of the reasons for the reported 25(OH)D deficiency of the ASD probands.

mtDNA "nomenclutter" and its consequences on the interpretation of genetic data.

Population-based studies of human mitochondrial genetic diversity often require the classification of mitochondrial DNA (mtDNA) haplotypes into more than 5400 described haplogroups, and further grouping those into hierarchically higher haplogroups. Such secondary haplogroup groupings (e.g., "macro-haplogroups") vary across studies, as they depend on the sample quality, technical factors of haplogroup calling, the aims of the study, and the researchers' understanding of the mtDNA haplogroup nomenclature. Retention of historical nomenclature coupled with a growing number of newly described mtDNA lineages results in increasingly complex and inconsistent nomenclature that does not reflect phylogeny well. This "clutter" leaves room for grouping errors and inconsistencies across scientific publications, especially when the haplogroup names are used as a proxy for secondary groupings, and represents a source for scientific misinterpretation. Here we explore the effects of phylogenetically insensitive secondary mtDNA haplogroup groupings, and the lack of standardized secondary haplogroup groupings on downstream analyses and interpretation of genetic data. We demonstrate that frequency-based analyses produce inconsistent results when different secondary mtDNA groupings are applied, and thus allow for vastly different interpretations of the same genetic data. The lack of guidelines and recommendations on how to choose appropriate secondary haplogroup groupings presents an issue for the interpretation of results, as well as their comparison and reproducibility across studies. To reduce biases originating from arbitrarily defined secondary nomenclature-based groupings, we suggest that future updates of mtDNA phylogenies aimed for the use in mtDNA haplogroup nomenclature should also provide well-defined and standardized sets of phylogenetically meaningful algorithm-based secondary haplogroup groupings such as "macro-haplogroups", "meso-haplogroups", and "micro-haplogroups". Ideally, each of the secondary haplogroup grouping levels should be informative about different human population history events. Those phylogenetically informative levels of haplogroup groupings can be easily defined using TreeCluster, and then implemented into haplogroup callers such as HaploGrep3. This would foster reproducibility across studies, provide a grouping standard for population-based studies, and reduce errors associated with haplogroup nomenclatures in future studies.

Lineage-informative microhaplotypes for recurrence classification and spatio-temporal surveillance of Plasmodium vivax malaria parasites.

Challenges in classifying recurrent Plasmodium vivax infections constrain surveillance of antimalarial efficacy and transmission. Recurrent infections may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or reinfection. Molecular inference of familial relatedness (identity-by-descent or IBD) can help resolve the probable origin of recurrences. As whole genome sequencing of P. vivax remains challenging, targeted genotyping methods are needed for scalability. We describe a P. vivax marker discovery framework to identify and select panels of microhaplotypes (multi-allelic markers within small, amplifiable segments of the genome) that can accurately capture IBD. We evaluate panels of 50-250 microhaplotypes discovered in a global set of 615 P. vivax genomes. A candidate global 100-microhaplotype panel exhibits high marker diversity in the Asia-Pacific, Latin America and horn of Africa (median HE = 0.70-0.81) and identifies 89% of the polyclonal infections detected with genome-wide datasets. Data simulations reveal lower error in estimating pairwise IBD using microhaplotypes relative to traditional biallelic SNP barcodes. The candidate global panel also exhibits high accuracy in predicting geographic origin and captures local infection outbreak and bottlenecking events. Our framework is open-source enabling customised microhaplotype discovery and selection, with potential for porting to other species or data resources.

Genetic structure of black soldier flies in northern Iran.

BACKGROUND: The black soldier fly (BSF), Hermetia illucens, is known for nutrient-recycling through the bioconversion of organic waste into protein-rich insect larvae that can be processed into an animal feed ingredient. However, information on species distribution and its genetic structure in Iran is scarce. METHODS AND RESULTS: We directed a survey on the Caspian Sea coast, with a reconstructing demographic relationships study using two parts of mitochondrial cytochrome C oxidase 1 (COI) gene (barcode and 3' end regions) and nuclear internal transcribed spacer 2 (ITS2) to identify BSF' genetic diversity in retrospect to the global diversity and the potential origin of the Iranian BSF population. Larvae and adults were recovered from highly decomposed poultry manure, in May 2020. Sequence analysis of both regions of COI gene (about 1500 bp) revealed a single haplotype, identical to that of haplotype C, a worldwide commercial strain originated from Nearctic, Palearctic, or African biogeographic regions. However, the ITS2 locus was confirmed to be invariable across samples from diverse biogeographic regions. CONCLUSION: The results proved the presence of BSF in north of Iran. However, it is not possible to determine with certainty when and where this species first established in Iran, and they have likely been released to nature due to the existence of companies importing and breeding such flies. Due to heavy international trading, the introduction and settlement of this fly in the southern coasts of the country is highly suggested.

Interbreeding between farmers and hunter-gatherers along the inland and Mediterranean routes of Neolithic spread in Europe.

The Neolithic (i.e., farming and stockbreeding) spread from the Near East across Europe since about 9000 years before the common era (BCE) until about 4000 yr BCE. It followed two main routes, namely a sea route along the northern Mediterranean coast and an inland one across the Balkans and central Europe. It is known that the dispersive behavior of farmers depended on geography, with longer movements along the Mediterranean coast than along the inland route. In sharp contrast, here we show that for both routes the percentage of farmers who interbred with hunter-gatherers and/or acculturated one of them was strikingly the same (about 3.6%). Therefore, whereas the dispersive behavior depended on the proximity to the Mediterranean sea, the interaction behavior (incorporation of hunter-gatherers) did not depend on geographical constraints but only on the transition in the subsistence economy (from hunting and gathering to farming) and its associated way of life. These conclusions are reached by analyzing the clines of haplogroup K, which was virtually absent in hunter-gatherers and the most frequent mitochondrial haplogroup in early farmers. Similarly, the most frequent Y-chromosome Neolithic haplogroup (G2a) displays an inland cline that agrees with the percentage of interbreeding reported above.

Linking watershed formation with the phylogenetic distribution of a soil microscopic fungus in Yunnan Province, China.

BACKGROUND: Phylogeographic studies have gained prominence in linking past geological events to the distribution patterns of biodiversity, primarily in mountainous regions. However, such studies often focus on plant taxa, neglecting the intricate biogeographical patterns of microbes, particularly soil microbial communities. This article explores the spatial distribution of the nematode-trapping fungus Arthrobotrys oligospora, a widespread microorganism, in a tectonically active region at the southeastern edge of the Qinghai-Tibetan Plateau. By analysing the genetic variation of this fungus alongside the historical structure of major river watersheds, we sought to uncover potential connections between the two. Our study involved sampling 149 strains from 116 sites across six major watersheds in the region. RESULTS: The resulting haplotype network revealed five distinct clusters, each corresponding closely to a specific watershed. These clusters exhibited high haplotype diversity and low nucleotide diversity, supporting the notion of watershed-based segregation. Further analysis of haplotypes shared across watersheds provided evidence for three proposed past river connections. In particular, we found numerous shared haplotypes between the Yangtze and Mekong basins, as well as between the Yangtze and the Red basins. Evidence for a Irrawaddy-Salween-Red and a Yangtze-Pearl-Red river connections were also portrayed in our mapping exercise. CONCLUSIONS: These findings emphasize the crucial role of historical geomorphological events in shaping the biogeography of microbial biodiversity, alongside contemporary biotic and abiotic factors. Watershed perimeters emerged as effective predictors of such patterns, suggesting their suitability as analytical units for regional-scale studies. Our study also demonstrates the potential of microorganisms and phylogeographic approaches to complement traditional geological analyses, providing a more comprehensive understanding of past landscape structure and its evolution.

Association analysis of HSP90AA1 polymorphism with thermotolerance in tropically adapted Indian crossbred cattle.

Raising cattle is a lucrative business that operates globally but is confronted by many obstacles, such as thermal stress, which results in substantial monetary losses. A vital role of heat shock proteins (HSPs) is to protect cells from cellular damage. HSP90 is a highly prevalent, extremely adaptable gene linked to physiological resilience in thermal stress. This study aimed to find genetic polymorphisms of the HSP90AA1 gene in Karan Fries cattle and explore their relationship to thermal tolerance and production traits. One SNP (g.3292 A > C) was found in the Intron 8 and three SNPs loci (g.4776 A > G, g.5218T > C and g.5224 A > C) were found in the exon 11 of 100 multiparous Karan Fries cattle. The association study demonstrated that the SNP1-g.3292 A > C was significantly (P < 0.01) linked to the variables respiratory rate (RR), heat tolerance coefficient (HTC) and total milk yield (TMY (kg)) attributes. There was no significant correlation identified between any of the other SNP sites (SNP2-g.4776 A > G; SNP3-g.5218T > C; SNP4-g.5224 A > C) with the heat tolerance and production attributes in Karan Fries cattle. Haploview 4.2 and SHEsis software programs were used to analyse pair linkage disequilibrium and construct haplotypes for HSP90AA1. Association studies indicated that the Hap3 (CATA) was beneficial for heat tolerance breeding in Karan Fries cattle. In conclusion, genetic polymorphisms and haplotypes in the HSP90AA1 were associated with thermal endurance attributes. This relationship can be utilized as a beneficial SNP or Hap marker for genetic heat resistance selection in cow breeding platforms.

Chromosome-level scaffolding of haplotype-resolved assemblies using Hi-C data without reference genomes.

Scaffolding is crucial for constructing most chromosome-level genomes. The high-throughput chromatin conformation capture (Hi-C) technology has become the primary scaffolding strategy due to its convenience and cost-effectiveness. As sequencing technologies and assembly algorithms advance, constructing haplotype-resolved genomes is increasingly preferred because haplotypes can provide additional genetic information on allelic and non-allelic variations. ALLHiC is a widely used allele-aware scaffolding tool designed for this purpose. However, its dependence on chromosome-level reference genomes and a higher chromosome misassignment rate still impede the unravelling of haplotype-resolved genomes. Here we present HapHiC, a reference-independent allele-aware scaffolding tool with superior performance on chromosome assignment as well as contig ordering and orientation. In addition, we provide new insights into the challenges in allele-aware scaffolding by conducting comprehensive analyses on various adverse factors. Finally, with the help of HapHiC, we constructed the haplotype-resolved allotriploid genome for Miscanthus × giganteus, an important lignocellulosic bioenergy crop.

Haplotype-based Analysis of 6q24-25 for Association with Rheumatoid Arthritis.

BACKGROUND: Rheumatoid arthritis (RA) is a multisystem inflammatory disorder. Family history of RA is an important risk factor as it is strongly linked with the inherited HLA-DR4 (most specifically DR0401 and 0404). The aim of this study is to conduct the haplotype-based analysis of 6q24-25 and evaluate its association with RA. MATERIALS AND METHODS: Case-control study which included all patients attending outpatient department (OPD) at Sardar Patel Medical College, Bikaner, Rajasthan and volunteers (only for blood samples). Blood samples of patients were collected. As per inclusion and exclusion criteria, a total of 103 subjects lacking history of disease were included under control group, while 48 cases were recruited as study group. Any significant departure of genotype distribution is evaluated using Hardy-Weinberg equilibrium by Chi-squared test. RESULTS: Case-control association was done using data from 151 genomic deoxyribonucleic acid (DNA) samples, which were allele typed. RA is significantly associated with >305bp (the longer allele of D6S1053 corresponding to 11 tetramer (GATA) repeats. Differences in individual allele frequency within the control population and RA cases were observed which shows the bimodal distribution of the 10 alleles of D6S1053 short tandem repeat (STR) marker observed in the cohort tested by us. No significant association with the risk for RA is shown by the allele for D6S1053 and the mutant allele 118G. Similarly, OPRM1 gene's haplotype frequency for rs1799972 for D6S1053 allele has not shown any added risk with the wild allele 17C compared to controls. The polymorphism showed that 17T depicted a higher odds ratio of 1.3 with an associated risk of 1.15 in the presence of longer allele of DS61053. Significance observed between short allele and RA was lost when haplotype analysis for the two genes were taken together. There was no difference observed between the short or long allele and 17T/C while there was a significant difference observed when haplotype frequencies were compared with alleles of A118G and the long and short allele of DS61053. CONCLUSION: We concluded that the short allele of µ-opioid receptor (MOR) gene offered a clear protection from a risk of getting RA.

HLA alleles associated to susceptibility to gliptin-associated bullous pemphigoid in Italian patients.

Bullous pemphigoid (BP), although a rare disease, is the most frequent subepidermal autoimmune disorder. Treatment with gliptins, used for type 2 diabetes, was reported as associated with BP onset. To identify HLA alleles that may reflect a higher susceptibility to BP in the Italian population, we analysed 30 patients affected by idiopathic bullous pemphigoid (IBP) and 86 gliptin-associated BP (GABP) patients. A significant association between HLA-DQB1*03:01 allele and IBP and GABP patients was found. Of note, both IBP and GABP were significantly associated with one of the following haplotypes: DRB1*11:01, DRB3*02:02, DQA1*05:05, DQB1*03:01 or DRB1*11:04, DRB3*02:02, DQA1*05:05 and DQB1*03:01. These data identify, for the first time, potential markers of susceptibility to BP in the Italian population, especially when associated with gliptin intake.

VDR Gene Polymorphisms (BsmI, FokI, TaqI, ApaI) in Total Hip Arthroplasty Outcome Patients.

A total hip arthroplasty (THA) can improve quality of life, but loosening of the hip prosthesis is a complex problem in which vitamin D may also play a role. The Vitamin D Receptor (VDR) is involved in the response of cells to the action of vitamin D, and its genetic variability raises the question of whether individual differences could influence the risk of prosthesis loosening. The aim of this study was to investigate the relationship between VDR single nucleotide polymorphisms (SNPs) (ApaI, BsmI, FokI and TaqI) and the serum VDR and 25(OH)D levels in three groups of patients: (1) arthroscopy patients after THA without loosening of the prosthesis (CA-Control Arthroplasty), (2) patients after THA with loosened hip prostheses (L-Loosening) and (3) the control group (C-Control). Our results suggest that the genotypes tt of TaqI, BB of BsmI, and FF of FokI may influence the VDR effect in patients with loosened protheses. Our results showed that the ACAC haplotype (AtBF) was over two times more frequent in the L group than in CA + C: OR =2.35 [95% CI 1.44-3.83; p = 0.001]. There was no significant correlation between the VDR and serum 25(OH)D levels, but there were differences between studied groups.

Genetic diversity of a recovering European roller (Coracias garrulus) population from Serbia.

The European Roller (Coracias garrulus), a long-distance migratory bird, faced a considerable decline in breeding pairs throughout Europe at the end of the 20th century. Due to conservation efforts and the installation of nesting boxes, the population of the European Roller in Serbia has made a remarkable recovery. Here, we used the variability of nucleotide sequences of the mitochondrial DNA (mtDNA) control region and 10 microsatellite loci to assess the genetic diversity and structuring, phylogeographic patterns and demographic history of this species using 224 individuals from Serbia. Our results showed moderate level of genetic diversity (HO = 0.392) and a slightly elevated level of inbreeding and homozygosity (FIS = 0.393). Genetic structuring based on microsatellite data indicated three genetic clusters, but without a clear spatial pattern. High haplotype diversity (Hd = 0.987) of the mtDNA control region sequences was detected, and neutrality tests indicated a recent demographic expansion. The phylogeographic analysis, which also included previously published sequences of the mtDNA control region, supported the subdivision into two distinct European and Asian haplogroups (ΦST = 0.712). However, the results of our study showed that a larger number of haplotypes sampled in Serbia are clustered in the Asian haplogroup as compared to previous studies, indicating a historically continuous distribution of this species and possibly a wider distribution of the subspecies Coracias garrulus semenovwi. Our results suggest that the European Roller population in Serbia is genetically stable, with no evidence of recent bottlenecks, and emphasize the importance of artificial nest boxes for promoting and maintaining population dynamics of European Rollers.

Combining distribution modelling and phylogeography to understand present, past and future of an endangered spider.

BACKGROUND: Understanding how endangered species respond to climatic changes is fundamental for their conservation. Due to its restricted geographic range, its sensitivity to the ongoing global warming and its continuing decline, the Southwestern-Alpine endemic wolf spider Vesubia jugorum is currently classified as Endangered in the IUCN Red List. Here, we combined species distribution modelling (SDM) and phylogeographic inference to describe the present, the past and the future of this species in light of the mtDNA genetic structure of extant populations. RESULTS: Phylogenetic and network analyses show a high level of genetic differentiation and a strong genetic structure of the populations, likely explicable by a long history of isolation and survival in separate refugia. The SDM projection into past climatic conditions supports these results by showing a smaller distribution range compared to present, mostly restricted to the Maritime and Ligurian Alps, which possibly served as main refugium. Future forecast shows a significant shift in the bioclimatic range towards higher altitudes and latitudes, with a drastic decrease of habitat suitability in the central and south-eastern parts of the range, with consequent general loss of haplotype diversity. CONCLUSION: SDM and phylogeographic inference support the hypothesis that the current distribution and the genetic structure of the extant populations mirror the survival in situ of Vesubia jugorum across repeated glacial and interglacial phases, in line with the 'long-term stability hypothesis'. Future predictions show a significant shift in the bioclimatic range that V. jugorum will be likely unable to track, with profound impact on its long-term survival and its genetic diversity. Our considerations have implication for conservation genetics, highlighting the pivotal role of the transboundary protected areas of the SW-Alps in promoting conservation efforts for this species.