Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

SET domain group encoding proteins function as histone lysine methyltransferases. These proteins are involved in various biological processes, including plant development and adaption to the environment by modifying the chromatin structures. So far, the SET domain genes (SDGs) have not been systematically investigated in Brassica napus (B. napus). In the current study, through genome-wide analysis, a total of 122 SDGs were identified in the B. napus genome. These BnSDGs were subdivided into seven (I–VII) classes based on phylogeny analysis, domain configurations, and motif distribution. Segmental duplication was involved in the evolution of this family, and the duplicated genes were under strong purifying selection. The promoter sequence of BnSDGs consisted of various growth, hormones, and stress-related cis-acting elements along with transcription factor binding sites (TFBSs) for 20 TF families in 59 of the 122 BnSDGs. The gene ontology (GO) analysis revealed that BnSDGs were closely associated with histone and non-histone methylation and metal binding capacity localized mostly in the nucleus. The in silico expression analysis at four developmental stages in leaf, stem root, floral organ, silique, and seed tissues showed a broad range of tissue and stage-specific expression pattern. The expression analysis under four abiotic stresses (dehydration, cold, ABA, and salinity) also provided evidence for the importance of BnSDGs in stress environments. Based on expression analysis, we performed reverse transcription-quantitative PCR for 15 target BnSDGs in eight tissues (young leaf, mature leaf, root, stem, carpel, stamen, sepal, and petals). Our results were in accordance with the in silico expression data, suggesting the importance of these genes in plant development. In conclusion, this study lays a foundation for future functional studies on SDGs in B. napus.

In plants, phloem is not only a vital structure that is used for nutrient transportation, but it is also the location of a response that defends against various stresses, named phloem-based defense (PBD). Phloem proteins (PP2s) are among the predominant proteins in phloem, indicating their potential functional role in PBD. Sclerotinia disease (SD), which is caused by the necrotrophic fungal pathogen S. sclerotiorum (Sclerotinia sclerotiorum), is a devastating disease that affects oil crops, especially Brassica napus (B. napus), mainly by blocking nutrition and water transportation through xylem and phloem. Presently, the role of PP2s in SD resistance is still largely estimated. Therefore, in this study, we identified 62 members of the PP2 gene family in the B. napus genome with an uneven distribution across the 19 chromosomes. A phylogenetic analysis classified the BnPP2s into four clusters (I–IV), with cluster I containing the most members (28 genes) as a consequence of its frequent genome segmental duplication. A comparison of the gene structures and conserved motifs suggested that BnPP2 genes were well conserved in clusters II to IV, but were variable in cluster I. Interestingly, the motifs in different clusters displayed unique features, such as motif 6 specifically existing in cluster III and motif 1 being excluded from cluster IV. These results indicated the possible functional specification of BnPP2s. A transcriptome data analysis showed that the genes in clusters II to IV exhibited dynamic expression alternation in tissues and the stimulation of S. sclerotiorum, suggesting that they could participate in SD resistance. A GWAS analysis of a rapeseed population comprising 324 accessions identified four BnPP2 genes that were potentially responsible for SD resistance and a transgenic study that was conducted by transiently expressing BnPP2-6 in tobacco (Nicotiana tabacum) leaves validated their positive role in regulating SD resistance in terms of reduced lesion size after inoculation with S. sclerotiorum hyphal plugs. This study provides useful information on PP2 gene functions in B. napus and could aid elaborated functional studies on the PP2 gene family.

The trans-Golgi network (TGN) acts as a central platform for sorting and secreting various cargoes to the cell surface, thus being essential for the full execution of plant immunity. However, the fine-tuned regulation of TGN components in plant defense and stress response has been not fully elucidated. Our study revealed that despite largely compromising penetration resistance, the loss-of-function mutation of the TGN component protein ECHIDNA (ECH) induced enhanced postinvasion resistance to powdery mildew in Arabidopsis thaliana. Genetic and transcriptome analyses and hormone profiling demonstrated that ECH loss resulted in salicylic acid (SA) hyperaccumulation via the ISOCHORISMATE SYNTHASE 1 biosynthesis pathway, thereby constitutively activating SA-dependent innate immunity that was largely responsible for the enhanced postinvasion resistance. Furthermore, the ech mutant displayed accelerated SA-independent spontaneous cell death and constitutive POWDERY MILDEW RESISTANCE 4-mediated callose depositions. In addition, ECH loss led to a chronically prolonged endoplasmic reticulum stress in the ech mutant. These results provide insights into understanding the role of TGN components in the regulation of plant immunity and stress responses.

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.