Additional file 1 of Extensive sequence duplication in Arabidopsis revealed by pseudo-heterozygosity
Creators: name: Jaegle, Benjamin
name: Pisupati, Rahul
name: Soto-Jiménez, Luz Mayela
name: Burns, Robin
name: Rabanal, Fernando A.
name: Nordborg, Magnus
Publisher: name: figshare
Year: 2023
Subjects: subject: Genetics
subject: FOS: Biological sciences

description: Additional file 1: Fig S1. Distribution of heterozygosity across individuals and the genome. Fig S2. Tract-length distribution. Fig S3. SNP filtering scheme. The SNPs matrix we started with contains 7 million SNPs. Of those, 3.3 million were detected as heterozygous in at least one line. We selected the 48,799 SNPs that we called heterozygous in at least 5 % of the lines, and focused on the 26,647 found in coding regions (according to the TAIR10 annotation). Fig S4. The distribution of distances between the “position” of a pseudo-SNPs and the corresponding GWAS peak. Out of the 411688 peaks we found that 236494 (57%) are on a different chromosome. The plots show the distribution of distances for the 175194 peaks that were on the same chromosome as the corresponding pseudo-SNPs. Roughly 50% of there are more distance than 50kb. Fig S5. Gene ontology analysis of the putatively duplicated genes. Fig S6. Pipeline to confirm GWAS peaks. Fig S7. Mapping reads tagging singleton pseudo-SNP. First, all reads overlapping the position of a specific pseudo-SNP were extracted based on mapping to the reference genome (TAIR10). This set of reads were then re-mapped to the appropriate Pacbio genome. Reads mapping to multiple regions indicate the presence of a duplicated segment. A decrease in coverage compared to the mapping to the reference genome is also a confirmation that reads map at different positions. An example is presented in Fig S9. Fig S8. Comparison of mapping coverage between reference genome and PacBio genomes for all regions surrounding pseudo-heterozygous doubleton positions. Fig S9. Mapping of reads overlapping singletons (example). Fig S10. Position of AT1G31390 in the reference (accession 6909) and each of the six newly assembled genomes. BLAST-thresholds of 70% identity were used, and only fragments of length greater than 50% of the original gene length are shown. Fig S11. Position of AT1G20400 in the reference (accession 6909) and each of the six newly assembled genomes. Cf. Fig S10. Fig S12. Dot-plots of the end of insertion B. LTR repeats can be detected on each side of the insertion. Cf. Fig S13. Fig S13. Dot plot of the ends of all insertions. Cf. Fig S12. Fig S14. Methylation profile across all copies of the TE+gene insertions. Bisulfite reads were mapped to the appropriate genomes and profiles extracted based on the inferred locations. Colors distinguish the different insertions found in each genome. Fig S15. Frequency of insertions of the new element in the 1001 Genomes. Fig S16. Mapping of the AT1G20400 region in multiple species. The rectangle corresponds to annotated genes around AT1G20400 in the A. thaliana reference genome. Fig S17. Chromosomal position of putatively duplicated genes that could not be confirmed using full genome sequence. The figure shows the distribution of these genes compared to the distribution of all genes, and RepeatMasker annotation. Fig S18. Illustration of GWAS peak calling. The p-values from GWAS were used to run the filterFFT function from the NucleR package in R, using 0.05 for pcKeepComp option. Fig S19. The relationship between the number of pseudo-SNPs and sequencing coverage. Each dot represents an accession.