Skip to main content

A seed germination transcriptomic study contrasting two soybean genotypes that differ in terms of their tolerance to the deleterious impacts of elevated temperatures during seed fill



Soybean seed development is negatively impacted by elevated temperatures during seed fill, which can decrease seed quality and economic value. Prior germplasm screens identified an exotic landrace able to maintain ~ 95% seed germination under stress conditions that reduce germination dramatically (> 50%) for typical soybean seeds. Seed transcriptomic analysis was performed for two soybean lines (a heat-tolerant landrace and a typical high-yielding adapted line) for dry, mature seed, 6-h imbibed seed and germinated seed. Seeds were produced in two environments: a typical Midwestern field and a heat stressed field located in the Midsouth soybean production region.


Transcriptomic analysis revealed 23–30K expressed genes in each seed tissue sample, and differentially expressed genes (DEGs) with ≥ twofold gene expression differences (at q-value < 0.05) comprised ~ 5–44% of expressed genes. Gene ontology (GO) enrichment analysis on DEGs revealed enrichment in heat-tolerant seeds for genes annotated for general and temperature-specific stress, as well as protein-refolding. DEGs were also clustered in modules using weighted co-expressed gene network analysis, which were examined for enrichment of GO biological process terms. Collectively, our results provide new and valuable insights into this unique form of genetic abiotic stress tolerance and to soybean seed physiological responses to elevated temperatures.


Soybean (Glycine max L. Merr.) is a major commodity crop, comprising ~ 34% (~ 36.5 million ha) of crop land in the United States in 2017 (, accessed 3-21-19). The value of the crop is principally derived from the high yield and quality of oil and protein in the seeds. The Midsouth soybean growing region of the United States experiences consistent late-season drought, which has resulted in historically reduced on-farm seed yield and economic return [1, 2]. Although irrigation can at least partially remedy these issues, fuel for pumping water is expensive and long term use of aquifers for agricultural irrigation may be unsustainable [3].

Traditionally, soybean maturity group (MG) 5–7 cultivars were planted in May and June with harvest in October and November. An alternative method, the Early Soybean Production System (ESPS) modifies soybean planting and harvest dates to avoid much of the late season endemic drought [2] by use of cultivars that flower and mature earlier (typically MG 3–4 and early 5 as compared with MG5-7) and by adjusting planting dates to early-to-mid-April with harvest typically occurring in September. The practices of the ESPS in the Midsouth region has increased seed yield and on-farm return on investment [1, 3] under both irrigated and non-irrigated conditions [1].

Soybean has traditionally been considered to be heat-tolerant, with a vegetative optimum temperature of ~ 30 °C [4]. However, the processes of pollination and seed growth/maturation are sensitive to elevated temperatures; the reproductive optimum is a relatively low 22–24 °C [5]. Despite economic and seed yield gains under the ESPS, seed produced in this system are exposed to much higher temperatures [6] during seed fill (≥ 32 °C maximum daytime temperature) than seeds of MG 5–7 cultivars produced in the traditional system. In typical MG 4 cultivars, exposure to such high temperatures during seed development reduces seed quality/germination, increases pathogen infection, and often results in economic loss through seed dockage [3, 7].

Soybean is a self-pollinating species, and modern high-yielding cultivars derive from an extremely limited genetic base; traditional breeding has exacerbated this problem [8, 9]. Exotic landraces may contain novel disease and stress resistance genes; a successful screen identified lines that can tolerate the high temperature associated with the ESPS [6]. An unimproved landrace (PI 587982A) has consistent and robust resistance (> 90% germination, near absence of Phomopsis longicolla infection). The first United States heat tolerant germplasm release, with tolerance derived from PI 587982A, was recently made by our group [10].

Transcriptomics, enabled by advances in DNA sequencing and computation, is a powerful tool to identify gene expression differences and correlations with genetic/developmental cues or environmental conditions. Detailed studies have generated “transcriptomic atlases” for soybean gene expression [11,12,13,14,15]. However, studies have ignored soybean seed germination, in favor of seed development or vegetative tissues (typically leaves or roots). In this study, we examined three soybean seed germination stages: (1) dry, mature seed; (2) imbibed seed; and (3) germinated seed and contrasted two soybean genotypes which differ in their tolerance to the impact of elevated temperature on seed quality, using seed produced in two environments differing in abiotic stress: (A) a lower temperature, Midwest location; and (B) the high temperature conditions of the ESPS.

