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.

Transcriptomic and metabolomic analyses reveal the differential accumulation of phenylpropanoids and terpenoids in hemp autotetraploid and its diploid progenitor.

BACKGROUND: Cannabis sativa, a dioecious plant that has been cultivated worldwide for thousands of years, is known for its secondary metabolites, especially cannabinoids, which possess several medicinal effects. In this study, we investigated the autopolyploidization effects on the biosynthesis and accumulation of these metabolites, transcriptomic and metabolomic analyses were performed to explore the gene expression and metabolic variations in industrial hemp autotetraploids and their diploid progenitors. RESULTS: Through these analyses, we obtained 1,663 differentially expressed metabolites and 1,103 differentially expressed genes. Integrative analysis revealed that phenylpropanoid and terpenoid biosynthesis were regulated by polyploidization. No substantial differences were found in the cannabidiol or tetrahydrocannabinol content between tetraploids and diploids. Following polyploidization, some transcription factors, including nine bHLH and eight MYB transcription factors, affected the metabolic biosynthesis as regulators. Additionally, several pivotal catalytic genes, such as flavonol synthase/flavanone 3-hydroxylase, related to the phenylpropanoid metabolic pathway, were identified as being modulated by polyploidization. CONCLUSIONS: This study enhances the overall understanding of the impact of autopolyploidization in C. sativa and the findings may encourage the application of polyploid breeding for increasing the content of important secondary metabolites in industrial hemp.

Differential Proteomic Analysis to Identify Proteins Associated with Apomeiosis in Boehmeria tricuspis (Hance) Makino Using an iTRAQ-Based Strategy.

Numerous candidate genes related to apomixis have been identified through transcriptomics; however, the molecular mechanism underlying apomixis remains unclear. Elucidation of the underlying mechanisms is essential to expand its application in crop breeding. Therefore, here, we employed the isobaric tags for a relative and absolute quantification labeling technology to investigate the protein expression in Boehmeria tricuspis generated through different reproductive modes at the functional megaspore stage. We identified 40 differential abundance proteins associated with apomeiosis, most of which were involved in "response to stress". Functional analysis suggested that lower levels of reactive oxygen species (ROS) play a role in inducing the development of apomeiosis. Proteins related to ROS regulation, cell wall modifications, and stability under heat stress play a crucial role in the development of diplosporic apomeiosis. Our results give evidence to the insight that stress can induce a switch from apomixis to sexuality by ROS content, and an increased composition of stress tolerance as well as secondary metabolites can buffer ROS effects. Precise coordination of these proteins involved in inter-related regulatory control mechanisms may act together in the transition from the sexual to apomixis development.

Establishment of an Agrobacterium-mediated genetic transformation and CRISPR/Cas9-mediated targeted mutagenesis in Hemp (Cannabis Sativa L.).

Hemp (Cannabis sativa L.) is an annual and typically dioecious crop. Due to the therapeutic potential for human diseases, phytocannabinoids as a medical therapy is getting more attention recently. Several candidate genes involved in cannabinoid biosynthesis have been elucidated using omics analysis. However, the gene function was not fully validated due to few reports of stable transformation for Cannabis tissues. In this study, we firstly report the successful generation of gene-edited plants using an Agrobacterium-mediated transformation method in C. sativa. DMG278 achieved the highest shoot induction rate, which was selected as the model strain for transformation. By overexpressing the cannabis developmental regulator chimera in the embryo hypocotyls of immature grains, the shoot regeneration efficiency was substantially increased. We used CRISPR/Cas9 technology to edit the phytoene desaturase gene and finally generated four edited cannabis seedlings with albino phenotype. Moreover, we propagated the transgenic plants and validated the stable integration of T-DNA in cannabis genome.

Gene coexpression network analysis and tissue-specific profiling of gene expression in jute (Corchorus capsularis L.).

