Dark jute (Corchorus olitorius L.) can absorb lead (Pb) in the roots and translocate it to the shoots, making it a potential candidate for phytoremediation of Pb-contaminated soils. However, the genetic and molecular mechanisms regulating this process remain unclear. In this study, VA-50, a dwarf dark jute accession that accumulates Pb only in roots, and 18BS-81, a tall, late-flowering accession capable of transferring Pb to its shoots, were selected as parents. An F2population from their cross was used to map QTL intervals and identify candidate genes related to Pb accumulation and tolerance. Thirty extreme sensitive and resistant F2individuals were subjected to Bulk Segregant Analysis (BSA) under 400 mu mol/L Pb stress, together with transcriptome analysis of the parents. BSA-seq identified 668,436 high-quality SNPs and overlapping SNP and InDel signals revealed six candidate regions on chromosome 5 spanning 0.46 Mb. Transcriptome data were also used to screen 142 genes within this region, from which six candidate genes were selected. COS05g_00848 and COS05g_00849 encode lipoxygenases (LOX), COS05g_00890 encodes a mitochondrial ADP/ATP carrier protein, COS05g_00801 encodes PTZ00169 protein, while COS05g_01275 and COS05g_01276 lack predicted protein families, suggesting novel functions. These genes are associated with lipid oxidation pathways linked to glycerolipid and glycerophospholipid metabolism, which support membrane stability and mitigate Pb toxicity in sensitive cellular components. This study integrates BSA-seq and RNA-seq to identify candidate genes regulating Pb accumulation and tolerance in dark jute, providing a foundation for genetic improvement of dark jute for phytoremediation applications.

The U6 promoter plays a pivotal role in the CRISPR/Cas9 system by driving the transcription of single guide RNA (sgRNA), which directs Cas9 to achieve precise genome editing. Endogenous U6 promoters typically exhibit superior transcriptional activation efficiency compared to exogenous counterparts, thereby enhancing the efficacy of genome editing. However, the endogenous U6 promoter in kenaf (Hibiscus cannabinus L.) remains uncharacterized. In this study, we conducted a homologous search of the kenaf genome using the Arabidopsis U6 (AtU6-26) RNA sequence as a reference, identifying two candidate promoters, HcU6-1 and HcU6-14. Promoter fragments were amplified from the genomic DNA of kenaf cultivar 'Fuhong 952' and subsequently cloned into a GUS fusion expression vector. Histochemical staining revealed transcriptional activity for both promoters, with HcU6-14 demonstrating significantly stronger activity. To evaluate editing efficiency, we constructed a CRISPR/Cas9 vector containing HcALS sgRNA, driven by either the kenaf U6-14P promoter or the cotton U6-9P (GbU6-9P) promoter. Kenaf hairy roots were regenerated via Agrobacterium rhizogenes K599-mediated transformation. Sequencing analysis of ALS gene fragments from these hairy roots confirmed successful targeted editing when using the kenaf U6-14P promoter, whereas no base mutations were detected with the cotton U6 promoter. These findings highlight the superior editing efficiency of the kenaf U6 promoter and provide a critical foundation for advancing functional genomics research in kenaf.

Dark jute (Corchorus olitorius L.) can absorb lead (Pb) in the roots and translocate it to the shoots, making it a potential candidate for phytoremediation of Pb-contaminated soils. However, the genetic and molecular mechanisms regulating this process remain unclear. In this study, VA-50, a dwarf dark jute accession that accumulates Pb only in roots, and 18BS-81, a tall, late-flowering accession capable of transferring Pb to its shoots, were selected as parents. An F₂ population from their cross was used to map QTL intervals and identify candidate genes related to Pb accumulation and tolerance. Thirty extreme sensitive and resistant F₂ individuals were subjected to Bulk Segregant Analysis (BSA) under 400 μmol/L Pb stress, together with transcriptome analysis of the parents. BSA-seq identified 668,436 high-quality SNPs and overlapping SNP and InDel signals revealed six candidate regions on chromosome 5 spanning 0.46 Mb. Transcriptome data were also used to screen 142 genes within this region, from which six candidate genes were selected. COS05g_00848 and COS05g_00849 encode lipoxygenases (LOX), COS05g_00890 encodes a mitochondrial ADP/ATP carrier protein, COS05g_00801 encodes PTZ00169 protein, while COS05g_01275 and COS05g_01276 lack predicted protein families, suggesting novel functions. These genes are associated with lipid oxidation pathways linked to glycerolipid and glycerophospholipid metabolism, which support membrane stability and mitigate Pb toxicity in sensitive cellular components. This study integrates BSA-seq and RNA-seq to identify candidate genes regulating Pb accumulation and tolerance in dark jute, providing a foundation for genetic improvement of dark jute for phytoremediation applications.

