genodive version 3.0: Easy-to-use software for the analysis of genetic data of diploids and polyploids

Diversity assessment of a lesser known buffalo population from Central India and its comparative evaluation reveals presence of sufficient genetic variation and absence of selection

Planned breeding and conservation strategies for a lesser-known population require an assessment of complete genetic diversity and population structure analysis in addition to its morphometric characteristics. In the present study, a comparative analysis of the genetic structure of a rare buffalo population, namely Chhattisgarhi, was extensively studied using a panel of FAO-recommended microsatellite markers along with well-established breeds namely Murrah, Nili-Ravi, Gojri, Kalahandi, and Nagpuri. Mode shift analysis indicated the absence of genetic bottleneck in the recent past. Assessment of genetic diversity indices across all loci indicated the presence of sufficient genetic variation within and between populations. Analysis of molecular variance between the six different buffalo populations attributed 19.05% of the variations to between-population differentiation. Cluster analyses using DAPC and Bayesian approach along with the phylogenetic tree based on UPGMA grouped six populations into three groups. The Chhattisgarhi population was revealed to be genetically closer to Nagpuri and Kalahandi populations. The study reveals the presence of sufficient genetic diversity within the Chhattisgarhi population and indicates the absence of a systematic selection program. We suggest improvement and conservation programs should be planned for this breed in the near future through short-term selection.

Severe hurricanes increase recruitment and gene flow in the clonal sponge Aplysina cauliformis

Upright branching sponges, such as Aplysina cauliformis, provide critical three-dimensional habitat for other organisms and assist in stabilizing coral reef substrata, but are highly susceptible to breakage during storms. Breakage can increase sponge fragmentation, contributing to population clonality and inbreeding. Conversely, storms could provide opportunities for new genotypes to enter populations via larval recruitment, resulting in greater genetic diversity in locations with frequent storms. The unprecedented occurrence of two Category 5 hurricanes in close succession during 2017 in the U.S. Virgin Islands (USVI) provided a unique opportunity to evaluate whether recolonization of newly available substrata on coral reefs was due to local (e.g. re-growth of remnants, fragmentation, larval recruitment) or remote (e.g. larval transport and immigration) sponge genotypes. We sampled A. cauliformis adults and juveniles from four reefs around St. Thomas and two in St. Croix (USVI). Using a 2bRAD protocol, all samples were genotyped for single-nucleotide polymorphisms (SNPs). Results showed that these major storm events favoured sponge larval recruitment but did not increase the genetic diversity of A. cauliformis populations. Recolonization of substratum post-storms via clonality was lower (15%) than expected and instead was mainly due to sexual reproduction (85%) via local larval recruitment. Storms did enhance gene flow among and within reef sites located south of St. Thomas and north of St. Croix. Therefore, populations of clonal marine species with low pelagic dispersion, such as A. cauliformis, may benefit from increased frequency and magnitude of hurricanes for the maintenance of genetic diversity and to combat inbreeding, enhancing the resilience of Caribbean sponge communities to extreme storm events.

Structure of Aspergillus flavus populations associated with maize in Greece, Spain, and Serbia: Implications for aflatoxin biocontrol on a regional scale

Aspergillus flavus is the most frequently identified producer of aflatoxins. Non-aflatoxigenic members of the A. flavus L strains are used in various continents as active ingredients of bioprotectants directed at preventing aflatoxin contamination by competitive displacement of aflatoxin producers. The current research examined the genetic diversity of A. flavus L strain across southern Europe to gain insights into the population structure and evolution of this species and to evaluate the prevalence of genotypes closely related to MUCL54911, the active ingredient of AF-X1. A total of 2173L strain isolates recovered from maize collected across Greece, Spain, and Serbia in 2020 and 2021 were subjected to simple sequence repeat (SSR) genotyping. The analysis revealed high diversity within and among countries and dozens of haplotypes shared. Linkage disequilibrium analysis indicated asexual reproduction and clonal evolution of A. flavus L strain resident in Europe. Moreover, haplotypes closely related to MUCL54911 were found to belong to the same vegetative compatibility group (VCG) IT006 and were relatively common in all three countries. The results indicate that IT006 is endemic to southern Europe and may be utilized as an aflatoxin mitigation tool for maize across the region without concern for potential adverse impacts associated with the introduction of an exotic microorganism.

Seascape genomics of the pink abalone (Haliotis corrugata): An insight into a cross-border species in the northeast Pacific coast

Seascape genomics gives insight into the geographic and environmental factors shaping local adaptations. It improves the understanding of the potential effects of climate change, which is relevant to provide the basis for the international management of fishery resources. The pink abalone (Haliotis corrugata) is distributed from California, United States to Baja California Sur, Mexico, exposed to a latitudinal environmental gradient in the California Current System. Management of the pink abalone contrasts between Mexico and the United States; Mexico has an active fishery organized in four administrative areas, while the United States has kept the fishery in permanent closure since 1996. However, the impact of environmental factors on genetic variation along the species distribution remains unknown, and understanding this relationship is crucial for effective spatial management strategies. This study aims to investigate the neutral and adaptive genomic structure of H. corrugata. A total of 203 samples from 13 locations were processed using ddRADseq, and covering the species' distribution. Overall, 2,231 neutral, nine potentially adaptive and three genomic-environmental association loci were detected. The neutral structure identified two groups: 1) California, United States and 2) Baja California Peninsula, México. In addition, the adaptive structure analysis also detected two groups with genetic divergence observed at Punta Eugenia. Notably, the seawater temperature significantly correlated with the northern group (temperate) and the southern (warmer) group. This study is a valuable foundation for future research and conservation initiatives, emphasizing the importance of considering neutral and adaptive genetic factors when developing management strategies for marine species.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

Genetic erosion in an endangered desert fish during a megadrought despite long-term supportive breeding

