Transcriptome brings variations of gene expression, alternative splicing, and structural variations into gene-scale trait dissection in soybean

Genome-wide association study (GWAS) identifies trait-associated loci, but due in part to slow decay of linkage disequilibrium (LD), identifying the causal genes can be a bottleneck. Transcriptome-wide association study (TWAS) addresses this by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci (eQTLs) with GWAS results. Here, we used self-pollinated soybean as a model to evaluate the application of TWAS in the genetic dissection of traits in plant species that exhibit slow LD decay. The first RNA-Seq analysis of a soybean diversity panel was conducted, which identified the genetic regulation of 29,286 genes. Different TWAS solutions were less affected by LD and robust with source of expression that identified known genes related to traits from different development stages and tissues. A novel gene named pod color L2 was identified via TWAS and functionally validated by genome editing. Our introduction of the new exon proportion feature significantly improves the capture of expression variations resulting from structural variations and alternative splicing. As a result, the genes identified by our TWAS approach exhibited a diverse range of causal variations, including SNP, insertion/deletion, gene fusion, copy number variation, and alternative splicing. Using our TWAS approach, we identified genes associated with flowering time, including both previously known candidates and novel genes that had not been linked to this trait before, providing complementary insights with GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.