Population Genomic Studies Research Articles

Overview of the Saccharomyces cerevisiae population structure through the lens of 3,034 genomes.

With the rise of high-throughput sequencing technologies, a holistic view of genetic variation within populations-through population genomics studies-appears feasible, although it remains an ongoing effort. Genetic variation arises from a diverse range of evolutionary forces, with mutation and recombination being key drivers in shaping genomes. Studying genetic variation within a population represents a crucial first step in understanding the relationship between genotype and phenotype and the evolutionary history of species. In this context, the budding yeast Saccharomyces cerevisiae has been at the forefront of population genomic studies. In addition, it has a complex history that involves adaptation to a wide range of wild and human-related ecological niches. Although to date more than 3,000 diverse isolates have been sequenced, there is currently a lack of a resource bringing together sequencing data and associated metadata for all sequenced isolates. To perform a comprehensive analysis of the population structure of S. cerevisiae, we collected genome sequencing data from 3,034 natural isolates and processed the data uniformly. We determined ploidy levels, identified single nucleotide polymorphisms (SNPs), small insertion-deletions (InDels), copy number variations (CNVs), and aneuploidies across the population, creating a publicly accessible resource for the yeast research community. Interestingly, we showed that this population captures ∼93% of the species diversity. Using neighbor-joining and Bayesian methods, we redefined the populations, revealing clustering patterns primarily based on ecological origin. This work represents a valuable resource for the community and efforts have been made to make it evolvable and integrable to future yeast population studies.

Transposable element accumulation drives genome size increase in Hylesia metabus (Lepidoptera: Saturniidae), an urticating moth species from South America

Abstract We present the first nuclear genome assembly and a complete mitogenome for Hylesia metabus (Arthropoda; Insecta; Lepidoptera; Saturniidae). The assembled nuclear genome sequence is 1,271 Mb long, which is among the 10 largest lepidopteran genome assemblies published to date. It is scaffolded in 31 pseudo chromosomes, has a BUSCO score of 99.5%, and has a highly conserved synteny compared to phylogenetically close species. Repetitive elements make up 67% of the nuclear genome and are mainly located in intergenic regions, among which LINEs were predominant, with CR1-Zenon being the most abundant. Phylogenetic and comparative analyses of H. metabus assembly and 17 additional Saturniidae and Sphingidae assemblies suggested that an accumulation of repetitive elements likely led to the increased size of H. metabus’ genome. Gene annotation using Helixer identified 26,122 transcripts. The Z scaffold was identified using both a synteny analysis and variations of coverage for two resequenced male and female H. metabus. The H. metabus nuclear genome and mitogenome assemblies can be found and browsed on the BIPAA website and constitute useful resources for future population and comparative genomics studies.

The molecular dialogue between Zymoseptoria tritici and wheat.

Zymoseptoria tritici is a highly damaging pathogen that causes high wheat yield losses in temperate climates. Z. tritici emerged during the domestication of wheat in the Fertile Crescent and has been extensively used as a model system for population genetic and genomic studies. New genetic tools and resources have provided a better understanding of the molecular components involved in the wheat-Z. tritici interaction, highlighted by the cloning of three wheat resistance genes and four Z. tritici avirulence genes. Despite the considerable progress made in the last few years, the mechanisms that mediate Z. tritici colonization remain largely unknown. In this review, we summarize the latest advances in understanding the molecular components mediating wheat-Z. tritici interactions and we discuss future research lines to close current knowledge gaps.

Improved genome assembly of Candida auris strain B8441 and annotation of B11205.

Candida auris is a fungal pathogen of significant worldwide concern, typically resistant to one or more antifungal drugs. We report a completed genome for clade Ia isolate B8441 and gene annotations of clade Ic isolate B11205. These resources will support public health investigations and population genomic studies of this pathogen.

Population connectivity and size reductions in the Anthropocene: the consequence of landscapes and historical bottlenecks in white forsythia fragmented habitats

