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.

Chromosome-Level Genome Assembly and Annotation of the Fiber Flax (Linum usitatissimum) Genome

Background: Drought severely limits sunflower production especially at the seedling stage. To investigate the response mechanism of sunflowers to drought stress, we utilized two genotypes of sunflower materials with different drought resistances as test materials. The physiological responses were investigated under well-watered (0 h) and drought-stressed conditions (24 h, 48 h, and 72 h). Results: ANOVA revealed the greatest differences in physiological indices between 72 h of drought stress and 0 h of drought stress. Transcriptome analysis was performed after 72 h of drought stress. At 0 h, there were 7482 and 5627 differentially expressed genes (DEGs) in the leaves of K55 and K58, respectively, and 2150 and 2527 DEGs in the roots of K55 and K58, respectively. A total of 870 transcription factors (TFs) were identified among theDEGs, among which the high-abundance TF families included AP2/ERF, MYB, bHLH,and WRKY. Five modules were screened using weighted gene coexpressionnetwork analysis (WGCNA), three and two of which were positively and negatively, respectively, related to physiological traits. KEGG analysis revealedthat under drought stress, “photosynthesis”, “carotenoid biosynthesis”, “starch and sucrose metabolism”, “ribosome”, “carotenoid biosynthesis”, “starch and sucrose metabolism”, “protein phosphorylation” and “phytohormone signaling” are six important metabolic pathways involved in the response of sunflower to drought stress. Cytoscape software was used to visualize the three key modules, and the hub genes were screened. Finally, a total of 99 important candidate genes that may be associated with the drought response in sunflower plants were obtained, and the homology of these genes was compared with that in Arabidopsis thaliana. Conclusions: Taken together, our findings could lead to a better understanding of drought tolerance in sunflowers and facilitate the selection of drought-tolerant sunflower varieties.

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for “Neiya No. 9”, a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

The MYB transcription factor family has numerous members, and is involved in biological activities, such as ABA signaling, which plays an important role in a plant’s resistance to abiotic stresses such as drought. However, the diversity of MYB members that respond to drought stress and their regulatory mechanisms in different flax varieties were unclear. In this study, we obtained 855.69 Gb of clean data from 120 flax root samples from 20 flax (Linum usitatissimum L.) varieties, assembled 92,861 transcripts, and identified 434 MYB family members in each variety. The expression profiles of the MYB transcription factor family from 20 flax varieties under drought stress were analyzed. The results indicated that there are four strategies by which the MYB family responds to drought stress in these 20 flax varieties, each of which has its own specific processes, such as development, reproduction, and localization processes. The four strategies also include common biological processes, such as stimulus responses, metabolic processes, and biological regulation. The WGCNA method was subsequently employed to identify key members of the MYB family involved in response strategies to drought stress. The results demonstrated that a 1R-MYB subfamily gene co-expression network is significantly related to the gibberellin response and cytokinin-activated signaling pathway processes in the ‘Strategy 4’ for MYB family response to drought, identifying core genes such as Lus.scaffold70.240. Our results showed a diversity of MYB family responses to drought stress within flax varieties, and these results contribute to deciphering the mechanisms of the MYB family regulation of drought resistance. This will promote the more accurate breeding development of flax to adapt to agricultural production under drought conditions.