Restriction Site-associated DNA Sequencing Research Articles

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease.

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Implications of headwater contact zones for the riverine barrier hypothesis: a case study of the Blue-capped Manakin (Lepidothrix coronata).

Rivers frequently delimit the geographic ranges of species in the Amazon Basin. These rivers also define the boundaries between genetic clusters within many species, yet river boundaries have been documented to break down in headwater regions where rivers are narrower. To explore the evolutionary implications of headwater contact zones in Amazonia, we examined genetic variation in the Blue-capped Manakin (Lepidothrix coronata), a species previously shown to contain several genetically and phenotypically distinct populations across the western Amazon Basin. We collected restriction site-associated DNA sequence data (RADcap) for 706 individuals and found that spatial patterns of genetic structure indicate several rivers, particularly the Amazon and Ucayali, are dispersal barriers for L. coronata. We also found evidence that genetic connectivity is elevated across several headwater regions, highlighting the importance of headwater gene flow for models of Amazonian diversification. The headwater region of the Ucayali River provided a notable exception to findings of headwater gene flow by harboring non-admixed populations of L. coronata on opposite sides of a < 1-km-wide river channel with a known dynamic history, suggesting that additional prezygotic barriers may be limiting gene flow in this region.

Phylogenomic insights into the polyphyletic nature of Altai falcons within eastern sakers (Falco cherrug) and the origins of gyrfalcons (Falco rusticolus)

The Altai falcon from Central Asia always attracted the attention of humans. Long considered a totemic bird in its native area, modern falconers still much appreciated this large-bodied and mighty bird of prey due to its rarity and unique look. The peculiar body characteristics halfway between the saker falcon (Falco cherrug) and the gyrfalcon (F. rusticolus) triggered debates about its contentious taxonomy. The weak phylogenetic signal associated with traditional genetic methods could not resolve this uncertainty. Here, we address the controversial evolutionary origin of Altai falcons by means of a genome-wide approach, Restriction-site Associated DNA sequencing, using sympatric eastern sakers falcons, allopatric western saker falcons and gyrfalcons as outgroup. This approach provided an unprecedented insight into the phylogenetic relationships of the studied populations by delivering 17,095 unlinked SNPs shedding light on the polyphyletic nature of Altai falcons within eastern sakers. Thus we concluded that the former must correspond to a low taxonomic rank, probably an ecotype or form of the latter. Also, we found that eastern sakers are paraphyletic without gyrfalcons, thus, these latter birds are best regarded as the direct sister lineage of the eastern sakers. This evolutionary relationship, corroborated also by re-analyzing the dataset with the inclusion of outgroup samples (F. biarmicus and F. peregrinus), put eastern sakers into a new light as the potential ancestral genetic source of high latitude and altitude adaptation in descendent populations. Finally, conservation genomic values hint at the stable genetic background of the studied saker populations.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene.

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Genomic insights into isolation of the threatened Florida crested caracara (Caracara plancus).

We conducted a population genomic study of the crested caracara (Caracara plancus) using samples (n = 290) collected from individuals in Florida, Texas, and Arizona, United States. Crested caracaras are non-migratory raptors ranging from the southern tip of South America to the southern United States, including a federally protected relict population in Florida long thought to have been isolated since the last ice age. Our objectives were to evaluate genetic diversity and population structure of Florida's apparently isolated population and to evaluate taxonomic relationships of crested caracaras at the northern edge of their range. Using DNA purified from blood samples, we conducted double-digest restriction site associated DNA sequencing and sequenced the mitochondrial ND2 gene. Analyses of population structure using over 9,000 SNPs suggest that two major clusters are best supported, one cluster including only Florida individuals and the other cluster including Arizona and Texas individuals. Both SNPs and mitochondrial haplotypes reveal the Florida population to be highly differentiated genetically from Arizona and Texas populations, whereas, Arizona and Texas populations are moderately differentiated from each other. The Florida population's mitochondrial haplotypes form a separate monophyletic group, while Arizona and Texas populations share mitochondrial haplotypes. Results of this study provide substantial genetic evidence that Florida's crested caracaras have experienced long-term isolation from caracaras in Arizona and Texas and thus, represent a distinct evolutionary lineage possibly warranting distinction as an Evolutionarily Significant Unit (ESU) or subspecies. This study will inform conservation strategies focused on long-term survival of Florida's distinct, panmictic population.

