Restriction Site-associated DNA Sequencing Research Articles

Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history.

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Neutral markers reveal complex population structure across the range of a widespread songbird.

Understanding how both contemporary and historical physical barriers influence gene flow is key to reconstructing evolutionary histories and can allow us to predict species' resilience to changing environmental conditions. During the last glacial maximum (LGM), many high latitude North American bird species were forced into glacial refugia, including mountain bluebirds (Silia currucoides). Within their current breeding range, mountain bluebirds still experience a wide variety of environmental conditions and barriers that may disrupt gene flow and isolate populations. Using single nucleotide polymorphisms (SNPs) obtained through restriction site-associated DNA sequencing, we detected at least four genetically distinct mountain bluebird populations. Based on this structure, we determined that isolation-by-distance, the northern Rocky Mountains, and discontinuous habitat are responsible for the low connectivity and the overall history of each population going back to the last glacial maximum. Finally, we identified five candidate genes under balancing selection and three loci under diversifying selection. This study provides the first look at connectivity and gene flow across the range of these high-altitude and high latitude songbirds.

Striking genetic homogeneity in the widespread South American bracken.

Bracken (Pteridium, Dennstaedtiaceae) is a cosmopolitan genus of aggressive disturbance colonizers that are toxic to agricultural livestock. The taxonomy of Pteridium has been treated in multiple schemes, ranging from one to six species worldwide, with numerous subspecies and varieties. Recent work has focused on the worldwide distribution and systematics of the bracken fern, but South America has been poorly represented. We present the first continent-wide sampling and analysis of Pteridium esculentum, a Southern Hemisphere diploid species. Within South America, P. esculentum has several morphotypes, distinguished into subspecies by variation in indument and lamina architecture. We used double digest restriction site-associated DNA sequencing (ddRADSeq) to assess the phylogenetic relationships of P. esculentum subspecies. We found a striking genetic homogeneity in the species, being able to support only two morphotypes from molecular data: P. e. arachnoideum and P. e. campestre. We had high confidence for shallow and deep phylogenetic relationships, but less support for relationships among crown groups. We describe an east-west geographic pattern that would explain the relationships between populations; and, in contrast to previous studies, we detected differences with P. esculentum from Australia. These results will lay the foundations for studying variations in this species' behavior as a weed, as well as its impact on the production of agricultural livestock in South America.

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers.

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

Analysis of genetic diversity by the SLAF-seq among the farmed Onychostoma macrolepis populations

ObjectiveThe objective of this study was to examine the genetic diversity within and between farmed populations of Onychostoma macrolepis, and to establish a foundation for enhancing the genetic resources of breeding groups through the introduction of new individuals and crossbreeding. A total of 49 individuals were subjected to sequencing using Specific-Locus Amplified Fragment Sequencing (SLAF-seq), one of the restriction site-associated DNA sequencing technologies. The single nucleotide polymorphisms(SNPs)were identified to conduct the analyzation of phylogeny population structure, principal component and genetic diversity.ResultsA total of 853,067 SNPs were identified. The results of the phylogenetic analysis revealed that each sample was genetically clustered into three distinct groups: ZhenPing (ZP), LanGao parents (LG), and their progeny population (LG-F1). Each population was observed to be clustered together. Analysis of population genetic diversity revealed that the observed heterozygosity (Ho) ranged from 0.200 to 0.230, the expected heterozygosity (He) ranged from 0.280 to 0.282, and the polymorphic information content (PIC) ranged from 0.228 to 0.230. These results indicate that the genetic diversity of the population is low and the signs of long-term interbreeding are obvious, but there are differences between the populations, and the genetic diversity of the population can be improved by hybridization in different regions.

Divergent ecological selection maintains species boundaries despite gene flow in a rare endemic tree, Quercus acerifolia (maple-leaf oak).

