Characterization of the Kenaf (Hibiscus cannabinus) Global Transcriptome Using Illumina Paired-End Sequencing and Development of EST-SSR Markers

Development of Roselle (Hibiscus sabdariffa L.) Transcriptome-Based Simple Sequence Repeat Markers and Their Application in Roselle

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.

Transcriptome profiling of kenaf (Hibiscus cannabinus L.) under plumbic stress conditions implies the involvement of NAC transcription factors regulating reactive oxygen species-dependent programmed cell death

Heavy metal contamination of soils has become a serious global issue, and bioremediation has been proposed as a potential solution. Kenaf (Hibiscus cannabinus L.) is a fast growing, non-woody multipurpose annual plant that is suitable for removing excess heavy metals from soils. However, there has been relatively little research on the kenaf molecular mechanisms induced in response to an exposure to heavy metal stress. Thus, whole kenaf seedlings grown under control (normal) and stress (plumbic treatment) conditions were sampled for transcriptome sequencing. Unigenes generated through the de novo assembly of clean reads were functionally annotated based on seven databases. Transcription factor (TF)-coding genes were predicted and the physiological traits of the seedlings were analyzed. A total of 44.57 Gb high-quality sequencing data were obtained, which were assembled into 136,854 unigenes. These unigenes included 1,697 that were regarded as differentially expressed genes (DEGs). A GO enrichment analysis of the DEGs indicated that many of them are related to catalytic activities. Moreover, the DEGs appeared to suggest that numerous KEGG pathways are suppressed (e.g., the photosynthesis-involving pathways) or enhanced (like the flavonoid metabolism pathways) in response to Pb stress. Of the 2,066 predicted TF-coding genes, only 55 were differentially expressed between the control and stressed samples. Further analyses suggested that the plumbic stress treatment induced reactive oxygen species-dependent programmed cell death in the kenaf plants via a process that may be regulated by the differentially expressed NAC TF genes.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

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.

Transcriptomic profiling of hemp bast fibres at different developmental stages

Bast fibres are long extraxylary cells which mechanically support the phloem and they are divided into xylan- and gelatinous-type, depending on the composition of their secondary cell walls. The former, typical of jute/kenaf bast fibres, are characterized by the presence of xylan and a high degree of lignification, while the latter, found in tension wood, as well as flax, ramie and hemp bast fibres, have a high abundance of crystalline cellulose. During their differentiation, bast fibres undergo specific developmental stages: the cells initially elongate rapidly by intrusive growth, subsequently they cease elongation and start to thicken. The goal of the present study is to provide a transcriptomic close-up of the key events accompanying bast fibre development in textile hemp (Cannabis sativa L.), a fibre crop of great importance. Bast fibres have been sampled from different stem regions. The developmental stages corresponding to active elongation and cell wall thickening have been studied using RNA-Seq. The results show that the fibres sampled at each stem region are characterized by a specific transcriptomic signature and that the major changes in cell wall-related processes take place at the internode containing the snap point. The data generated also identify several interesting candidates for future functional analysis.