Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally‐inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole‐genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5 Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S‐locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold‐responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.     

Genomewide investigation of adaptation to harmful algal blooms in common bottlenose dolphins (Tursiops truncatus)

Harmful algal blooms (HABs), which can be lethal in marine species and cause illness in humans, are increasing worldwide. In the Gulf of Mexico, HABs of Karenia brevis produce neurotoxic brevetoxins that cause large‐scale marine mortality events. The long history of such blooms, combined with the potentially severe effects of exposure, may have produced a strong selective pressure for evolved resistance. Advances in next‐generation sequencing, in particular genotyping‐by‐sequencing, greatly enable the genomic study of such adaptation in natural populations. We used restriction site‐associated DNA (RAD) sequencing to investigate brevetoxicosis resistance in common bottlenose dolphins (Tursiops truncatus). To improve our understanding of the epidemiology and aetiology of brevetoxicosis and the potential for evolved resistance in an upper trophic level predator, we sequenced pools of genomic DNA from dolphins sampled from both coastal and estuarine populations in Florida and during multiple HAB‐associated mortality events. We sequenced 129 594 RAD loci and analysed 7431 single nucleotide polymorphisms (SNPs). The allele frequencies of many of these polymorphic loci differed significantly between live and dead dolphins. Some loci associated with survival showed patterns suggesting a common genetic‐based mechanism of resistance to brevetoxins in bottlenose dolphins along the Gulf coast of Florida, but others suggested regionally specific mechanisms of resistance or reflected differences among HABs. We identified candidate genes that may be the evolutionary target for brevetoxin resistance by searching the dolphin genome for genes adjacent to survival‐associated SNPs.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Lineage‐specific adaptation to climate involves flowering time in North American Arabidopsis lyrata

Adaptation to local climatic conditions is commonly found within species, but whether it involves the same intraspecific genomic variants is unknown. We studied this question in North American Arabidopsis lyrata, whose current distribution is shaped by post‐glacial range expansion from two refugia, resulting in two distinct genetic clusters covering comparable climatic gradients. Using pooled whole‐genome sequence data of 41 outcrossing populations, we identified loci associated with three niche‐determining climatic variables in the two clusters and compared these outliers. Little evidence was found for parallelism in climate adaptation for single nucleotide polymorphisms (SNPs) and for genes with an accumulation of outlier SNPs. Significantly increased selection coefficients supported them as candidates of climate adaptation. However, the fraction of gene ontology (GO) terms shared between clusters was higher compared to outlier SNPs and outlier genes, suggesting that selection acts on similar pathways but not necessarily the same genes. Enriched GO terms involved responses to abiotic and biotic stress, circadian rhythm and development, with flower development and reproduction being among the most frequently detected. In line with GO enrichment, regulators of flowering time were detected as outlier genes. Our results suggest that while adaptation to environmental gradients on the genomic level are lineage‐specific in A. lyrata, similar biological processes seem to be involved. Differential loss of standing genetic variation, probably driven by genetic drift, can in part account for the lack of parallel evolution on the genomic level.     

Genomic variation at the tips of the adaptive radiation of Darwin's finches

Adaptive radiation unfolds as selection acts on the genetic variation underlying functional traits. The nature of this variation can be revealed by studying the tips of an ongoing adaptive radiation. We studied genomic variation at the tips of the Darwin's finch radiation; specifically focusing on polymorphism within, and variation among, three sympatric species of the genus Geospiza. Using restriction site‐associated DNA (RAD‐seq), we characterized 32 569 single‐nucleotide polymorphisms (SNPs), from which 11 outlier SNPs for beak and body size were uncovered by a genomewide association study (GWAS). Principal component analysis revealed that these 11 SNPs formed four statistically linked groups. Stepwise regression then revealed that the first PC score, which included 6 of the 11 top SNPs, explained over 80% of the variation in beak size, suggesting that selection on these traits influences multiple correlated loci. The two SNPs most strongly associated with beak size were near genes associated with beak morphology across deeper branches of the radiation: delta‐like 1 homologue (DLK1) and high‐mobility group AT‐hook 2 (HMGA2). Our results suggest that (i) key adaptive traits are associated with a small fraction of the genome (11 of 32 569 SNPs), (ii) SNPs linked to the candidate genes are dispersed throughout the genome (on several chromosomes), and (iii) micro‐ and macro‐evolutionary variation (roots and tips of the radiation) involve some shared and some unique genomic regions.     