Strong gene flow from outcrossing relatives tends to blur species boundaries, while divergent ecological selection can counteract gene flow. To better understand how these two forces affect the maintenance of species boundaries, we focused on a species complex including a rare species, maple-leaf oak (Quercus acerifolia), which is found in only four disjunct ridges in Arkansas. Its limited range and geographic proximity to co-occurring close relatives create the possibility for genetic swamping. In this study, we gathered genome-wide SNPs using restriction-site associated DNA sequencing (RADseq) from 190 samples of Q. acerifolia and three of its close relatives, Q. shumardii, Q. buckleyi, and Q. rubra. We found that Q. shumardii and Q. acerifolia are reciprocally monophyletic with low support, suggesting incomplete lineage sorting, introgression between Q. shumardii and Q. acerifolia, or both. Analyses that model allele distributions demonstrate that admixture contributes strongly to this pattern. Populations of Q. acerifolia experience gene flow from Q. shumardii and Q. rubra, but we found evidence that divergent selection is likely maintaining species boundaries: 1) ex situ collections of Q. acerifolia have a higher proportion of hybrids compared to the mature trees of the wild populations, suggesting ecological selection against hybrids at the seed/seedling stage; 2) ecological traits co-vary with genomic composition; and 3) Q. acerifolia shows genetic differentiation at loci hypothesized to influence tolerance of radiation, drought, and high temperature. Our findings strongly suggest that in maple-leaf oak, selection results in higher divergence at regions of the genome despite gene flow from close relatives.

Genetic structure of a narrowly distributed species Chimonobambusa tumidissinoda in the Yunnan-Guizhou Plateau, and its implications for conservation

Chimonobambusa tumidissinoda is a narrowly distributed and endangered bamboo species, which is endemic to the Yunnan-Guizhou Plateau and possesses significant ecological and economic value, particularly due to its delectable shoots. In this study, we employed the restriction site-associated DNA sequencing (RAD-seq) method to elucidate the population structure, genetic diversity, and demographic history of C. tumidissinoda across its entire distribution range. Unlike the genetic pattern revealed previously, three distinct genetic groups within the C. tumidissinoda samples were consistently discriminated based on multiple analyses, and low levels of genetic diversity (HO = 0.085, HE = 0.115, π = 0.116) were found. Furthermore, moderate genetic differentiation (FST = 0.070 - 0.095; PhiPT = 0.129 - 0.156) was detected among these groups with Mantel tests suggesting that isolation by distance (IBD) could potentially drive population genetic differentiation in this species. Moreover, demographic analyses incicated that C. tumidissinoda likely experienced a bottleneck event since the Pliocene, which may be attributed to drier and colder winters in Yunnan resulting from intensified Asian monsoon activity. Based on these findings, we suggest that the endangered status of C. tumidissinoda has been underestimated. Conservation measures should be implemented both in situ and ex situ to enhance its adaptability towards future climate fluctuations. Additionally, seed collection for germplasm storage is recommended.

Conservation genomics assessment of Tharp's bluestar (Amsonia tharpii) with comparisons to widespread (A. longilora) and narrowly endemic (A. fugatei) congeners.

Land-use change and habitat fragmentation are threats to biodiversity. The decrease in available habitat, increase in isolation, and mating within populations can lead to elevated inbreeding, lower genetic diversity, and poor fitness. Here we investigate the genetics of two rare and threatened plant species, Amsonia tharpii and A. fugatei, and we compare them to a widespread congener A. longiflora. We also report the first phylogenetic study of the genus Amsonia (Apocynaceae), including 10 of the 17 taxa and multiple sampling locations, to understand species relationships. We used a double digest restriction-site associated DNA sequencing (ddRADseq) approach to investigate the genetic diversity and gene flow of each species and to create a maximum likelihood phylogeny. The ddRADseq data was mapped to a reference genome to separate out the chloroplast and nuclear markers for population genetic analysis. Our results show that genetic diversity and inbreeding were low across all three species. The chloroplast and nuclear dataset in A. tharpii were highly structured, whereas they showed no structure for A. fugatei, while A. longiflora lacked structure for nuclear data but not chloroplast. Phylogenetic results revealed that A. tharpii is distinct and sister to A. fugatei, and together they are distantly related to A. longiflora. Our results demonstrated that evolutionary history and contemporary ecological processes largely influences genetic diversity within Amsonia. Interestingly, we show that in A. tharpii there was significant structure despite being pollinated by large, bodied hawkmoths that are known to be able to carry pollen long distances, suggesting that other factors are contributing to the structure observed among A. tharpii populations. Conservation efforts should focus on protecting all of the A. tharpii populations, as they contain unique genetic diversity, and a protection plan for A. fugatei needs to be established due to its limited distribution.

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense

BackgroundRelict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun.ResultsA total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM).ConclusionOur study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

Whole-genome selection signature differences between Chaohu and Ji'an red ducks.

Assessing the genetic structure of local varieties and understanding their genetic data are crucial for effective management and preservation. However, the genetic differences among local breeds require further explanation. To enhance our understanding of their population structure and genetic diversity, we conducted a genome-wide comparative study of Chaohu and Ji'an Red ducks using genome sequence and restriction site-associated DNA sequencing technology. Our analysis revealed a distinct genetic distinction between the two breeds, leading to divided groups. The phylogenetic tree for Chaohu duck displayed two branches, potentially indicating minimal impact from artificial selection. Additionally, our ROH (runs of homozygosity) analysis revealed that Chaohu ducks had a lower average inbreeding coefficient than Ji'an Red ducks. We identified several genomic regions with high genetic similarity in these indigenous duck breeds. By conducting a selective sweep analysis, we identified 574 candidate genes associated with muscle growth (BMP2, ITGA8, MYLK, and PTCH1), fat deposits (ELOVL1 and HACD2), and pigmentation (ASIP and LOC101797494). These results offer valuable insights for the further enhancement and conservation of Chinese indigenous duck breeds.

Low-coverage sequencing and Wahlund effect severely bias estimates of inbreeding, heterozygosity and effective population size in North American wolves.

vonHoldt etal. ((2024), Molecular Ecology, 33, e17231) (vH24) used low-coverage (average ~ 7X read depth) restriction site-associated DNA sequence data to estimate individual inbreeding and heterozygosity, and recent effective population size (Ne), in Great Lakes (GL) and Northern Rocky Mountain (RM) wolves. They concluded that RM heterozygosity rapidly declined between 1991 and 2020, and that Ne declined substantially in GL and RM over the last 50 generations. Here, we evaluate the effects of low sequence coverage and sampling strategy on vH24's findings and provide general recommendations for using sequence data to evaluate inbreeding, heterozygosity and Ne. Low-coverage sequencing resulted in downwardly biased estimates of individual inbreeding and heterozygosity, and an erroneous large temporal decline in RM heterozygosity due to declining read depth through time. Additionally, vH24's sampling strategy-which combined individuals from several genetically differentiated populations and across at least eight wolf generations-is expected to have resulted in severe downward bias in estimates of recent Ne for RM. We recommend using high-coverage sequence data ( 15-20X) to estimate inbreeding and heterozygosity. Carefully filtering individuals, loci and genotypes, and using genotype imputation or likelihoods can help to minimise bias when low-coverage sequence data must be used. For estimation of contemporary Ne, the marginal benefits of using more than 103-104 loci are small, so aggressive filtering of loci with low average read depth potentially can retain most individuals without sacrificing much precision. Individuals are relatively more valuable than loci because analyses of contemporary Ne should focus on roughly single-generation samples from local breeding populations.

Development of 51 K liquid-phased probe array for Loblolly and Slash pines and its application to GWAS of Slash pine breeding population