BACKGROUND: Jute (Corchorus spp.), belonging to the Malvaceae family, is an important natural fiber crop, second only to cotton, and a multipurpose economic crop. Corchorus capsularis L. is one of the only two commercially cultivated species of jute. Gene expression is spatiotemporal and is influenced by many factors. Therefore, to understand the molecular mechanisms of tissue development, it is necessary to study tissue-specific gene expression and regulation. We used weighted gene coexpression network analysis, to predict the functional roles of gene coexpression modules and individual genes, including those underlying the development of different tissue types. Although several transcriptome studies have been conducted on C. capsularis, there have not yet been any systematic and comprehensive transcriptome analyses for this species. RESULTS: There was significant variation in gene expression between plant tissues. Comparative transcriptome analysis and weighted gene coexpression network analysis were performed for different C. capsularis tissues at different developmental stages. We identified numerous tissue-specific differentially expressed genes for each tissue, and 12 coexpression modules, comprising 126 to 4203 genes, associated with the development of various tissues. There was high consistency between the genes in modules related to tissues, and the candidate upregulated genes for each tissue. Further, a gene network including 21 genes directly regulated by transcription factor OMO55970.1 was discovered. Some of the genes, such as OMO55970.1, OMO51203.1, OMO50871.1, and OMO87663.1, directly involved in the development of stem bast tissue. CONCLUSION: We identified genes that were differentially expressed between tissues of the same developmental stage. Some genes were consistently up- or downregulated, depending on the developmental stage of each tissue. Further, we identified numerous coexpression modules and genes associated with the development of various tissues. These findings elucidate the molecular mechanisms underlying the development of each tissue, and will promote multipurpose molecular breeding in jute and other fiber crops.

Construction of a high-resolution genetic map and identification of quantitative trait loci for salt tolerance in jute (Corchous spp.).

BACKGROUND: Jute (Corchorus spp.) is the most important natural fiber crop after cotton in terms of cultivation area and production. Salt stress greatly restricts plant development and growth. A high-density genetic linkage map is the basis of quantitative trait locus (QTLs) mapping. Several high-density genetic maps and QTLs mapping related to salt tolerance have been developed through next-generation sequencing in many crop species. However, such studies are rare for jute. Only several low-density genetic maps have been constructed and no salt tolerance-related QTL has been mapped in jute to date. RESULTS: We developed a high-density genetic map with 4839 single nucleotide polymorphism markers spanning 1375.41 cM and an average distance of 0.28 cM between adjacent markers on seven linkage groups (LGs) using an F2 jute population, LGs ranged from LG2 with 299 markers spanning 113.66 cM to LG7 with 1542 markers spanning 350.18 cM. In addition, 99.57% of gaps between adjacent markers were less than 5 cM. Three obvious and 13 minor QTLs involved in salt tolerance were identified on four LGs explaining 0.58-19.61% of the phenotypic variance. The interval length of QTL mapping varied from 1.3 to 20.2 cM. The major QTL, qJST-1, was detected under two salt stress conditions that explained 11.81 and 19.61% of the phenotypic variation, respectively, and peaked at 19.3 cM on LG4. CONCLUSIONS: We developed the first high-density and the most complete genetic map of jute to date using a genotyping-by-sequencing approach. The first QTL mapping related to salt tolerance was also carried out in jute. These results should provide useful resources for marker-assisted selection and transgenic breeding for salt tolerance at the germination stage in jute.

Development of an InDel polymorphism database for jute via comparative transcriptome analysis.

Jute (Corchorus spp.) is one of the most commercially important bast fiber crops in the world. However, molecular markers and high-density genetic maps are still lacking on jute compared with other crops. Insertion/deletion (InDel) markers, one of the most abundant sources of DNA/RNA variations in plant genomes, can easily be distinguished among different accessions using high-throughput sequencing. Using three transcriptome datasets, we identified and developed InDel markers. Altogether, 51 172 InDel sites in 18 800 unigenes were discovered, and the number of InDel loci per unigene varied from 1 to 31. Further, we found 94 InDel types, varying from 1 to 159 bp; the most common were single-nucleotide (23 028), binucleotide (9824), and trinucleotide (9182). In total, 49 563 InDels in 18 445 transcripts were discovered in the comparison between TC and YG, followed by 48 934 InDels in 18 408 transcripts between NY and YG, and 3570 InDels in 2701 unigenes between NY and TC. Additionally, there were 1273 InDel sites in 1129 unigenes with polymorphisms between any two of the three accessions. Twenty-nine (58%) primer pairs represented polymorphisms when compared to the jute accessions, and PIC varied from 0.340 to 0.680, with an average of 0.491.