BackgroundWhite forsythia (Abeliophyllum distichum) is an endangered Korean Peninsula endemic that has been subjected to recent population genomics studies using SNPs via RAD sequencing. Here, we primarily employed the often underutilized haplotype information from RAD loci to further describe the species’ previously uninvestigated haplotype-based genomic variation and structure, and genetic-geographic characteristics and gene flow patterns among its five earlier identified genetic groups. We also inferred the time of past events that may have impacted the effective population size of these groups, as well as the species’ potential future distribution amidst the warming climate and anthropogenic threats.ResultsOur findings emphasized the recognition of the species’ regional patterns of genetic structure, and the role of topography and its associated gene flow patterns as some of the possible factors that may have influenced the species’ present-day fragmented population distribution. The inferred bottleneck events during the Anthropocene, some of which aligned with the time of historical catastrophic events on the Peninsula (e.g., the Korean War), were revealed to have contributed to the generally low effective population size of its five lineages, particularly those with marginal distributional range. Future distribution under both optimistic and pessimistic climatic scenarios suggests unlikely suitable habitats for these populations to expand from their current range limits, at least in the next 80 years.ConclusionsThe small effective population size and landscape-driven limited gene flow among white forsythia populations will remain a big challenge for the conservation management of the species’ already fragmented population distribution. To help mitigate these impacts, the merging of various research approaches and the use of genomic data to their full potential is recommended to provide the optimized knowledge-based tools for the conservation of this endangered species, and other similar plants under pressure.

Screening Familial Risk for Hereditary Breast and Ovarian Cancer.

Most patients with pathogenic or likely pathogenic (P/LP) variants for breast cancer have not undergone genetic testing. To identify patients meeting family history criteria for genetic testing in the electronic health record (EHR). This study included both cross-sectional (observation date, February 1, 2024) and retrospective cohort (observation period, January 1, 2018, to February 1, 2024) analyses. Participants included patients aged 18 to 79 years enrolled in Renown Health, a large health system in Northern Nevada. Genotype was known for 38 003 patients enrolled in Healthy Nevada Project (HNP), a population genomics study. An EHR indicating that a patient is positive for criteria according to the Seven-Question Family History Questionnaire (hereafter, FHS7 positive) assessing familial risk for hereditary breast and ovarian cancer (HBOC). The primary outcomes were the presence of P/LP variants in the ATM, BRCA1, BRCA2, CHEK2, or PALB2 genes (cross-sectional analysis) or a diagnosis of cancer (cohort analysis). Age-adjusted cancer incidence rates per 100 000 patients per year were calculated using the 2020 US population as the standard. Hazard ratios (HRs) for cancer attributable to FHS7-positive status were estimated using cause-specific hazard models. Among 835 727 patients, 423 393 (50.7%) were female and 29 913 (3.6%) were FHS7 positive. Among those who were FHS7 positive, 24 535 (82.0%) had no evidence of prior genetic testing for HBOC in their EHR. Being FHS7 positive was associated with increased prevalence of P/LP variants in BRCA1/BRCA2 (odds ratio [OR], 3.34; 95% CI, 2.48-4.47), CHEK2 (OR, 1.62; 95% CI, 1.05-2.43), and PALB2 (OR, 2.84; 95% CI, 1.23-6.16) among HNP female individuals, and in BRCA1/BRCA2 (OR, 3.35; 95% CI, 1.93-5.56) among HNP male individuals. Being FHS7 positive was also associated with significantly increased risk of cancer among 131 622 non-HNP female individuals (HR, 1.44; 95% CI, 1.22-1.70) but not among 114 982 non-HNP male individuals (HR, 1.11; 95% CI, 0.87-1.42). Among 1527 HNP survey respondents, 352 of 383 EHR-FHS7 positive patients (91.9%) were survey-FHS7 positive, but only 352 of 883 survey-FHS7 positive patients (39.9%) were EHR-FHS7 positive. Of the 29 913 FHS7-positive patients, 19 764 (66.1%) were identified only after parsing free-text family history comments. Socioeconomic differences were also observed between EHR-FHS7-negative and EHR-FHS7-positive patients, suggesting disparities in recording family history. In this cross-sectional study, EHR-derived FHS7 identified thousands of patients with familial risk for breast cancer, indicating a substantial gap in genetic testing. However, limitations in EHR family history data suggested that other identification methods, such as direct-to-patient questionnaires, are required to fully address this gap.

The repertoire of G-protein-coupled receptor variations in the Japanese population 54KJPN.