Dark jute (Corchorus olitorius L.) is a high-biomass fiber crop with potential value for phytoremediation, but the genetic basis of whole-plant cadmium (Cd) accumulation and tolerance remains poorly understood. Identifying loci and candidate genes associated with Cd-related traits can support molecular improvement of dark jute for cultivation in Cd-affected environments. In this study, we integrated bulked segregant analysis sequencing (BSA-seq) with time-course RNA sequencing (RNA-seq), Cd phenotyping, and qRT-PCR validation to investigate Cd accumulation and tolerance in dark jute. A Cd-tolerant, low-accumulating accession (CO57) and a Cd-sensitive, high-accumulating accession (CO83) were used to generate an F₂ population for extreme-phenotype bulking. Time-course RNA-seq of the tolerant parent revealed substantial Cd-responsive transcriptional reprogramming, with 3,784-4,167 differentially expressed genes (DEGs) detected across root control-versus-Cd comparisons and 2,216-3,532 DEGs detected across aboveground tissue comparisons. Overlap analysis showed 1,207 shared DEGs between root and aboveground DEG sets. GO classification and KEGG pathway annotation provided functional context for the Cd-responsive gene set, indicating the involvement of cellular regulation, metabolic responses, stress response, transport, detoxification, sulfur metabolism, glutathione-related metabolism, and protein degradation pathways. For BSA-seq, 50 extremely low-accumulating and 50 extremely high-accumulating F₂ individuals were pooled, yielding 625,354 high-quality SNPs. ΔSNP-index analysis detected five candidate intervals spanning 7.13 Mb and containing 662 annotated genes. By integrating positional evidence, Cd-responsive expression patterns, and functional annotations, seven candidate genes were prioritized, including genes encoding ABC transporter-like proteins, a NAC-domain transcription factor, and two PLAC8-domain proteins annotated as Cd-resistance-related proteins. qRT-PCR analysis supported the expression patterns of the prioritized candidates. These results identify candidate genomic intervals and prioritized genes associated with whole-plant Cd accumulation and tolerance in dark jute. The study provides a genetic and transcriptomic framework for future functional validation, marker development, and phytoremediation-oriented improvement of dark jute.

Cannabinoids, the bioactive compounds unique to Cannabis plants, demonstrate significant therapeutic potential in treating epilepsy, pain management, and inflammatory conditions. Despite their clinical importance, the genetic mechanisms underlying cannabinoid biosynthesis in seed hemp remain poorly understood. Our study systematically identified 25 key genes involved in the cannabinoid synthesis pathway within the seed hemp genome. Comprehensive bioinformatics analyses were conducted to characterize their physicochemical properties, evolutionary relationships, gene structures, promoter cis-elements, and expression patterns. Transcriptomic profiling of ‘Yushe’ cultivar revealed tissue-specific expression patterns across roots, stems, leaves, and male/female flowers. These genes were organized into four distinct biosynthetic pathways: MEP (methylerythritol phosphate), GPP (geranyl diphosphate), acetic acid, and cannabinoid-specific pathways. Promoter analysis identified critical cis-elements associated with light responsiveness, phytohormone signaling, and stress adaptation. Notably, seven genes (CsaDXS, CsaHDS, CsaHDR, CsaGPP(lsu), CsaHPL, CsaTHCAS, CsaCBDAS) exhibited higher expression in male leaves, while four genes (CsaDXR, CsaGPP(lsu), CsaHPL, CsaCBDAS) were upregulated in female flowers. Comparative genomic analysis revealed that medicinal cannabis types possessed higher gene copy numbers for CsaDXS, CsaHDS, CsaTHCAS, and CsaCBDAS, which may account for their elevated cannabinoid production. Functional validation through VIGS-mediated silencing of PT gene demonstrated significant downregulation of THCAS and CBDAS expression. This study provides crucial genetic resources for functional characterization of cannabinoid synthase genes and facilitates the rational utilization of seed hemp germplasm.