Insights into the mechanism of multi-walled carbon nanotubes phytotoxicity in Arabidopsis through transcriptome and m6A methylome analysis.

With the increasing production and wide application of carbon nanotubes (CNTs), they are inevitably released into the natural environment and ecosystems, where plants are the main primary producers. Hence, it is imperative to understand the toxic effects of CNTs on plants. The molecular mechanisms underlying the toxic effects of CNTs on plants are still unclear. Therefore, in the present study, we investigated the effects of high concentrations of multi-walled CNTs (MWCNTs) on Arabidopsis. Root elongation and leaf development were severely inhibited after MWCNT exposure. Excess production of H2O2, O2-, and malondialdehyde was observed, indicating that MWCNTs induced oxidative stress. The antioxidant system was activated to counter MWCNTs-induced oxidative stress. Combinatorial transcriptome and m6A methylome analysis revealed that MWCNTs suppressed auxin signaling and photosynthesis. Reactive oxygen species metabolism, toxin metabolism, and plant responses to pathogens were enhanced to cope with the phytotoxicity of MWCNTs. Our results provide new insights into the molecular mechanisms of CNT phytotoxicity and plant defense responses to CNTs.

A Full-Length Reference Floral Transcriptome of Boehmeria tricuspis Provides Insights into Apomeiosis and Polyploidy.

Boehmeria tricuspis (Hance) Makino constitutes a hardy herbaceous or shrubby perennial native to East Asia that includes different ploidy levels and reproductive modes (diplosporous to sexual). Although several apomeiosis-associated genes have been described, the genetic control and molecular mechanisms underlying apomeiosis remain poorly understood. Moreover, the basis of the correlation between polyploidy and apomixis has not yet been clarified. We utilized long-read sequencing to produce a full-length reference floral transcriptome of B. tricuspis. Based on the generated database, gene expression of the female flowers of different ploidy levels and reproductive mode cytotypes was compared. Overall, 1,387 genes related to apomeiosis, 217 genes related to ploidy, and 9 genes associated with both apomixis and ploidy were identified. Gene Ontology analyses of this set of transcripts indicated reproductive genes, especially those related to "cell differentiation" and "cell cycle process," as significant factors regulating apomeiosis. Furthermore, our results suggested that different expressions of stress response genes might be important in the preparation for apomeiosis transition. In addition, our observations indicated that the expression of apomeiosis may not depend on polyploidy but rather on deregulation of the sexual pathway in B. tricuspis.

Genome-Wide Expression Profiles of Hemp (Cannabis sativa L.) in Response to Drought Stress.

Drought is the main environmental factor impairing hemp growth and yield. In order to decipher the molecular responses of hemp to drought stress, transcriptome changes of drought-stressed hemp (DS1 and DS2), compared to well-watered control hemp (CK1 and CK2), were studied with RNA-Seq technology. RNA-Seq generated 9.83, 11.30, 11.66, and 11.31 M clean reads in the CK1, CK2, DS1, and DS2 libraries, respectively. A total of 1292 differentially expressed genes (DEGs), including 409 (31.66%) upregulated and 883 (68.34%) downregulated genes, were identified. The expression patterns of 12 selected genes were validated by qRT-PCR, and the results were accordant with Illumina analysis. Gene Ontology (GO) and KEGG analysis illuminated particular important biological processes and pathways, which enriched many candidate genes such as NAC, B3, peroxidase, expansin, and inositol oxygenase that may play important roles in hemp tolerance to drought. Eleven KEGG pathways were significantly influenced, the most influenced being the plant hormone signal transduction pathway with 15 differentially expressed genes. A similar expression pattern of genes involved in the abscisic acid (ABA) pathway under drought, and ABA induction, suggested that ABA is important in the drought stress response of hemp. These findings provide useful insights into the drought stress regulatory mechanism in hemp.

Salt-Stress Response Mechanisms Using de Novo Transcriptome Sequencing of Salt-Tolerant and Sensitive Corchorus spp. Genotypes.

