Reference genomes of the two cultivated jute species

Whole genome resequencing unveils low-temperature stress tolerance specific genomic variations in jute (Corchorus sp.)

Jute (Corchorus sp.), a commercially important and eco-friendly crop, is widely cultivated in Bangladesh, India, and China. Some varieties of this tropical plant such as the Corchorus olitorius. Variety accession no. 2015 (acc. 2015) has been found to be low-temperature tolerant. The current study was designed to explore the genome-wide variations present in the tolerant plant acc. 2015 in comparison to the sensitive farmer popular variety Corchorus olitorius var. O9897 using the whole genome resequencing technique. Among different variations, intergenic Single Nucleotide Polymorphism (SNPs) and Insertion-Deletion (InDels) were found in the highest percentage whereas approximately 3% SNPs and 2% InDels were found in exonic regions in both plants. Gene enrichment analysis indicated the presence of acc. 2015 specific SNPs in the genes encoding peroxidase, ER lumen protein retaining receptor, and hexosyltransferase involved in stress response (GO:0006950) which were not present in sensitive variety O9897. Besides, distinctive copy number variation regions (CNVRs) comprising 120 gene loci were found in acc. 2015 with a gain of function from multiple copy numbers but absent in O9897. Gene ontology analysis revealed these gene loci to possess different receptors like kinases, helicases, phosphatases, transcription factors especially Myb transcription factors, regulatory proteins containing different binding domains, annexin, laccase, acyl carrier protein, potassium transporter, and vesicular transporter proteins that are responsible for low temperature induced adaptation pathways in plants. This work of identifying genomic variations linked to cold stress tolerance traits will help to develop successful markers that will pave the way to develop genetically modified cold-resistant jute lines for year-round cultivation to meet the demand for a sustainable fiber crop economy.

Cloning and functional mechanism of the dwarf gene gba affecting stem elongation and cellulose biosynthesis in jute (Corchorus olitorius)

Plant height (PH) is an important factor affecting bast fiber yield in jute. Here, we report the mechanism of dwarfism in the 'Guangbaai' (gba) of jute. The mutant gba had shorter internode length and cell length compared to the standard cultivar 'TaiZi 4' (TZ4). Exogenous GA3 treatment indicated that gba is a GA-insensitive dwarf mutant. Quantitative trait locus (QTL) analysis of three PH-related traits via a high-density genetic linkage map according to re-seq showed that a total of 25 QTLs were identified, including 13 QTLs for PH, with phenotypic variation explained ranging from 2.42 to 74.16%. Notably, the functional mechanism of the candidate gene CoGID1a, the gibberellic acid receptor, of the major locus qPHIL5 was evaluated by transgenic analysis and virus-induced gene silencing. A dwarf phenotype-related single nucleotide mutation in CoGID1a was identified in gba, which was also unique to the dwarf phenotype of gba among 57 cultivars. Cogid1a was unable to interact with the growth-repressor DELLA even in the presence of highly accumulated gibberellins in gba. Differentially expressed genes between transcriptomes of gba and TZ4 after GA3 treatment indicated up-regulation of genes involved in gibberellin and cellulose synthesis in gba. Interestingly, it was found that up-regulation of CoMYB46, a key transcription factor in the secondary cell wall, by the highly accumulated gibberellins in gba promoted the expression of cellulose synthase genes CoCesA4 and CoCesA7. These findings provide valuable insights into fiber development affected by endogenous gibberellin accumulation in plants.

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Proteo-metabolomic Dissection of Extracellular Matrix Reveals Alterations in Cell Wall Integrity and Calcium Signaling Governs Wall-Associated Susceptibility during Stem Rot Disease in Jute