G-protein-coupled receptors (GPCRs) are the largest superfamily in the human genome and the major targets for the market drugs. Recent massive genomics studies revealed numerous natural variations in the general population. 54KJPN is the most extensive Japanese population genomics study, curating the whole genome sequences from about 54,000 individuals. Here, by analyzing 390 non-olfactory GPCR genes in the 54KJPN dataset, we annotated 25,443 missense single-nucleotide variations. Among them, we found 120 major variations that appear with an allele frequency greater than 0.5, including variations that occurred on posttranslational modification sites. Structural alignment of GPCRs using the generic numbering system in the GPCRdb reveals enrichment of alterations in the conserved arginine residue within the DRY motif, which contributes to downstream G-protein signaling. A comparison with the worldwide 1000 Genomes Project (1KGP) dataset found 23 variations that were present exclusively in the 54KJPN dataset. This study will be the basis for future pharmacogenomics studies for the Japanese population.

Population genomics of Plasmodium ovale species in sub-Saharan Africa.

Plasmodium ovale curtisi (Poc) and Plasmodium ovale wallikeri (Pow) are relapsing malaria parasites endemic to Africa and Asia that were previously thought to represent a single species. Amid increasing detection of ovale malaria in sub-Saharan Africa, we performed a population genomic study of both species across the continent. We conducted whole-genome sequencing of 25 isolates from Central and East Africa and analyzed them alongside 20 previously published African genomes. Isolates were predominantly monoclonal (43/45), with their genetic similarity aligning with geography. Pow showed lower average nucleotide diversity (1.8×10-4) across the genome compared to Poc (3.0×10-4) (p < 0.0001). Signatures of selective sweeps involving the dihydrofolate reductase gene were found in both species, as were signs of balancing selection at the merozoite surface protein 1 gene. Differences in the nucleotide diversity of Poc and Pow may reflect unique demographic history, even as similar selective forces facilitate their resilience to malaria control interventions.

Evolution of the correlated genomic variation landscape across a divergence continuum in the genus Castanopsis.

The heterogeneous landscape of genomic variation has been well documented in population genomic studies. However, disentangling the intricate interplay of evolutionary forces influencing the genetic variation landscape over time remains challenging. In this study, we assembled a chromosome-level genome for Castanopsis eyrei and sequenced the whole genomes of 276 individuals from 12 Castanopsis species, spanning a broad divergence continuum. We found highly correlated genomic variation landscapes across these species. Furthermore, variations in genetic diversity and differentiation along the genome were strongly associated with recombination rates and gene density. These results suggest that long-term linked selection and conserved genomic features have contributed to the formation of a common genomic variation landscape. By examining how correlations between population summary statistics change throughout the species divergence continuum, we determined that background selection alone does not fully explain the observed patterns of genomic variation; the effects of recurrent selective sweeps must be considered. We further revealed that extensive gene flow has significantly influenced patterns of genomic variation in Castanopsis species. The estimated admixture proportion correlated positively with recombination rate and negatively with gene density, supporting a scenario of selection against gene flow. Additionally, putative introgression regions exhibited strong signals of positive selection, an enrichment of functional genes, and reduced genetic burdens, indicating that adaptive introgression has played a role in shaping the genomes of hybridizing species. This study provides insights into how different evolutionary forces have interacted in driving the evolution of the genomic variation landscape.

Population genomics and demographic modeling enhance our understanding of trophic level interactions in biological control

In this paper, we focus on the application of population genomics, including next-generation sequencing and demographic modeling, to enhance the understanding of historical and contemporary trophic level interactions that are a fundamental basis of biological control. We outline how the integration of ecological approaches with population genomics and demographic modeling can provide additional insights into existing trophic interactions and how natural enemies, their parasitoids, and pest species evolve through these interactions, potentially affecting future levels of biological control. We briefly review and summarize our recent population genomic studies of predatory ladybird beetles and discuss how our results provide guidelines for the use of population genomics techniques and population demographic modeling to enhance our understanding of trophic level interactions. Our intent is to demonstrate as a proof of concept, how these methods can supplement currently used techniques and add a complementary genomic aspect to ecological understandings of trophic level interactions, including estimates of effective population sizes of natural enemies and pest species, examination of correlations between these effective population sizes, and evolutionary changes in natural enemies and pest species which may influence levels of biological control. This in combination with experiments and simulations provide a well-rounded framework for using population genomics to inform multi-tier trophic level interactions.

Cover Image:

Grapes, first written about in The Book of Songs, have been cultivated for millennia and have become globally valued as a fresh fruit and for making wine. The cover depicts grapevines growing along the coast in Shenzhen, China, where they produce a bountiful harvest brought by the gentle breeze. However, these grapevines are not delicate—they can withstand the sea wind and the splashing of seawater, symbolizing the tenacious vitality of grapes that are adapted to the coastal environment. The population genomics study by Zhang et al. (pages 1408–1426) offers valuable insights for breeding salt‐tolerant grape rootstocks and identifies key structural variations, advancing horticulture of this economically and culturally significant crop.

PfHDAC1 is an essential regulator of P. falciparum asexual proliferation and host cell invasion genes with a dynamic genomic occupancy responsive to artemisinin stress.

Malaria is a major public health problem, with the parasite Plasmodium falciparum causing most of the malaria-associated mortality. It is spread by the bite of infected mosquitoes and results in symptoms such as cyclic fever, chills, and headache. However, if left untreated, it can quickly progress to a more severe and life-threatening form. The World Health Organization currently recommends the use of artemisinin combination therapy, and it has worked as a gold standard for many years. Unfortunately, certain countries in southeast Asia and Africa, burdened with a high prevalence of malaria, have reported cases of drug-resistant infections. One of the major problems in controlling malaria is the emergence of artemisinin resistance. Population genomic studies have identified mutations in the Kelch13 gene as a molecular marker for artemisinin resistance. However, several reports thereafter indicated that Kelch13 is not the main mediator but rather hinted at transcriptional deregulation as a major determinant of drug resistance. Earlier, we identified PfGCN5 as a global regulator of stress-responsive genes, which are known to play a central role in artemisinin resistance generation. In this study, we have identified PfHDAC1, a histone deacetylase as a cell cycle regulator, playing an important role in artemisinin resistance generation. Taken together, our study identified key transcriptional regulators that play an important role in artemisinin resistance generation.

Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

Genomic insights into the critically endangered King Island scrubtit.

Small, fragmented or isolated populations are at risk of population decline due to fitness costs associated with inbreeding and genetic drift. The King Island scrubtit Acanthornis magna greeniana is a critically endangered subspecies of the nominate Tasmanian scrubtit A. m. magna, with an estimated population of < 100 individuals persisting in three patches of swamp forest. The Tasmanian scrubtit is widespread in wet forests on mainland Tasmania. We sequenced the scrubtit genome using PacBio HiFi and undertook a population genomic study of the King Island and Tasmanian scrubtits using a double-digest restriction site-associated DNA (ddRAD) dataset of 5,239 SNP loci. The genome was 1.48 Gb long, comprising 1,518 contigs with an N50 of 7.715Mb. King Island scrubtits formed one of four overall genetic clusters, but separated into three distinct subpopulations when analysed independently of the Tasmanian scrubtit. Pairwise FST values were greater among the King Island scrubtit subpopulations than among most Tasmanian scrubtit subpopulations. Genetic diversity was lower and inbreeding coefficients were higher in the King Island scrubtit than all except one of the Tasmanian scrubtit subpopulations. We observed crown baldness in 8/15 King Island scrubtits, but 0/55 Tasmanian scrubtits. Six loci were significantly associated with baldness, including one within the DOCK11 gene which is linked to early feather development. Contemporary gene flow between King Island scrubtit subpopulations is unlikely, with further field monitoring required to quantify the fitness consequences of its small population size, low genetic diversity and high inbreeding. Evidence-based conservation actions can then be implemented before the taxon goes extinct.

Species-specific dynamics may cause deviations from general biogeographical predictions - evidence from a population genomics study of a New Guinean endemic passerine bird family (Melampittidae).

The family Melampittidae is endemic to New Guinea and consists of two monotypic genera: Melampitta lugubris (Lesser Melampitta) and Megalampitta gigantea (Greater Melampitta). Both Melampitta species have scattered and disconnected distributions across New Guinea in the central mountain range and in some of the outlying ranges. While M. lugubris is common and found in most montane regions of the island, M. gigantaea is elusive and known from only six localities in isolated pockets on New Guinea with very specific habitats of limestone and sinkholes. In this project, we apply museomics to determine the population structure and demographic history of these two species. We re-sequenced the genomes of all seven known M. gigantaea samples housed in museum collections as well as 24 M. lugubris samples from across its distribution. By comparing population structure between the two species, we investigate to what extent habitat dependence, such as in M. gigantaea, may affect population connectivity. Phylogenetic and population genomic analyses, as well as acoustic variation revealed that M. gigantaea consists of a single population in contrast to M. lugubris that shows much stronger population structure across the island. We suggest a recent collapse of M. gigantaea into its fragmented habitats as an explanation to its unexpected low diversity and lack of population structure. The deep genetic divergences between the M. lugubris populations on the Vogelkop region, in the western central range and the eastern central range, respectively, suggests that these three populations should be elevated to full species level. This work sheds new light on the mechanisms that have shaped the intriguing distribution of the two species within this family and is a prime example of the importance of museum collections for genomic studies of poorly known and rare species.

