Hello darkness, my old friend: 3-KETOACYL-COENZYME A SYNTHASE4 is a branch point in the regulation of triacylglycerol synthesis in Arabidopsis thaliana

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

A genome-wide association study identifies genes associated with cuticular wax metabolism in maize

The plant cuticle is essential in plant defense against biotic and abiotic stresses. To systematically elucidate the genetic architecture of maize (Zea mays L.) cuticular wax metabolism, 2 cuticular wax-related traits, the chlorophyll extraction rate (CER) and water loss rate (WLR) of 389 maize inbred lines, were investigated and a genome-wide association study (GWAS) was performed using 1.25 million single nucleotide polymorphisms (SNPs). In total, 57 nonredundant quantitative trait loci (QTL) explaining 5.57% to 15.07% of the phenotypic variation for each QTL were identified. These QTLs contained 183 genes, among which 21 strong candidates were identified based on functional annotations and previous publications. Remarkably, 3 candidate genes that express differentially during cuticle development encode β-ketoacyl-CoA synthase (KCS). While ZmKCS19 was known to be involved in cuticle wax metabolism, ZmKCS12 and ZmKCS3 functions were not reported. The association between ZmKCS12 and WLR was confirmed by resequencing 106 inbred lines, and the variation of WLR was significant between different haplotypes of ZmKCS12. In this study, the loss-of-function mutant of ZmKCS12 exhibited wrinkled leaf morphology, altered wax crystal morphology, and decreased C32 wax monomer levels, causing an increased WLR and sensitivity to drought. These results confirm that ZmKCS12 plays a vital role in maize C32 wax monomer synthesis and is critical for drought tolerance. In sum, through GWAS of 2 cuticular wax-associated traits, this study reveals comprehensively the genetic architecture in maize cuticular wax metabolism and provides a valuable reference for the genetic improvement of stress tolerance in maize.

Tetracosanoic acids produced by 3-ketoacyl-CoA synthase 17 are required for synthesizing seed coat suberin in Arabidopsis

Very long-chain fatty acids (VLCFAs) are precursors for the synthesis of membrane lipids, cuticular waxes, suberins, and storage oils in plants. 3-Ketoacyl CoA synthase (KCS) catalyzes the condensation of C2 units from malonyl-CoA to acyl-CoA, the first rate-limiting step in VLCFA synthesis. In this study, we revealed that Arabidopsis KCS17 catalyzes the elongation of C22-C24 VLCFAs required for synthesizing seed coat suberin. Histochemical analysis of Arabidopsis plants expressing GUS (β-glucuronidase) under the control of the KCS17 promoter revealed predominant GUS expression in seed coats, petals, stigma, and developing pollen. The expression of KCS17:eYFP (enhanced yellow fluorescent protein) driven by the KCS17 promoter was observed in the outer integument1 of Arabidopsis seed coats. The KCS17:eYFP signal was detected in the endoplasmic reticulum of tobacco epidermal cells. The levels of C22 VLCFAs and their derivatives, primary alcohols, α,ω-alkane diols, ω-hydroxy fatty acids, and α,ω-dicarboxylic acids increased by ~2-fold, but those of C24 VLCFAs, ω-hydroxy fatty acids, and α,ω-dicarboxylic acids were reduced by half in kcs17-1 and kcs17-2 seed coats relative to the wild type (WT). The seed coat of kcs17 displayed decreased autofluorescence under UV and increased permeability to tetrazolium salt compared with the WT. Seed germination and seedling establishment of kcs17 were more delayed by salt and osmotic stress treatments than the WT. KCS17 formed homo- and hetero-interactions with KCR1, PAS2, and ECR, but not with PAS1. Therefore, KCS17-mediated VLCFA synthesis is required for suberin layer formation in Arabidopsis seed coats.

Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

Genome-wide identification and expression analysis of the KCS gene family in soybean (Glycine max) reveal their potential roles in response to abiotic stress

Very long chain fatty acids (VLCFAs) are fatty acids with chain lengths of 20 or more carbon atoms, which are the building blocks of various lipids that regulate developmental processes and plant stress responses. 3-ketoacyl-CoA synthase encoded by the KCS gene is the key rate-limiting enzyme in VLCFA biosynthesis, but the KCS gene family in soybean (Glycine max) has not been adequately studied thus far. In this study, 31 KCS genes (namely GmKCS1 - GmKCS31) were identified in the soybean genome, which are unevenly distributed on 14 chromosomes. These GmKCS genes could be phylogenetically classified into seven groups. A total of 27 paralogous GmKCS gene pairs were identified with their Ka/Ks ratios indicating that they had undergone purifying selection during soybean genome expansion. Cis-acting element analysis revealed that GmKCS promoters contained multiple hormone- and stress-responsive elements, indicating that GmKCS gene expression levels may be regulated by various developmental and environmental stimuli. Expression profiles derived from RNA-seq data and qRT-PCR experiments indicated that GmKCS genes were diversely expressed in different organs/tissues, and many GmKCS genes were found to be differentially expressed in the leaves under cold, heat, salt, and drought stresses, suggesting their critical role in soybean resistance to abiotic stress. These results provide fundamental information about the soybean KCS genes and will aid in their further functional elucidation and exploitation.

Genome-wide identification and expression analysis of 3-ketoacyl-CoA synthase gene family in rice (Oryza sativa L.) under cadmium stress

3-Ketoacyl-CoA synthase (KCS) is the key rate-limiting enzyme for the synthesis of very long-chain fatty acids (VLCFAs) in plants, which determines the carbon chain length of VLCFAs. However, a comprehensive study of KCSs in Oryza sativa has not been reported yet. In this study, we identified 22 OsKCS genes in rice, which are unevenly distributed on nine chromosomes. The OsKCS gene family is divided into six subclasses. Many cis-acting elements related to plant growth, light, hormone, and stress response were enriched in the promoters of OsKCS genes. Gene duplication played a crucial role in the expansion of the OsKCS gene family and underwent a strong purifying selection. Quantitative Real-time polymerase chain reaction (qRT-PCR) results revealed that most KCS genes are constitutively expressed. We also revealed that KCS genes responded differently to exogenous cadmium stress in japonica and indica background, and the KCS genes with higher expression in leaves and seeds may have functions under cadmium stress. This study provides a basis for further understanding the functions of KCS genes and the biosynthesis of VLCFA in rice.