High salinity is a major environmental stressor for crops. To understand the regulatory mechanisms underlying salt tolerance, we conducted a comparative transcriptome analysis between salt-tolerant and salt-sensitive jute (Corchorus spp.) genotypes in leaf and root tissues under salt stress and control conditions. In total, 68,961 unigenes were identified. Additionally, 11,100 unigenes (including 385 transcription factors (TFs)) exhibited significant differential expression in salt-tolerant or salt-sensitive genotypes. Numerous common and unique differentially expressed unigenes (DEGs) between the two genotypes were discovered. Fewer DEGs were observed in salt-tolerant jute genotypes whether in root or leaf tissues. These DEGs were involved in various pathways, such as ABA signaling, amino acid metabolism, etc. Among the enriched pathways, plant hormone signal transduction (ko04075) and cysteine/methionine metabolism (ko00270) were the most notable. Eight common DEGs across both tissues and genotypes with similar expression profiles were part of the PYL-ABA-PP2C (pyrabactin resistant-like/regulatory components of ABA receptors-abscisic acid-protein phosphatase 2C). The methionine metabolism pathway was only enriched in salt-tolerant jute root tissue. Twenty-three DEGs were involved in methionine metabolism. Overall, numerous common and unique salt-stress response DEGs and pathways between salt-tolerant and salt-sensitive jute have been discovered, which will provide valuable information regarding salt-stress response mechanisms and help improve salt-resistance molecular breeding in jute.

De novo transcriptome sequencing of two cultivated jute species under salinity stress.

Soil salinity, a major environmental stress, reduces agricultural productivity by restricting plant development and growth. Jute (Corchorus spp.), a commercially important bast fiber crop, includes two commercially cultivated species, Corchorus capsularis and Corchorus olitorius. We conducted high-throughput transcriptome sequencing of 24 C. capsularis and C. olitorius samples under salt stress and found 127 common differentially expressed genes (DEGs); additionally, 4489 and 492 common DEGs were identified in the root and leaf tissues, respectively, of both Corchorus species. Further, 32, 196, and 11 common differentially expressed transcription factors (DTFs) were detected in the leaf, root, or both tissues, respectively. Several Gene Ontology (GO) terms were enriched in NY and YY. A Kyoto Encyclopedia of Genes and Genomes analysis revealed numerous DEGs in both species. Abscisic acid and cytokinin signal pathways enriched respectively about 20 DEGs in leaves and roots of both NY and YY. The Ca2+, mitogen-activated protein kinase signaling and oxidative phosphorylation pathways were also found to be related to the plant response to salt stress, as evidenced by the DEGs in the roots of both species. These results provide insight into salt stress response mechanisms in plants as well as a basis for future breeding of salt-tolerant cultivars.

Transcriptome Analysis of Two Species of Jute in Response to Polyethylene Glycol (PEG)- induced Drought Stress.

Drought stress results in significant crop yield losses. Comparative transcriptome analysis between tolerant and sensitive species can provide insights into drought tolerance mechanisms in jute. We present a comprehensive study on drought tolerance in two jute species-a drought tolerant species (Corchorus olitorius L., GF) and a drought sensitive species (Corchorus capsularis L., YY). In total, 45,831 non-redundant unigenes with average sequence length of 1421 bp were identified. Higher numbers of differentially expressed genes (DEGs) were discovered in YY (794) than in GF (39), implying that YY was relatively more vulnerable or hyper-responsive to drought stress at the molecular level; the two main pathways, phenylpropanoid biosynthesis and peroxisome pathway, significantly involved in scavenging of reactive oxygen species (ROS) and 14 unigenes in the two pathways presented a significant differential expression in response to increase of superoxide. Our classification analysis showed that 1769 transcription factors can be grouped into 81 families and 948 protein kinases (PKs) into 122 families. In YY, we identified 34 TF DEGs from and 23 PK DEGs, including 19 receptor-like kinases (RLKs). Most of these RLKs were downregulated during drought stress, implying their role as negative regulators of the drought tolerance mechanism in jute.

Diplosporous development in Boehmeria tricuspis: Insights from de novo transcriptome assembly and comprehensive expression profiling.

