Flax is an important crop for oil and fiber, however, no high-density genetic maps have been reported for this species. Specific length amplified fragment sequencing (SLAF-seq) is a high-resolution strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms. In this study, SLAF-seq was employed to develop SNP markers in an F2 population to construct a high-density genetic map for flax. In total, 196.29 million paired-end reads were obtained. The average sequencing depth was 25.08 in male parent, 32.17 in the female parent, and 9.64 in each F2 progeny. In total, 389,288 polymorphic SLAFs were detected, from which 260,380 polymorphic SNPs were developed. After filtering, 4,638 SNPs were found suitable for genetic map construction. The final genetic map included 4,145 SNP markers on 15 linkage groups and was 2,632.94 cM in length, with an average distance of 0.64 cM between adjacent markers. To our knowledge, this map is the densest SNP-based genetic map for flax. The SNP markers and genetic map reported in here will serve as a foundation for the fine mapping of quantitative trait loci (QTLs), map-based gene cloning and marker assisted selection (MAS) for flax.

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−6). 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.

No abstract

The yellow nutsedge (Cyperus esculentus L. 1753) is an unconventional oil plant with oil-rich tubers, and a potential alternative for traditional oil crops. Here, we reported the first high-quality and chromosome-level genome assembly of the yellow nutsedge generated by combining PacBio HiFi long reads, Novaseq short reads, and Hi-C data. The final genome size is 225.6 Mb with an N50 of 4.3 Mb. More than 222.9 Mb scaffolds were anchored to 54 pseudochromosomes with a BUSCO score of 96.0%. We identified 76.5 Mb (33.9%) repetitive sequences across the genome. A total of 23,613 protein-coding genes were predicted in this genome, of which 22,847 (96.8%) were functionally annotated. A whole-genome duplication event was found after the divergence of Carex littledalei and Rhynchospora breviuscula, indicating the rich genetic resources of this species for adaptive evolution. Several significantly enriched GO terms were related to invasiveness of the yellow nutsedge, which may explain its plastic adaptability. In addition, several enriched Kyoto Encyclopedia of Genes and Genomes pathways and expanded gene families were closely related with substances in tubers, partially explaining the genomic basis of characteristics of this oil-rich tuber.