Genetic erosion in a tropical tree species demonstrates the need to conserve wide-ranging germplasm amid extreme habitat fragmentation

Comparing the genetic diversity across different generations within tropical tree populations is an understudied topic. To assess the potential genetic consequences and conservation implications of contemporary disturbances, a population genomic study of Palaquium obovatum across age classes was undertaken. Trees and juveniles were sampled from ten different localities (eight in Singapore, two in Peninsular Malaysia) and subjected to double digest restriction-site associated DNA-sequencing (ddRAD-seq) to assess intergenerational genetic differences and investigate population structure in a hexaploid lineage. Genetic erosion, characterised by reduced heterozygosity, was found to have occurred in almost all wild populations over time, the exceptions being in one isolated coastal population and some areas with cultivated occurences. Population structure was highly localised with the number of genetically distinct populations usually following geographically separated districts, which indicates limitations in pollen and seed dispersal between fragments, possibly due to declines in the associated assemblage of dispersers. For this reason, the germplasm for conserving species diversity in degraded habitats and forest fragments should be selected from a wide range of wild populations across the landscape.

The chromosome-scale reference genome for the pinfish (Lagodon rhomboides) provides insights into their evolutionary and demographic history.

The pinfish (Lagodon rhomboides) is an ecologically, economically, and culturally relevant member of the family Sparidae, playing crucial roles in the marine food webs of the western Atlantic Ocean and Gulf of Mexico. Despite their high abundance and ecological importance, there is a scarcity of genomic resources for this species. We assembled and annotated a chromosome-scale genome for the pinfish, resulting in a highly contiguous 785 Mb assembly of 24 scaffolded chromosomes. The high-quality assembly contains 98.9% complete BUSCOs and shows strong synteny to other chromosome-scale genomes of fish in the family Sparidae, with a limited number of large-scale genomic rearrangements. Leveraging this new genomic resource, we found evidence of significant expansions of dietary gene families over the evolutionary history of the pinfish, which may be associated with an ontogenetic shift from carnivory to herbivory seen in this species. Estimates of historical patterns of population demography using this new reference genome identified several periods of population growth and contraction which were associated with ancient climatic shifts and sea level changes. This genome serves as a valuable reference for future studies of population genomics and differentiation and provides a much-needed genomic resource for this western Atlantic sparid.

The structural landscape and diversity of Pyricularia oryzae MAX effectors revisited.

Magnaporthe AVRs and ToxB-like (MAX) effectors constitute a family of secreted virulence proteins in the fungus Pyricularia oryzae (syn. Magnaporthe oryzae), which causes blast disease on numerous cereals and grasses. In spite of high sequence divergence, MAX effectors share a common fold characterized by a ß-sandwich core stabilized by a conserved disulfide bond. In this study, we investigated the structural landscape and diversity within the MAX effector repertoire of P. oryzae. Combining experimental protein structure determination and in silico structure modeling we validated the presence of the conserved MAX effector core domain in 77 out of 94 groups of orthologs (OG) identified in a previous population genomic study. Four novel MAX effector structures determined by NMR were in remarkably good agreement with AlphaFold2 (AF2) predictions. Based on the comparison of the AF2-generated 3D models we propose a classification of the MAX effectors superfamily in 20 structural groups that vary in the canonical MAX fold, disulfide bond patterns, and additional secondary structures in N- and C-terminal extensions. About one-third of the MAX family members remain singletons, without strong structural relationship to other MAX effectors. Analysis of the surface properties of the AF2 MAX models also highlights the high variability within the MAX family at the structural level, potentially reflecting the wide diversity of their virulence functions and host targets.

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants.

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BackgroundGenomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.ResultsHere, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.ConclusionsOur results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.