Boehmeria tricuspis includes sexually reproducing diploid and apomictic triploid individuals. Previously, we established that triploid B. tricuspis reproduces through obligate diplospory. To understand the molecular basis of apomictic development in B. tricuspis, we sequenced and compared transcriptomic profiles of the flowers of sexual and apomictic plants at four key developmental stages. A total of 283,341 unique transcripts were obtained from 1,463 million high-quality paired-end reads. In total, 18,899 unigenes were differentially expressed between the reproductive types at the four stages. By classifying the transcripts into gene ontology categories of differentially expressed genes, we showed that differential plant hormone signal transduction, cell cycle regulation, and transcription factor regulation are possibly involved in apomictic development and/or a polyploidization response in B. tricuspis. Furthermore, we suggest that specific gene families are possibly related to apomixis and might have important effects on diplosporous floral development. These results make a notable contribution to our understanding of the molecular basis of diplosporous development in B. tricuspis.

Diversity analysis in Cannabis sativa based on large-scale development of expressed sequence tag-derived simple sequence repeat markers.

Cannabis sativa L. is an important economic plant for the production of food, fiber, oils, and intoxicants. However, lack of sufficient simple sequence repeat (SSR) markers has limited the development of cannabis genetic research. Here, large-scale development of expressed sequence tag simple sequence repeat (EST-SSR) markers was performed to obtain more informative genetic markers, and to assess genetic diversity in cannabis (Cannabis sativa L.). Based on the cannabis transcriptome, 4,577 SSRs were identified from 3,624 ESTs. From there, a total of 3,442 complementary primer pairs were designed as SSR markers. Among these markers, trinucleotide repeat motifs (50.99%) were the most abundant, followed by hexanucleotide (25.13%), dinucleotide (16.34%), tetranucloetide (3.8%), and pentanucleotide (3.74%) repeat motifs, respectively. The AAG/CTT trinucleotide repeat (17.96%) was the most abundant motif detected in the SSRs. One hundred and seventeen EST-SSR markers were randomly selected to evaluate primer quality in 24 cannabis varieties. Among these 117 markers, 108 (92.31%) were successfully amplified and 87 (74.36%) were polymorphic. Forty-five polymorphic primer pairs were selected to evaluate genetic diversity and relatedness among the 115 cannabis genotypes. The results showed that 115 varieties could be divided into 4 groups primarily based on geography: Northern China, Europe, Central China, and Southern China. Moreover, the coefficient of similarity when comparing cannabis from Northern China with the European group cannabis was higher than that when comparing with cannabis from the other two groups, owing to a similar climate. This study outlines the first large-scale development of SSR markers for cannabis. These data may serve as a foundation for the development of genetic linkage, quantitative trait loci mapping, and marker-assisted breeding of cannabis.

Phytochemical profiles and antioxidant activities in six species of ramie leaves.

Increased consumption of vegetables or plant food has been associated with decreased risk of developing major chronic diseases, such as cancers, diabetes, cardiovascular diseases, and age-related functional decline. Ramie leaves are rich in phenolics and flavonoids, which have been suggested for human health benefits. Phenolic contents, flavonoid contents, phenolic compounds, and anti-cancer properties in six species of ramie leaves were analyzed by Folin-reagent method, sodium borohydride/chloranil-based assay (SBC), HPLC method and antiproliferation, cytoxicity, respectively. Antioxidant activities were measured through peroxyl radical scavenging capacity (PSC) method, oxygen radical absorbance capacity (ORAC) method, and cellular antioxidant activity (CAA). Research indicated that Boehmeria penduliflora contained the highest total phenolic content (2313.7±27.28 mg GAE/100 g FW), and flavonoid content (1682.4±27.70 mg CAE/100 g FW). Boehmeria tricuspis showed the highest PSC value (9574.8±117.63 µM vit. C equiv./100 g FW), while Boehmeria penduliflora indicated the highest ORAC value (330.44±16.88 µmol Trolox equiv./g FW). The antioxidant activities were correlated with phenolic contents and flavonoid contents. Boehmeria tricuspis had the highest antiproliferative capacity with the lowest EC50 (4.11±0.19 mg/mL). The results for the analyzed ramie for CAA were significantly different from each other (p<0.05), Boehmeria tricuspis had the highest CAA value (133.63±7.10 µmol QE/100 g). Benzoic acid, 4-coumaric acid, caffeic acid, and ferulic acid were the dominant phenolic ingredients in the ramie leaves according to HPLC analysis. Our research is the first report to study the phytochemical profiles and antioxidant activities in different species of ramie leaves for their health benefit.