Abundant recent activity of retrovirus‐like retrotransposons within and among flycatcher species implies a rich source of structural variation in songbird genomes

Transposable elements (TEs) are genomic parasites capable of inserting virtually anywhere in the host genome, with manifold consequences for gene expression, DNA methylation and genomic stability. Notably, they can contribute to phenotypic variation and hence be associated with, for example, local adaptation and speciation. However, some organisms such as birds have been widely noted for the low densities of TEs in their genomes and this has been attributed to a potential dearth in transposition during their evolution. Here, we show that avian evolution witnessed diverse and abundant transposition on very recent timescales. First, we made an in‐depth repeat annotation of the collared flycatcher genome, including identification of 23 new, retrovirus‐like LTR retrotransposon families. Then, using whole‐genome resequencing data from 200 Ficedula flycatchers, we detected 11,888 polymorphic TE insertions (TE presence/absence variations, TEVs) that segregated within and among species. The density of TEVs was one every 1.5–2.5 Mb per individual, with heterozygosities of 0.12–0.16. The majority of TEVs belonged to some 10 different LTR families, most of which are specific to the flycatcher lineage. TEVs were validated by tracing the segregation of hundreds of TEVs across a three‐generation pedigree of collared flycatchers and also by their utility as markers recapitulating the phylogenetic relationships among flycatcher species. Our results suggest frequent germline invasions of songbird genomes by novel retroviruses as a rich source of structural variation, which may have had underappreciated phenotypic consequences for the diversification of this species‐rich group of birds.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

珍稀濒危飘带兜兰叶绿体全基因组种内变异研究(附录)

飘带兜兰(Paphiopedilum parishii)分布范围狭窄,仅在中国、缅甸、泰国以及老挝有少量分布。近年来,因生境破坏和人为滥采而导致飘带兜兰野生种群极度缩减。为开发种内多态性的分子标记用于保护生物学研究,该研究对飘带兜兰4个野生个体经测序、组装、注释获得的叶绿体基因组序列,与已公布的飘带兜兰2个个体的叶绿体全基因组序列进行比对,分析飘带兜兰叶绿体基因组的种内差异。结果表明:(1)飘带兜兰叶绿体基因组具有典型被子植物叶绿体基因组环状四分体结构,基因组长度为154 403～154 809 bp,共编码129个基因,包括78个蛋白质编码基因、39个tRNA基因、8个rRNA基因,以及4个假基因。(2)在飘带兜兰6个个体叶绿体基因组中检测到103～107个SSRs(simple sequence repeats)位点,其中21个SSR位点具有多态性。此外,在6个个体叶绿体基因组中还检测到60个长序列重复,包括17～21个正向重复、18～29个反向重复、9～16个回文重复、4～9个互补重复。(3)通过比较6个个体叶绿体基因组序列的核苷酸多样性,共发现70处变异,包括10个SNPs(single nucleotide polymorphism)、60个插入缺失(InDels)。其中,有3个SNP位点发生了非同义替换,导致编码功能基因的氨基酸发生改变; 19个插入缺失多态性较高,具有开发为分子标记的潜力。(4)通过计算核苷酸多样性值(P_i)共发现8个有变异的区域,P_i值为0～0.006 32,其中变异度较大的是rps3-rpl22、trnL-UAC-rpl32、rpoB-trnC-GCA以及ycf4,这些高变区可开发为分子标记用于评估飘带兜兰遗传多样性。(5)系统发生分析结果表明,飘带兜兰6个个体叶绿体基因组序列聚在一起,与长瓣兜兰互为姐妹群。综上表明,飘带兜兰叶绿体基因组的SSRs、长序列重复、SNPs、InDels以及核苷酸序列呈现了足够的种内多样性,可开发成分子标记用于该种的系统演化及保护生物学研究。     

