Genome-wide transcriptome analyses of developing seeds from low and normal phytic acid soybean lines

Background: Low phytic acid (lpa) crops are potentially eco-friendly alternative to conventional normal phytic acid (PA) crops, improving mineral bioavailability in monogastric animals as well as decreasing phosphate pollution. The lpa crops developed to date carry mutations that are directly or indirectly associated with PA biosynthesis and accumulation during seed development. These lpa crops typically exhibit altered carbohydrate profiles, increased free phosphate, and lower seedling emergence, the latter of which reduces overall crop yield, hence limiting their large-scale cultivation. Improving lpa crop yield requires an understanding of the downstream effects of the lpa genotype on seed development. Towards that end, we present a comprehensive comparison of gene-expression profiles between lpa and normal PA soybean lines (Glycine max) at five stages of seed development using RNA-Seq approaches. The lpa line used in this study carries single point mutations in a myo-inositol phosphate synthase gene along with two multidrug-resistance protein ABC transporter genes. Results: RNA sequencing data of lpa and normal PA soybean lines from five seed-developmental stages (total of 30 libraries) were used for differential expression and functional enrichment analyses. A total of 4235 differentially expressed genes, including 512-transcription factor genes were identified. Eighteen biological processes such as apoptosis, glucan metabolism, cellular transport, photosynthesis and 9 transcription factor families including WRKY, CAMTA3 and SNF2 were enriched during seed development. Genes associated with apoptosis, glucan metabolism, and cellular transport showed enhanced expression in early stages of lpa seed development, while those associated with photosynthesis showed decreased expression in late developmental stages. The results suggest that lpa-causing mutations play a role in inducing and suppressing plant defense responses during early and late stages of seed development, respectively. Conclusions: This study provides a global perspective of transcriptomal changes during soybean seed development in an lpa mutant. The mutants are characterized by earlier expression of genes associated with cell wall biosynthesis and a decrease in photosynthetic genes in late stages. The biological processes and transcription factors identified in this study are signatures of lpa-causing mutations.