<figure><img src="image_1.png" /></figure>

HAL
open science

# Transcriptomic profiling of sesame during waterlogging and recovery

Komivi Dossa, Jun You, Linhai Wang, Yanxin Zhang, Donghua Li, Rong

Zhou, Jingyin Yu, Xin Wei, Xiaodong Zhu, Shiyang Jiang, et al.

► To cite this version:

Komivi Dossa, Jun You, Linhai Wang, Yanxin Zhang, Donghua Li, et al. Transcriptomic profiling of sesame during waterlogging and recovery. Scientific Data, 2019, 6, pp.204. 10.1038/s41597-019-0226-z. hal-05181641

HAL Id: hal-05181641

https://hal.science/hal-05181641v1

Submitted on 23 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC DATA 110110
0111101
11011110
011101101

OPEN

DATA DESCRIPTOR

Received: 21 February 2019

Accepted: 21 August 2019

Published online: 15 October 2019

Transcriptomic profiling of sesame
during waterlogging and recovery

Komivi Dossa<fnref n="1" />, Jun You<fnref n="1" />, Linhai Wang<fnref n="1" />, Yanxin Zhang<fnref n="1" />, Donghua Li<fnref n="1" />, Rong Zhou<fnref n="1" />, Jingyin Yu<fnref n="1" />, Xin Wei<fnref n="1" />,², Xiaodong Zhu<fnref n="1" />, Shiyang Jiang<fnref n="1" />, Yuan Gao<fnref n="1" />, Marie A. Mmadi<fnref n="1" /> & Xiurong Zhang<fnref n="1" />

Sesame is naturally adapted to arid environments but highly susceptible to waterlogging stress. A few hours of waterlogging (lasting over 36 h) are detrimental to the crop growth, yield and survival. To better understand the molecular mechanisms underlying sesame responses to waterlogging and recovery, it is essential to design a high-resolution time-series experiment. In this study, we reported the RNA-seq profiling of two contrasting genotypes under waterlogging and recovery. The plants were grown in pots and subjected to waterlogging treatment at the flowering stage for 36 h and subsequently, 12 h drainage. Root samples were collected in triplicate at 22 time points under waterlogging/drainage treatments and at 10 time points in the control condition. This represents a total of 195 biological samples and the RNA-seq yielded over eight billion reads. Basic data analyses demonstrated a clear separation of transcriptomes from control, waterlogging and drainage treatments. Overall, the generated high-quality and comprehensive RNA-seq resources will undoubtedly advance our understanding of waterlogging/drainage responses in a non-model sensitive crop.

Background & Summary

With the recent surge in flooding events in many regions of the world, waterlogging has become a serious problem for agricultural production<fnref n="1" />. Soil waterlogging leads to hypoxia, high CO₂ in the root zone and slows down photosynthesis, which ultimately impairs normal crop growth and yields. Like most of typical dryland crops, sesame (*Sesamum indicum* L.) is highly sensitive to waterlogging stress². Waterlogging occurring in several sesame growing areas in Asia and Africa engenders catastrophic economic loss for smallholders. According to Sun *et al*.³, most of the sesame cultivars hardly survive over 36 h of waterlogging in the fields. Previous studies that analyzed sesame transcriptional responses to waterlogging and drainage were conducted with limited temporal resolution⁴,⁵, hence, the architecture and dynamics of the waterlogging/recovery gene regulatory network are yet to be elucidated. The present study generated high-quality and high-resolution time series RNA sequencing data (195 RNA-seq in total) from root of two contrasting sesame cultivars (ZZM2541 and Ezhi-2) during the waterlogging and recovery stages. We believe that these precious resources from a non-model crop will help unlock novel genes-pathways-mechanisms modulating waterlogging/drainage responses and assist in crop improvement strategies.

Methods

**Plant materials and stress treatment.** Two genotypes of sesame (*Sesamum indicum* L.) were obtained from the China National Genebank, Oil Crops Research Institute, Chinese Academy of Agricultural Sciences and used in this experiment. The genotype ZZM2541 (R2G) displays a strong tolerance to waterlogging stress while Ezhi No. 2 (EG) is highly susceptible as demonstrated by Wei *et al*.². The experiment was conducted in a greenhouse as described by Wang *et al*.⁵ and Dossa *et al*.⁶. Plants were grown in pots (25 cm diameter and 30 cm depth) containing 7 Kg of loam soil mixed with 10% compound fertilizer. The plants were irrigated every 3 days and the soil volumetric water content (vwc) was maintained at ~35%. A pot tray was placed under each pot to avoid water loss. A completely randomized blocking design with 3 replicates was employed. 15 days after the initiation of flowering, half of the pots were waterlogged by standing in a plastic bucket filled with tap water up to 3 cm above the soil surface. Each pot contains 3 seedlings, which were maintained waterlogged for 36 h and afterwards, pots were drained to allow plants to recover for 12 h (vwc = ~35%). In parallel, half of the pots were kept under normal

