Exploring the Transcriptome of Non-model Oleaginous Microalga Dunaliella Tertiolecta Through High-throughput Sequencing and High Performance Computing

Overview

Journal BMC Bioinformatics

Publisher Biomed Central

Specialty Biology

Date 2017 Feb 24

PMID 28228091

Citations 3

Authors

Lina Yao

Kenneth Wei Min Tan

Tin Wee Tan

Yuan Kun Lee

Affiliations

Soon will be listed here.

Abstract

Background: RNA-Seq technology has received a lot of attention in recent years for microalgal global transcriptomic profiling. It is widely used in transcriptome-wide analysis of gene expression., particularly for microalgal strains with potential as biofuel sources. However, insufficient genomic or transcriptomic information of non-model microalgae has limited the understanding of their regulatory mechanisms and hampered genetic manipulation to enhance biofuel production. As such, an optimal microalgal transcriptomic database construction is a subject of urgent investigation.

Results: Dunaliella tertiolecta, a non-model oleaginous microalgal species, was sequenced via Illumina MISEQ and HISEQ 4000 in RNA-Seq studies. The high quality high-throughout sequencing data were explored using high performance computing (HPC) in a petascale data center and subjected to de novo assembly and parallelized mpiBLASTX search with multiple species. As a result, a transcriptome database of 17,845 was constructed (~95% completeness). This enlarged database constructed fueled the RNA-Seq data analysis, which was validated by a nitrogen deprivation (ND) study that induces triacylglycerol (TAG) production.

Conclusions: The new paralleled assembly and annotation method under HPC presented here allows the solution of large-scale data processing problems in acceptable computation time. There is significant increase in the number of transcriptomic data achieved and observable heterogeneity in the performance to identify differentially expressed genes in the ND treatment paradigm. The results provide new insights as to how response to ND treatment in microalgae is regulated. ND analyses highlight the advantages of this database generated in this study that could also serve as a useful resource for future gene manipulation and transcriptome-wide analysis. We thus demonstrate the usefulness of exploring the transcriptome as an informative platform for functional studies and genetic manipulations in similar species.

Citing Articles

The Time-Resolved Salt Stress Response of -A Comprehensive System Biology Perspective.

Keil L, Mehlmer N, Cavelius P, Garbe D, Haack M, Ritz M Int J Mol Sci. 2023; 24(20).

PMID: 37895054 PMC: 10607294. DOI: 10.3390/ijms242015374.

Transcriptomic profiles of Dunaliella salina in response to hypersaline stress.

He Q, Lin Y, Tan H, Zhou Y, Wen Y, Gan J BMC Genomics. 2020; 21(1):115.

PMID: 32013861 PMC: 6998148. DOI: 10.1186/s12864-020-6507-2.

Evolutionary and genomic analysis of the caleosin/peroxygenase (CLO/PXG) gene/protein families in the Viridiplantae.

Rahman F, Hassan M, Rosli R, Almousally I, Hanano A, Murphy D PLoS One. 2018; 13(5):e0196669.

PMID: 29771926 PMC: 5957377. DOI: 10.1371/journal.pone.0196669.

References

Thompson J, Gibson T, Higgins D . Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2008; Chapter 2:Unit 2.3. DOI: 10.1002/0471250953.bi0203s00. View

Wang Z, Gerstein M, Snyder M . RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2008; 10(1):57-63. PMC: 2949280. DOI: 10.1038/nrg2484. View

White O, Soderlund C, Shanmugan P, Fields C . Information contents and dinucleotide compositions of plant intron sequences vary with evolutionary origin. Plant Mol Biol. 1992; 19(6):1057-64. DOI: 10.1007/BF00040537. View

Li R, Yu C, Li Y, Lam T, Yiu S, Kristiansen K . SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009; 25(15):1966-7. DOI: 10.1093/bioinformatics/btp336. View