Main text


Field seed production, seed imbibition and germination measurement, and RNA isolation and RNA sequencing, mapping and statistical analysis

Full details are provided in Additional file 1.

GO term enrichment and venn diagrams

GO term enrichment was performed using the tool present on Soybase ( using DEGs identified through Cuffdiff analysis. Venn diagrams were generated using the Venny tool at and the Venn diagram tool at

Whole genome comparative network analysis and gene ontology enrichment of co-expressed gene modules

Modules of genes with highly correlated expression patterns were described using weighted gene co-expression network analysis (WGCNA). We expect these modules to correspond to networks of genes that are co-expressed and thus interact and share biological processes. We constructed unsigned weighted gene co-expression modules using the WGCNA [16] package in R. The blockwiseModules function was run with the Pearson correlation coefficient and a soft thresholding power of 18. The resulting genes modules were named by assigning them different colors arbitrarily. Additionally, we further analyzed each module by conducting significant associations for Gene Ontology (GO) function annotations enrichment analysis (Additional file 5) and used hierarchical clustering to group differentially expressed genes across samples (Additional file 2).

g:GOSt ( was used to examine modules detected by WGCNA, in order to detect statistically significant enriched GO terms within specific modules, using the Benjamin–Hochberg FDR method at α = 0.05 as significant.

qRT-PCR analyses

qRT-PCR analysis was performed as described [17], using the ΔΔCt method [18]. FPKM output was normalized to the cons14 [19] gene (Glyma16g32510) and expressed as log2 ratio for comparison to CuffDuff output.


Germination assays

We examined germination kinetics for two soybean genotypes: (1) a heat-tolerant soybean plant introduction line (PI 587982A) henceforth referred to as “PI”; and (2) S99-11986, a conventional high yielding improved line [20], comparable to cultivars commonly grown in the Midsouth and Midwest regions, henceforth referred to as “SG”. Seed to be germinated were produced (Fig. 1a) in one of two environments: (1) a location with endemic high temperature stress associated with the Early Soybean Production System (henceforth refered to as ESPS—Stoneville, MS; Fig. 1a); or (2) a less stressful Conventional Soybean Production System (CSPS—Columbia, MO; Fig. 1a).

Fig. 1

Details on production of soybean seed and RNAseq analysis of seed germination. a Weather parameters, planting/harvest dates and approximate period of soybean seed fill (stage R5-R8) for genotypes PI 598982A and S99-11986 produced in Columbia, MO (Conventional Soybean Production System) and Stoneville, MS (Early Soybean Production System). b Kinetics of seed germination in this study. Seed of two genotypes, PI 587982A-heat tolerant and S99-11986-heat sensitive, produced in either Columbia, MO (Conventional Soybean Production System) and Stoneville, MS (Early Soybean Production System) were imbibed and germination observed for 72 h (n = 3, 25 seeds per replicated). The mean of three replicates is plotted and error bars indicate standard deviation. c Total genes expressed in each sample with an average FPKM ≥ 0.3. d Distribution of differential gene expression amongst 20 different comparisons. Upregulation indicates higher level of expression in the top sample as compared to the bottom sample. The names indicated in the horizontal axis for c and d are three-letter codes which indicates genotype (P = PI 587982A, S = S99-11986), location seed were produced (C = CSPS, E = ESPS), and tissue (M = dry, mature seed; I = 6-h imbibed seed; G = germinated seed)

Seed of the PI line were found to germinate much more rapidly than those of the SG line in both environments (Fig. 1b), and PI seed from both unstressed and heat-stressed locations germinated with very high efficiency (> 80%, Fig. 1a). In contrast, only 75% of CSPS-produced seed from SG germinated by the end of 72 h. A dramatic reduction in germination was noted for SG seed produced under the heat-stress of the ESPS (~ 30% germination at 72 h, Fig. 1b). Our germination results are concordant with our previous metabolic study [21].

We then selected three stages (Fig. 1b, Table 1) to obtain transcriptomic data: (1) mature, dry seed; (2) 6-h imbibed seed; and (3) germinated seed with emerged radicle for each genotype grown in both environments (Fig. 1b, Table 1). It is important to note that the time from imbibition to germination varied between genotype/environments (Fig. 1b). Three biological replicates (each consisting of 5 seed) per genotype/condition/timepoint were used for analysis to quantify gene expression. The number of genes expressed (FPKM > 0.3) in each sample ranged from 23,560 to 30,349 (Fig. 1c, Table 1).

Table 1 RNAseq sample details

A core set of genes expressed was identified: (A) 21,082 in all mature seed tissues; (B) 26,372 genes expressed in all 6 h imbibed seed tissues; and (C) 21,843 genes in all germinated seed tissues (Fig. 2a, Additional file 3).

Fig. 2

Venn diagrams of genes expression and GO term enrichment of soybean seeds during germination in response to genotype and environment. a Venn diagram showing total numbers of genes expressed in mature, 6-h imbibed or germinated soybean seeds. Two-letter combination indicates genotype (P = PI 587982A, S = S99-11986) and environment (C = CSPS, E = ESPS). b Venn diagram indicating overlap of gene ontology (GO) terms in upregulated differential genes; upregulation indicates higher level of expression in the first code listed

Differential expressed gene analysis

We utilized a Tuxedo RNAseq analysis pipeline to make 20 distinct comparisons, which can be divided into four general categories: (1) environmental effects; (2) genotypic effects; (3) the transition between mature seeds to 6-h imbibed seeds; and (4) the transition between imbibed seeds to germinated seeds (Table 1, Additional file 4).

An average of 7385 differentially expressed genes (DEGs) were detected between environments (threshold for all comparisons was q-value < 0.05). An average of 7789 DEGs were detected between genotypes. An average of 11,833 DEGs were detected between mature and 6-h imbibed seeds, across genotypes and environments (Fig. 1d). Lastly, an average of 13,344 DEGs were detected between imbibed and germinated seeds (Fig. 1d, full details are presented in Additional file 4).

Gene ontology enrichment

We utilized a gene ontology (GO) term enrichment tool ( to examine lists of differentially expressed genes for the 20 comparisons (Figs. 1d, 2b, c, Additional file 5). A larger number of DEGs were found in comparisons of developmental transitions (from mature-to-imbibed or imbibed-to-germinated seeds) than within the same developmental stage (either between environments or between genotypes). Seed development/maturation under the high temperature conditions of the ESPS (Figs. 1d, 2b, Additional file 5) was associated with significant enrichment (at a threshold p-value < 0.05) for gene annotations for heat stress, response to oxidative stress and protein folding; “response to hydrogen peroxide”, “response to high light intensity”, and “response to heat” were both overrepresented and GO-term enriched in all 6 comparisons for environmental effects, with “response to wounding” overrepresented and GO term enriched in 4/6 comparisons (excepting comparisons #2, #4; comparison numbers specified in Additional file 5). The enrichment of numerous GO terms associated with abiotic stress response gives a clear indication that the mRNA pools of both genotypes are responsive to the higher temperatures of the ESPS as compared to the less stressful CSPS.

Despite this environmental response, seed mRNA pools of the heat-tolerant line were further enriched (Fig. 2b, Additional file 5) for genes with GO-terms associated with abiotic stress response [e.g. “response to high light intensity”, “response to hydrogen peroxide” and “response to heat” in 5/6 comparisons (excepting #9)” and “response to water deprivation” in 4/6 (excepting #7, #9); “response to cadmium ion” in 4/6 comparisons (excepting #9, #11); “response to salt stress” in 4/6 comparisons (excepting #9, #10)]. In addition we observed enrichment in the tolerant PI mRNA pools for protein refolding-associated GO terms: “nucleosome assembly” in 4/6 comparisons (excepting #7, #11) and “response to endoplasmic reticulum stress. The GO term “l-ascorbic acid biosynthesis” was also observed to be enriched in seed of the stress tolerant PI under the ESPS; these results are concordant our previous metabolomics study [21], which conclusively demonstrated that higher levels of ascorbate precursors were found in seeds of a heat-tolerant soybean line. Collectively, these results suggest fundamental differences exist between seed mRNA pools between the two genotypes; the more stress tolerant PI genotype is effectively “genetically primed” to more effectively manage abiotic stress as well as for higher levels of seed antioxidant compounds. This mRNA priming trend persists through seed germination and ultimately biologically translates to more efficient and effective seed germination (and in field conditions seedling emergence).

Weighted gene co-expression network analysis

Weighted gene co-expression network analysis (WGCNA) is a systems biology method for describing the correlation patterns among genes across samples [16]. We utilized the WGCNA package in R on FPKM data of all samples to find modules (clusters) of highly correlated expressed differential genes (≥ twofold) and a total 16 clusters were detected (Additional files 6 and 7; clusters are color coded).

Co-expressed gene clusters were then examined for overabundance of GO Biological Process terms (Additional file 8, GO Molecular Function are also provided in Additional file 9). For brevity, only Biological Process results will be discussed here. For 6/16 GO:BP clusters no significantly enriched GO terms were found (salmon, pink, purple, midnightblue, magenta, black). Several gene clusters were enriched for gene expression/chromatin remodeling (Yellow, GreenYellow, Brown) for translation/ribosome components (Yellow, Blue), mRNA splicing (Cyan, Brown) Actin/Cytoskeleton (Red, Grey) for Cell wall/Carbohydrate metabolism (Turquoise, Red). Of particular interest is the Green gene expression module, which displayed enrichment for numerous GO:BP terms annotated for abiotic stress responses (e.g. response to temperature stimulus, response to reactive oxygen species, response to salt stress, response to heat, response to stress, etc.).

Validation of RNAseq data using qRT-PCR

Four differentially expressed genes (Additional file 10) were selected for qRT-PCR validation of the RNAseq data. Two genes were highly expressed (KTI-1, average 86 FPKM; HSP20, average 856 FPKM) and two were lower expressed genes (SAM-methyltransferase, average 31 FPKM; UDP-glycosyl transferase, average 4.0 FPKM). qRT-PCR were tested via the ΔΔCt method and expressed as log2 ratios (Additional file 11). Correlations between qRT-PCR and FPKM results were robust for mature (r2 = 0.9729) and imbibed (r2 = 0.9919), but less robust for germinated samples (r2 = 0.6844). The high concordance between RNAseq and qRT-PCR highlights the high quality of our RNAseq dataset.


In this study we provide substantial new mRNA sequencing data that defines the very early stages of soybean seed germination (mature seed > imbibed seed > germinated seed). We also contrasted two genotypes which differ in terms of tolerance to high temperature stress during seed development, which were produced under two distinct temperature stress field locations. We demonstrate that the more temperature-tolerant PI genotype is primed at the mRNA level to handle higher levels of temperature stress. In addition, we demonstrate that the PI line has faster, more efficient and more effective seed germination regardless of seed production location/environmental stress. These results highlight some of the genetic gains possible by leveraging exotic soybean germplasm as sources of novel traits in soybean breeding programs.


The experiment mandated a need to visually rate seeds (exposure to light) during germination on prewetted filter paper. Therefore, the transcriptomes may not completely reflect how germination of seeds in soil would proceed.

We observed poor clustering of RNAseq data for germinated seeds of the PI produced in the ESPS with other samples (PEG, Additional file 2), which is most evident in the large number of significant DEGs detected (Additional file 4).

Availability of data and materials

All sequence data obtained have been deposited in the NCBI Sequence Read Archive under project SRP090036. Analyzed datasets have also been uploaded to the SoyKB community resource ( [22,23,24] and is freely available to all researchers for visualization and interactive data analysis purposes, within the “Differential Expression Suite of Tools” and gene card pages in SoyKB. PI 587982A and S99-11986 are available from the USDA-GRIN germplasm repository ( Seed from public germplasm release DS25-1, which has heat tolerance from PI 587982A in an agronomically and yield improved line, is available by contacting Dr. Rusty Smith (Rusty.Smith@ARS.USDA.GOV). DS25-1 is also available from the USDA-GRIN germplasm repository where and can be located under the identifier PI 684675.



conventional soybean production system


differentially expressed genes


early soybean production system


fragment per kilobase million


gene ontology


maturity group


messenger ribonucleic acid


plant introduction


quantitative real time polymerase chain reaction


ribonucleic acid sequencing


susceptible genotype S99-11986


tolerant genotype PI587982A


United States Department of Agriculture-Germplasm Resources Information Network


weighted co-expressed gene network analysis


  1. 1.

    Heatherly LG, Spurlock SR. Yield and economics of traditional and early soybean production system (ESPS) seedings in the midsouthern United States. Field Crops Res. 1999;63(1):35–45.

    Article  Google Scholar 

  2. 2.

    Heatherly LG. Early soybean production system (ESPS). In: Heatherly LG, Hodges HF, editors. Soybean production system in the midsouth. Boca Raton: CRC Press; 1999. p. 103–15.

    Google Scholar 

  3. 3.

    Heatherly LG. Yield and germinability of seed from irrigated and nonirrigated early- and late-planted MG IV and V soybean. Crop Sci. 1996;36(4):1000–6.

    Article  Google Scholar 

  4. 4.

    Hesketh JD, Myhre DL, Willey CR. Temperature control of time intervals between vegetative and reproductive events in soybeans. Crop Sci. 1973;13(2):250–4.

    Article  Google Scholar 

  5. 5.

    Hatfield JL, Boote KJ, Kimball BA, Ziska LH, Izaurralde RC, Ort D, Thomson AM, Wolfe D. Climate impacts on agriculture: implications for crop production. Agron J. 2011;103(2):351–70.

    Article  Google Scholar 

  6. 6.

    Smith JR, Mengistu A, Nelson RL, Paris RL. Identification of soybean accessions with high germinability in high-temperature environments. Crop Sci. 2008;48(6):2279–88.

    Article  Google Scholar 

  7. 7.

    Mengistu A, Heatherly LG. Planting date, irrigation, maturity group, year, and environment effects on Phomopsis longicolla, seed germination, and seed health rating of soybean in the early soybean production system of the midsouthern USA. Crop Prot. 2006;25(4):310–7.

    Article  Google Scholar 

  8. 8.

    Hyten DL, Song Q, Zhu Y, Choi IY, Nelson RL, Costa JM, Specht JE, Shoemaker RC, Cregan PB. Impacts of genetic bottlenecks on soybean genome diversity. Proc Natl Acad Sci USA. 2006;103(45):16666–71.

    CAS  Article  Google Scholar 

  9. 9.

    Gizlice Z, Carter TE, Burton JW. Genetic base for north american public soybean cultivars released between 1947 and 1988. Crop Sci. 1994;34(5):1143–51.

    Article  Google Scholar 

  10. 10.

    Smith JR. Soybean germplasm line DS25-1 with heat tolerance and competitive yield under heat stress. Columbia: USDA-ARS; 2017.

    Google Scholar 

  11. 11.

    Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, He J, Xu D, May G, Stacey G. An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010;63(1):86–99.

    CAS  PubMed  Google Scholar 

  12. 12.

    Severin AJ, Woody JL, Bolon Y-T, Joseph B, Diers BW, Farmer AD, Muehlbauer GJ, Nelson RT, Grant D, Specht JE, Graham MA, Cannon SB, May GD, Vance CP, Shoemaker RC. RNA-Seq atlas of glycine max: a guide to the soybean transcriptome. BMC Plant Biol. 2010;10(1):160.

    Article  Google Scholar 

  13. 13.

    Woody JL, Severin AJ, Bolon Y-T, Joseph B, Diers BW, Farmer AD, Weeks N, Muehlbauer GJ, Nelson RT, Grant D, Specht JE, Graham MA, Cannon SB, May GD, Vance CP, Shoemaker RC. Gene expression patterns are correlated with genomic and genic structure in soybean. Genome. 2010;54(1):10–8.

    Article  Google Scholar 

  14. 14.

    Komatsu S, Yamamoto R, Nanjo Y, Mikami Y, Yunokawa H, Sakata K. A comprehensive analysis of the soybean genes and proteins expressed under flooding stress using transcriptome and proteome techniques. J Proteome Res. 2009;8(10):4766–78.

    CAS  Article  Google Scholar 

  15. 15.

    Le DT, Nishiyama R, Watanabe Y, Tanaka M, Seki M, Ham LH, Yamaguchi-Shinozaki K, Shinozaki K. Tran L-SP: differential gene expression in soybean leaf tissues at late developmental stages under drought stress revealed by genome-wide transcriptome analysis. PLoS ONE. 2012;7(11):e49522.

    CAS  Article  Google Scholar 

  16. 16.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9(1):559.

    Article  Google Scholar 

  17. 17.

    Gillman JD, Kim W-S, Krishnan HB. Identification of a new soybean Kunitz trypsin inhibitor mutation and its effect on Bowman–Birk protease inhibitor content in soybean seed. J Agric Food Chem. 2015;63(5):1352–9.

    CAS  Article  Google Scholar 

  18. 18.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    CAS  Article  Google Scholar 

  19. 19.

    Libault M, Thibivilliers S, Bilgin DD, Radwan O, Benitez M, Clough SJ, Stacey G. Identification of four soybean reference genes for gene expression normalization. Plant Genome. 2008;1:44–54.

    CAS  Article  Google Scholar 

  20. 20.

    Shannon JG, Nelson RL, Wrather JA. Registration of S99-11509 and S99-11986 improved soybean germplasm with diverse pedigree registration by CSSA. Crop Sci. 2005;45(4):1672–3.

    Article  Google Scholar 

  21. 21.

    Chebrolu KK, Fritschi FB, Ye S, Krishnan HB, Smith JR, Gillman JD. Impact of heat stress during seed development on soybean seed metabolome. Metabolomics. 2016;12(2):28.

    Article  Google Scholar 

  22. 22.

    Joshi T, Patil K, Fitzpatrick MR, Franklin LD, Yao Q, Cook JR, Wang Z, Libault M, Brechenmacher L, Valliyodan B, Wu X, Cheng J, Stacey G, Nguyen HT, Xu D. Soybean knowledge base (SoyKB): a web resource for soybean translational genomics. BMC Genom. 2012;13(Suppl 1):S15.

    CAS  Article  Google Scholar 

  23. 23.

    Joshi T, Wang J, Zhang H, Chen S, Zeng S, Xu B, Xu D. The evolution of soybean knowledge base (SoyKB). Methods Mol Biol. 2017;1533:149–59.

    CAS  Article  Google Scholar 

  24. 24.

    Joshi T, Fitzpatrick MR, Chen S, Liu Y, Zhang H, Endacott RZ, Gaudiello EC, Stacey G, Nguyen HT, Xu D. Soybean knowledge base (SoyKB): a web resource for integration of soybean translational genomics and molecular breeding. Nucleic Acids Res. 2014;42(Database issue):D1245–52.

    CAS  Article  Google Scholar 

Download references


Not applicable.


Funds for sequencing and bioinformatics analysis (to WS, SG) was drawn from internal USDA-ARS project funds (5070-21000-038-00D). The salary for SY was from a United Soybean Board Grant (1420-532-5613), and salary for JB from a Missouri Soybean Merchandising Council Grant (#14-359). Funding agencies provided financial support to certain authors, but played no direct role in study design, data collection/analysis/interpretation nor in writing the manuscript

Author information




JDG conceived of the study, performed germination experiments, isolated RNA, prepared libraries, and main-authored the paper. SY and WS carried out bioinformatics analyses under the supervision of SG. JB performed statistical analyses, and created the first draft of the manuscript. JRS was involved in the conception of the study, provided all soybean lines, and co-authored the manuscript. FB was involved in conception of the study, supervision of SY and JB and in co-authoring the manuscript. ZL performed the WGCNA gene module analysis. TJ was responsible for bioinformatics analyses related to the SoyKB ( data depository and informatics tools hosted therein. All authors have reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jason D. Gillman.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1.

Additional field, sample, RNA sequencing and mapping methods.

Additional file 2.

Hierarchical clustering and heatmap of differentially expressed genes.

Additional file 3.

Cufflinks fragment per kilobase million (FPKM) results for all samples.

Additional file 4.

Cuffdiff differentially expressed gene results.

Additional file 5.

Gene ontology (GO) enrichment results for 20 DEG comparisons.

Additional file 6.

Heatmap of genes analyzed by weighted gene co-expression network analysis (WGCNA, ≥ twofold difference on a log2 scale).

Additional file 7.

Weighted gene co-expression network analysis (WGCNA) for gene function for differentially expressed genes (≥ twofold difference on a log2 scale).

Additional file 8.

Gene ontology biological process term enrichment analysis of modules identified by weighted gene co-expression network analysis (WGCNA).

Additional file 9.

Gene ontology molecular function term enrichment analysis of modules identified by weighted gene co-expression network analysis (WGCNA).

Additional file 10.

Primer sequences used for qRT-PCR.

Additional file 11.

Figure displaying correlation of qRT-PCR and RNAseq data.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gillman, J.D., Biever, J.J., Ye, S. et al. A seed germination transcriptomic study contrasting two soybean genotypes that differ in terms of their tolerance to the deleterious impacts of elevated temperatures during seed fill. BMC Res Notes 12, 522 (2019).

Download citation


  • Soybean
  • Transcriptomic analysis
  • Seed germination
  • Temperature stress
  • Seed development