Loblolly pine (Pinus taeda) and Slash pine (Pinus elliottii) are important economically and ecologically tree species worldwide. To implement genomic based breeding strategy, genotyping large number of samples is needed. However, huge and complex pines genome make it difficult to obtain abundant and cost-effectively markers. In this study, we designed a SNP array containing 51 K liquid-phased probes based on the published SNPs and probes of the two species, novel Slash pine-specific long-read transcript, and SNPs obtained from a double digest Restriction-Site Associated DNA sequencing (ddRADseq). We performed targeted capture sequencing using the developed SNP array on three pine species of Loblolly, Slash, and Caribbean pine (Pinus caribaea). The probe capture had an efficiency of 68.26% on average. Principal Component Analysis (PCA), phylogenetic analysis, and genetic structure analysis of 28,670 SNPs genotyped for the three pine species demonstrated the transferability of the SNP array for the three pines. Furthermore, we genotyped 221 range-wide Slash pine genotypes, analyzed their genetic diversity and performed a genome-wide association study (GWAS). We identified 95 SNPs significantly associated with growth and wood quality traits of Slash pine. We have identified potential roles of the MYB2 and MYB46 transcription factors in regulating cell cycle and differentiation in Slash pine, while also suggesting that AP2/ERF and Glycosyl hydrolase may play significant roles in wood formation processes. This new SNP array can be widely applied in future molecular breeding research of conifers.

Genomic evaluation of mud sunfish in the atlantic slope drainage

Mud Sunfish (Acantharchus pomotis, Baird), a cryptic species native to the Atlantic and Gulf slope drainages of North America, have experienced a documented reduction in their geographic range. Atlantic Slope populations of Mud Sunfish were evaluated to determine the population structure, assess genetic diversity, and to identify presence of evolutionary significant units (ESU) in order to develop fundamental understandings of its occurrence. Tissue samples were collected from across the Atlantic Slope for genomic comparison. Restriction site Associated DNA Sequencing (RAD-Seq) was performed on 95 individuals from 14 populations. Mantel test provided evidence of a correlation between genetic and geographical distances, suggesting a pattern of isolation-by-distance. Differences among populations of Mud Sunfish largely occur along a north-south latitudinal gradient. Bayesian clustering of putatively outlier SNP loci, suggests that regional adaptive pressure is occurring. Despite a large geographic distance and significant differences between the majority of the populations of Mud Sunfish, Bayesian clustering analysis suggests that two unresolved clades are persisting within the Atlantic Slope Drainage: a northern clade of Mud Sunfish ranging from New Jersey to North Carolina and a southern clade ranging from South Carolina to Florida. A population of Mud Sunfish in Florida appears distinct from all others in the Atlantic Coastal plain and likely represents a distinct management unit.

Genome-wide association study reveals candidate genes critical for skin pigmentation in common carp (Cyprinus carpio) strains including koi

The differentiation and variation of body color in fish is an important component of biological evolution and phylogenetic research, and is influenced by several factors such as heredity, nutrition, physiological hormones and rearing conditions, with heredity being the most critical factor. An improved comprehension of the genetic regulation in different carp populations will facilitate future applications of precision breeding and pigment cell development-related genes editing. In this study, we used restriction site-associated DNA sequencing (2d-RAD) to determine the distribution characteristics of single nucleotide polymorphisms (SNPs) in eight common carp strains including six koi strains, i.e., Kohaku (Hb), Taisho (Dz), Showa (Zh), Golden (Hj), Matsukawabake (Sch) and Five-color (Ws), and two genetically improved breeding strains, i.e., FFRC No. 2 strain (Fr) and red FFRC strain (Hfr), and then SNP markers to resolve genetic information on the population structure and skin color formation of these strains. We found that most individuals of eight strains were clearly divided into eight clades. Genetic diversity within Fr, Sch and Ws is relatively high, while it is low within Dz, Hb and Hfr. A comparison between the populations revealed that Hfr and Ws have the highest degree of genetic differences, while Hb and Dz have the lowest. The GWAS analysis for two strategies identified several SNPs that are significantly related to skin color variation and affect melanin synthesis (ephrinb2, slc22a23, herc2, ephb2), carotenoid metabolism (oca2, bcmo1, bco2) and circadian rhythm(cipc). To summarize, our study provides some information on the genetic structure and diversity of Japanese koi and Chinese FFRC carp. Furthermore, based on the data obtained on relevant genes related to skin color differentiation and pigmentation, it can provide a theoretical basis for the subsequent integration of dominant traits and precision breeding in carps.