Urzica E, Vieler A, Hong-Hermesdorf A, Page M, Casero D, Gallaher S . Remodeling of membrane lipids in iron-starved Chlamydomonas. J Biol Chem. 2013; 288(42):30246-30258. PMC: 3798491. DOI: 10.1074/jbc.M113.490425. View

Leng N, Dawson J, Thomson J, Ruotti V, Rissman A, Smits B . EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics. 2013; 29(8):1035-43. PMC: 3624807. DOI: 10.1093/bioinformatics/btt087. View

Cloonan N, Forrest A, Kolle G, Gardiner B, Faulkner G, Brown M . Stem cell transcriptome profiling via massive-scale mRNA sequencing. Nat Methods. 2008; 5(7):613-9. DOI: 10.1038/nmeth.1223. View

Yang Z, Niu Y, Ma Y, Xue J, Zhang M, Yang W . Molecular and cellular mechanisms of neutral lipid accumulation in diatom following nitrogen deprivation. Biotechnol Biofuels. 2013; 6(1):67. PMC: 3662598. DOI: 10.1186/1754-6834-6-67. View

Li B, Dewey C . RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011; 12:323. PMC: 3163565. DOI: 10.1186/1471-2105-12-323. View

10.

Huang T, Shi X, Wang P, He Z, Feng K, Hu L . Analysis and prediction of the metabolic stability of proteins based on their sequential features, subcellular locations and interaction networks. PLoS One. 2010; 5(6):e10972. PMC: 2881046. DOI: 10.1371/journal.pone.0010972. View

11.

Mouchiroud D, DOnofrio G, Aissani B, Macaya G, Gautier C, Bernardi G . The distribution of genes in the human genome. Gene. 1991; 100:181-7. DOI: 10.1016/0378-1119(91)90364-h. View

12.

Merchant S, Prochnik S, Vallon O, Harris E, Karpowicz S, Witman G . The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science. 2007; 318(5848):245-50. PMC: 2875087. DOI: 10.1126/science.1143609. View

13.

Martin J, Wang Z . Next-generation transcriptome assembly. Nat Rev Genet. 2011; 12(10):671-82. DOI: 10.1038/nrg3068. View

14.

Zerbino D, Birney E . Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008; 18(5):821-9. PMC: 2336801. DOI: 10.1101/gr.074492.107. View

15.

Yao L, Tan T, Ng Y, Ban K, Shen H, Lin H . RNA-Seq transcriptomic analysis with Bag2D software identifies key pathways enhancing lipid yield in a high lipid-producing mutant of the non-model green alga Dunaliella tertiolecta. Biotechnol Biofuels. 2015; 8:191. PMC: 4660794. DOI: 10.1186/s13068-015-0382-0. View

16.

Jabbari K, Bernardi G . CpG doublets, CpG islands and Alu repeats in long human DNA sequences from different isochore families. Gene. 1999; 224(1-2):123-7. DOI: 10.1016/s0378-1119(98)00474-0. View

17.

Goodall G, Filipowicz W . Different effects of intron nucleotide composition and secondary structure on pre-mRNA splicing in monocot and dicot plants. EMBO J. 1991; 10(9):2635-44. PMC: 452964. DOI: 10.1002/j.1460-2075.1991.tb07806.x. View

18.

Blanc G, Agarkova I, Grimwood J, Kuo A, Brueggeman A, Dunigan D . The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation. Genome Biol. 2012; 13(5):R39. PMC: 3446292. DOI: 10.1186/gb-2012-13-5-r39. View

19.

Reddy A . Alternative splicing of pre-messenger RNAs in plants in the genomic era. Annu Rev Plant Biol. 2007; 58:267-94. DOI: 10.1146/annurev.arplant.58.032806.103754. View

20.

Rismani-Yazdi H, Haznedaroglu B, Bibby K, Peccia J . Transcriptome sequencing and annotation of the microalgae Dunaliella tertiolecta: pathway description and gene discovery for production of next-generation biofuels. BMC Genomics. 2011; 12:148. PMC: 3061936. DOI: 10.1186/1471-2164-12-148. View