<footnote n="1">¹Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, No.2 Xudong 2nd Road, Wuhan, 430062, China. ²College of Life Sciences, Shanghai Normal University, Shanghai, 200234, China. These authors contributed equally: Komivi Dossa, Jun You and Linhai Wang. Correspondence and requests for materials should be addressed to Xiurong Zhang (email: zhangxr@oilcrops.cn)</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Flowchart of experimental design of this study. A high-resolution temporal waterlogging/drainage treatment was used to capture the gene expression changes in roots of two *S. indicum* genotypes (EG and R2G). The untreated plants were used as controls. Three biological replicates per condition were sampled for transcriptome sequencing. All raw reads were quality controlled prior to aligning to *S. indicum* reference genome (v1.0). The uniquely aligned reads were employed for expression profile analyses.</figcaption></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 2** RNA sequencing quality details. An example of the FastQC report illustrating the average quality scores across all bases of the paired-end datasets for the sample CK_EG_12_rep1 (SRR8490161). (a,b) Phred quality scores for each nucleotide position are represented as a box and whisker plot. The central red line is the median value. The yellow box represents the inter-quartile range (25–75%). The upper and lower whiskers represent the 10 and 90% points. The blue line represents the mean quality. (c,d) Average quality per read along the reads.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 3** Heatmap clustering of the transcriptome data. Correlations were estimated using the Pearson correlation coefficient based on TPM values. In total two main groups of samples were obtained, including control and stressed samples.</figcaption></figure>

growth conditions (vwc = ~35%) during the whole experiment. Root samples were collected from a single plant/ pot from the three replicated pots (3 biological replicates) in the stress and control treatments at the different time points following the flowchart presented in Fig. 1. After 48 h treatment, several EG plants were dead while few survived. We therefore sampled separately both dead (EG-48-W) and survived (EG-48) EG plants. In total, 195 root samples were collected and snap frozen in liquid nitrogen for follow-up analyses (Sample information for the study) available at figshare⁷.

**RNA extraction, library preparation and RNA sequencing.** Total RNA was extracted from 195 root samples using the TRIzol reagent (Invitrogen), and treated with DNase I and Oligo (dT) to isolate mRNAs. The concentration and quality were determined using an ultraviolet spectrophotometer and 2% denaturing agarose gels. The cDNA was synthesized using the mRNA fragments as templates. The short fragments (~300 bp) were ligated with adapters and the suitable fragments were selected for PCR amplification. The libraries were paired-end sequenced using the Illumina platform Hiseq 2500⁸.