Genetic diversity and population structure of four Anatolian sheep revealed by genome-wide ddRADseq data

This is the first wide-scale study to reveal next-generation sequencing-based and comprehensive study to reveal genetic diversity, population structure, and phylogenetic relationships of some of Türkiye's indigenous fat-tailed sheep breeds. A total of 80 animals belonging to Akkaraman (AKR), Güney Karaman (GKR), Morkaraman (MKR), and Karakaş (KRK) sheep breeds were genotyped with genome-wide of 296.097 single nucleotide polymorphisms (SNPs) data recovered from double digest restriction site-associated DNA sequencing (ddRADseq) libraries processed by Illumina HiSeq X Ten instrument. Minor allele frequency (MAF) ranged from 0.311 (MKR) to 0.316 (AKR and GKR). Observed heterozygosity (HO) values were higher than expected values (HE) across all populations with a mean of 0.307. Negative inbreeding coefficient (FIS) values were observed in four sheep populations varying between −0.034 (AKR) to −0.060 (MKR). Cross-entropy criterion-based ADMIXTURE analysis demonstrated that four Anatolian sheep breeds were represented by two ancestral populations in which GKR and AKR were distinct, while KRK and MKR were mixed populations. Genetic distance-based tree, discriminant analysis of principal components (DAPC), and ADMIXTURE analyses confirmed that GKR, which is believed to be derived from the AKR breed, has become genetically different from other sheep populations. Moreover, no migration from and into the GKR breed was detected, while a migration edge was drawn from the KRK and MKR clade into the AKR breed by the TreeMix algorithm. This study confirmed that high-density SNP data obtained from advanced next-generation sequencing (NGS) platforms are efficient in genetically distinguishing animal breeds, even varieties. It is recommended that further NGS-based studies could be carried out to better characterize different indigenous livestock populations raised in Türkiye.

Genetic divergence and adaptation of an isolated European lobster population in the Netherlands

Abstract Identifying isolated populations is a key step towards enacting effective conservation management. European lobsters (Homarus gammarus) from Oosterschelde in the Netherlands are subject to fishery pressure and have previously been reported as genetically differentiated. They are also putatively of transplanted origin and have subsequently endured recent bottlenecking and environmental change. We assessed Oosterschelde lobsters to evaluate their demographic independence and appraise potential founder effects and evolutionary responses to isolation. Using restriction-site associated DNA sequencing, we genotyped 6185 single nucleotide polymorphisms across 188 individuals from 27 sites across the Atlantic range of H. gammarus to investigate population genetic diversity, structure, and potential adaptation. Our results show that Oosterschelde lobsters are genetically divergent from other stocks. We evidence extensive differentiation via both neutral and outlier loci, indicative of strong biophysical and demographic isolation, and detect signatures of reduced genetic diversity that may reflect weak founder effects or subsequent population contractions. Among outlier loci, we identify candidates for range-wide local adaptation via variants in genes of important biological functionality and link a private allele of Oosterschelde to a locus potentially conveying adaptive tolerance to environmental hypoxia. Given our findings, we advise proactive monitoring of Oosterschelde lobsters to explore whether existing management measures effectively conserve this discrete, self-recruiting population.

A Diverging Species within the Stewartia gemmata (Theaceae) Complex Revealed by RAD-Seq Data.