A resource of genome‐wide single‐nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel

Reduced representation genome sequencing such as restriction‐site‐associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single‐nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome‐wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome‐wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19 703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long‐term effective population size was estimated to range between 132 000 and 1 320 000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82 425 loci and 376 918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.     

Specific alleles at immune genes,rather than genome‐wide heterozygosity,are related to immunity and survival in the critically endangered Attwater's prairie‐chicken

The negative effects of inbreeding on fitness are serious concerns for populations of endangered species. Reduced fitness has been associated with lower genome‐wide heterozygosity and immune gene diversity in the wild; however, it is rare that both types of genetic measures are included in the same study. Thus, it is often unclear whether the variation in fitness is due to the general effects of inbreeding, immunity‐related genes or both. Here, we tested whether genome‐wide heterozygosity (20 990 SNPs) and diversity at nine immune genes were better predictors of two measures of fitness (immune response and survival) in the endangered Attwater's prairie‐chicken (Tympanuchus cupido attwateri). We found that postrelease survival of captive‐bred birds was related to alleles of the innate (Toll‐like receptors, TLRs) and adaptive (major histocompatibility complex, MHC) immune systems, but not to genome‐wide heterozygosity. Likewise, we found that the immune response at the time of release was related to TLR and MHC alleles, and not to genome‐wide heterozygosity. Overall, this study demonstrates that immune genes may serve as important genetic markers when monitoring fitness in inbred populations and that in some populations specific functional genes may be better predictors of fitness than genome‐wide heterozygosity.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula

Bottom‐up evolutionary approaches, including geographically explicit population genomic analyses, have the power to reveal the mechanistic basis of adaptation. Here, we conduct a population genomic analysis in the model legume, Medicago truncatula, to characterize population genetic structure and identify symbiosis‐related genes showing evidence of spatially variable selection. Using RAD‐seq, we generated over 26,000 SNPs from 191 accessions from within three regions of the native range in Europe. Results from STRUCTURE analysis identify five distinct genetic clusters with divisions that separate east and west regions in the Mediterranean basin. Much of the genetic variation is maintained within sampling sites, and there is evidence for isolation by distance. Extensive linkage disequilibrium was identified, particularly within populations. We conducted genetic outlier analysis with F_ST‐based genome scans and a Bayesian modeling approach (PCAdapt). There were 70 core outlier loci shared between these distinct methods with one clear candidate symbiosis related gene, DMI1. This work sets that stage for functional experiments to determine the important phenotypes that selection has acted upon and complementary efforts in rhizobium populations.     

Chironomus riparius (Diptera) genome sequencing reveals the impact of minisatellite transposable elements on population divergence

Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky–Muller incompatibilities (DMI), thus contributing to the evolution of endogenous genetic barriers and eventually causing population divergence. We investigated differential TE abundance among conspecific populations of the nonbiting midge Chironomus riparius and evaluated their potential role in causing endogenous genetic incompatibilities between these populations. We focussed on a Chironomus‐specific TE, the minisatellite‐like Cla‐element, whose activity is associated with speciation in the genus. Using a newly generated and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population‐specific TE insertion patterns with many private insertions. A significant correlation of the pairwise F_ST estimated from genomewide single‐nucleotide polymorphisms (SNPs) and the F_ST estimated from TEs is consistent with drift as the major force driving TE population differentiation. However, the significantly higher Cla‐element F_ST level due to a high proportion of differentially fixed Cla‐element insertions also indicates selection against segregating (i.e. heterozygous) insertions. With reciprocal crossing experiments and fluorescent in situ hybridization of Cla‐elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings. We propose that the inferred negative selection on heterozygous Cla‐element insertions may cause endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs may thus have a direct impact on population divergence that is operationally different from drift and local adaptation.     

High‐throughput genomics in sorghum: from whole‐genome resequencing to a SNP screening array