The plant surveillance system confers specificity to disease and immune states by activating distinct molecular pathways linked to cellular functionality. The extracellular matrix (ECM), a preformed passive barrier, is dynamically remodeled at sites of interaction with pathogenic microbes. Stem rot, caused by Macrophomina phaseolina, adversely affects fiber production in jute. However, how wall related susceptibility affects the ECM proteome and metabolome remains undetermined in bast fiber crops. Here, stem rot responsive quantitative temporal ECM proteome and metabolome were developed in jute upon M. phaseolina infection. Morpho-histological examination revealed that leaf shredding was accompanied by reactive oxygen species production in patho-stressed jute. Electron microscopy showed disease progression and ECM architecture remodeling due to necrosis in the later phase of fungal attack. Using isobaric tags for relative and absolute quantitative proteomics and liquid chromatography-tandem mass spectrometry, we identified 415 disease-responsive proteins involved in wall integrity, acidification, proteostasis, hydration, and redox homeostasis. The disease-related correlation network identified functional hubs centered on α-galactosidase, pectinesterase, and thaumatin. Gas chromatography-mass spectrometry analysis pointed toward enrichment of disease-responsive metabolites associated with the glutathione pathway, TCA cycle, and cutin, suberin, and wax metabolism. Data demonstrated that wall-degrading enzymes, structural carbohydrates, and calcium signaling govern rot responsive wall-susceptibility. Proteomics data were deposited in Pride (PXD046937; PXD046939).

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata)

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

Genome-wide association study reveals loci and candidate genes of flowering time in jute (Corchorus L.)

Suitable flowering time can improve fiber yield and quality, which is of great significance for jute biological breeding. In this study, 242 jute accessions were planted in Fujian for 2 consecutive years, and 244,593 SNPs distributed in jute genome were used for genome-wide association analysis of flowering time. A total of 19 candidate intervals (P < 0.0001) were identified by using GLM and FaST-LMM and were significantly associated with flowering time, with phenotypic variation explained (PVE) ranging from 5.8 to 18.61%. Six stable intervals that were repeatedly detected in different environments were further identified by the linkage disequilibrium heatmap. The most likely 7 candidate genes involved to flowering time were further predicted according to the gene functional annotations. Notably, functional analysis of the candidate gene CcPRR7 of the major loci qFT-3-1, a key factor in circadian rhythm in the photoperiodic pathway, was evaluated by linkage, haplotype, and transgenic analysis. β-glucuronidase (GUS) and luciferase (LUC) activity assay of the promoters with two specific haplotypes confirmed that the flowering time can be controlled by regulating the expression of CcPRR7. The model of CcPRR7 involved in the photoperiod regulation pathway under different photoperiods was proposed. These findings provide insights into genetic loci and genes for molecular marker-assisted selection in jute and valuable information for genetically engineering PRR7 homologs in plants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01435-8.

Systematic evolution of bZIP transcription factors in Malvales and functional exploration of AsbZIP14 and AsbZIP41 in Aquilaria sinensis

Introduction

Agarwood, the dark-brown resin produced by Aquilaria trees, has been widely used as incense, spice, perfume or traditional medicine and 2-(2-phenethyl) chromones (PECs) are the key markers responsible for agarwood formation. But the biosynthesis and regulatory mechanism of PECs were still not illuminated. The transcription factor of basic leucine zipper (bZIP) presented the pivotal regulatory roles in various secondary metabolites biosynthesis in plants, which might also contribute to regulate PECs biosynthesis. However, molecular evolution and function of bZIP are rarely reported in Malvales plants, especially in Aquilaria trees.

Methods and results

Here, 1,150 bZIPs were comprehensively identified from twelve Malvales and model species genomes and the evolutionary process were subsequently analyzed. Duplication types and collinearity indicated that bZIP is an ancient or conserved TF family and recent whole genome duplication drove its evolution. Interesting is that fewer bZIPs in A. sinensis than that species also experienced two genome duplication events in Malvales. 62 AsbZIPs were divided into 13 subfamilies and gene structures, conservative domains, motifs, cis-elements, and nearby genes of AsbZIPs were further characterized. Seven AsbZIPs in subfamily D were significantly regulated by ethylene and agarwood inducer. As the typical representation of subfamily D, AsbZIP14 and AsbZIP41 were localized in nuclear and potentially regulated PECs biosynthesis by activating or suppressing type III polyketide synthases (PKSs) genes expression via interaction with the AsPKS promoters.

Discussion

Our results provide a basis for molecular evolution of bZIP gene family in Malvales and facilitate the understanding the potential functions of AsbZIP in regulating 2-(2-phenethyl) chromone biosynthesis and agarwood formation.

Identification of AREB/ABF Gene Family Involved in the Response of ABA under Salt and Drought Stresses in Jute (Corchorus olitorius L.)