**RNA-seq data processing.** The program FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was employed to determine the base quality of the raw reads (in FASTQ format) and we removed the paired-end reads containing more than 5% ambiguous residues (Ns) and those containing >10% bases with a Phred quality score of 10 (Fig. 2). Then, the raw reads were trimmed using Trimmomatic, version 0.32⁹. After cleaning and quality reassessment with FastQC, approximately 31.6–56.8 million high-quality reads of 90-bp length remained in each sample (Quality check report of 195 sesame transcriptomes under control, waterlogging and drainage treatments) available at figshare⁷. The high-quality reads were mapped to the sesame (*Sesamum indicum* L.) reference genome v1.0 (http://ocri-genomics.org/Sinbase/login.htm)¹⁰ using the STAR software¹¹, allowing no more than one mis-match in the alignment. Approximately, 88.9–97.4% of the clean reads were uniquely
<figure><img src="image_9.png" /><figcaption>**Fig. 4** Tsne scatter plot of the samples. The plot depicts the clustering patterns of the samples according to the genotypes and stress treatments. In the diagram, the transverse ordinates represent the first and second principal components; the symbols in the graph represent the samples, and the different colors represent the time points of sampling.</figcaption></figure>

mapped to the reference genome, with 94.3–98.6% of them uniquely mapped to the genic regions (Statistics of the clean read mapping) available at figshare<sup>7</sup>. Using the featureCount package<sup>12</sup>, the gene expression levels were calculated based on the number of unique matched reads to the sesame genome v1.0<sup>10</sup> and were normalized to Transcripts Per Million (TPM).

## Data Records

The RNA-seq raw data of the 195 samples are deposited to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession SRP181800<sup>13</sup>. The gene expression data of the 195 RNA-seq samples are deposited to the Gene Expression Omnibus under the accession GSE133186<sup>14</sup>. Supplementary files accompanying this manuscript have been deposited to figshare<sup>7</sup>.

## Technical Validation

**Quality control.** In this project, a total of 195 RNA libraries were prepared and sequenced generating over 8 billion raw reads (Quality check report of 195 sesame transcriptomes under control, waterlogging and drainage treatments) available at figshare<sup>7</sup>. By applying FastQC to the whole dataset, we successfully obtained high quality clean data with 97% of bases scoring Q30 and above (Fig. 2). Approximately, 98% of reads could be uniquely mapped to the reference genome of *S. indicum*.

**Basic analysis of RNA-seq data.** A heatmap for cluster relationships among the samples representing Pearson distance was generated with the TPM values using the R package 'Plothmap' v.1.0.12 (https://cran.r-project.org/web/packages/plothmap/index.html) (Fig. 3). It could be observed that samples from the control conditions (CK) and those from the stress treatments formed distinct groups. To further assess sample relationships and the time-course expression patterns, we implemented the tsne reduction scatter plot, which is a non-linear dimensionality reduction method for embedding high dimensional data into a low-dimensional space<sup>15</sup>. The analysis was performed with the R package 'tsne' v.0.1–3. The results showed a clear separation between samples from control (CK), waterlogging treatment (0-36 h) and drainage phase (36–48 h) (Fig. 4). In addition, the difference between the two contrasting genotypes used in this study could be observed through their differential gene expression during the waterlogging treatment.

## Usage Notes

This study generated dense and high-resolution gene expression data under waterlogging stress-drainage in a non-model crop. In contrast to many experimental designs relating to temporal transcriptome profiling under stress which set the beginning of the stress application as the control, in the present project, we collected samples at different time points under the control condition. Comparing expression data from samples harvested at the same time point under stress and control conditions will provide more accurate differential gene expression records. Concerning the time points for which samples were uniquely harvested under stress treatments, the users can still use the previous time point in the control condition for gene differential expression analysis. The whole datasets were publicly deposited at NCBI SRA and we anticipate that our RNA-seq data would provide new insights into the molecular basis of waterlogging stress responses in plants.

## Code Availability

Codes that were used for the RNA-seq data processing are available at figshare<sup>7</sup>. Software and their versions were described in Methods.
## References

1. Hirabayashi, Y. et al. Global flood risk under climate change. *Nat. Climate Change* **3**, 816–821 (2013).

2. Wei, W. et al. Morpho-anatomical and physiological responses to waterlogging of sesame (*Sesamum indicum* L.). *Plant Sci.* **208**, 102–111 (2013).

3. Sun, J., Zhang, X., Zhang, Y., Wang, L. & Huang, B. Effects of waterlogging on leaf protective enzyme activities and seed yield of sesame at different growth stages. *Chin. J. Appl. Environ. Biol.* **15**, 790–795 (2009).

4. Wang, L. et al. Global gene expression responses to waterlogging in roots of sesame (*Sesamum indicum* L.). *Acta Physiol. Plant.* **34**, 2241 (2012).

5. Wang, L. et al. Tolerant and susceptible sesame genotypes reveal waterlogging stress response patterns. *PloS ONE* **11**, e0149912 (2016).

6. Dossa, K. et al. Transcriptomic biochemical and physio-anatomical investigations shed more light on responses to drought stress in two contrasting sesame genotypes. *Sci. Rep.* **7**, 8755 (2017).

7. Dossa, K. et al. A new arsenal for deciphering waterlogging and recovery response, 195 high-resolution RNA-seq data from sesame. figshare. https://doi.org/10.6084/m9.fishshare.c.4551407 (2019).

8. Dossa, K. et al. The contrasting response to drought and waterlogging is underpinned by divergent DNA methylation programs associated with transcript accumulation in sesame. *Plant Sci.* **277**, 207–217 (2018).

9. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).

10. Wang, L. et al. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. *Genome Biol.* **15**, R39 (2014).

11. Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. *Bioinformatics* **29**, 15–21 (2013).

12. Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. *Bioinformatics* **30**, 923–930 (2014).

13. NCBI Sequence Read Archive, https://identifiers.org/ncbi/insdc.sra:SRP181800 (2019).

14. Dossa, K. et al. A new arsenal for deciphering waterlogging/drainage responses: 195 dense time-series RNA-seq dataset from sesame. *Gene Expression Omnibus*, https://identifiers.org/geo:GSE133186 (2019).

15. Van der Maaten, L. J. P. & Hinton, G. E. Visualizing data using t-SNE. *J. Mach. Learn. Res.* **9**, 2579–2605 (2008).

## Acknowledgements

The study was supported by the China Agriculture Research System (CARS-14), the Agricultural Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2013-OCRI), Wuhan Cutting-edge Application Technology Fund (2018020401011303), and the Fundamental Research Funds for Central Non-profit Scientific Institution (1610172017003).

## Author Contributions

X.Z. and L.W. designed the project. K.D., J.Y., L.W., Y.Z., D.L., R.Z., J.Y., X.W., X.Z., S.J., Y.G., M.A.M., X.Z. executed the experiment. K.D., J.Y., L.W. analyzed the results. K.D. and J.Y. wrote the manuscript. All authors contributed to reading and editing the manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

© The Author(s) 2019