Informed species delimitation is crucial in diverse biological fields; however, it can be problematic for species complexes. Showing a peripatric distribution pattern, Stewartia gemmata and S. acutisepala (the S. gemmata complex) provide us with an opportunity to study species boundaries among taxa undergoing nascent speciation. Here, we generated genomic data from representative individuals across the natural distribution ranges of the S. gemmata complex using restriction site-associated DNA sequencing (RAD-seq). Based on the DNA sequence of assembled loci containing 41,436 single-nucleotide polymorphisms (SNPs) and invariant sites, the phylogenetic analysis suggested strong monophyly of both the S. gemmata complex and S. acutisepala, and the latter was nested within the former. Among S. gemmata individuals, the one sampled from Mt. Tianmu (Zhejiang) showed the closest evolutionary affinity with S. acutisepala (which is endemic to southern Zhejiang). Estimated from 2996 high-quality SNPs, the genetic divergence between S. gemmata and S. acutisepala was relatively low (an Fst of 0.073 on a per-site basis). Nevertheless, we observed a proportion of genomic regions showing relatively high genetic differentiation on a windowed basis. Up to 1037 genomic bins showed an Fst value greater than 0.25, accounting for 8.31% of the total. After SNPs subject to linkage disequilibrium were pruned, the principal component analysis (PCA) showed that S. acutisepala diverged from S. gemmata along the first and the second PCs to some extent. By applying phylogenomic analysis, the present study determines that S. acutisepala is a variety of S. gemmata and is diverging from S. gemmata, providing empirical insights into the nascent speciation within a species complex.

Intergeneric hybrid origin of the invasive tetraploid Cirsium vulgare.

The invasive tetraploid Cirsium vulgare hybridizes with both Cirsium and Lophiolepis. Its conflicted position in molecular phylogenies, and its peculiar combination of morphological, anatomical, and genomic features that are alternatively shared with representatives of Cirsium or Lophiolepis, strongly suggest its intergeneric hybrid origin. Genetic relationships of C. vulgare (8 samples) with genus Lophiolepis (11 species) and other representatives of genus Cirsium (12 species) were evaluated using restriction site-associated DNA sequencing (RADseq) and examined using analytical and imaging approaches, such as NeighborNet, Heatmap, and STRUCTURE, to identify nuclear genomes admixture. Estimation of the intensity of spontaneous hybridization within and between Cirsium and Lophiolepis was based on herbarium revisions and published data for all reported hybrids pertinent to taxa currently included in Cirsium or Lophiolepis. The genome of any examined Cirsium species is more similar to C. vulgare than to any Lophiolepis species, and vice versa. The nuclear genome of the tetraploid C. vulgare is composed of two equivalent parts, each attributable either to Lophiolepis or to Cirsium; the organellar RADseq data clustered C. vulgare with the genus Cirsium. Spontaneous hybridization between Cirsium and Lophiolepis is significantly less intensive than within these genera. Our analyses provide compelling evidence that the invasive species C. vulgare has an allotetraploid intergeneric origin, with the maternal parent from Cirsium and the paternal from Lophiolepis. For the purpose of delimiting monophyletic genera, we propose keeping Lophiolepis separate from Cirsium and segregating C. vulgare into the hybridogenous genus Ascalea.

Different reference genomes determine different results: Comparing SNP calling in RAD-seq of Engelhardia roxburghiana using different reference genomes

Advances in next-generation sequencing (NGS) have significantly reduced the cost and improved the efficiency of obtaining single nucleotide polymorphism (SNP) markers, particularly through restriction site-associated DNA sequencing (RAD-seq). Meanwhile, the progression in whole genome sequencing has led to the utilization of an increasing number of reference genomes in SNP calling processes. This study utilized RAD-seq data from 242 individuals of Engelhardia roxburghiana, a tropical tree of the walnut family (Juglandaceae), with SNP calling conducted using the STACKS pipeline. We aimed to compare both reference-based approaches, namely, employing a closely related species as the reference genome versus the species itself as the reference genome, to evaluate their respective merits and limitations. Our findings indicate a substantial discrepancy in the number of obtained SNPs between using a closely related species as opposed to the species itself as reference genomes, the former yielded approximately an order of magnitude fewer SNPs compared to the latter. While the missing rate of individuals and sites of the final SNPs obtained in the two scenarios showed no significant difference. The results showed that using the reference genome of the species itself tends to be prioritized in RAD-seq studies. However, if this is unavailable, considering closely related genomes is feasible due to their wide applicability and low missing rate as alternatives. This study contributes to enrich the understanding of the impact of SNP acquisition when utilizing different reference genomes.

Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data.

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.