With its small, diploid and completely sequenced genome, sorghum (Sorghum bicolor L. Moench) is highly amenable to genomics‐based breeding approaches. Here, we describe the development and testing of a robust single‐nucleotide polymorphism (SNP) array platform that enables polymorphism screening for genome‐wide and trait‐linked polymorphisms in genetically diverse S. bicolor populations. Whole‐genome sequences with 6× to 12× coverage from five genetically diverse S. bicolor genotypes, including three sweet sorghums and two grain sorghums, were aligned to the sorghum reference genome. From over 1 million high‐quality SNPs, we selected 2124 Infinium Type II SNPs that were informative in all six source genomes, gave an optimal Assay Design Tool (ADT) score, had allele frequencies of 50% in the six genotypes and were evenly spaced throughout the S. bicolor genome. Furthermore, by phenotype‐based pool sequencing, we selected an additional 876 SNPs with a phenotypic association to early‐stage chilling tolerance, a key trait for European sorghum breeding. The 3000 attempted bead types were used to populate half of a dual‐species Illumina iSelect SNP array. The array was tested using 564 Sorghum spp. genotypes, including offspring from four unrelated recombinant inbred line (RIL) and F₂ populations and a genetic diversity collection. A high call rate of over 80% enabled validation of 2620 robust and polymorphic sorghum SNPs, underlining the efficiency of the array development scheme for whole‐genome SNP selection and screening, with diverse applications including genetic mapping, genome‐wide association studies and genomic selection.     

Population‐specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome‐wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population‐specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.     

Validating genome‐wide associated signals for clinical mastitis in German Holstein cattle

A validation study for six genomic regions previously identified by a genome‐wide association study for somatic cell score was conducted with data of clinical mastitis in German Holstein cattle. Out of 10 tested SNPs, five on chromosomes 6, 13 and 19 were significantly associated with clinical mastitis (P < 0.05). Three SNPs on chromosomes 6 and 19 had the same direction of effect as those previously reported in the initial genome‐wide association study for somatic cell score. The other two SNPs on chromosome 13 had opposite effects. As well as validating associations within known QTL from previous studies, e.g. chromosomes 6 and 19, novel loci on chromosome 13 were confirmed. Promising candidate genes are, for example: deoxycytidine kinase, immunoglobulin J chain, vitamin D binding protein, forkhead box K2, sodium/hydrogen exchanger 8 and cytoplasmic nuclear factor of activated T‐cells 2. Our confirmation study provides additional evidence for the functional role of the linked genomic regions to immune response. This information can be used as a basis for further functional studies for those potential genes.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

A survey of single nucleotide polymorphisms identified from whole‐genome sequencing and their functional effect in the porcine genome,

Genetic variants detected from sequence have been used to successfully identify causal variants and map complex traits in several organisms. High and moderate impact variants, those expected to alter or disrupt the protein coded by a gene and those that regulate protein production, likely have a more significant effect on phenotypic variation than do other types of genetic variants. Hence, a comprehensive list of these functional variants would be of considerable interest in swine genomic studies, particularly those targeting fertility and production traits. Whole‐genome sequence was obtained from 72 of the founders of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the founding boars (12 Duroc and 12 Landrace) and 48 Yorkshire–Landrace composite sows. Sequence reads were mapped to the Sscrofa10.2 genome build, resulting in a mean of 6.1 fold (×) coverage per genome. A total of 22 342 915 high confidence SNPs were identified from the sequenced genomes. These included 21 million previously reported SNPs and 79% of the 62 163 SNPs on the PorcineSNP60 BeadChip assay. Variation was detected in the coding sequence or untranslated regions (UTRs) of 87.8% of the genes in the porcine genome: loss‐of‐function variants were predicted in 504 genes, 10 202 genes contained nonsynonymous variants, 10 773 had variation in UTRs and 13 010 genes contained synonymous variants. Approximately 139 000 SNPs were classified as loss‐of‐function, nonsynonymous or regulatory, which suggests that over 99% of the variation detected in our pigs could potentially be ignored, allowing us to focus on a much smaller number of functional SNPs during future analyses.