Presence-absence Variation Research Articles

Wheat Panache: A pangenome graph database representing presence-absence variation across sixteen bread wheat genomes.

Bread wheat (Triticum aestivum L.) is one of humanity's most important staple crops, characterized by a large and complex genome with a high level of gene presence-absence variation (PAV) between cultivars, hampering genomic approaches for crop improvement. With the growing global population and the increasing impact of climate change on crop yield, there is an urgent need to apply genomic approaches to accelerate wheat breeding. With recent advances in DNA sequencing technology, a growing number of high-quality reference genomes are becoming available, reflecting the genetic content of a diverse range of cultivars. However, information on the presence or absence of genomic regions has been hard to visualize and interrogate because of the size of these genomes and the lack of suitable bioinformatics tools. To address this limitation, we have produced a wheat pangenome graph maintained within an online database to facilitate interrogation and comparison of wheat cultivar genomes. The database allows users to visualize regions of the pangenome to assess PAV between bread wheat genomes.

The genome sequencing and comparative analysis of a wild kiwifruit Actinidia eriantha

The current kiwifruit industry is mainly based on the cultivars derived from the species Actinidia chinensis (Ac) which may bring risks such as canker disease. Introgression of desired traits from wild relatives is an important method for improving kiwifruit cultivars. Actinidia eriantha (Ae) is a particularly important taxon used for hybridization or introgressive breeding of new kiwifruit cultivars because of its valued species-specific traits. Here, we assembled a chromosome-scale high-quality genome of a Ae sample which was directly collected from its wild populations. Our analysis revealed that 41.3% of the genome consists of repetitive elements, comparable to the percentage in Ac and Ae cultivar “White” genomes. The genomic structural variation, including the presence/absence-variation (PAV) of genes, is distinct between Ae and Ac, despite both sharing the same two kiwifruit-specific whole genome duplication (WGD) events. This suggests that a post-WGD divergence mechanism occurred during their evolution. We further investigated genes involved in ascorbic acid biosynthesis and disease-resistance of Ae, and we found introgressive genome could contribute to the complex relationship between Ae and other representative kiwifruit taxa. Collectively, the Ae genome offers valuable genetic resource to accelerate kiwifruit breeding applications.

Expression complementation of gene presence/absence polymorphisms in hybrids contributes importantly to heterosis in sunflower

IntroductionNumerous crops have transitioned to hybrid seed production to increase yields and yield stability through heterosis. However, the molecular mechanisms underlying heterosis and its stability across environments are not yet fully understood. ObjectivesThis study aimed to (1) elucidate the genetic and molecular mechanisms underlying heterosis in sunflower, and (2) determine how heterosis is maintained under different environments. MethodsGenome-wide association (GWA) analyses were employed to assess the effects of presence/absence variants (PAVs) and stop codons on 16 traits phenotyped in the sunflower association mapping population at three locations. To link the GWA results to transcriptomic variation, we sequenced the transcriptomes of two sunflower cultivars and their F1 hybrid (INEDI) under both control and drought conditions and analyzed patterns of gene expression and alternative splicing. ResultsThousands of PAVs were found to affect phenotypic variation using a relaxed significance threshold, and at most such loci the “absence” allele reduced values of heterotic traits, but not those of non-heterotic traits. This pattern was strengthened for PAVs that showed expression complementation in INEDI. Stop codons were much rarer than PAVs and less likely to reduce heterotic trait values. Hybrid expression patterns were enriched for the GO category, sensitivity to stimulus, but all genotypes responded to drought similarily – by up-regulating water stress response pathways and down-regulating metabolic pathways. Changes in alternative splicing were strongly negatively correlated with expression variation, implying that alternative splicing in this system largely acts to reinforce expression responses. ConclusionOur results imply that complementation of expression of PAVs in hybrids is a major contributor to heterosis in sunflower, consistent with the dominance model of heterosis. This mechanism can account for yield stability across different environments. Moreover, given the much larger numbers of PAVs in plant vs. animal genomes, it also offers an explanation for the stronger heterotic responses seen in the former.

Intraspecific de novo gene birth revealedby presence-absence variant genes in Caenorhabditis elegans.

