Comparative sequence analysis across Brassicaceae, regulatory diversity in KCS5 and KCS6 homologs from Arabidopsis thaliana and Brassica juncea, and intronic fragment as a negative transcriptional regulator

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

BrWAX3, Encoding a β-ketoacyl-CoA Synthase, Plays an Essential Role in Cuticular Wax Biosynthesis in Chinese Cabbage

In this study, we identified a novel glossy mutant from Chinese cabbage, named SD369, and all wax monomers longer than 26 carbons were significantly decreased. Inheritance analysis revealed that the glossy trait of SD369 was controlled by a single recessive locus, BrWAX3. We fine-mapped the BrWAX3 locus to an interval of 161.82 kb on chromosome A09. According to the annotated genome of Brassica rapa, Bra024749 (BrCER60.A09), encoding a β-ketoacyl-CoA synthase, was identified as the candidate gene. Expression analysis showed that BrCER60.A09 was significantly downregulated in all aerial organs of glossy plants. Subcellular localization indicated that the BrCER60.A09 protein functions in the endoplasmic reticulum. A 5567-bp insertion was identified in exon 1 of BrCER60.A09 in SD369, which lead to a premature stop codon, thus causing a loss of function of the BrCER60.A09 enzyme. Moreover, comparative transcriptome analysis revealed that the ‘cutin, suberine, and wax biosynthesis’ pathway was significantly enriched, and genes involved in this pathway were almost upregulated in glossy plants. Further, two functional markers, BrWAX3-InDel and BrWAX3-KASP1, were developed and validated. Overall, these results provide a new information for the cuticular wax biosynthesis and provide applicable markers for marker-assisted selection (MAS)-based breeding of Brassica rapa.

AtMYB31 is a wax regulator associated with reproductive development in Arabidopsis

AtMYB31, a R2R3-MYB transcription factor that modulates wax biosynthesis in reproductive tissues, is involved in seed development in Arabidopsis. R2R3-MYB transcription factors play important roles in plant development; yet, the exact role of each of them remains to be resolved. Here we report that the Arabidopsis AtMYB31 is required for wax biosynthesis in epidermis of reproductive tissues, and is involved in seed development. AtMYB31 was ubiquitously expressed in both vegetative and reproductive tissues with higher expression levels in siliques and seeds, while AtMYB31 was localized to the nucleus and cytoplasm. Loss of function of AtMYB31 reduced wax accumulation in the epidermis of silique and flower tissues, disrupted seed coat epidermal wall development and mucilage production, altered seed proanthocyanidin and polyester content. AtMYB31 could direct activate expressions of several wax biosynthetic target genes. Altogether, AtMYB31, a R2R3-MYB transcription factor, regulates seed development in Arabidopsis.

ArabidopsisKCS5 and KCS6 Play Redundant Roles in Wax Synthesis

3-ketoacyl-CoA synthases (KCSs), as components of a fatty acid elongase (FAE) complex, play key roles in determining the chain length of very-long-chain fatty acids (VLCFAs). KCS6, taking a predominate role during the elongation from C26 to C28, is well known to play an important role in wax synthesis. KCS5 is one paralog of KCS6 and its role in wax synthesis remains unknown. Wax phenotype analysis showed that in kcs5 mutants, the total amounts of wax components derived from carbon 32 (C32) and C34 were apparently decreased in leaves, and those of C26 to C32 derivatives were obviously decreased in flowers. Heterologous yeast expression analysis showed that KCS5 alone displayed specificity towards C24 to C28 acids, and its coordination with CER2 and CER26 catalyzed the elongation of acids exceeding C28, especially displaying higher activity towards C28 acids than KCS6. BiLC experiments identified that KCS5 physically interacts with CER2 and CER26. Wax phenotype analysis of different organs in kcs5 and kcs6 single or double mutants showed that KCS6 mutation causes greater effects on the wax synthesis than KCS5 mutation in the tested organs, and simultaneous repression of both protein activities caused additive effects, suggesting that during the wax biosynthesis process, KCS5 and KCS6 play redundant roles, among which KCS6 plays a major role. In addition, simultaneous mutations of two genes nearly block drought-induced wax production, indicating that the reactions catalyzed by KCS5 and KCS6 play a critical role in the wax biosynthesis in response to drought.