The abscisic acid (ABA)-responsive element binding protein/ABRE-binding factor (AREB/ABF) subfamily members are essential to ABA signaling pathways and plant adaptation to various environmental stresses. Nevertheless, there are no reports on AREB/ABF in jute (Corchorus L.). Here, eight AREB/ABF genes were identified in the C. olitorius genome and classified into four groups (A-D) based on their phylogenetic relationships. A cis-elements analysis showed that CoABFs were widely involved in hormone response elements, followed by light and stress responses. Furthermore, the ABRE response element was involved in four CoABFs, playing an essential role in the ABA reaction. A genetic evolutionary analysis indicated that clear purification selection affects jute CoABFs and demonstrated that the divergence time was more ancient in cotton than in cacao. A quantitative real-time PCR revealed that the expression levels of CoABFs were upregulated and downregulated under ABA treatment, indicating that CoABF3 and CoABF7 are positively correlated with ABA concentration. Moreover, CoABF3 and CoABF7 were significantly upregulated in response to salt and drought stress, especially with the application of exogenous ABA, which showed higher intensities. These findings provide a complete analysis of the jute AREB/ABF gene family, which could be valuable for creating novel jute germplasms with a high resistance to abiotic stresses.

Caffeoyl-CoA 3-O-methyltransferase gene family in jute: Genome-wide identification, evolutionary progression and transcript profiling under different quandaries

Jute (Corchorus sp.), is a versatile, naturally occurring, biodegradable material that holds the promising possibility of diminishing the extensive use of plastic bags. One of the major components of the cell wall, lignin plays both positive and negative roles in fiber fineness and quality. Although it gives mechanical strength to plants, an excess amount of it is responsible for the diminution of fiber quality. Among various gene families involved in the lignin biosynthesis, Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is the most significant and has remained mostly unexplored. In this study, an extensive in-silico characterization of the CCoAOMT gene family was carried out in two jute species (C. capsularis L. and C. olitoroius L.) by analyzing their structural, functional, molecular and evolutionary characteristics. A total of 6 CCoAOMT gene members were identified in each of the two species using published reference genomes. These two jute species showed high syntenic conservation and the identified CCoAOMT genes formed four clusters in the phylogenetic tree. Histochemical assay of lignin in both jute species could shed light on the deposition pattern in stems and how it changes in response to abiotic stresses. Furthermore, expression profiling using qPCR showed considerable alteration of CCoAOMT transcripts under various abiotic stresses and hormonal treatment. This study will lay a base for further analysis and exploration of target candidates for overexpression of gene silencing using modern biotechnological techniques to enhance the quality of this economically important fiber crop.

Multi-omics provides new insights into the domestication and improvement of dark jute (Corchorus olitorius)

Jute (Corchorus sp.) is the most important bast fiber crop worldwide; however, the mechanisms underlying domestication and improvement remain largely unknown. We performed multi-omics analysis by integrating de novo sequencing, resequencing, and transcriptomic and epigenetic sequencing to clarify the domestication and improvement of dark jute Corchorus olitorius. We demonstrated that dark jute underwent early domestication and a relatively moderate genetic bottleneck during improvement breeding. A genome-wide association study of 11 important agronomic traits identified abundant candidate loci. We characterized the selective sweeps in the two breeding stages of jute, prominently, soil salinity differences played an important role in environmental adaptation during domestication, and the strongly selected genes for improvement had an increased frequency of favorable haplotypes. Furthermore, we speculated that an encoding auxin/indole-3-acetic acid protein COS07g_00652 could enhance the flexibility and strength of the stem to improve fiber yield. Our study not only provides valuable genetic resources for future fiber breeding in jute, but also is of great significance for reviewing the genetic basis of early crop breeding.

Intelligent Identification of Jute Pests Based on Transfer Learning and Deep Convolutional Neural Networks

Pest attacks pose a substantial threat to jute production and other significant crop plants. Jute farmers in Bangladesh generally distinguish between different pests that appear to be the same using their eyes and expertise, which isn't always accurate. We developed an intelligent model for jute pests identification based on transfer learning (TL) and deep convolutional neural networks (DCNN) to solve this practical problem. The proposed DCNN model can realize fast and accurate automatic identification of jute pests based on photographs. Specifically, the VGG19 CNN model was trained by TL on the ImageNet database. A well-structured image dataset of four dominant jute pests is also established. Our model shows a final accuracy of 95.86% on the four most vital jute pest classes. The model's performance is further demonstrated by the precision, recall, F1-score, and confusion matrix results. The proposed model is integrated into Android and IOS applications for practical uses.