Resource Status of Wild Hawthorn in Xinjiang and Research Progress of Genetic Diversity in Hawthorn

Wild hawthorn is a rare wild plant resource endemic to the western Tianshan Mountain in Xinjiang, containing three species: Crataegus songorica, C. alorocarpa and C. sanguirnta. However, the area of wild hawthorn forest has shrunken under the disturbance of human activities, genetic diversity has also suffered serious damage. Therefore, it is urgent to carry out the present status investigation of genetic diversity and assess the genetic differentiation of geographic populations of wild hawthorns. This paper summarizes the resource status of wild hawthorn in Xinjiang and research progress of genetic diversity in hawthorn, and developed an approach to screen molecular markers of single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) using Restriction-site Associated DNA sequencing (RAD-seq). SNP and SSR markers will be used to assess genetic diversity, genetic differentiations and genetic relationships of different populations of wild hawthorns.

Optimizing ddRAD sequencing for population genomic studies with ddgRADer.

Double-digest Restriction-site Associated DNA sequencing (ddRADseq) is widely used to generate genomic data for non-model organisms in evolutionary and ecological studies. Along with affordable paired-end sequencing, this method makes population genomic analyses more accessible. However, multiple factors should be considered when designing a ddRADseq experiment, which can be challenging for new users. The generated data often suffer from substantial read overlaps and adaptor contamination, severely reducing sequencing efficiency and affecting data quality. Here, we analyse diverse datasets from the literature and carry out controlled experiments to understand the effects of enzyme choice and size selection on sequencing efficiency. The empirical data reveal that size selection is imprecise and has limited efficacy. In certain scenarios, a substantial proportion of short fragments pass below the lower size-selection cut-off resulting in low sequencing efficiency. However, enzyme choice can considerably mitigate inadvertent inclusion of these shorter fragments. A simple model based on these experiments is implemented to predict the number of genomic fragments generated after digestion and size selection, number of SNPs genotyped, number of samples that can be multiplexed and the expected sequencing efficiency. We developed ddgRADer - http://ddgrader.haifa.ac.il/ - a user-friendly webtool and incorporated these calculations to aid in ddRADseq experimental design while optimizing sequencing efficiency. This tool can also be used for single enzyme protocols such as Genotyping-by-Sequencing. Given user-defined study goals, ddgRADer recommends enzyme pairs and allows users to compare and choose enzymes and size-selection criteria. ddgRADer improves the accessibility and ease of designing ddRADseq experiments and increases the probability of success of the first population genomic study conducted in labs with no prior experience in genomics.

Isolation and characterization of 32 SNP markers in Spinibarbus hollandi based on restriction site-associated DNA sequencing (RADseq)

Isolation and characterization of 32 SNP markers in Spinibarbus hollandi based on restriction site-associated DNA sequencing (RADseq)

RAD-seq data reveals robust phylogeny and morphological evolutionary history of Rhododendron

Rhododendron is famous for its high ornamental value. However, the genus is taxonomically difficult and the relationships within Rhododendron remain unresolved. In addition, the origin of key morphological characters with high horticulture value need to be explored. Both problems largely hinder utilization of germplasm resources. Most studies attempted to disentangle the phylogeny of Rhododendron, but only used a few genomic markers and lacked large-scale sampling, resulting in low clade support and contradictory phylogenetic signals. Here, we used restriction-site associated DNA sequencing (RAD-seq) data and morphological traits for 144 species of Rhododendron, representing all subgenera and most sections and subsections of this species-rich genus, to decipher its intricate evolutionary history and reconstruct ancestral state. Our results revealed high resolutions at subgenera and section levels of Rhododendron based on RAD-seq data. Both optimal phylogenetic tree and split tree recovered five lineages among Rhododendron. Subg. Therorhodion (clade I) formed the basal lineage. Subg. Tsutsusi and Azaleastrum formed clade II and had sister relationships. Clade III included all scaly rhododendron species. Subg. Pentanthera (clade IV) formed a sister group to Subg. Hymenanthes (clade V). The results of ancestral state reconstruction showed that Rhododendron ancestor was a deciduous woody plant with terminal inflorescence, ten stamens, leaf blade without scales and broadly funnelform corolla with pink or purple color. This study shows significant distinguishability to resolve the evolutionary history of Rhododendron based on high clade support of phylogenetic tree constructed by RAD-seq data. It also provides an example to resolve discordant signals in phylogenetic trees and demonstrates the application feasibility of RAD-seq with large amounts of missing data in deciphering intricate evolutionary relationships. Additionally, the reconstructed ancestral state of six important characters provides insights into the innovation of key characters in Rhododendron.