Genes embed their evolutionary history in the form of various alleles. Presence–absence variants (PAVs) are extreme cases of such alleles, where a gene present in one haplotype does not exist in another. Because PAVs may result from either birth or death of a gene, PAV genes and their alternative alleles, if available, can represent a basis for rapid intraspecific gene evolution. Using long-read sequencing technologies, this study traced the possible evolution of PAV genes in the PD1074 and CB4856 C. elegans strains as well as their alternative alleles in 14 other wild strains. We updated the CB4856 genome by filling 18 gaps and identified 46 genes and 7,460 isoforms from both strains not annotated previously. We verified 328 PAV genes, out of which 46 were C. elegans-specific. Among these possible newly born genes, 12 had alternative alleles in other wild strains; in particular, the alternative alleles of three genes showed signatures of active transposons. Alternative alleles of three other genes showed another type of signature reflected in accumulation of small insertions or deletions. Research on gene evolution using both species-specific PAV genes and their alternative alleles may provide new insights into the process of gene evolution.

Pan-Genome Analysis Reveals the Abundant Gene Presence/Absence Variations Among Different Varieties of Melon and Their Influence on Traits.

Melon (Cucumismelo L.) is an important vegetable crop that has been subjected to domestication and improvement. Several varieties of melons with diverse phenotypes have been produced. In this study, we constructed a melon pan-genome based on 297 accessions comprising 168 Mb novel sequences and 4,325 novel genes. Based on the results, there were abundant genetic variations among different melon groups, including 364 unfavorable genes in the IMP_A vs. LDR_A group, 46 favorable genes, and 295 unfavorable genes in the IMP_M vs. LDR_M group. The distribution of 709 resistance gene analogs (RGAs) was also characterized across 297 melon lines, of which 603 were core genes. Further, 106 genes were found to be variable, 55 of which were absent in the reference melon genome. Using gene presence/absence variation (PAV)-based genome-wide association analysis (GWAS), 13 gene PAVs associated with fruit length, fruit shape, and fruit width were identified, four of which were located in pan-genome additional contigs.

Fine mapping of Pm58 from Aegilops tauschii conferring powdery mildew resistance.

The powdery mildew resistance gene Pm58 was traced to a 141.3-kb interval with the co-segregating marker Xkasp68500 in wheat breeding. Pm58 is a powdery mildew resistance gene identified in Aegilops tauschii accession TA1662 and effective in a common wheat background. To finely map Pm58, an F2 population of 676 plants derived from the cross T093 × TA1662 was used for recombinant screening. We obtained 13 recombinants that occurred between the flanking markers Xhnu670 and Xhnu186. Genotyping and phenotyping these recombinant F2:3 families delimited Pm58 to a 0.22-cM interval (Xsts20220-Xkasp61553) on chromosome arm 2DS. The region carrying the Pm58 locus was approximately 141.3-kb, which contained eight annotated genes according to the reference genome sequence of Ae. tauschii AL8/78. Haplotype analysis of 178 Ae. tauschii accessions using the candidate gene-specific markers identified a disease resistance gene AET2Gv20068500 as a candidate for Pm58. Comparative mapping of the Pm58-containing interval revealed two presence/absence variations (PAVs) between AL8/78 and common wheat Chinese Spring. PAV-1 resides in the 3'-end of AET2Gv20068500. The majority of 158 common wheat cultivars (84.8%) displayed the absence of a 14.1-kb fragment in the PAV-1 region, which was confirmed by aligning the targeted genome sequences of the other sequenced Ae. tauschii accessions and common wheat cultivars. A co-segregating marker Xkasp68500 developed from AET2Gv20068500 can distinguish TA1662 from all randomly selected common wheat cultivars and will be instrumental for tracking Pm58 in breeding programs.

The power of chromosome-scale, haplotype-resolved genomes