Human water use combined with a recent megadrought have reduced river and stream flow through the southwest United States and led to periodic drying of formerly perennial river segments. Reductions in snowmelt runoff and increased extent of drying collectively threaten short-lived, obligate aquatic species, including the endangered Rio Grande silvery minnow (Hybognathus amarus). This species is subject to boom-and-bust population dynamics, under which large fluctuations in abundance are expected to lower estimates of effective population size and erode genetic diversity over time. Rates of diversity loss are also affected by additions of hatchery-origin fish used to supplement the wild population. We used demographic and genetic data from wild and hatchery individuals to examine the relationship of genetic diversity and effective population size to abundance over the last two decades. Genetic diversity was low during the early 2000s, but diversity and demographic metrics stabilized after the hatchery program was initiated and environmental conditions improved. Yet, from 2017 onward, allelic diversity declined (Cohen's d = 1.34) and remained low despite hatchery stocking and brief wild population recovery. Across the time series, single-sample estimates of effective population size based on linkage disequilibrium (LD Ne ) were positively associated (r = 0.53) with wild abundance and total abundance, but as the proportion of hatchery-origin spawners increased, LD Ne declined (r = -0.55). Megadrought limited wild spawner abundance and precluded refreshment of hatchery brood stocks with wild fish; hence, we predict a riverine population increasingly dominated by hatchery-origin individuals and accelerated loss of genetic diversity despite supplementation. We recommend an adaptive and accelerated management plan that integrates river flow management and hatchery operations to slow the pace of genetic diversity loss exacerbated by megadrought.

Habitat diversification associated with urban development has a little effect on genetic structure in the annual native plant Commelina communis in an East Asian megacity

Urban development greatly alters the natural and semi-natural habitats of native plants. Urbanisation results in a range of diverse habitats including remnant agricultural lands, urban parks, and roadside habitats. This habitat diversity often promotes trait divergence within urban areas. However, the mechanisms by which diverse urban habitats influence the population genetic structure of individual plant species remain poorly understood. We investigated the effects of urbanisation on genetic diversity and structure within 24 Commelina communis populations across diverse habitat types (rural agricultural land, urban agricultural land, urban park land, and urban roadsides) within the Kyoto-Osaka-Kobe megacity in Japan. We conducted multiplexed inter-simple sequence repeat genotyping to compare genetic diversity among populations in different habitats. We also examined the correlation between Nei's genetic distance and geographic and environmental distances and performed principal coordinate analysis (PCoA) to evaluate genetic differentiation among urban habitats. There were no significant differences in genetic diversity indices between urban and rural populations and among urban habitat types. Although we detected no isolation-by-distance structure in population pairs of the same habitat type and in those of different habitats, the difference in surrounding landscape facilitated genetic differentiation not only between urban and rural habitats but also between different urban habitats. PCoA revealed no clear genetic differentiation among rural and urban habitat populations. Our findings indicate that the establishment of diverse habitat types through urbanisation has no and little impact on genetic diversity and structure, respectively, in C. communis, likely due to its high selfing rate and ability to adapt to urban conditions.

How long have we been mistaken? Multi-tools shedding light into the systematics of the widespread deep-water genus Madrepora Linnaeus, 1758 (Scleractinia)

Deep-water coral reefs are found worldwide and harbor biodiversity levels that are comparable to their shallow-water counterparts. However, the genetic diversity and population structure of deep-water species remain poorly explored, and historical taxonomical issues still need to be resolved. Here we used microsatellite markers as well as ultraconserved elements (UCE) and exons to shed light on the population structure, genetic diversity, and phylogenetic position of the genus Madrepora, which contains M. oculata, one of the most widespread scleractinian species. Population structure of 107 samples from three Southwestern Atlantic sedimentary basins revealed the occurrence of a cryptic species, herein named M. piresae sp. nov. (authored by Kitahara, Capel and Zilberberg), which can be found in sympatry with M. oculata. Phylogeny reconstructions based on 134 UCEs and exon regions corroborated the population genetic data, with the recovery of two well-supported groups, and reinforced the polyphyly of the family Oculinidae. In order to better accommodate the genus Madrepora, while reducing taxonomical confusion associated with the name Madreporidae, we propose the monogeneric family Bathyporidae fam. nov. (authored by Kitahara, Capel, Zilberberg and Cairns). Our findings advance the knowledge on the widespread deep-water genus Madrepora, resolve a long-standing question regarding the phylogenetic position of the genus, and highlight the need of a worldwide review of the genus.

Discovery of ex situ individuals of Andrias sligoi, an extremely endangered species and one of the largest amphibians worldwide

The South China giant salamander, Andrias sligoi, is one of the largest extant amphibian species worldwide. It was recently distinguished from another Chinese species, the Chinese giant salamander, Andrias davidianus, which is considered Critically Endangered according to the International Union for Conservation of Nature (IUCN) Red List. It appears too late to save this extremely rare and large amphibian in situ. Another extant species of the same genus, Andrias japonicus, inhabits Japan. However, the introduction of Chinese giant salamanders into some areas of Japan has resulted in hybridization between the Japanese and Chinese species. During our genetic screening of giant salamanders in Japan, we unexpectedly discovered four individuals of the South China giant salamander: two were adult males in captivity, and one had recently died. The last individual was a preserved specimen. In this study, we report these extremely rare individuals of A. sligoi in Japan and discuss the taxonomic and conservational implications of these introduced individuals.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Single-cell genomics of a bloom-forming phytoplankton species reveals population genetic structure across continents

The study of microbial diversity over time and space is fundamental to the understanding of their ecology and evolution. The underlying processes driving these patterns are not fully resolved but can be studied using population genetic approaches. Here we investigated the population genetic structure of Gonyostomum semen, a bloom-forming phytoplankton species, across two continents. The species appears to be expanding in Europe, whereas similar trends are not observed in the USA. Our aim was to investigate if populations of Gonyostomum semen in Europe and in the USA are genetically differentiated, if there is population genetic structure within the continents, and what the potential drivers of differentiation are. To this end, we used a novel method based on single-amplified genomes combined with Restriction-site Associated DNA sequencing that allows de novo genotyping of natural single-cell isolates without the need for culturing. We amplified over 900 single-cell genomes from 25 lake populations across Europe and the USA and identified two distinct population clusters, one in Europe and another in the USA. Low genetic diversity in European populations supports the hypothesized recent expansion of Gonyostomum semen on this continent. Geographic population structure within each continent was associated with differences in environmental variables that may have led to ecological divergence of population clusters. Overall, our results show that single-amplified genomes combined with Restriction-site Associated DNA sequencing can be used to analyze microalgal population structure and differentiation based on single-cell isolates from natural, uncultured samples.

Genetic approaches reveal a healthy population and an unexpectedly recent origin for an isolated desert spring fish

Foskett Spring in Oregon's desert harbors a historically threatened population of Western Speckled Dace (Rhinichthys klamathensis). Though recently delisted, the dace's recruitment depends upon regular removal of encroaching vegetation. Previous studies assumed that Foskett Dace separated from others in the Warner Valley about 10,000 years ago, thereby framing an enigma about the population's surprising ability to persist for so long in a tiny habitat easily overrun by plants. To investigate that persistence and the effectiveness of interventions to augment population size, we assessed genetic diversity among daces inhabiting Foskett Spring, a refuge at Dace Spring, and three nearby streams. Analysis revealed a robust effective population size (Ne) of nearly 5000 within Foskett Spring, though Ne in the Dace Spring refuge is just 10% of that value. Heterozygosity is slightly lower than expected based on random mating at all five sites, indicating mild inbreeding, but not at a level of concern. These results confirm the genetic health of Foskett Dace. Unexpectedly, genetic differentiation reveals closer similarity between Foskett Dace and a newly discovered population from Nevada's Coleman Creek than between Foskett Dace and dace elsewhere in Oregon. Demographic modeling inferred Coleman Creek as the ancestral source of Foskett Dace fewer than 1000 years ago, much more recently than previously suspected and possibly coincident with the arrival of large herbivores whose grazing may have maintained open water suitable for reproduction. These results solve the enigma of persistence by greatly shortening the duration over which Foskett Dace have inhabited their isolated spring.

Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

Isolated Grauer's gorilla populations differ in diet and gut microbiome

The animal gut microbiome has been implicated in a number of key biological processes, ranging from digestion to behaviour, and has also been suggested to facilitate local adaptation. Yet studies in wild animals rarely compare multiple populations that differ ecologically, which is the level at which local adaptation may occur. Further, few studies simultaneously characterize diet and gut microbiome from the same sample, despite their probable interdependence. Here, we investigate the interplay between diet and gut microbiome in three geographically isolated populations of the critically endangered Grauer's gorilla (Gorilla beringei graueri), which we show to be genetically differentiated. We find population- and social group-specific dietary and gut microbial profiles and covariation between diet and gut microbiome, despite the presence of core microbial taxa. There was no detectable effect of age, and only marginal effects of sex and genetic relatedness on the microbiome. Diet differed considerably across populations, with the high-altitude population consuming a lower diversity of plants compared to low-altitude populations, consistent with plant availability constraining dietary choices. The observed pattern of covariation between diet and gut microbiome is probably a result of long-term social and environmental factors. Our study suggests that the gut microbiome is sufficiently plastic to support flexible food selection and hence contribute to local adaptation.

Molecular phylogeography and historical demography of a widespread herbaceous species from eastern North America, Podophyllum peltatum

PREMISE

Glacial/interglacial cycles and topographic complexity are both considered to have shaped today's diverse phylogeographic patterns of taxa from unglaciated eastern North America (ENA). However, few studies have focused on the phylogeography and population dynamics of wide-ranging ENA herbaceous species occurring in forest understory habitat. We examined the phylogeographic pattern and evolutionary history of Podophyllum peltatum L., a widely distributed herb inhabiting deciduous forests of ENA.

METHODS

Using chloroplast DNA (cpDNA) sequences and nuclear microsatellite loci, we investigated the population structure and genetic diversity of the species. Molecular dating, demographic history analyses, and ecological niche modeling were also performed to illustrate the phylogeographic patterns.

RESULTS

Our cpDNA results identified three main groups that are largely congruent with boundaries along the Appalachian Mountains and the Mississippi River, two major geographic barriers in ENA. Populations located to the east of the Appalachians and along the central Appalachians exhibited relatively higher levels of genetic diversity. Extant lineages may have diverged during the late Miocene, and range expansions of different groups may have happened during the Pleistocene glacial/interglacial cycles.

CONCLUSIONS

Our findings indicate that geographic barriers may have started to facilitate the population divergence in P. peltatum before the Pleistocene. Persistence in multiple refugia, including areas around the central Appalachians during the Quaternary glacial period, and subsequent expansions under hospitable climatic condition, especially westward expansion, are likely responsible for the species' contemporary genetic structure and phylogeographic pattern.

Functional and evolutionary analysis of host Synaptogyrin-2 in porcine circovirus type 2 susceptibility

Mammalian evolution has been influenced by viruses for millions of years, leaving signatures of adaptive evolution within genes encoding for viral interacting proteins. Synaptogyrin-2 (SYNGR2) is a transmembrane protein implicated in promoting bacterial and viral infections. A genome-wide association study of pigs experimentally infected with porcine circovirus type 2b (PCV2b) uncovered a missense mutation (SYNGR2 p.Arg63Cys) associated with viral load. In this study, CRISPR/Cas9-mediated gene editing of the porcine kidney 15 (PK15, wtSYNGR2+p.63Arg) cell line generated clones homozygous for the favorable SYNGR2 p.63Cys allele (emSYNGR2+p.63Cys). Infection of edited clones resulted in decreased PCV2 replication compared to wildtype PK15 (P<0.05), with consistent effects across genetically distinct PCV2b and PCV2d isolates. Sequence analyses of wild and domestic pigs (n>700) revealed the favorable SYNGR2 p.63Cys allele is unique to domestic pigs and more predominant in European than Asian breeds. A haplotype defined by the SYNGR2 p.63Cys allele was likely derived from an ancestral haplotype nearly fixed within European (0.977) but absent from Asian wild boar. We hypothesize that the SYNGR2 p.63Cys allele arose post-domestication in ancestral European swine. Decreased genetic diversity in homozygotes for the SYNGR2 p.63Cys allele compared to SYNGR2 p.63Arg, corroborates a rapid increase in frequency of SYGNR2 p.63Cys via positive selection. Signatures of adaptive evolution across mammalian species were also identified within SYNGR2 intraluminal loop domains, coinciding with the location of SYNGR2 p.Arg63Cys. Therefore, SYNGR2 may reflect a novel component of the host-virus evolutionary arms race across mammals with SYNGR2 p.Arg63Cys representing a species-specific example of putative adaptive evolution.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

Identification of local camel populations in Turkiye using morphological, genetic, and breeding info

The aim of this study was to investigate the demographic, morphological, and genetic characteristics of local camel populations reared in the Turkiye provinces of Aydin, Denizli, and Antalya, which have a long history of camel breeding. Although Turkiye has an old history of camel breeding in its historical process, the number of scientific studies aimed at identifying camel populations in Turkiye is almost negligible. In this study, local camel populations in Aydin, Denizli, and Antalya cities of Turkiye were examined in three dimensions as demographic, morphological, and genetic. A face-to-face survey of 117 breeders was used to determine demographic definitions. While the region where local camels were detected the most was determined as Antalya Region with 78.6%, it was determined that 82.6% of the breeders participating in the survey preferred to breed camels due to their docile temperament. Body measurements were made on 45 camels for morphological identification. Moreover, DNA were sampled with oral swabs from 57 camels for phylogenetic analyses using 16 SSR microsatellite loci to identify the genetic structure of local camel populations. The genetic analyses using SSR markers revealed that the camel populations in the Antalya region had a considerably more isolated genetic structure than the Aydin and Denizli populations, and consequently, these populations may be regarded the native camel population.

Genetic imprints of grafting in wild iron walnut populations in southwestern China

BACKGROUND

Anthropogenic activities are causing unprecedented loss of genetic diversity in many species. However, the effects on genetic diversity from large-scale grafting onto wild plants of crop species are largely undetermined. Iron walnut (Juglans sigillata Dode) is a deciduous nut tree crop endemic to southwestern China with a long history of cultivation. Due to the rapid expansion of the walnut industry, many natural populations are now being replaced by cultivars grafted onto wild rootstocks. However, little is known about the potential genetic consequences of such action on natural populations.

RESULTS

We sampled the scion and the rootstock from each of 149 grafted individuals within nine wild populations of J. sigillata from Yunnan Province which is the center of walnut diversity and cultivation in China, and examined their genetic diversity and population structure using 31 microsatellite loci. Scions had lower genetic diversity than rootstocks, and this pattern was repeated in seven of the nine examined populations. Among those seven populations, AMOVA and clustering analyses showed a clear genetic separation between all rootstocks and all scions. However, the two remaining populations, both from northern Yunnan, showed genetic similarity between scions and rootstocks, possibly indicating that wild populations here are derived from feralized local cultivars. Moreover, our data indicated probable crop-to-wild gene flow between scions and rootstocks, across all populations.

CONCLUSIONS

Our results indicate that large-scale grafting has been causing genetic diversity erosion and genetic structure breakdown in the wild material of J. sigillata within Yunnan. To mitigate these effects, we caution against the overuse of grafting in wild populations of iron walnut and other crop species and recommend the preservation of natural genotypes through in situ  and ex situ conservation.

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.

A Comparative Phylogeography of Three Marine Species with Different PLD Modes Reveals Two Genetic Breaks across the Southern Caribbean Sea

The comparative phylogeography of marine species with contrasting dispersal potential across the southern Caribbean Sea was evaluated by the presence of two putative barriers: the Magdalena River plume (MRP) and the combination of the absence of a rocky bottom and the almost permanent upwelling in the La Guajira Peninsula (ARB + PUG). Three species with varying biological and ecological characteristics (i.e., dispersal potentials) that inhabit shallow rocky bottoms were selected: Cittarium pica (PLD < 6 days), Acanthemblemaria rivasi (PLD < 22 days), and Nerita tessellata (PLD > 60 days). We generated a set of SNPs for the three species using the ddRad-seq technique. Samples of each species were collected in five locations from Capurganá to La Guajira. For the first time, evidence of a phylogeographic break caused by the MRP is provided, mainly for A. rivasi (AMOVA: ΦCT = 0.420). The ARB + PUG barrier causes another break for A. rivasi (ΦCT = 0.406) and C. pica (ΦCT = 0.224). Three populations (K = 3) were identified for A. rivasi and C. pica, while N. tessellata presented one population (K = 1). The Mantel correlogram indicated that A. rivasi and C. pica fit the hierarchical population model, and only the A. rivasi and C. pica comparisons showed phylogeographic congruence. Our results demonstrate how the biological traits of these three species and the biogeographic barriers have influenced their phylogeographic structure.

Clones on the run: The genomics of a recently expanded partially clonal species

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Exceptional population genomic homogeneity in the black brittle star Ophiocomina nigra (Ophiuroidea, Echinodermata) along the Atlantic-Mediterranean coast

The Atlantic-Mediterranean marine transition is characterised by strong oceanographic barriers and steep environmental gradients that generally result in connectivity breaks between populations from both basins and may lead to local adaptation. Here, we performed a population genomic study of the black brittle star, Ophiocomina nigra, covering most of its distribution range along the Atlantic-Mediterranean region. Interestingly, O. nigra is extremely variable in its coloration, with individuals ranging from black to yellow-orange, and different colour morphs inhabiting different depths and habitats. In this work, we used a fragment of the mitochondrial COI gene and 2,374 genome-wide ddRADseq-derived SNPs to explore: (a) whether the different colour morphs of O. nigra represent different evolutionary units; (b) the disruptive effects of major oceanographic fronts on its population structure; and (c) genomic signals of local adaptation to divergent environments. Our results revealed exceptional population homogeneity, barely affected by oceanographic fronts, with no signals of local adaptation nor genetic differentiation between colour morphs. This remarkable panmixia likely results from a long pelagic larval duration, a large effective population size and recent demographic expansions. Our study unveils an extraordinary phenotypic plasticity in O. nigra, opening further research questions on the ecological and molecular mechanisms underpinning coloration in Ophiuroidea.

Climate and Soil Microsite Conditions Determine Local Adaptation in Declining Silver Fir Forests

Ongoing climatic change is threatening the survival of drought-sensitive tree species, such as silver fir (Abies alba). Drought-induced dieback had been previously explored in this conifer, although the role played by tree-level genetic diversity and its relationship with growth patterns and soil microsite conditions remained elusive. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to describe different genetic characteristics of five silver fir forests in the Spanish Pyrenees, including declining and non-declining trees. Single nucleotide polymorphisms (SNPs) were used to investigate the relationships between genetics, dieback, intraspecific trait variation (functional dendrophenotypic traits and leaf traits), local bioclimatic conditions, and rhizosphere soil properties. While there were no noticeable genetic differences between declining and non-declining trees, genome-environment associations with selection signatures were abundant, suggesting a strong influence of climate, soil physicochemical properties, and soil microbial diversity on local adaptation. These results provide novel insights into how genetics and diverse environmental factors are interrelated and highlight the need to incorporate genetic data into silver fir forest dieback studies to gain a better understanding of local adaptation.

Genetic diversity and forensic parameters of autosomal STR markers included PowerPlex® CS7 panel in the general population of the Russian Federation

PowerPlex® CS7 multiplex is commonly used as a source of supplementary markers in parentage and kinship studies. We analysed a total of 687 unrelated individuals from 94 geographically different localities across all Federal Districts of the Russian Federation and provide forensically important parameters and allele frequencies. The paper also presents results of an intra-population genetic diversity study between the populations of the Federal Districts and comparison with populations from various regions of the world.

Non-random mating within an Island rookery of Hawaiian hawksbill turtles: demographic discontinuity at a small coastline scale

Hawksbill sea turtles (Eretmochelys imbricata) from the Hawaiian archipelago form a small and genetically isolated population, consisting of only a few tens of individuals breeding annually. Most females nest on the island of Hawai'i, but little is known about the demographics of this rookery. This study used genetic relatedness, inferred from 135 microhaplotype markers, to determine breeding sex-ratios, estimate female nesting frequency and assess relationships between individuals nesting on different beaches. Samples were collected during the 2017 nesting season and final data included 13 nesting females and 1002 unhatched embryos, salvaged from 41 nests, of which 13 had no observed mother. Results show that most females used a single nesting beach laying 1-5 nests each. From female and offspring alleles, the paternal genotypes of 12 breeding males were reconstructed and many showed high relatedness to their mates. Pairwise relatedness of offspring revealed one instance of polygyny but otherwise suggested a 1 : 1 breeding-sex ratio. Relatedness analysis and spatial-autocorrelation of genotypes indicate that turtles from different nesting areas do not regularly interbreed, suggesting that strong natal homing tendencies in both sexes result in non-random mating across the study area. Complexes of nearby nesting beaches also showed unique patterns of inbreeding across loci, further indicating that Hawaiian hawksbill turtles have demographically discontinuous nesting populations separated by only tens of km.

Drivers of strong isolation and small effective population size at a leading range edge of a widespread plant

Climate change has influenced species distributions worldwide with upward elevational shifts observed in many systems. Leading range edge populations, like those at upper elevation limits, are crucial for climate change responses but can exhibit low genetic diversity due to founder effects, isolation, or limited outbreeding. These factors can hamper local adaptation at range limits. Using the widespread herb, Argentina anserina, we measured ecological attributes (population density on the landscape, area of population occupancy, and plant and flower density) spanning a 1000 m elevation gradient, with high elevation populations at the range limit. We measured vegetative clonal potential in the greenhouse for populations spanning the gradient. We combined these data with a ddRAD-seq dataset to test the hypotheses that high elevation populations would exhibit ecological and genomic signatures of leading range edge populations. We found that population density on the landscape declined towards the high elevation limit, as is expected towards range edges. However, plant density was elevated within edge populations. In the greenhouse, high elevation plants exhibited stronger clonal potential than low elevation plants, likely explaining increased plant density in the field. Phylogeographic analysis supported more recent colonization of high elevation populations which were also more genetically isolated, had more extreme heterozygote excess and had smaller effective population size than low. Results support that colonization of high elevations was likely accompanied by increased asexuality, contributing to a decline in effective population size. Despite high plant density in leading edge populations, their small effective size, isolation and clonality could constrain adaptive potential.

Population Structure and Association Mapping for Agronomical and Biochemical Traits of a Large Spanish Apple Germplasm

A basic knowledge of linkage disequilibrium and population structure is necessary in order to determine the genetic control and identify significant associations with agronomical and phytochemical compounds in apple (Malus × domestica Borkh). In this study, 186 apple accessions (Pop1), representing both Spanish native accessions (94) and non-Spanish cultivars (92) from the EEAD-CSIC apple core collection, were assessed using 23 SSRs markers. Four populations were considered: Pop1, Pop2, Pop3, and Pop4. The initial Pop1 was divided into 150 diploid (Pop2) and 36 triploid accessions (Pop3), while for the inter-chromosomal linkage disequilibrium and the association mapping analysis, 118 phenotype diploid accessions were considered Pop4. Thus, the average number of alleles per locus and observed heterozygosity for the overall sample set (Pop1) were 15.65 and 0.75, respectively. The population structure analysis identified two subpopulations in the diploid accessions (Pop2 and Pop4) and four in the triploids (Pop3). Regarding the Pop4, the population structure with K = 2 subpopulations segregation was in agreement with the UPGMA cluster analysis according to the genetic pairwise distances. Moreover, the accessions seemed to be segregated by their origin (Spanish/non-Spanish) in the clustering analysis. One of the two subpopulations encountered was quite-exclusively formed by non-Spanish accessions (30 out of 33). Furthermore, agronomical and basic fruit quality parameters, antioxidant traits, individual sugars, and organic acids were assessed for the association mapping analysis. A high level of biodiversity was exhibited in the phenotypic characterization of Pop4, and a total of 126 significant associations were found between the 23 SSR markers and the 21 phenotypic traits evaluated. This study also identified many new marker-locus trait associations for the first time, such as in the antioxidant traits or in sugars and organic acids, which may be useful for predictions and for a better understanding of the apple genome.

Multiple introductions, polyploidy and mixed reproductive strategies are linked to genetic diversity and structure in the most widespread invasive plant across Southern Ocean archipelagos

Biological invasions in remote areas that experience low human activity provide unique opportunities to elucidate processes responsible for invasion success. Here we study the most widespread invasive plant species across the isolated islands of the Southern Ocean, the annual bluegrass, Poa annua. To analyse geographical variation in genome size, genetic diversity and reproductive strategies, we sampled all major sub-Antarctic archipelagos in this region and generated microsatellite data for 470 individual plants representing 31 populations. We also estimated genome sizes for a subset of individuals using flow cytometry. Occasional events of island colonization are expected to result in high genetic structure among islands, overall low genetic diversity and increased self-fertilization, but we show that this is not the case for P. annua. Microsatellite data indicated low population genetic structure and lack of isolation by distance among the sub-Antarctic archipelagos we sampled, but high population structure within each archipelago. We identified high levels of genetic diversity, low clonality and low selfing rates in sub-Antarctic P. annua populations (contrary to rates typical of continental populations). In turn, estimates of selfing declined in populations as genetic diversity increased. Additionally, we found that most P. annua individuals are probably tetraploid and that only slight variation exists in genome size across the Southern Ocean. Our findings suggest multiple independent introductions of P. annua into the sub-Antarctic, which promoted the establishment of genetically diverse populations. Despite multiple introductions, the adoption of convergent reproductive strategies (outcrossing) happened independently in each major archipelago. The combination of polyploidy and a mixed reproductive strategy probably benefited P. annua in the Southern Ocean by increasing genetic diversity and its ability to cope with the novel environmental conditions.

Chloroplast genome sequence of triploid Toxicodendron vernicifluum and comparative analyses with other lacquer chloroplast genomes

BACKGROUND

Toxicodendron vernicifluum, belonging to the family Anacardiaceae, is an important commercial arbor species, which can provide us with the raw lacquer, an excellent adhesive and painting material used to make lacquer ware. Compared with diploid, triploid lacquer tree has a higher yield of raw lacquer and stronger resistance to stress. Triploid T. vernicifluum was a newly discovered natural triploid lacquer tree. However, the taxonomy of triploid T. vernicifluum has remained uncertain. Here, we sequenced and analyzed the complete chloroplast (cp) genome of triploid T. vernicifluum and compared it with related species of Toxicodendron genus based on chloroplast genome and SSR markers.

RESULTS

The plastome of triploid T. vernicifluum is 158,221 bp in length, including a pair of inverted repeats (IRs) of 26,462 bp, separated by a large single-copy region of 86,951 bp and a small single-copy region of 18,346 bp. In total, 132 genes including 87 protein-coding genes, 37 tRNA genes and 8 rRNA genes were identified in the triploid T. vernicifluum. Among these, 16 genes were duplicated in the IR regions, 14 genes contain one intron, while three genes contain two introns. After nucleotide substitutions, seven small inversions were analyzed in the chloroplast genomes, eight hotspot regions were found, which could be useful molecular genetic markers for future population genetics. Phylogenetic analyses showed that triploid T. vernicifluum was a sister to T. vernicifluum cv. Dahongpao and T. vernicifluum cv. Hongpigaobachi. Moreover, phylogenetic clustering based on the SSR markers showed that all the samples of triploid T. vernicifluum, T. vernicifluum cv. Dahongpao and T. vernicifluum cv. Hongpigaobachi in one group, while the samples of T. vernicifluum and T. succedaneum in another group, which is consistent with the cp genome and morphological analysis.

CONCLUSIONS

The current genomic datasets provide pivotal genetic resources to determine the phylogenetic relationships, variety identification, breeding and resource exploitation, and future genetic diversity-related studies of T. vernicifluum.

Genetic characteristics influence the phenotype of marine macroalga Fucus vesiculosus (Phaeophyceae)

Intraspecific variation is an important component of heterogeneity in biological systems that can manifest at the genotypic and phenotypic level. This study investigates the influence of genetic characteristics on the phenotype of free-living Fucus vesiculosus using traditional morphological measures and microsatellite genotyping. Two sympatric morphotypes were observed to be significantly genetically and morphologically differentiated despite experiencing analogous local environmental conditions; indicating a genetic element to F. vesiculosus morphology. Additionally, the observed intraclonal variation established divergent morphology within some genets. This demonstrated that clonal lineages have the ability to alter morphological traits by either a plastic response or somatic mutations. We provide support for the potential occurrence of the Gigas effect (cellular/organ enlargement through genome duplication) in the Fucus genus, with polyploidization appearing to correlate with a general increase in the size of morphological features. Phenotypic traits, as designated by morphology within the study, of F. vesiculosus are partially controlled by the genetic characteristics of the thalli. This study suggests that largely asexually reproducing algal populations may have the potential to adapt to changing environmental conditions through genome changes or phenotypic plasticity.

Characterization of 35 new microsatellite markers for the blacktip reef shark (Carcharhinus melanopterus) and cross-species amplification in eight other shark species

BACKGROUND

Shark species are overfished at a global scale, as they are poached for the finning industry or are caught as bycatch. Efficient conservation measures require fine-scale spatial and temporal studies to characterize shark habitat use, infer migratory habits, analyze relatedness, and detect population genetic differentiation. Gathering these types of data is costly and time-consuming, especially when it requires collection of shark tissue samples.

METHODS AND RESULTS

Genetic tools, such as microsatellite markers, are the most economical sampling method for collecting genetic data, as they enable the estimation of genetic diversity, population structure and parentage relationships and are thus an efficient way to inform conservation strategies. Here, a set of 45 microsatellite loci was tested on three blacktip reef shark (Carcharhinus melanopterus) populations from three Polynesian islands: Moorea, Morane and Tenararo. The set was composed of 10 previously published microsatellite markers and 35 microsatellite markers that were developed specifically for C. melanopterus as part of the present study. The 35 novel and 10 existing loci were cross-amplified on eight additional shark species (Carcharhinus amblyrhynchos, C. longimanus, C. sorrah, Galeocerdo cuvier, Negaprion acutidens, Prionacea glauca, Rhincodon typus and Sphyrna lewini). These species had an average of 69% of successful amplification, considered if at least 50% of the individual samples being successfully amplified per species and per locus.

CONCLUSIONS

This novel microsatellite marker set will help address numerous knowledge gaps that remain, concerning genetic stock identification, shark behavior and reproduction via parentage analysis.

The contrary conservation situations of two local critically endangered species, Vaccinium emarginatum (Ericaceae) and Elatostema platyphyllum (Urticaceae), growing on the eastern edge of the distribution

As biodiversity loss continues, there is an urgent need to develop efficient conservation measures to protect diversity with limited conservation resources. Conservation targets have generally been selected based on their population size, but more detailed assessments clarifying the phylogenetic genetic status, history, and phylogenetic uniqueness of rare species is crucial to set more appropriate and effective conservation measures. In Japan, the Act on Conservation of Endangered Species of Wild Fauna and Flora designated endangered plants with high conservation priority, but &gt;40% of these species also grow overseas. We conducted comparative analyses based on ddRADseq and MIG-seq to evaluate the population conservation status and value of Vaccinium emarginatum and Elatostema platyphyllum which are growing across national borders at the eastern edge of their species distribution range. The analyses revealed contrasting conservation status between the two species; the Japanese population of V. emarginatum had lower genetic diversity at the individual level and phylogenetically differentiated from Taiwanese populations, while that of E. platyphyllum had higher diversity at the individual level and is a relatively recent migrant with little phylogenetical differentiation from Taiwanese populations. The two species, which share the common feature of being critically rare in Japan, showed contrasting genetic/phylogenetic characteristics. This study provided useful information for appropriate conservation measures based on species’ phylogenetic traits and genetic diversity.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

Genomic evidence for three distinct species in the Erebia manto complex in Central Europe (Lepidoptera, Nymphalidae)

A problem to implement conservation strategies is that in many cases recognized taxa are in fact complexes of several cryptic species. Failure to properly delineate species may lead to misplaced priorities or to inadequate conservation measures. One such species complex is the yellow-spotted ringlet Erebia manto, which comprises several phenotypically distinct lineages, whose degree of genomic isolation has so far not been assessed. Some of these lineages are geographically restricted and thus possibly represent distinct units with conservation priorities. Using several thousand nuclear genomic markers, we evaluated to which degree the bubastis lineage from the Alps and the vogesiaca lineage from the Vosges, are genetically isolated from the widespread manto lineage. Our results suggest that both lineages are genetically as strongly differentiated from manto as other taxonomically well separated sibling species in this genus from each other, supporting a delineation of bubastis and vogesiaca as independent species. Given the restricted and isolated range of vogesiaca as well as the disjunct distribution of bubastis, our findings have significant implication for future conservation efforts on these formerly cryptic species and highlight the need to investigate the genomic identity within species complexes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01501-w.

A low-density single nucleotide polymorphism panel for brown trout (Salmo trutta L.) suitable for exploring genetic diversity at a range of spatial scales

The rivers of southern England and northern France which drain into the English Channel contain several genetically unique groups of trout (Salmo trutta L.) that have suffered dramatic declines in numbers over the past 40 years. Knowledge of levels and patterns of genetic diversity is essential for effective management of these vulnerable populations. Using restriction site-associated DNA sequencing (RADseq) data, we describe the development and characterisation of a panel of 95 single nucleotide polymorphism (SNP) loci for trout from this region and investigate their applicability and variability in both target (i.e., southern English) and non-target trout populations from northern Britain and Ireland. In addition, we present three case studies which demonstrate the utility and resolution of these genetic markers at three levels of spatial separation:(a) between closely related populations in nearby rivers, (b) within a catchment and (c) when determining parentage and familial relationships between fish sampled from a single site, using both empirical and simulated data. The SNP loci will be useful for population genetic and assignment studies on brown trout within the UK and beyond.

Temporal Monitoring of the Floreana Island Galapagos Giant Tortoise Captive Breeding Program

Captive breeding programs benefit from genetic analyses that identify relatedness between individuals, assign parentage to offspring, and track levels of genetic diversity. Monitoring these parameters across breeding cycles is critical to the success of a captive breeding program as it allows conservation managers to iteratively evaluate and adjust program structure. However, in practice, genetic tracking of breeding outcomes is rarely conducted. Here, we examined the first three offspring cohorts (2017-2020) of the genetically informed captive breeding program for the Floreana Island Galapagos giant tortoise, Chelonoidis niger. This captive breeding program is unique as the Floreana tortoise has been extinct since the 1800s, but its genome persists, in part, in the form of living hybrids with the extant Volcano Wolf tortoise, C. becki. Breeding over the study period took place at the Galapagos National Park Directorate breeding facility in four corrals, each containing three females and two males. Using 17 microsatellite markers, we were able to assign parentage to 94 of the 98 offspring produced over the study period. We observe that despite the addition of more founders since the pilot breeding program, the effective population size remains low, and changes to the arrangements of breeding corrals may be necessary to encourage more equal reproductive output from the males. This study demonstrates the value of hybrids for species restoration and the importance of continually reassessing the outcomes of captive breeding.

A new GTSeq resource to facilitate multijurisdictional research and management of walleye Sander vitreus

Conservation and management professionals often work across jurisdictional boundaries to identify broad ecological patterns. These collaborations help to protect populations whose distributions span political borders. One common limitation to multijurisdictional collaboration is consistency in data recording and reporting. This limitation can impact genetic research, which relies on data about specific markers in an organism's genome. Incomplete overlap of markers between separate studies can prevent direct comparisons of results. Standardized marker panels can reduce the impact of this issue and provide a common starting place for new research. Genotyping-in-thousands (GTSeq) is one approach used to create standardized marker panels for nonmodel organisms. Here, we describe the development, optimization, and early assessments of a new GTSeq panel for use with walleye (Sander vitreus) from the Great Lakes region of North America. High genome-coverage sequencing conducted using RAD capture provided genotypes for thousands of single nucleotide polymorphisms (SNPs). From these markers, SNP and microhaplotype markers were chosen, which were informative for genetic stock identification (GSI) and kinship analysis. The final GTSeq panel contained 500 markers, including 197 microhaplotypes and 303 SNPs. Leave-one-out GSI simulations indicated that GSI accuracy should be greater than 80% in most jurisdictions. The false-positive rates of parent-offspring and full-sibling kinship identification were found to be low. Finally, genotypes could be consistently scored among separate sequencing runs >94% of the time. Results indicate that the GTSeq panel that we developed should perform well for multijurisdictional walleye research throughout the Great Lakes region.

Cryptic population structure at the northern range margin of the service tree Sorbus domestica

Climate change has aroused interest in planting warm- and drought-adapted trees in managed forests and urban areas. An option is to focus on tree species that occur naturally, but have centers of distribution in warmer and drier areas. However, in order to protect the species pool of genetic diversity, efforts of planting and promotion should be informed by knowledge on the local genetic diversity. Here, we studied the macro- and micro-scale population genetic structure of the rare European fruit tree Sorbus domestica at its northern range margin, in western Switzerland. New microsatellite data were combined with published data from across the European distribution of the species. Analyses revealed the presence of mainly one of two species-wide ancestral clusters, i.e., the western European cluster, with evidence that it consists of two cryptic sub-clusters. Average pairwise F _ST of 0.118 was low across the range, and only allelic richness was reduced in the northern margin compared to more southern and southeastern areas of Europe. Based on our finding of considerable genetic diversity of the species in western and northern Switzerland, we suggest that a national propagation program should focus on collecting seeds from natural, high-density tree stands and propagate locally. More generally, our study shows that rare tree species in marginal areas of their distributions do not necessarily have low genetic diversity or heightened levels of inbreeding, and in those cases probably need no assisted migration in efforts to propagate them.

Molecular and Phytochemical Variability of Endemic Juniperus sabina var. balkanensis from Its Natural Range

Juniperus sabina L. var. balkanensis R.P. Adams & Tashev is a recently described endemic variety from the Balkan Peninsula. Its strong sprouting ability and fast vegetative propagation, on one hand, and fragmented distribution, on the other, can lead to lower genetic diversity in local populations and to the differentiation of populations. As there has been no detailed investigation of this variety, we studied Balkan natural populations using phytochemical and molecular markers. Leaf essential oils (EOs) were chosen based on their proven usability in the population studies of Juniperus taxa, while ISSRs (Inter Simple Sequence Repeats) have been used due to their high resolution. In addition, since this variety is best described using molecular markers, the chloroplast trnS-trnG region was amplified from individuals from different populations having different chemotypes. Based on the essential oil profile, three chemotypes could be identified with a difference in their distribution. The analysis of molecular variance showed moderate differentiation of populations and regions, attesting to the start of the separation of three regions in the Balkans: west, east and south. The bioclimatic and environmental parameters and sex of the individual did not influence the EO profile, although some of the compounds present in low-to-medium concentrations showed strong correlation with several bioclimatic parameters.

Genomic Analyses Implicate the Amazon-Orinoco Plume as the Driver of Cryptic Speciation in a Swimming Crab

The Amazon-Orinoco plume (AOP) is the world's largest freshwater and sediment discharge into the ocean. Previous studies limited to mtDNA suggest that the swimming crab Callinectes ornatus Ordway, 1863 exists as two distinct genetic clusters separated by the AOP. However, questions concerning migration, diversification time, and species delimitation are unresolved. Densely sampling markers across the genome (SNPs) could elucidate the evolutionary processes within this species. Here, we combined mtDNA data and ddRAD-seq to explore the diversification patterns and processes within the swimming crab C. ornatus. We show great genetic differentiation between groups on the north and south sides of the plume but also signs of hybridization. Demographic modeling indicates the divergence between groups starting around 8 Mya following the AOP's formation. After a period of isolation, we detect two incidences of secondary contact with stronger migration in concordance with the North Brazil Current flow. Our results suggest speciation with gene flow explained by the interplay among the AOP, oceanographic currents, and long larval dispersal. This work represents the first investigation employing ddRAD-seq in a marine invertebrate species with distribution encompassing the north and south Atlantic and sheds light on the role of the AOP in the diversification of a marine species.

Population Structure of White Sturgeon (Acipenser transmontanus) in the Columbia River Inferred from Single-Nucleotide Polymorphisms

White sturgeon (Acipenser transmontanus) are the largest freshwater fish in North America, with reproducing populations in the Sacramento-San Joaquin, Fraser, and Columbia River Basins. Of these, the Columbia River is the largest, but it is also highly fragmented by hydroelectric dams, and many segments are characterized by declining abundance and persistent recruitment failure. Efforts to conserve and supplement these fish requires an understanding of their spatial genetic structure. Here, we assembled a large set of samples from throughout the Columbia River Basin, along with representative collections from adjacent basins, and genotyped them using a panel of 325 single-nucleotide markers. Results from individual- and group-based analyses of these data indicate that white sturgeon in the uppermost Columbia River Basin, in the Kootenai and upper Snake Rivers, are the most distinct, while the remaining populations downstream in the basin can be described as a genetic gradient consistent with an isolation-by-distance effect. Notably, the population in the lowest reaches of the Columbia River is more distinct from the middle or upper reaches than from outside basins, and suggests historically a higher or more recent gene exchange through coastal routes than with populations in the interior Columbia Basin. Nonetheless, proximal reaches were generally only marginally or non-significantly divergent, suggesting that transplanting larvae or juveniles from nearby sources poses relatively little risk of outbreeding depression. Indeed, we inferred examples of dispersal between reaches via close-kin mark-recapture and genetic mark-recapture that indicate movement between nearby reaches is not unusual. Samples from the Kootenai and upper Snake Rivers exhibited notably lower genetic diversity than the remaining samples as a result of population bottlenecks, genetic drift, and/or historical divergence. Conservation actions, such as supplementation, are underway to maintain population viability and will require balanced efforts to increase demographic abundance while maintaining genetic diversity.

Population genetics of an invasive mosquito vector, Aedes albopictus in the Northeastern USA

The Asian tiger mosquito (Aedes albopictus) arrived in the USA in the 1980’s and rapidly spread throughout eastern USA within a decade. The predicted northern edge of its overwintering distribution on the East Coast of the USA roughly falls across New York, Connecticut, and Massachusetts, where the species has been recorded as early as 2000. It is unclear whether Ae. albopictus populations have become established and survive the cold winters in these areas or are recolonized every year. We genotyped and analyzed populations of Ae. albopictus from the northeast USA using 15 microsatellite markers and compared them with other populations across the country and to representatives of the major global genetic clades to investigate their connectivity and stability. Founder effects or bottlenecks were rare at the northern range of the Ae. albopictus distribution in the northeastern USA, with populations displaying high levels of genetic diversity and connectivity along the East Coast. There is no evidence of population turnover in Connecticut during the course of three consecutive years, with consistent genetic structure throughout this period. Overall, these results support the presence of established populations of Ae. albopictus in New York, Connecticut, and Massachusetts, successfully overwintering and migrating in large numbers. Given the stability and interconnectedness of these populations, Ae. albopictus has the potential to continue to proliferate and expand its range northward under mean warming conditions of climate change. Efforts to control Ae. albopictus in these areas should thus focus on vector suppression rather than eradication strategies, as local populations have become firmly established and are expected to reemerge every summer.

Assessment of the Genetic Distinctiveness and Uniformity of Pre-Basic Seed Stocks of Italian Ryegrass Varieties

Lolium multiflorum Lam., commonly known as Italian ryegrass, is a forage grass mostly valued for its high palatability and digestibility, along with its high productivity. However, Italian ryegrass has an outbreeding nature and therefore has high genetic heterogeneity within each variety. Consequently, the exclusive use of morphological descriptors in the existing varietal identification and registration process based on the Distinctness, Uniformity, and Stability (DUS) test results in an inadequately precise assessment. The primary objective of this work was to effectively test whether the uniformity observed at the phenological level within each population of Italian ryegrass was confirmed at the genetic level through an SSR marker analysis. In this research, using 12 polymorphic SSR loci, we analyzed 672 samples belonging to 14 different Italian ryegrass commercial varieties to determine the pairwise genetic similarity (GS), verified the distribution of genetic diversity within and among varieties, and investigated the population structure. Although the fourteen commercial varieties did not show elevated genetic differentiation, with only 13% of the total variation attributable to among-cultivar genetic variation, when analyzed as a core, each variety constitutes a genetic cluster on its own, resulting in distinct characteristics from the others, except for two varieties. In this way, by combining a genetic tool with the traditional morphological approach, we were able to limit biases linked to the environmental effect of field trials, assessing the real source of diversity among varieties and concretely answering the key requisites of the Plant Variety Protection (PVP) system.

Genetic analysis of walnut cultivars from southwest China: Implications for germplasm improvement

Walnuts are highly valued for their rich nutritional profile and wide medicinal applications. This demand has led to the intensification of breeding activities in major walnut production areas such as southwest China, in order to develop more superior cultivars. With the increasing number of cultivars, accurate identification becomes fundamental to selecting the right cultivar for grafting, industrial processing or development of new cultivars. To ensure proper identification of cultivars and understand the genetic structure of wild and cultivated material, we genotyped 362 cultivated and wild individuals of walnut trees from southwest China (with two additional populations from Xinjiang, plus three cultivars from Canada, France and Belgium) using 36 polymorphic microsatellite loci. We found relatively low indices of genetic diversity (H _O = 0.570, H _E = 0.404, N _A = 2.345) as well as a high level of clonality (>85% of cultivars), indicating reliance on genetically narrow sources of parental material for breeding. Our STRUCTURE and PCoA analyses generally delineated the two species, though considerable levels of introgression were also evident. More significantly, we detected a distinct genetic group of cultivated Juglans sigillata, which mainly comprised individuals of the popular 'Yangbidapao' landrace. Finally, a core set of 18 SSR loci was selected, which was capable of identifying 32 cultivars. In a nutshell, our results call for more utilization of genetically disparate material, including wild walnut trees, as parental sources to breed for more cultivars. The data reported herein will significantly contribute towards the genetic improvement and conservation of the walnut germplasm in southwest China.

Cryptic genetic divergence of the red dwarf rasbora, Microrasbora rubescens, in and around Inle Lake: implications for the origin of endemicity in the ancient lake in Myanmar

Inle Lake, an ancient lake located in the Shan Plateau of Myanmar, is a biogeographically attractive region with high fish endemism. Some endemic species inhabit the lake as well as the surrounding areas. The genetic and ecological relationships between populations in the lake and surrounding areas provide important insights into the process underlying ichthyofaunal formation in Inle Lake. In this study, the authors focused on red dwarf rasbora Microrasbora rubescens, an endemic genus and species in this region, and estimated its population structure and evolutionary scenario based on genome-wide polymorphism, mtDNA and geometric morphometric analyses using samples from Inle Lake and three areas surrounding the lake. The results showed that M. rubescens comprises at least three genetically divergent lineages (Inle, Heho and Hopong) with distinct geographic structures consistent with nuclear and mtDNA data. In contrast, there was no clear regional differentiation in morphology. The divergence time estimation based on mtDNA suggests that the Hopong lineage diverged at 2.7 Ma and the Inle and Heho lineages diverged at 1.9 Ma - consistent with the nuclear DNA results. The deep divergence observed in the endemic species supports the ancient history of ichthyofaunal development in this region. The distinct regional differentiation and morphological conservatism of this species might have been shaped by niche conservatism in stagnant water environments that limit dispersal and morphological diversification. Future comprehensive genetic and morphological analyses and comparisons for other native species should reveal the geographic and ecological processes that shaped the ichthyofauna in this region.

Pre-domestication bottlenecks of the cultivated seaweed Gracilaria chilensis

Gracilaria chilensis is the main cultivated seaweed in Chile. The low genetic diversity observed in the Chilean populations has been associated with the over-exploitation of natural beds and/or the founder effect that occurred during post-glacial colonization from New Zealand. How these processes have affected its evolutionary trajectory before farming and incipient domestication is poorly understood. In this study, we used 2232 single nucleotide polymorphisms (SNPs) to assess how the species' evolutionary history in New Zealand (its region of origin), the founder effect linked to transoceanic dispersion and colonization of South America, and the recent over-exploitation of natural populations have influenced the genetic architecture of G. chilensis in Chile. The contrasting patterns of genetic diversity and structure observed between the two main islands in New Zealand attest to the important effects of Quaternary glacial cycles on G. chilensis. Approximate Bayesian Computation (ABC) analyses indicated that Chatham Island and South America were colonized independently near the end of the Last Glacial Maximum and emphasized the importance of coastal and oceanic currents during that period. Furthermore, ABC analyses inferred the existence of a recent and strong genetic bottleneck in Chile, matching the period of over-exploitation of the natural beds during the 1970s, followed by rapid demographic expansion linked to active clonal propagation used in farming. Recurrent genetic bottlenecks strongly eroded the genetic diversity of G. chilensis prior to its cultivation, raising important challenges for the management of genetic resources in this incipiently domesticated species.

Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.