Phylogenomic data reveal cryptic diversity and deep phylogeographical structure within the common chuckwalla, Sauromalus ater (Squamata: Iguanidae)

Abstract Understanding how historical geological processes drive diversification and shape the contemporary distribution of species is fundamental to phylogeography. We take a genomic approach in order to elucidate the deep phylogeographical history and species limits of chuckwallas (Sauromalus), a conspicuous group of lizards of the arid lands of south-western North America. Phylogenetic and population genomic analyses of double digest restriction site-associated DNA sequencing data confirm the presence of at least two major lineages, peninsular and continental groups, within the widespread and morphologically variable common chuckwalla (Sauromalus ater). These lineages diversified in the vicinity of the head of the Gulf of California in north-eastern Baja California in the early Pliocene to late Miocene, during the formation of the northern gulf. The peninsular lineage of S. ater subsequently gave rise to the four insular endemic species of Sauromalus associated with the Baja California peninsula. Genomic analyses strongly support the continued recognition of the insular gigantics Sauromalus varius and Sauromalus hispidus as distinct species, although their relationship as sister species remains unresolved. Weaker phylogenetic signal for the insular species Sauromalus slevini and Sauromalus klauberi is provided by the genomic data; thus, it is advocated to continue recognizing these species until additional data can be analysed to evaluate their distinctiveness.

Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree.

In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado. The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated. High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction. The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Limited population genetic structure in the littoral gastropod Nerita tessellata (Neritimorpha, Neritidae) suggests high dispersal across the Caribbean Sea

The Caribbean Sea is characterized by an incredible biodiversity, including several endemic mollusk species, most of which have larval dispersal as their main mechanism of gene flow. It is known that the present-day population structure of a species reflects the combination of oceanographic currents, life-history traits, and historical events. Nerita tessellata is a common gastropod species in coastal, rocky, and exposed shores across the Caribbean Sea; however, little is known about the mechanisms that have shaped its distribution. Here, we tested potential barriers to gene flow and the associated genetic structure of N. tessellata in 12 populations across the Caribbean using digest restriction-site associated DNA sequencing (ddRAD-seq). Likewise, we evaluated demographic history and migration patterns through two coalescence approaches using Approximate Bayesian Computation (ABC) combined with supervised machine learning. The global analysis, based on 26,421 single nucleotide polymorphism (SNP) markers derived from ddRAD-seq, revealed low nucleotide diversity with high levels of gene flow and a lack of significant population differentiation within the survey area, indicating that N. tessellata is panmictic across the Caribbean Sea. Our results reveal that a pattern of north to south movement along Caribbean currents could potentially explain the limited genetic structure and high gene flow among N. tessellata populations. The scenarios simulated with ABC showed a constant effective population size for the ancestral population. Additionally, the four populations tested experienced demographic expansions that date to the Pleistocene. An extended pelagic larval stage could explain the large-scale genetic uniformity of this species across its distribution. Our findings contribute to understanding genetic diversity and population structure in a common gastropod within the Caribbean Sea.

Using genomics and morphometrics to monitor data-poor and commercially exploited octopod populations

Over 150 species of benthic octopods have been described within the ‘catch-all’ Octopus genus (Family: Octopodidae) and yet, many Octopus species harvested by fisheries remain unidentified to species-level due to a lack of distinguishing traits. Within species, there is also limited information on how populations differ genetically and the level of connectivity between populations. Therefore, we sampled octopods from commercial fisheries in southeast Australia, in order to identify the species, examine the phylogeographic relationships among species and the level of population genetic structuring within species, as well as to look for any adaptive genetic variation. The mitochondrial gene, cytochrome oxidase subunit III (COIII), was sequenced in 346 octopods along with single nucleotide polymorphisms using double digest restriction site-associated DNA sequencing (ddRADseq). Morphometric traits were also measured in mature specimens. The southern keeled octopus (‘Octopus’berrima) and pale octopus (‘Octopus’pallidus) were identified using COIII data. For ‘Octopus’berrima, we found that some populations whilst being morphologically similar were genetically distinct. In contrast, ‘Octopus’pallidus populations were both morphologically and genetically distinct across the studied regions. Our results provide key information to better inform conservation and management decisions for developing octopod fisheries in southeast Australia and highlight the importance of genomics tools in the conservation management of commercially and recreationally important species.