The ARF (Auxin response factor), GH3 (Gretchen Hagen 3), and Aux/IAA (Auxin/indole-3-acetic acid) gene families are key components in auxin signaling pathway and function as regulators of growth in plants. However, this research is rarely reported in jute, which severely limits the understandings of mechanisms involved in fiber development. In this study, 13 ARF, 12 GH3 and 20 Aux/IAA putative genes were identified in the whole genome of jute. Exon-intron structures revealed the high conservation among these auxin-related gene family members. Chromosomal localization and synteny analysis showed that segmental duplication contributed to the expansion of CcARF, CcGH3 and CcIAA gene families. Phylogenetic and conserved motif analysis revealed that they have distinct functional and CcARF, CcGH3 and CcIAA-specific domains, respectively. The expression pattern analysis based on RNA-seq and qRT-PCR indicated that 7 CcARF, 5 CcGH3, and 14 CcIAA genes showed higher expression in stem barks than leaves at the vigorous vegetative growth stage in an elite cultivar Huangma 179 with normal plant height, respectively, suggesting they might regulate the development of bast fiber. Moreover, the expression of 5 CcARF, 4 CcGH3, and 12 CcIAA genes was differentially expressed in stem barks of a typical GA3 sensitive dwarf germplasm in comparison to Huangma 179. The cis-element analysis showed that promoters of 4 CcARF, 3 CcGH3, and 7 CcIAA genes had 1 to 3 cis-elements involved in gibberellin-responsiveness, giving a hint that they could respond to endogenous gibberellin accumulation in Huangma 179 and form a complicated network to regulate hormone regulatory network and plant height. This study provides useful information for functional analysis of ARF, GH3, and Aux/IAA genes, which would be taken as candidates for genetic improvement of bast fiber quality in jute.

Jute (Corchorus capsularis L.) is the second most important natural plant fiber source after cotton. However, developing an efficient gene editing system for jute remains a challenge. In this study, the transgenic hairy root system mediated by Agrobacterium rhizogenes strain K599 was developed for Meifeng 4, an elite jute variety widely cultivated in China. The transgenic hairy root system for jute was verified by subcellular localization and bimolecular fluorescence complementation (BiFC) assays. The CHLOROPLASTOS ALTERADOS 1 (CcCLA1) gene, which is involved in the development of chloroplasts, was targeted for editing at two sites in Meifeng 4. Based on this hairy root transformation, the gRNA scaffold was placed under the control of cotton ubiquitin GhU6.7 and -GhU6.9 promoters, respectively, to assess the efficiency of gene editing. Results indicated the 50.0% (GhU6.7) and 38.5% (GhU6.9) editing events in the target 2 alleles (gRNA2), but no mutation was detected in the target 1 allele (gRNA1) in transgenic-positive hairy roots. CcCLA1 gene editing at gRNA2 under the control of GhU6.7 in Meifeng 4 was also carried out by Agrobacterium tumefaciens-mediated transformation. Two CcCLA1 mutants were albinic, with a gene editing efficiency of 5.3%. These findings confirm that the CRISPR/Cas9 system, incorporating promoter GhU6.7, can be used as a gene editing tool for jute.

Protein phosphatase 2 C (PP2C) gene family is known for regulating a variety of abiotic stresses in plants. However, little is currently understood about this gene family in jute. In this study, a total of 60 CcPP2C genes were identified in jute, which were named from CcPP2C1 to CcPP2C60 according to the order of their physical coordinates. These CcPP2C proteins contained about 162-1324 amimo acids. Their isoelectric points varied in the range of 4.51-9.45, and their molecular weights varied from 18.307 to 146.906 kDa. The evolutionary tree of PP2Cs in jute and Arabidopsis indicated that PP2Cs can be divided into 12 branches. Sixty-seven collinear gene pairs were identified between 42 CcPP2Cs and 62 AtPP2Cs by analyzing the whole genome collinearity data of jute and Arabidopsis. The analyses of conserved motifs in these CcPP2C proteins and cis-elements in the promoters revealed that CcPP2Cs might play an extensive role in responses to various abiotic stresses in jute. Particularly, the expression patterns of CcPP2Cs in group A under salt treatments showed that CcPP2C51 was up-regulated under various stress treatments. And CcPP2C51 was localized in the nucleus and cell membrane. In Arabidopsis, the overexpression of CcPP2C51 positively regulated the salt-tolerance genes. Also, the overexpression of CcPP2C51 improved the salt tolerance of jute regenerated roots. This study revealed the origin and evolution of PP2C genes in jute, and provides potential possibilities for molecular breeding of jute.