The assembly of chromosome-scale and haplotype-resolved reference genomes is now more easily attainable, largely due to various improvements in both assembly algorithms and long-read sequencing technologies (see recent review by Michael and VanBuren, 2020 Michael T.P. VanBuren R. Building near-complete plant genomes. Curr. Opin. Plant Biol. 2020; 54: 26-33 Google Scholar ). Due to these technological advancements, there has been a shift away from using highly inbred accessions, in some instances double haploids, to instead assembling the genomes of important reference genotypes. For example, a haplotype-resolved genome was recently assembled for an important cultivated garden strawberry (Fragaria x ananassa) cultivar Royal Royce, which is not only highly heterozygous but also octoploid (8×) ( Hardigan et al., 2021 Hardigan M.A. Feldmann M.J. Pincot D.D.A. Famula R.A. Vachev M.V. Madera M.A. Zerbe P. Mars K. Peluso P. Rank D. et al. Blueprint for phasing and assembling the genomes of heterozygous polyploids: application to the octoploid genome of strawberry. Preprint at bioRxiv. 2021; https://doi.org/10.1101/2021.11.03.467115 Google Scholar ). Similarly, haplotype-phased genomes have been assembled for other important crops, including cassava (Manihot esculenta) ( Hu et al., 2021 Hu W. Ji C. Shi H. Liang Z. Ding Z. Ye J. Ou W. Zhou G. Tie W. Yan Y. et al. Allele-defined genome reveals biallelic differentiation during cassava evolution. Mol. Plant. 2021; 14: 851-854 Google Scholar ) and potato (Solanum tuberosum) ( Hoopes et al., 2022 Hoopes G. Meng X. Hamilton J.P. Achakkagari S.R. de Alves Freitas Guesdes F. Bolger M.E. Coombs J.J. Esselink D. Kaiser N.R. Kodde L. et al. Phased, chromosome-scale genome assemblies of tetraploid potato reveals a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity. Mol. Plant. 2022; (S1674-2052, 00003-X)https://doi.org/10.1016/j.molp.2022.01.003 Google Scholar ). This has resulted in the development of new genomic resources that can be directly used in guiding molecular breeding efforts, without worrying about the gene presence-absence variation that exists within a crop or species ( Golicz et al., 2016 Golicz A.A. Batley J. Edwards D. Towards plant pangenomics. Plant Biotechnol. J. 2016; 14: 1099-1105 Google Scholar ; Bayer et al., 2021 Bayer P.E. Scheben A. Golicz A.A. Yuan Y. Faure S. Lee H. Chawla H.S. Anderson R. Bancroft I. Raman H. et al. Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids. Plant Biotechnol. J. 2021; 19: 2488-2500 Google Scholar ). An additional benefit of having a reference genome of an important cultivar and/or genotype, versus from a highly inbred line, is that valuable insights that are potentially gained from the analysis of genetic variants present among haplotypes are not lost, allowing for additional insights.

Genome assembly of the JD17 soybean provides a new reference genome for comparative genomics.

Cultivated soybean (Glycine max) is an important source for protein and oil. Many elite cultivars with different traits have been developed for different conditions. Each soybean strain has its own genetic diversity, and the availability of more high-quality soybean genomes can enhance comparative genomic analysis for identifying genetic underpinnings for its unique traits. In this study, we constructed a high-quality de novo assembly of an elite soybean cultivar Jidou 17 (JD17) with chromosome contiguity and high accuracy. We annotated 52,840 gene models and reconstructed 74,054 high-quality full-length transcripts. We performed a genome-wide comparative analysis based on the reference genome of JD17 with 3 published soybeans (WM82, ZH13, and W05), which identified 5 large inversions and 2 large translocations specific to JD17, 20,984–46,912 presence–absence variations spanning 13.1–46.9 Mb in size. A total of 1,695,741–3,664,629 SNPs and 446,689–800,489 Indels were identified and annotated between JD17 and them. Symbiotic nitrogen fixation genes were identified and the effects from these variants were further evaluated. It was found that the coding sequences of 9 nitrogen fixation-related genes were greatly affected. The high-quality genome assembly of JD17 can serve as a valuable reference for soybean functional genomics research.

Pan-genome-based rice terpene synthase (OsTPS) gene family analysis reveals the gene presence-absence variations and atypical gene structures caused by structure variations

Pan-genome-based rice terpene synthase (OsTPS) gene family analysis reveals the gene presence-absence variations and atypical gene structures caused by structure variations

Differential Regulation and Production of Secondary Metabolites among Isolates of the Fungal Wheat Pathogen Zymoseptoria tritici.

ABSTRACTThe genome of the wheat-pathogenic fungus Zymoseptoria tritici represents extensive presence-absence variation in gene content. Here, we addressed variation in biosynthetic gene cluster (BGC) content and biochemical profiles among three isolates. We analyzed secondary metabolite properties based on genome, transcriptome, and metabolome data. The isolates represent highly distinct genome architecture but harbor similar repertoires of BGCs. Expression profiles for most BGCs show comparable patterns of regulation among the isolates, suggesting a conserved biochemical infection program. For all three isolates, we observed a strong upregulation of a putative abscisic acid (ABA) gene cluster during biotrophic host colonization, indicating that Z. tritici interferes with host defenses by the biosynthesis of this phytohormone. Further, during in vitro growth, the isolates show similar metabolomes congruent with the predicted BGC content. We assessed if secondary metabolite production is regulated by histone methylation using a mutant impaired in formation of facultative heterochromatin (H3K27me3). In contrast to other ascomycete fungi, chromatin modifications play a less prominent role in regulation of secondary metabolites. In summary, we show that Z. tritici has a conserved program of secondary metabolite production, contrasting with the immense variation in effector expression, and some of these metabolites might play a key role during host colonization.IMPORTANCE Zymoseptoria tritici is one of the most devastating pathogens of wheat. So far the molecular determinants of virulence and their regulation are poorly understood. Previous studies have focused on proteinaceous virulence factors and their extensive diversity. In this study, we focus on secondary metabolites produced by Z. tritici. Using a comparative framework, we characterize core and noncore metabolites produced by Z. tritici by combining genome, transcriptome, and metabolome data sets. Our findings indicate highly conserved biochemical profiles with contrasting genetic and phenotypic diversity of the field isolates investigated here. This discovery has relevance for future crop protection strategies.

GC content of plant genes is linked to past gene duplications.

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a “biased gene duplication / biased mutation” (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers.

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

Detecting Structural Variants and Associated Gene Presence-Absence Variation Phenomena in the Genomes of Marine Organisms.

As complete genomes become easier to attain, even from previously difficult-to-sequence species, and as genomic resequencing becomes more routine, it is becoming obvious that genomic structural variation is more widespread than originally thought and plays an important role in maintaining genetic variation in populations. Structural variants (SVs) and associated gene presence-absence variation (PAV) can be important players in local adaptation, allowing the maintenance of genetic variation and taking part in other evolutionarily relevant phenomena. While recent studies have highlighted the importance of structural variation in Mollusca, the prevalence of this phenomenon in the broader context of marine organisms remains to be fully investigated.Here, we describe a straightforward and broadly applicable method for the identification of SVs in fully assembled diploid genomes, leveraging the same reads used for assembly. We also explain a gene PAV analysis protocol, which could be broadly applied to any species with a fully sequenced reference genome available. Although the strength of these approaches have been tested and proven in marine invertebrates, which tend to have high levels of heterozygosity, possibly due to their lifestyle traits, they are also applicable to other species across the tree of life, providing a ready means to begin investigations into this potentially widespread phenomena.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice.

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

Sugarcane Mosaic Virus Resistance in the Wisconsin Sweet Corn Diversity Panel

Sugarcane mosaic virus [SCMV (Potyvirus sugarcane mosaic virus)] is an ssRNA virus that negatively affects yield in maize (Zea mays) worldwide. Resistance to SCMV is controlled primarily by a single dominant gene (Scm1). The goal of this study was to identify sweet corn (Z. mays) inbreds that demonstrate resistance to SCMV, confirm the presence of genomic regions previously identified in maize associated with resistance, and identify other resistant loci in sweet corn. Eight plants from each of 563 primarily sweet corn inbred lines were tested for SCMV resistance. Plants were inoculated 14 d after planting and observed for signs of infection 24 d after planting. A subset of 420 inbred lines were genotyped using 7504 high-quality genotyping-by-sequencing single-nucleotide polymorphism markers. Population structure of the panel was observed, and a genome-wide association study was conducted to identify loci associated with SCMV resistance. Forty-six of the inbreds were found to be resistant to SCMV 10 d after inoculation. The Scm1 locus was confirmed with the presence of two significant loci on chromosome 6 (P = 2.5 × 10−8 and 1.6 × 10−8), 5 Mb downstream of the Scm1 gene previously located at Chr6: 14194429.14198587 and the surrounding 2.7-Mb presence–absence variation. We did not identify other loci associated with resistance. This research has increased information on publicly available SCMV-resistant germplasm useful to future breeding projects and demonstrated that SCMV resistance in this sweet corn panel is driven by the Scm1 gene.

Investigating the Transcriptomic and Expression Presence-Absence Variation Exist in Japanese Eel (Anguilla japonica), a Primitive Teleost.

The pan-genome was defined as the complete gene set across strains, and it is built upon genes displaying presence-absence variations (PAVs); the pan-transcriptome is defined by recalling the pan-genome. Indeed, a PAV is reflected from the expression presence-absence variation (ePAV). In this study, treated with androgen, eels, which are a primitive fish from the basal lineage of Teleost, with different ovarian developments were chosen and submitted to RAN-sequencing. Transcriptomes were the assembly against eel genome scaffolds; a pair was the unit (the same eel before and after treatment) to analyze DEGs (differentially expressed genes); the core, unique, or accessory genes were identified, and the list of DEGs was analyzed to investigate ePAV. The results suggest that there was ePAV in Japanese eel, and the ePAV of eel was analyzed by pathway enrichment. These results signify the importance of genetic differential expression on the variations of phenotypes by androgen, and a transcriptomic approach appears to enable extracting multiple layers of genomic data.

Panache: a web browser-based viewer for linearized pangenomes.

MotivationPangenomics evolved since its first applications on bacteria, extending from the study of genes for a given population to the study of all of its sequences available. While multiple methods are being developed to construct pangenomes in eukaryotic species there is still a gap for efficient and user-friendly visualization tools. Emerging graph representations come with their own challenges, and linearity remains a suitable option for user-friendliness.ResultsWe introduce Panache, a tool for the visualization and exploration of linear representations of gene-based and sequence-based pangenomes. It uses a layout similar to genome browsers to display presence absence variations and additional tracks along a linear axis with a pangenomics perspective.Availability and implementationPanache is available at github.com/SouthGreenPlatform/panache under the MIT License.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids.

SummaryPlant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.

Megabase-scale presence-absence variation with Tripsacum origin was under selection during maize domestication and adaptation

BackgroundStructural variants (SVs) significantly drive genome diversity and environmental adaptation for diverse species. Unlike the prevalent small SVs (< kilobase-scale) in higher eukaryotes, large-size SVs rarely exist in the genome, but they function as one of the key evolutionary forces for speciation and adaptation.ResultsIn this study, we discover and characterize several megabase-scale presence-absence variations (PAVs) in the maize genome. Surprisingly, we identify a 3.2 Mb PAV fragment that shows high integrity and is present as complete presence or absence in the natural diversity panel. This PAV is embedded within the nucleolus organizer region (NOR), where the suppressed recombination is found to maintain the PAV against the evolutionary variation. Interestingly, by analyzing the sequence of this PAV, we not only reveal the domestication trace from teosinte to modern maize, but also the footprints of its origin from Tripsacum, shedding light on a previously unknown contribution from Tripsacum to the speciation of Zea species. The functional consequence of the Tripsacum segment migration is also investigated, and environmental fitness conferred by the PAV may explain the whole segment as a selection target during maize domestication and improvement.ConclusionsThese findings provide a novel perspective that Tripsacum contributes to Zea speciation, and also instantiate a strategy for evolutionary and functional analysis of the “fossil” structure variations during genome evolution and speciation.

The Chicken Pan-Genome Reveals Gene Content Variation and a Promoter Region Deletion in IGF2BP1 Affecting Body Size.

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional approximately 66.5-Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression levels are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based genome-wide association studies identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size quantitative trait locus located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.