The pan-genome and local adaptation of Arabidopsis thaliana

Arabidopsis thaliana has been used as a model species for research in a diverse collection of plant species. However, previous studies based on single reference genomes and short-read sequencing data are restricted to detecting variable genes and large structural variation (SV) underlying local adaptation. Here we de novo assemble high-quality chromosomal genomes of 38 A. thaliana ecotypes (with 6 relict ones) using PacBio-HiFi long-read sequencing. From these newly assembled genomes, we annotate several thousand new genes through pan-genomic analysis in comparison to the previous reference genome. The identified variable genes are mainly enriched in and associated with ecological adaptation and this species substantially expands its gene repertoire for local adaptation. We construct a graph-based pan-genome and identify 62,525 SVs which overlap with 14,243 genes. These genes are enriched in multiple ecological adaptation functions, including secondary metabolic processes, enzyme regulation, and biotic/abiotic stimulus. For example, a 566 bp insertion in the promoter of the light-adaptation KNAT3 gene was specific to the high-altitude relict Tibet-0 ecotype. This SV reduces the expression level of KNAT3 and promotes A. thaliana adaptation to habitats high in light radiation. In addition, compared with the SNPs, the SVs identified in this study captured the missing heritability and we detected novel SV associations with environmental variables in their native range, highlighting the value of SVs in environmental adaptation. The genome resources presented here will help pinpoint genetic changes that include both SVs and the ecotype-specific genes for local adaptation of A. thaliana and increase our understanding of the molecular mechanisms in this model species to respond to varied habitats.