Genetic analysis of heat tolerance in hot pepper: insights from comprehensive phenotyping and QTL mapping.

High temperatures present a formidable challenge to the cultivation of hot pepper, profoundly impacting not only vegetative growth but also leading to flower and fruit abscission, thereby causing a significant reduction in yield. To unravel the intricate genetic mechanisms governing heat tolerance in hot pepper, an F2 population was developed through the crossing of two distinct genotypes exhibiting contrasting heat tolerance characteristics: DLS-161-1 (heat tolerant) and DChBL-240 (heat susceptible). The F2 population, along with the parental lines, was subjected to comprehensive phenotyping encompassing diverse morphological, physiological, and biochemical heat-related traits under high temperature conditions (with maximum temperature ranging from 31 to 46.5°C and minimum temperature from 15.4 to 30.5°C). Leveraging the Illumina Nova Seq-6000 platform, Double digest restriction-site associated DNA sequencing (ddRAD-seq) was employed to generate 67.215 Gb data, with subsequent alignment of 218.93 million processed reads against the reference genome of Capsicum annuum. Subsequent variant calling and ordering resulted in 5806 polymorphic SNP markers grouped into 12 LGs. Further QTL analysis identified 64 QTLs with LOD values ranging from 2.517 to 11.170 and explained phenotypic variance ranging from 4.05 to 19.39%. Among them, 21 QTLs, explaining more than 10% phenotypic variance, were identified as major QTLs controlling 9 morphological, 3 physiological, and 2 biochemical traits. Interestingly, several QTLs governing distinct parameters were found to be colocalized, suggesting either a profound correlation between the QTLs regulating these traits or their significant genomic proximity. In addition to the QTLs, we also identified 368380 SSR loci within the identified QTL regions, dinucleotides being the most abundant type (211,381). These findings provide valuable insights into the genetics of heat tolerance in hot peppers. The identified QTLs and SSR markers offer opportunities to develop heat-tolerant varieties, ensuring better crop performance under high-temperature conditions.

A new genomic resource to enable standardized surveys of SNPs across the native range of brook trout (Salvelinus fontinalis).

Understanding how genetic diversity is distributed across spatiotemporal scales in species of conservation or management concern is critical for identifying large-scale mechanisms affecting local conservation status and implementing large-scale biodiversity monitoring programmes. However, cross-scale surveys of genetic diversity are often impractical within single studies, and combining datasets to increase spatiotemporal coverage is frequently impeded by using different sets of molecular markers. Recently developed molecular tools make surveys based on standardized single-nucleotide polymorphism (SNP) panels more feasible than ever, but require existing genomic information. Here, we conduct the first survey of genome-wide SNPs across the native range of brook trout (Salvelinus fontinalis), a cold-adapted species that has been the focus of considerable conservation and management effort across eastern North America. Our dataset can be leveraged to easily design SNP panels that allow datasets to be combined for large-scale analyses. We performed restriction site-associated DNA sequencing for wild brook trout from 82 locations spanning much of the native range and domestic brook trout from 24 hatchery strains used in stocking efforts. We identified over 24,000 SNPs distributed throughout the brook trout genome. We explored the ability of these SNPs to resolve relationships across spatial scales, including population structure and hatchery admixture. Our dataset captures a wide spectrum of genetic diversity in native brook trout, offering a valuable resource for developing SNP panels. We highlight potential applications of this resource with the goal of increasing the integration of genomic information into decision-making for brook trout and other species of conservation or management concern.

A genome-wide association study of coat color in Chinese Rex rabbits.

