Integration of eQTL and GWAS analysis uncovers a genetic regulation of natural ionomic variation in Arabidopsis

Abstract The functional interpretation of traits associated variants by expression quantitative trait loci (eQTL) analysis is usually performed in bulk tissue samples. While the regulation of gene expression is context-dependent, such as cell-type-specific manner. In this study, we estimated cell type abundances from 728 bulk tissue samples using single-cell RNA-sequencing dataset, and performed cis-eQTL mapping to identify cell-type interaction eQTL (cis-eQTLs(ci)) in A. thaliana. Also, we performed Genome-wide association studies (GWAS) analyses for 999 accessions to identify the genetic basis of variations in A. thaliana leaf ionome. As a result, a total of 5,664 unique eQTL genes and 15,038 unique cis-eQTLs(ci) were significant. The majority (62.83%) of cis-eQTLs(ci) were cell-type-specific eQTLs. Using colocalization, we uncovered one interested gene AT2G25590 in Phloem cell, encoding a kind of plant Tudor-like protein with possible chromatin-associated functions, which colocalized with the most significant cis-eQTL(ci) of a Mo-related locus (Chr2:10908806:A:C; P = 3.27×10-27). Furthermore, we prioritized eight target genes associated with AT2G25590, which were previously reported in regulating the concentration of Mo element in A. thaliana. This study revealed the genetic regulation of ionomic variations and provided a foundation for further studies on molecular mechanisms of genetic variants controlling the A. thaliana ionome.