Background Satellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes. Results In this study, we investigated the satellite repeat signature in five okra ( Abelmoschus esculentus ) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA. Conclusions Our findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.

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 GA(3) 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 GA(3) 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.

Roselle (Hibiscus sabdariffa L.) simple sequence repeat (SSR) markers were developed using RNA sequencing technology, providing a foundation for genetic analysis and the identification of roselle varieties. In this study, 10 785 unigenes containing 12 994 SSR loci with an average of one SSR locus per 6.87 Kb were identified, and the occurrence frequency of the SSR loci was 11.36%. Trinucleotide repeat motifs were the most abundant, followed by dinucleotide repeats, with AAG/CTT and AT/AT being the predominant types, respectively. After screening 100 primer pairs with a polymorphic ratio of 32.0%, we obtained 32 primer pairs, resulting in clear and stable polymorphic bands. Twenty-seven primer pairs were highly or moderately polymorphic, and seven primer pairs were highly polymorphic. Genetic relationship analysis based on the selected SSR primers showed that 38 roselle accessions were classified into different clades, with those from the same regions clustered into the same subgroups. In contrast, individuals with unique morphological traits were separated. DNA fingerprints of 38 roselle varieties were constructed using five SSR primers, providing an effective method for identifying roselle varieties at a molecular level. Our data provide novel insights into the genetics of H. sabdariffa and may be used in SSR-assisted roselle breeding.

Jute is considered one of the most important crops for fiber production and multipurpose usages. Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is a crucial enzyme involved in lignin biosynthesis in plants. The potential functions of CCoAOMT in lignin biosynthesis of jute have been reported in several studies. However, little is known about the evolution of the CCoAOMT gene family, and either their expression level at different developing stages in different jute cultivars, as well as under abiotic stresses including salt and drought stress. In the present study, 66 CCoAOMT genes from 12 species including 12 and eight CCoAOMTs in Corchorus olitorius and C. capsularis were identified. Phylogenetic analysis revealed that CCoAOMTs could be divided into six groups, and gene expansion was observed in C. olitorius. Furthermore, gene expression analysis of developing jute fibers was conducted at different developmental stages (15, 30, 45, 60, and 90 days after sowing [DAS]) in six varieties (Jute-179 [J179], Lubinyuanguo [LB], and Qiongyueqing [QY] for C. capsularis; Funong No.5 [F5], Kuanyechangguo [KY], and Cvlv [CL] for C. olitorius). The results showed that CCoAOMT1 and CCoAOMT2 were the dominant genes in the CCoAOMT family. Of these two dominant CCoAOMTs, CCoAOMT2 showed a constitutive expression level during the entire growth stages, while CCoAOMT1 exhibited differential expression patterns. These two genes showed higher expression levels in C. olitorius than in C. capsularis. The correlation between lignin content and CCoAOMT gene expression levels indicated that this gene family influences the lignin content of jute. Using real-time quantitative reverse transcription PCR (qRT-PCR), a substantial up-regulation of CCoAOMTs was detected in stem tissues of jute 24 h after drought treatment, with an up to 17-fold increase in expression compared to that of untreated plants. This study provides a basis for comprehensive genomic studies of the entire CCoAOMT gene family in C. capsularis and C. olitorius. Comparative genomics analysis among the CCoAOMT gene families of 12 species revealed the close evolutionary relationship among Corchorus, Theobroma cacao and Gossypium raimondii. This study also shows that CCoAOMTs are not only involved in lignin biosynthesis, but also are associated with the abiotic stress response in jute, and suggests the potential use of these lignin-related genes to genetically improve the fiber quality of jute.