Coat color is an important phenotypic characteristic of the domestic rabbit (Oryctolagus cuniculus) and has specific economic importance in the Rex rabbit industry. Coat color varies considerably among different populations of rabbits, and several causal genes for this variation have been thoroughly studied. Nevertheless, the candidate genes affecting coat color variation in Chinese Rex rabbits remained to be investigated. In this study, we collected blood samples from 250 Chinese Rex rabbits with six different coat colors. We performed genome sequencing using a restriction site-associated DNA sequencing approach. A total of 91,546 single nucleotide polymorphisms (SNPs), evenly distributed among 21 autosomes, were identified. Genome-wide association studies (GWAS) were performed using a mixed linear model, in which the individual polygenic effect was fitted as a random effect. We detected a total of 24 significant SNPs that were located within a genomic region on chromosome 4 (OCU4). After re-fitting the most significant SNP (OCU4:13,434,448, p = 1.31e-12) as a covariate, another near-significant SNP (OCU4:11,344,946, p = 7.03e-07) was still present. Hence, we conclude that the 2.1-Mb genomic region located between these two significant SNPs is significantly associated with coat color in Chinese Rex rabbits. The well-studied coat-color-associated agouti signaling protein (ASIP) gene is located within this region. Furthermore, low genetic differentiation was also observed among the six coat color varieties. In conclusion, our results confirmed that ASIP is a putative causal gene affecting coat color variation in Chinese Rex rabbits.

Ice holes microrefugia harbor genetically and functionally distinct populations of Vaccinium vitis-idaea (Ericaceae)

In the mountain terrain, ice holes are little depressions between rock boulders that are characterized by the exit of cold air able to cool down the rock surface even in summer. This cold air creates cold microrefugia in warmer surroundings that preserve plant species probably over thousands of years under extra-zonal climatic conditions. We hypothesized that ice hole populations of the model species Vaccinium vitis-idaea (Ericaceae) show genetic differentiation from nearby zonal subalpine populations, and high functional trait distinctiveness, in agreement with genetic patterns. We genotyped almost 30,000 single nucleotide polymorphisms using restriction site-associated DNA sequencing and measured eight functional traits indicative of individual performance and ecological strategies. Genetic results showed high differentiation among the six populations suggesting isolation. On siliceous bedrock, ice hole individuals exhibited higher levels of admixture than those from subalpine populations which could have experienced more bottlenecks during demographic fluctuations related to glacial cycles. Ice hole and subalpine calcareous populations clearly separated from siliceous populations, indicating a possible effect of bedrock in shaping genetic patterns. Trait analysis reflected the bedrock effect on populations’ differentiation. The significant correlation between trait and genetic distances suggests the genetic contribution in shaping intraspecific functional differentiation. In conclusion, extra-zonal populations reveal a prominent genetic and phenotypic differentiation determined by history and ecological contingency. Therefore, microrefugia populations can contribute to the overall variability of the species and lead to intraspecific-driven responses to upcoming environmental changes.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios.

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study.

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Source tracing of Anthonomus grandis captured in areas of the USA where the species had previously been eradicated

Anthonomus grandis Boheman (Coleoptera: Curculionidae) is comprised of two subspecies: the boll weevil (A. g. grandis) and the thurberia weevil (A. g. thurberiae). The thurberia weevil is generally regarded as a non-pest variant of the species, whereas the boll weevil is an important economic pest of commercial cotton. In the USA, a nationwide eradication program has been largely successful against the boll weevil, but reintroductions of eradicated zones remain a threat. In 2015, A. grandis was trapped in the functionally eradicated Winter Garden growing region in Texas, where it had not been found since 2012. Weevils established in the area and triggered a multi-year, multi-million-dollar mitigation effort. In 2017, the species was trapped in Hidalgo County, New Mexico, another area where it had been previously eradicated. In this study, we obtained specimens from both trapping events and used double digest restriction site-associated DNA sequencing to obtain genetic data in homology with previously published A. grandis population genetic data. We conducted source tracing to match the captured weevils to previously described populations to determine a likely source population. We also compared the accuracy of the approach to other diagnostic methods currently used by the United States Department of Agriculture. The results suggested that a small founder population of weevils dispersed from the Lower Rio Grande Valley population in southern Texas and northern Tamaulipas and established in the Winter Garden. In Hidalgo Co., NM, the weevils were found to be genetically indistinct from non-pest populations of thurberia weevil associated with Arizona wild cotton (Gossypium thurberi). Thurberia weevils in Arizona remain unlikely to threaten commercial growers. Improved rapid diagnostic tools may be desirable to enable better assessment of trap captures and allow for a reduction in the application of mitigation efforts when the non-pest variant is identified.