Drought is one of the most serious abiotic stress factors limiting crop yields. Although sunflower is considered a moderate drought-tolerant plant, drought stress still has a negative impact on sunflower yield as cultivation expands into arid regions. The extent of drought stress is varieties and time-dependent, however, the molecular response mechanisms of drought tolerance in sunflower with different varieties are still unclear. Here, we performed comparative physiological and transcriptome analyses on two sunflower inbred lines with different drought tolerance at the seedling stage. The analysis of nine physiological and biochemical indicators showed that the leaf surface area, leaf relative water content, and cell membrane integrity of drought tolerance inbred line were higher than those of drought-sensitive inbred line under drought stress, indicating that DT had stronger drought resistance. Transcriptome analyses identified 24,234 differentially expressed genes (DEGs). Gene ontology (GO) analysis showed the up-regulated genes were mainly enriched in gibberellin metabolism and rRNA processing, while the down-regulated genes were mainly enriched in cell-wall, photosynthesis, and terpene metabolism. Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway analysis showed genes related to GABAergic synapse, ribosome biogenesis were up-regulated, while genes related with amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, photosynthesis were down-regulated. Mapman analysis revealed differences in plant hormone-signaling genes over time and between samples. A total of 1,311 unique putative transcription factors (TFs) were identified from all DEGs by iTAK, among which the high abundance of transcription factor families include bHLH, AP2/ERF, MYB, C2H2, etc. Weighted gene co-expression network analysis (WGCNA) revealed a total of 2,251 genes belonging to two modules(blue 4, lightslateblue), respectively, which were significantly associated with six traits. GO and KEGG enrichment analysis of these genes was performed, followed by visualization with Cytoscape software, and the top 20 Hub genes were screened using the CytoHubba plugin.

Sunflower is one of the most important oil crops in the world, and drought stress can severely limit its production and quality. To understand the underlying mechanism of drought tolerance, and identify candidate genes for drought tolerance breeding, we conducted a combined genome-wide association studies (GWAS) and RNA-seq analysis. A total of 226 sunflower inbred lines were collected from different regions of China and other countries. Eight phenotypic traits were evaluated under control and drought stress conditions. Genotyping was performed using a Specific-Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 934.08 M paired-end reads were generated, with an average Q30 of 91.97%. Based on the 243,291 polymorphic SLAF tags, a total of 94,162 high-quality SNPs were identified. Subsequent analysis of linkage disequilibrium (LD) and population structure in the 226 accessions was carried out based on the 94,162 high-quality SNPs. The average LD decay across the genome was 20 kb. Admixture analysis indicated that the entire population most likely originated from 11 ancestors. GWAS was performed using three methods (MLM, FarmCPU, and BLINK) simultaneously. A total of 80 SNPs showed significant associations with the 8 traits (p < 1.062 × 10⁻⁶). Next, a total of 118 candidate genes were found. To obtain more reliable candidate genes, RNA-seq analysis was subsequently performed. An inbred line with the highest drought tolerance was selected according to phenotypic traits. RNA was extracted from leaves at 0, 7, and 14 days of drought treatment. A total of 18,922 differentially expressed genes were obtained. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed up-regulated genes were mainly enriched in the branched-chain amino acid catabolic process, while the down-regulated genes were mainly enriched in the photosynthesis-related process. Six DEGs were randomly selected from all DEGs for validation; these genes showed similar patterns in RNA-seq and RT-qPCR analysis, with a correlation coefficient of 0.8167. Through the integration of the genome-wide association study and the RNA-sequencing, 14 candidate genes were identified. Four of them (LOC110885273, LOC110872899, LOC110891369, LOC110920644) were abscisic acid related protein kinases and transcription factors. These genes may play an important role in sunflower drought response and will be used for further study. Our findings provide new insights into the response mechanisms of sunflowers against drought stress and contribute to further genetic breeding.

Sunflower is an important oilseed crop across the world. It is considered as a moderately drought tolerant plant, however, its yield is still negatively affected by drought stress. Improving drought tolerance is of the outmost important for breeding. Although several studies have documented the relationship between the sunflower phenotype and genotype under drought stress, but relatively few studies have simultaneously investigated the molecular mechanisms of drought tolerance in the sunflower at different growth stages. In this study, we conducted quantitative trait locus (QTL) analysis for different sunflower traits during the germination and seedling stages. Eighteen phenotypic traits were evaluated under well-watered and drought stress conditions. We determined that the germination rate, germination potential, germination index, and root-to-shoot ratio can be used as effective indexes for drought tolerance selection and breeding. A total of 33 QTLs were identified on eight chromosomes (PVE: 0.016%–10.712% with LOD: 2.017–7.439). Within the confidence interval of the QTL, we identified 60 putative drought-related genes. Four genes located on chromosome 13 may function in both germination and seedling stages for drought response. Genes LOC110898128, LOC110898092, LOC110898071, and LOC110898072 were annotated as aquaporin SIP1-2-like, cytochrome P450 94C1, GABA transporter 1-like, and GABA transporter 1-like isoform X2, respectively. These genes will be used for further functional validation. This study provides insight into the molecular mechanisms of the sunflower’s in response to drought stress. At the same time, it lays a foundation for sunflower drought tolerance breeding and genetic improvement.