[目的]丰富玫瑰茄(Hibiscus sabdariffa L.)SSR分子标记的种类和数量,验证锦葵科作物SSR标记在玫瑰茄中的通用性,进而分析玫瑰茄的遗传多样性.[方法]选取锦葵科红麻、棉花、黄秋葵和木槿等4种作物共62对SSR引物(木槿17对,其他各15对),对8份不同来源的玫瑰茄种质进行多态性分析.[结果]62对SSR引物中有58对具有通用性,占比为93.5％;通用性引物中有48对具有多态性,占比为82.8％;这些引物共扩增出246条条带,其中多态性条带239条,占比97.2％.同时,棉花有14对SSR引物具有多态性,多态性引物比率最高,达86.7％;其次是木槿,有82.4％的引物具有多态性,且其引物多态信息含量值(PIC)最高,为0.685;而黄秋葵和红麻多态性引物比率较低,分别为73.3％和66.7％.另外,聚类分析结果表明,在相似系数0.49处可将8份玫瑰茄种质分为3个各具特点的亚类,其中19FZ-76被单独聚为一类,表明这一种质与其他7份种质亲缘关系最远,遗传差异最大,聚类结果与其表型性状表现一致;而来自缅甸的19FZ-74和漳州本地种遗传相似系数最高,说明这两份种质的亲缘关系较近,遗传差异较小,表明我国福建种植的玫瑰茄品种与东南亚的品种有较近的血缘关系,有可能从东南亚引种而来.[结论]锦葵科作物的SSR引物在玫瑰茄中具有较好的通用性,筛选出的SSR标记可为玫瑰茄品种鉴定、亲缘关系分析、遗传图谱构建等遗传学研究提供可靠方法.

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.

Virus-induced gene silencing (VIGS) is a potent way to study gene function in plants. In this study, a tobacco rattle virus (TRV)-mediated gene silencing system (TRV-VIGS) was established with the phytoene desaturase (HcPDS) as the marker gene in kenaf. Firstly, the seed soak agroinoculation method of the VIGS system was explored. HcPDS gene silencing led to photobleaching phenotype in kenaf leaves. According to qRT-PCR analysis, photobleached leaves had significantly lower gene expression levels compared with the control. Secondly, an optimizing approach with leaf infiltration was used to improve silencing efficiency. After the kenaf leaves were injected with Agrobacterium cultures with OD600 0.8, the plants were cultivated in chamber growth at 24 °C with approximately 60

Dear Editor, Hemp(Cannabis sativa L.)(2n=20)is an annual herb of the Can-nabaceae family,mostly dioecious and occasionally monoecious(Van et al.,2011;Long et al.,2017).A major source of food,textiles,and oilseed,it is believed to have originated in central Asia and spread rapidly throughout Asia and Europe(Andre et al.,2016;Braich et al.,2019;Gao et al.,2020).Depending on its usage and cannabinoidcontent,cannabis is classified as fiber type(hemp or industrial hemp)or drug type(medicinal cannabis or marijuana)(Struik et al.,2000;Hurgobin et al.,2021).As a type of industrial hemp,seed hemp is widely used to produce seed oil,but we still know little about the genetic basis of its fatty acid and vitamin E synthesis.

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

MYB transcription factors (TFs) are vital players in various plant biochemical and physiological processes. To date, MYB TFs have been comprehensively investigated in numerous plant species. Corchorus capsularis is one of the most important fiber crops worldwide. However, systematic investigation of the MYB TFs in C. capsularis remains elusive, which limits our insights into the contribution of MYB TFs regulating numerous biochemical processes such as anthocyanin accumulation and physiological processes (Color, growth, etc.) in this crop. Therefore, this study systematically investigated MYB genes comprising the gene structure, phylogeny, chromosomal locations, protein motifs, expression patterns, protein-protein interaction, and collinearity analysis. A total of 88 Corchorus capsularis R2R3-MYB (CcMYB) genes were characterized and distributed on seven chromosomes. The collinearity study revealed that the segmental duplication events were responsible for the enlargement of CcMYB genes. The highest syntenic relationship was observed between CcMYB genes (70/88) and Gossypium hirsutism genes. Further, phylogenetic tree analysis classified the 88 CcMYBs into 25 subgroups based on domain similarity and phylogenetic topology. Furthermore, the phylogenetic investigation of CcMYBs and Arabidopsis MYBs displayed the functional discrepancy of the CcMYB genes throughout evolution. Changes in anthocyanin content were observed among white jute's stem bark tissues, which may be due to their different colors. Expression study of 11 CcMYB genes in total by real-time quantitative PCR revealed the diverse expression profiles in various tissues organs and functional predictions found that some CcMYBs might be associated with stem bark color and anthocyanin content, cellulose, and lignin biosynthesis. Co-expression interaction analysis showed that CcMYB genes could be involved in physiological, biological, and molecular processes. Our study provides comprehensive details on MYB genes and contributes to a better understanding of their functions in Corchorus capsularis.