Next Article in Journal
Monitoring the Emergence of Algal Toxins in Shellfish: First Report on Detection of Brevetoxins in French Mediterranean Mussels
Next Article in Special Issue
Screening of Neutrophil Activating Factors from a Metagenome Library of Sponge-Associated Bacteria
Previous Article in Journal
Blue Biotechnology: Computational Screening of Sarcophyton Cembranoid Diterpenes for SARS-CoV-2 Main Protease Inhibition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Marine Drugs
Volume 19
Issue 7
10.3390/md19070392
marinedrugs-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes

by
Haomiao Cheng
1,2,3,
Chris Bowler
4,
Xiaohui Xing
5,6,7,
Vincent Bulone
5,6,7,
Zhanru Shao
1,2,* and
Delin Duan
1,2,8,*
1
CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
Laboratory for Marine Biology and Biotechnology, Pilot Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
Institut de Biologie de l’ENS (IBENS), Département de Biologie, École Normale Supérieure, CNRS, INSERM, Université PSL, 75005 Paris, France
5
Division of Glycoscience, Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology (KTH), AlbaNova University Centre, 10691 Stockholm, Sweden
6
Australian Research Council Centre of Excellence in Plant Cell Walls, School of Agriculture, Food and Wine, University of Adelaide, Waite Campus, Urrbrae 5064, Australia
7
Adelaide Glycomics, School of Agriculture Food and Wine, University of Adelaide, Waite Campus, Urrbrae 5064, Australia
8
State Key Laboratory of Bioactive Seaweed Substances, Qingdao Bright Moon Seaweed Group Co., Ltd., Qingdao 266400, China
*
Authors to whom correspondence should be addressed.
Mar. Drugs 2021, 19(7), 392; https://doi.org/10.3390/md19070392
Submission received: 10 June 2021 / Revised: 6 July 2021 / Accepted: 8 July 2021 / Published: 13 July 2021
(This article belongs to the Special Issue Genetic Approach and Data Mining in Discovery of Marine Natural Products)
Browse Figures
Versions Notes

Abstract

:
β-Chitin produced by diatoms is expected to have significant economic and ecological value due to its structure, which consists of parallel chains of chitin, its properties and the high abundance of diatoms. Nevertheless, few studies have functionally characterised chitin-related genes in diatoms owing to the lack of omics-based information. In this study, we first compared the chitin content of three representative Thalassiosira species. Cell wall glycosidic linkage analysis and chitin/chitosan staining assays showed that Thalassiosira weissflogii was an appropriate candidate chitin producer. A full-length (FL) transcriptome of T. weissflogii was obtained via PacBio sequencing. In total, the FL transcriptome comprised 23,362 annotated unigenes, 710 long non-coding RNAs (lncRNAs), 363 transcription factors (TFs), 3113 alternative splicing (AS) events and 3295 simple sequence repeats (SSRs). More specifically, 234 genes related to chitin metabolism were identified and the complete biosynthetic pathways of chitin and chitosan were explored. The information presented here will facilitate T. weissflogii molecular research and the exploitation of β-chitin-derived high-value enzymes and products.

Graphical Abstract

1. Introduction

Chitin, a polymer of 1,4-linked β-D-N-acetylglucosaminyl residues, is the second most abundant natural biopolymer after cellulose, and is widely distributed across taxa [1,2,3]. Together with its partially de-N-acetylated derivative chitosan, chitin has many biomedical applications, such as wound healing, artificial organs and drug delivery [4]. The parallel arrangement of chitin chains in the β-chitin allomorph confers specific properties to the polymer, such as higher solubility, reactivity and swelling compared to the most abundant form of chitin, α-chitin (Figure 1) [5,6]. There are very few classes of organisms that produce β-chitin, such as the diatom Thalassiosira sp. [7,8,9]. Various diatoms have been previously reported to produce chitin, mainly in the genera Thalassiosira and Cyclotella [10,11]. The first report of the occurrence of chitin in diatoms was in T. weissflogii (fluviatilis), which showed that T. weissflogii chitin represented up to 34% of the total cell mass (including the silica) [12]. The content and structure of chitin have been widely studied in diatoms, but its biosynthetic pathway based on the sequencing and gene annotation is incomplete [13,14]. In this study, we first measured the chitin content in the cell walls of three diatom species, Thalassiosira rotula, T. pseudonana and T. weissflogii, and found that T. weissflogii had the highest chitin content. However, apart from the identification of chitinous fibres [9], our knowledge of chitin in T. weissflogii is limited, which impedes the exploitation of chitin and chitin-related enzymes from this class of organisms.
Diatoms are a major group of phytoplankton and contribute approximately 20% of global primary productivity [15,16]. The last two decades have witnessed a surge in diatom gene information, with complete genomes [17,18,19] and a number of transcriptomes sequenced [20,21,22]. By contrast, T. weissflogii genes can only be referred to the transcripts released by the Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) [23], and no full-length reference database is openly available so far. Genes encoding enzymes metabolising high-value products such as β-chitin and the full potential of T. weissflogii remain to be unveiled. Pacific BioSciences (PacBio) single-molecule real-time (SMRT) sequencing enables the discovery of novel genes and gene product isoforms. It provides longer reads, making them suitable for biological problems that are poorly solved by second-generation sequencing (SGS) [24]. A full-length (FL) transcriptome obtained by PacBio sequencing is an alternative complete transcript assembly for a non-model species. In the present study, the FL transcriptome of T. weissflogii was sequenced to elucidate the genetic profile and uncover the chitin-related genes of this diatom species. Our study highlights information for β-chitin metabolism in T. weissflogii, provides a reference genetic dataset for future studies and assists in evolutionary studies to a broad taxonomy. The FL transcriptome will enable the construction of chitin metabolic pathways and facilitate the in vitro application of chitin-related enzymes.

2. Results and Discussion

2.1. Glycosidic Linkage Analysis

Three representative species of the Thalassiosira genus were selected and subjected to glycosidic linkage analysis of cell wall polysaccharides (Table 1). Glucosyl (Glc) and N-acetylglucosaminyl (GlcNAc) residues were the two most abundant monosaccharides in all three Thalassiosira species analysed. Notably, T. weissflogii contained more GlcNAc than the other two species, representing 11.5 and 4.9 times that in T. pseudonana and T. rotula, respectively (Table 1). This indicates that T. weissflogii might be the most appropriate candidate of all three species analysed for exploring chitin-related enzymes in diatoms.

2.2. Staining of Chitin and Chitosan

In order to observe the localisation of chitin and chitosan, chitin-binding protein (CBP) and chitosan-affinity protein (CAP) tagged with enhanced green fluorescence protein (eGFP) were used to stain T. weissflogii live cells. The results showed that the fluorescence of chitin was present in the cytoplasm and at the cell boundary (Figure 2A). Notably, strong, clumped fluorescent signals were detected in the cell suspension (Figure 2B), which indicated that a large quantity of chitin detached from the cell and aggregated into clusters. This is the evidence that T. weissflogii could produce a large amount of external chitin microfibrils, which was consistent with the scanning electron microscopy-derived results of Ogawa et al. (2011) [9]. After CAP-eGFP staining, continuous fluorescence appeared around the cells (Figure 2C), which is a new discovery indicating that Thalassiosira can synthesise chitosan in the cell wall.

2.3. Sequencing and Data Processing

As a source of high chitin content, T. weissflogii was subjected to PacBio sequencing. A total of 70,038,125,361 basepairs (bp) containing 44,233,932 subreads were yielded (Table 2). In total, 1,021,310 circular consensus sequences (CCSs) with an average length of 2127 bp were generated after merging transcripts with at least two full passes. The full-length non-chimeric (FLNC) sequences of the CCSs were then clustered and polished, producing 110,527 high-quality isoforms and 338 low-quality isoforms. After removing the sequence redundancy, the polished high-quality isoforms were trimmed to 25,412 unigenes for further analyses (Table 2). Only 25,412 (2.5%) of the 1,021,310 CCSs were retained as unigenes after redundancy removal, which was less than the transcript numbers of the two T. weissflogii strains released in the MMETSP project (55,443 of strain CCMP1336 and 282,372 of strain CCMP1010) [23]. The significant level of redundancy showed the good depth and integrity of the T. weissflogii FL transcriptome achieved via PacBio sequencing technology. The GC content of the T. weissflogii FL transcriptome was comparable to those of T. rotula and three Pseudo-nitzschia species [25,26], indicating high homogeneity within diatoms. A relatively high percentage of complete Benchmarking Universal Single-Copy Orthologs (BUSCOs), 200 out of 303 total BUSCOs (66.01%), showed the high-quality assembly completeness of our transcriptomes (Figure 3).
Transcriptomes reflect gene expression potentially affecting the physiological and biochemical processes from a molecular perspective. In diatoms, SGS has revealed genetic information in terms of underwater adhesion, biofuel accumulation and nitric oxide synthase genes in Amphora coffeaeformis, Fistulifera solaris and Pseudo-nitzschia [26,27,28]. The PacBio sequencing captures full-length transcripts without assembly and overcomes the limitations of SGS, such as complex genomic region assembly and determination, isoform and methylation detection [24]. For a non-model diatom without a reference genome, SGS sequencing is inadequate for producing an FL transcriptome. PacBio sequencing has been used to analyse the FL transcriptome of animals and higher plants [29,30]. The T. weissflogii FL transcriptome here represents the first general transcription encyclopaedia of the species and FL transcriptome from a diatom. The average read length of the T. weissflogii FL transcriptome obtained by PacBio sequencing was longer than those of other diatom species acquired by SGS [20,27]. Furthermore, compared with the SGS transcriptome of T. rotula, the T. weissflogii FL transcriptome contained higher N50 and longer transcripts, and retained much fewer genes [25].

2.4. Analyses of Coding Sequence, Long Non-Coding RNAs and Transcription Factors

The coding sequence (CDS) of a gene is a singular section of DNA or RNA that encodes the corresponding protein. In 25,412 unigenes with an average length of 2045.04 bp (Figure 4A), a total of 24,500 (96.4%) CDSs were predicted, with the most represented length range being 401–600 bp (Figure 4B). LncRNAs are defined as transcripts longer than 200 nucleotides that are not translated into protein. They are fundamentally involved in biological processes such as transcription, translation, protein localisation, cellular structure integrity, reprogramming and other cellular activities [31]. Altogether, 710 unigenes were characterised as lncRNAs by the joint prediction of CNCI and CPC (Figure 4C). TFs are proteins that recognise and bind specific nucleotide sequences, mediating the primary expression of nearby genes during transcription [32]. A total number of 363 TF unigenes were identified and categorised into 17 families, the top ten of which are shown in Figure 4D. Over half of the TFs were members of the heat shock transcription factor (HSF) family (187, 51.5%). The proportion of HSF in the TFs of T. weissflogii was even higher than that reported in the Phaeodactylum tricornutum and T. pseudonana TFs, where the HSF family constituted the most abundant TF class (33.0% and 36.4%, respectively) [33]. The number of TF genes identified in this FL transcriptome was higher than those reported in the P. tricornutum and Nitzschia SGS transcriptomes, although the TFs of Nitzschia were sorted into 38 families [20,34].

2.5. Annotation Analyses

A total of 23,362 unigenes were annotated (91.9%), of which 7564 could be annotated in all four of the Non-Redundant Protein Sequence Database (NR), Swiss-Prot, euKaryotic Ortholog Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The number of unigenes assigned annotation terms in the above four databases were 23,340, 12,489, 10,913 and 8624, respectively (Figure 5A). According to the prediction by the NR database, 23,340 unigenes were annotated in 189 homologous species, where the top ten species are shown in Figure 5B. T. pseudonana, whose genome was the first to be described for a eukaryotic marine phytoplankton [17], was found to contain the highest number of homologous sequences (9587, 41.1%) with T. weissflogii. The second was a pennate diatom Fragilariopsis cylindrus and the third was the centric diatom Thalassiosira oceanica. This suggested that T. weissflogii might be phylogenetically more closely related to F. cylindrus than its centric companion. Two macroalgal species, Ectocarpus siliculosus and Klebsormidium flaccidum, were ranked among the top ten species.
From the KOG analysis, 10,913 unigenes were categorised into 25 classes; the most prominent class was “posttranslational modification, protein turnover, chaperones” (Figure 5C). A total of 8624 unigenes were annotated in 127 pathways in the KEGG database, with five A classes and nineteen B classes (Figure 5D). The B classes “carbohydrate metabolism” and “amino acid metabolism” were at the second and third places, which demonstrated the pivotal roles of these relevant genes for T. weissflogii. Within marine phytoplankton, diatoms outcompeted others for ample NO3, which fuels CO2 fixation and the proliferation of diatoms in nutrient-rich environments such as upwellings [35,36]. These relevant genes might enable diatoms to rapidly respond to variations in carbon and nutrient sources and contribute to their ecological success.

2.6. Analyses of Alternative Splicing and Simple Sequence Repeats

The AS process provides eukaryotes with peculiarly versatile means of genetic regulation. Splicing site alterations from a gene generates multiple mRNA and protein products [37]. In the present study, 3113 AS events were detected, with the genes containing two isoforms ranked the highest (1105) (Figure 6A). Only 189 events were classified into four AS types (Figure 6B), and the major AS type was the retained intron, which was previously reported in P. tricornutum [38]. No exonskipping, mutually exclusive exons and alternative first exon AS events were detected. However, the detection of AS events in this study was limited due to the lack of a reference genome of T. weissflogii. This would result in the missing and underestimation of some types of AS [29]. An SSR or microsatellite is a repetitive DNA sequence where certain motifs are repeated. A total of 3295 SSRs were detected (Figure 6C), exceeding the SSRs identified in P. tricornutum on numbers (1135 in Phatr2 and 255 in Phatr3) [38]. In these SSRs, the largest group was trinucleotide repeats (2754, 83.6%), most of which were 4–7 repeats (Figure 6D). SSRs are highly polymorphic genetic markers [39], and the information here will be of convenience for phylogenetic studies of diatoms.

2.7. Chitin-Related Gene Mining

Transcriptome profiling has revealed information of chitin-related genes in many organisms, e.g., a chitin elicitor receptor kinase gene in barley [40], chitin utilisation-related genes in Vibrio coralliilyticus and Photobacterium galatheae [41] and chitin metabolism-related genes in Glyphodes pyloalis Walker [42]. To date, reports identifying chitin-related genes on whole-genome or -transcriptome scales in diatoms are scarce. Traller et al. (2016) compiled chitin metabolism pathway genes in Cyclotella cryptica [14]. Cheng et al. (2021) characterised a gene family of 24 members encoding chitinases in T. pseudonana [43]. In contrast, the identification of these classes of genes in T. weissflogii at the whole-genome or -transcriptome level has been lacking.
From the FL transcriptome dataset of T. weissflogii, 234 unigenes potentially related to chitin metabolism were identified, including 25 glutamine-fructose-6-phosphate transaminases (isomerising), 2 N-acetylglucosamine-6-phosphate deacetylases, 5 phosphoacetylglucosamine mutases, 5 UDP-N-acetylglucosamine diphosphorylases, 30 chitin synthases, 124 chitinases, 17 beta-N-acetylhexosaminidases, 4 chitin deacetylases and 22 chitin-binding proteins (Figure 7). We constructed the chitin metabolism pathway in T. weissflogii and compared it with pathways in P. tricornutum, T. pseudonana and C. cryptica. We found that all four diatom species harbour rather complete chitin metabolism pathways (Figure 7). The total number of identified chitin-related genes in T. weissflogii, T. pseudonana, C. cryptica and P. tricornutum were 234, 141, 50 and 48, respectively. The higher abundance of these genes in T. weissflogii, particularly the chitin synthase, chitinase and chitin deacetylase genes that are directly related to chitin, implies a more active chitin metabolism in T. weissflogii. As P. tricornutum does not generate chitin fibrils, it is not surprising that it harbours the lowest number and fewest types of chitin-related genes. The N-acetylglucosamine kinase that phosphorylates GlcNAc to GlcNAc-6-P was absent in all four species. We speculate that this may be because diatoms lack the GlcNAc phosphorylation process, rather than due to inconclusive annotations or distinction from characterised N-acetylglucosamine kinases, as Traller et al. (2016) suggested [14]. In T. weissflogii, chitinases accounted for over half of the sum (53.0%), indicating the presence of an active chitin degradation process. However, these genes identified in the FL transcriptome require further sequence analyses for each isoform, which is needed for a better and more accurate understanding of the entire chitin metabolic pathway.

3. Materials and Methods

3.1. Glycosidic Linkage Analysis

T. weissflogii, T. pseudonana and T. rotula cultures were grown at 19 °C with 12 h:12 h light:dark cycles (100 μmol m−2 s−1). Approximately 100 mg lyophilised samples were collected from large-scale cultures of each Thalassiosira species for glycosidic linkage composition analysis. The cell wall sample preparation and linkage analysis (methylation–GC–MS analysis) was performed as described by Shao et al. (2019) [6]. Experiments were conducted in duplicate.

3.2. Staining of Chitin and Chitosan

T. weissflogii cell cultures at exponential period were centrifuged at 13,000× g, prior to a suspension of cell pellets in fresh f/2 liquid medium. CBP-eGFP and CAP-eGFP were used to visualise the presence and localisation of chitin and chitosan, respectively, according to the methods by Hardt and Laine (2004) and Nampally et al. (2012) [44,45]. GFP signal was observed using a Zeiss Axio Imager M2 fluorescent microscope (Zeiss, Germany). ImageJ software was used to analyse the epifluoresent photos.

3.3. T. weissflogii Collection and RNA Extraction

T. weissflogii (9021) for sequencing was acquired from the Microalgae Collection Center at Ningbo University, Ningbo, China, and grown in optimised f/2 liquid medium provided by Shanghai Guangyu Biological Technology Co., Ltd., Shanghai, China. Cells were cultured at 19 °C with 12 h:12 h light:dark cycles (100 μmol m−2 s−1) and swirled at 100 rpm. Cells reaching the exponential phase were harvested and frozen in liquid nitrogen and stored at −80 °C for further experiments. Total RNA was extracted by grinding the T. weissflogii frozen sample in TRIzol reagent (Life Technologies, Carlsbad, CA, USA) and processed following the protocol provided by the manufacturer.

3.4. Library Construction and Sequencing

RNA concentration was checked using a Nanodrop micro-spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and RNA integrity was verified with a Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). mRNA was enriched using Oligo (dT) magnetic beads and then reverse-transcribed into cDNA using the Clontech SMARTer PCR cDNA Synthesis Kit (Clontech, Palo Alto, CA, USA). The cDNA was then amplified by PCR. All cDNAs were DNA damage-repaired, end-repaired and connected with adaptors. PacBio sequencing was then conducted on a PacBio Sequel platform (Gene Denovo Biotechnology Co., Guangzhou, China).

3.5. Data Processing

The raw data were initially processed following the SMRT Link v6.0 standard pipeline [46]. The offline transcripts with full passes ≥2 were merged into CCSs. FLNC sequences of CCSs were then clustered and corrected by the interactive clustering and error correction (ICE) algorithm. The corrected isoforms were aligned with non-full-length transcripts using Quiver algorithm, generating polished consensus isoforms of high- and low-quality. The software CD-HIT-v4.6.7 was used to remove the redundancy of the polished high-quality isoforms by merging the sequences with a threshold of 99% identities, ultimately obtaining the FL transcriptome (unigenes). The assembly completeness was assessed by BUSCO software using Eukaryota dataset of OrthoDB [47].

3.6. Functional Annotation

All of the unigenes obtained were functionally annotated by BLASTx analysing against the NR, Swiss-Prot, KEGG and KOG databases with an E-value < 1 × 10−5. Each unigene was annotated with the information of the protein with the highest sequence similarity. GO annotation was performed using the Blast2GO software with the NR annotation results [48]. Unigenes with the first 20 highest scores and ≥33 high-scoring segment pair (HSPs) hits were selected to undergo Blast2GO analysis. Subsequently, the WEGO software was used to classify the functional annotation of the unigenes [49].

3.7. Gene Type Analyses

The unigenes were blastx searched against the databases with the E-values < 1 × 10−5 to retrieve a protein sequence for each unigene from either of the four databases in the order of NR, Swiss-Prot, KEGG and KOG, which then located the CDS of the unigene. The unigenes that failed to retrieve a protein sequence were subjected to ANGEL for CDS prediction [50]. LncRNA analysis was conducted on the unigenes that were not annotated to the four databases by combining the prediction of Coding-Non-Coding Index (CNCI) and Coding Potential Calculator (CPC) software [51,52]. Unigenes predicted as non-coding by both CNCI and CPC were considered lncRNAs. TF analysis was performed using hmmscan against the Plant TFdb database [53].

3.8. Detection of Alternative Splicing and Simple Sequence Repeats

The AS events were identified using the software Cogent and SUPPA [54,55]. Cogent partitioned high-quality FLNC sequences into gene families by K-mer = 30 and K-mer similarity >95%, and built each family a reference sequence by De Bruijn graph methods. Then, the AS events were detected using SUPPA with references. The software MISA was employed to identify SSRs within the FL transcriptome with the parameters of length-minimum number of repetitions = 2–6 or 3–5 or 4–4 or 5–4 or 6–4 and interruptions of 100 bp [56].

3.9. Mining of Chitin Metabolism Genes in Diatoms

Chitin-related genes in the T. weissflogii FL transcriptome were annotated by the four databases of NR, Swiss-Prot, KEGG and KOG. These genes in the genomes of P. tricornutum and T. pseudonana were retrieved from the Joint Genomics Institute PhycoCosm database (JGI PhycoCosm, https://jgi.doe.gov/data-and-tools/phycocosm/, accessed on 20 April 2021). Gene retrieval was performed by searching with enzyme names on the ENZYME database (https://enzyme.expasy.org/, accessed on 20 April 2021) and with genes published in previous literature as a supplement. Numbers of the corresponding genes in C. cryptica were obtained from Traller et al. (2016) [14].

4. Conclusions

In this study, we constructed a full-length transcriptome of T. weissflogii using PacBio sequencing. The transcriptome consists of 25,412 unigenes, 23,362 annotated unigenes, 710 lncRNAs, 363 TFs, 3113 AS events and 3295 SSRs. Furthermore, we identified 234 genes related to chitin metabolism. The whole metabolic pathway of chitin biosynthesis and degradation was explored. The information published here will pave the way for T. weissflogii molecular research in the future, expand the resource of β-chitin and promote the development of high-value enzymes.

Author Contributions

Conceptualisation, Z.S., C.B. and D.D.; methodology, H.C., X.X. and V.B.; writing—original draft preparation, H.C.; writing—review and editing, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (41806175), Youth Project of Marine Biology and Biotechnology Laboratory in Pilot National Laboratory for Marine Science and Technology (YQ2018NO06).

Data Availability Statement

The obtained raw data of the T. weissflogii FL transcriptome were deposited into the NCBI SRA database with the accession number PRJNA717330. Transcripts of T. weissflogii in MMETSP are openly available on the iMicrobe website with the accession number CAM_P_0001000 (https://www.imicrobe.us/, accessed on 15 February 2021).

Acknowledgments

We thank Guangzhou Gene Denovo Biotechnology Co., Ltd. for the assistance with sequencing and bioinformatics analyses. We would also like to thank Xiaohui Li from Ningbo University, China, for providing the T. weissflogii 9021 strain.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wagner, G.P. Evolution and multi-functionality of the chitin system. In Molecular Ecology and Evolution: Approaches and Applications; Schierwater, B., Streit, B., Wagner, G.P., DeSalle, R., Eds.; Birkhäuser: Basel, Switzerland, 1994; Volume 69, pp. 559–577. [Google Scholar]
  2. Muzzarelli, R.A.A.; Boudrant, J.; Meyer, D.; Manno, N.; DeMarchis, M.; Paoletti, M.G. Current views on fungal chitin/chitosan, human chitinases, food preservation, glucans, pectins and inulin: A tribute to Henri Braconnot, precursor of the carbohydrate polymers science, on the chitin bicentennial. Carbohydr. Polym. 2012, 87, 995–1012. [Google Scholar] [CrossRef]
  3. Tang, W.J.; Fernandez, J.; Sohn, J.J.; Amemiya, C.T. Chitin is endogenously produced in vertebrates. Curr. Biol. 2015, 25, 897–900. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Shigemasa, Y.; Minami, S. Applications of chitin and chitosan for biomaterials. Biotechnol. Genet. Eng. Rev. 1996, 13, 383–420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cuong, H.N.; Minh, N.C.; Van Hoa, N.; Trung, T.S. Preparation and characterization of high purity beta-chitin from squid pens (Loligo chenisis). Int. J. Biol. Macromol. 2016, 93, 442–447. [Google Scholar] [CrossRef]
  6. Shao, Z.; Thomas, Y.; Hembach, L.; Xing, X.; Duan, D.; Moerschbacher, B.M.; Bulone, V.; Tirichine, L.; Bowler, C. Comparative characterization of putative chitin deacetylases from Phaeodactylum tricornutum and Thalassiosira pseudonana highlights the potential for distinct chitin-based metabolic processes in diatoms. New Phytol. 2019, 221, 1890–1905. [Google Scholar] [CrossRef] [PubMed]
  7. Revol, J.F.; Chanzy, H. High-Resolution Electron Microscopy of β-Chitin Microfibrils. Biopolymers 1986, 25, 1599–1601. [Google Scholar] [CrossRef]
  8. Brunner, E.; Richthammer, P.; Ehrlich, H.; Paasch, S.; Simon, P.; Ueberlein, S.; van Pee, K.H. Chitin-based organic networks: An integral part of cell wall biosilica in the diatom Thalassiosira pseudonana. Angew. Chem. Int. Ed. Engl. 2009, 48, 9724–9727. [Google Scholar] [CrossRef] [PubMed]
  9. Ogawa, Y.; Kimura, S.; Wada, M. Electron diffraction and high-resolution imaging on highly-crystalline beta-chitin microfibril. J. Struct. Biol. 2011, 176, 83–90. [Google Scholar] [CrossRef]
  10. McLachlan, J.; Craigie, J.S. Chitan fibres in Cyclotella cryptica and growth of C. cryptica and Thalassiosira fluviatilis. In Some Contemporary Studies in Marine Science; Barnes, H., Ed.; George Allen and Unwin Ltd.: London, UK, 1966; pp. 511–517. [Google Scholar]
  11. Smucker, R.A. Chitin primary production. Biochem. Syst. Ecol. 1991, 19, 357–369. [Google Scholar] [CrossRef]
  12. McLachlan, J.; McInnes, A.G.; Falk, M. Studies on chitan (chitinpoly-n-acetylglucosamine) fibers of diatom Thalassiosira fluviatilis Hustedt. 1. Production and isolation of chitan fibers. Can. J. Bot. 1965, 43, 707. [Google Scholar] [CrossRef]
  13. Durkin, C.; Mock, T.; Armbrust, E. Chitin in Diatoms and Its Association with the Cell Wall. Eukaryot. Cell 2009, 8, 1038–1050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Traller, J.C.; Cokus, S.J.; Lopez, D.A.; Gaidarenko, O.; Smith, S.R.; McCrow, J.P.; Gallaher, S.D.; Podell, S.; Thompson, M.; Cook, O.; et al. Genome and methylome of the oleaginous diatom Cyclotella cryptica reveal genetic flexibility toward a high lipid phenotype. Biotechnol. Biofuels 2016, 9, 258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Field, C.; Behrenfeld, M.; Randerson, J.; Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 1998, 281, 237–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Nelson, D.M.; Tréguer, P.; Brzezinski, M.A.; Leynaert, A.; Quéguiner, B. Production and dissolution of biogenic silica in the ocean: Revised global estimates, comparison with regional data and relationship to biogenic sedimentation. Glob. Biogeochem. Cycles 1995, 9, 359–372. [Google Scholar] [CrossRef]
  17. Armbrust, E.; Berges, J.; Bowler, C.; Green, B.; Martinez, D.; Putnam, N.; Zhou, S.; Allen, A.; Apt, K.; Bechner, M.; et al. The Genome of the Diatom Thalassiosira Pseudonana: Ecology, Evolution, and Metabolism. Science 2004, 306, 79–86. [Google Scholar] [CrossRef] [Green Version]
  18. Bowler, C.; Allen, A.E.; Badger, J.H.; Grimwood, J.; Jabbari, K.; Kuo, A.; Maheswari, U.; Martens, C.; Maumus, F.; Otillar, R.P.; et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature 2008, 456, 239–244. [Google Scholar] [CrossRef]
  19. Paajanen, P.; Strauss, J.; van Oosterhout, C.; McMullan, M.; Clark, M.D.; Mock, T. Building a locally diploid genome and transcriptome of the diatom Fragilariopsis cylindrus. Sci. Data 2017, 4, 170149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Cheng, R.; Feng, J.; Zhang, B.; Huang, Y.; Cheng, J.; Zhang, C. Transcriptome and Gene Expression Analysis of an Oleaginous Diatom Under Different Salinity Conditions. BioEnergy Res. 2013, 7, 192–205. [Google Scholar] [CrossRef]
  21. Nanjappa, D.; Sanges, R.; Ferrante, M.I.; Zingone, A. Diatom flagellar genes and their expression during sexual reproduction in Leptocylindrus danicus. BMC Genom. 2017, 18, 813. [Google Scholar] [CrossRef] [Green Version]
  22. Galachyants, Y.P.; Zakharova, Y.R.; Volokitina, N.A.; Morozov, A.A.; Likhoshway, Y.V.; Grachev, M.A. De novo transcriptome assembly and analysis of the freshwater araphid diatom Fragilaria radians, Lake Baikal. Sci. Data 2019, 6, 183. [Google Scholar] [CrossRef] [Green Version]
  23. Keeling, P.J.; Burki, F.; Wilcox, H.M.; Allam, B.; Allen, E.E.; Amaral-Zettler, L.A.; Armbrust, E.V.; Archibald, J.M.; Bharti, A.K.; Bell, C.J.; et al. The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): Illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing. PLoS Biol. 2014, 12, e1001889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rhoads, A.; Au, K.F. PacBio Sequencing and Its Applications. Genom. Proteom. Bioinf. 2015, 13, 278–289. [Google Scholar] [CrossRef] [Green Version]
  25. Di Dato, V.; Di Costanzo, F.; Barbarinaldi, R.; Perna, A.; Ianora, A.; Romano, G. Unveiling the presence of biosynthetic pathways for bioactive compounds in the Thalassiosira rotula transcriptome. Sci. Rep. 2019, 9, 9893. [Google Scholar] [CrossRef] [PubMed]
  26. Di Dato, V.; Musacchia, F.; Petrosino, G.; Patil, S.; Montresor, M.; Sanges, R.; Ferrante, M.I. Transcriptome sequencing of three Pseudo-nitzschia species reveals comparable gene sets and the presence of Nitric Oxide Synthase genes in diatoms. Sci. Rep. 2015, 5, 12329. [Google Scholar] [CrossRef] [PubMed]
  27. Buhmann, M.T.; Poulsen, N.; Klemm, J.; Kennedy, M.R.; Sherrill, C.D.; Kroger, N. A tyrosine-rich cell surface protein in the diatom Amphora coffeaeformis identified through transcriptome analysis and genetic transformation. PLoS ONE 2014, 9, e110369. [Google Scholar] [CrossRef]
  28. Tanaka, T.; Maeda, Y.; Veluchamy, A.; Tanaka, M.; Abida, H.; Marechal, E.; Bowler, C.; Muto, M.; Sunaga, Y.; Tanaka, M.; et al. Oil accumulation by the oleaginous diatom Fistulifera solaris as revealed by the genome and transcriptome. Plant Cell 2015, 27, 162–176. [Google Scholar] [CrossRef] [Green Version]
  29. Hong, F.; Mo, S.H.; Lin, X.Y.; Niu, J.; Yin, J.; Wei, D. The PacBio Full-Length Transcriptome of the Tea Aphid as a Reference Resource. Front. Genet. 2020, 11, 558394. [Google Scholar] [CrossRef]
  30. Chen, C.; Shi, X.; Zhou, T.; Li, W.; Li, S.; Bai, G. Full-length transcriptome analysis and identification of genes involved in asarinin and aristolochic acid biosynthesis in medicinal plant Asarum sieboldii. Genome 2020. [Google Scholar] [CrossRef]
  31. Ma, L.; Bajic, V.B.; Zhang, Z. On the classification of long non-coding RNAs. RNA Biol. 2013, 10, 925–933. [Google Scholar] [CrossRef]
  32. Babu, M.M. Structure, evolution and dynamics of transcriptional regulatory networks. Biochem. Soc. Trans. 2010, 38, 1155–1178. [Google Scholar] [CrossRef]
  33. Rayko, E.; Maumus, F.; Maheswari, U.; Jabbari, K.; Bowler, C. Transcription factor families inferred from genome sequences of photosynthetic stramenopiles. New Phytol. 2010, 188, 52–66. [Google Scholar] [CrossRef]
  34. Cruz de Carvalho, M.H.; Sun, H.X.; Bowler, C.; Chua, N.H. Noncoding and coding transcriptome responses of a marine diatom to phosphate fluctuations. New Phytol. 2016, 210, 497–510. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sarthou, G.; Timmermans, K.R.; Blain, S.; Tréguer, P. Growth physiology and fate of diatoms in the ocean: A review. J. Sea Res. 2005, 53, 25–42. [Google Scholar] [CrossRef]
  36. Allen, A.E.; Vardi, A.; Bowler, C. An ecological and evolutionary context for integrated nitrogen metabolism and related signaling pathways in marine diatoms. Curr. Opin. Plant Biol. 2006, 9, 264–273. [Google Scholar] [CrossRef]
  37. Black, D.L. Mechanisms of Alternative Pre-Messenger RNA Splicing. Annu. Rev. Biochem. 2003, 72, 291–336. [Google Scholar] [CrossRef] [Green Version]
  38. Rastogi, A.; Maheswari, U.; Dorrell, R.G.; Vieira, F.R.J.; Maumus, F.; Kustka, A.; McCarthy, J.; Allen, A.E.; Kersey, P.; Bowler, C.; et al. Integrative analysis of large scale transcriptome data draws a comprehensive landscape of Phaeodactylum tricornutum genome and evolutionary origin of diatoms. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Sharopova, N.; McMullen, M.D.; Schultz, L.; Schroeder, S.; Sanchez-Villeda, H.; Gardiner, J.; Bergstrom, D.; Houchins, K.; MeliaHancock, S.; Musket, T.; et al. Development and mapping of SSR markers for maize. Plant Mol. Biol. 2002, 48, 463–481. [Google Scholar] [CrossRef]
  40. Karre, S.; Kumar, A.; Dhokane, D.; Kushalappa, A.C. Metabolo-transcriptome profiling of barley reveals induction of chitin elicitor receptor kinase gene (HvCERK1) conferring resistance against Fusarium graminearum. Plant Mol. Biol. 2017, 93, 247–267. [Google Scholar] [CrossRef] [PubMed]
  41. Giubergia, S.; Phippen, C.; Nielsen, K.F.; Gram, L. Growth on Chitin Impacts the Transcriptome and Metabolite Profiles of Antibiotic-Producing Vibrio coralliilyticus S2052 and Photobacterium galatheae S2753. mSystems 2017, 2. [Google Scholar] [CrossRef] [Green Version]
  42. Shao, Z.M.; Li, Y.J.; Zhang, X.R.; Chu, J.; Ma, J.H.; Liu, Z.X.; Wang, J.; Sheng, S.; Wu, F.A. Identification and Functional Study of Chitin Metabolism and Detoxification-Related Genes in Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) Based on Transcriptome Analysis. Int. J. Mol. Sci. 2020, 21, 1904. [Google Scholar] [CrossRef] [Green Version]
  43. Cheng, H.; Shao, Z.; Lu, C.; Duan, D. Genome-wide identification of chitinase genes in Thalassiosira pseudonana and analysis of their expression under abiotic stresses. BMC Plant Biol. 2021, 21, 87. [Google Scholar] [CrossRef] [PubMed]
  44. Hardt, M.; Laine, R.A. Mutation of active site residues in the chitin-binding domain ChBDChiA1 from chitinase A1 of Bacillus circulans alters substrate specificity: Use of a green fluorescent protein binding assay. Arch Biochem. Biophys. 2004, 426, 286–297. [Google Scholar] [CrossRef] [PubMed]
  45. Nampally, M.; Moerschbacher, B.M.; Kolkenbrock, S. Fusion of a novel genetically engineered chitosan affinity protein and green fluorescent protein for specific detection of chitosan in vitro and in situ. Appl. Environ. Microbiol. 2012, 78, 3114–3119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gordon, S.P.; Tseng, E.; Salamov, A.; Zhang, J.; Meng, X.; Zhao, Z.; Kang, D.; Underwood, J.; Grigoriev, I.V.; Figueroa, M.; et al. Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing. PLoS ONE 2015, 10, e0132628. [Google Scholar] [CrossRef] [Green Version]
  47. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  50. Shimizu, K.; Adachi, J.; Muraoka, Y. ANGLE: A sequencing errors resistant program for predicting protein coding regions in unfinished cDNA. J. Bioinf. Comput. Biol. 2006, 4, 649–664. [Google Scholar] [CrossRef]
  51. Sun, L.; Luo, H.; Bu, D.; Zhao, G.; Yu, K.; Zhang, C.; Liu, Y.; Chen, R.; Zhao, Y. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. [Google Scholar] [CrossRef]
  52. Kong, L.; Zhang, Y.; Ye, Z.Q.; Liu, X.Q.; Zhao, S.Q.; Wei, L.; Gao, G. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345–W349. [Google Scholar] [CrossRef]
  53. Tian, F.; Yang, D.C.; Meng, Y.Q.; Jin, J.; Gao, G. PlantRegMap: Charting functional regulatory maps in plants. Nucleic Acids Res. 2020, 48, D1104–D1113. [Google Scholar] [CrossRef] [PubMed]
  54. Li, J.; Harata-Lee, Y.; Denton, M.D.; Feng, Q.; Rathjen, J.R.; Qu, Z.; Adelson, D.L. Long read reference genome-free reconstruction of a full-length transcriptome from Astragalus membranaceus reveals transcript variants involved in bioactive compound biosynthesis. Cell Discov. 2017, 3, 17031. [Google Scholar] [CrossRef] [PubMed]
  55. Alamancos, G.P.; Pages, A.; Trincado, J.L.; Bellora, N.; Eyras, E. Leveraging transcript quantification for fast computation of alternative splicing profiles. RNA 2015, 21, 1521–1531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Beier, S.; Thiel, T.; Munch, T.; Scholz, U.; Mascher, M. MISA-web: A web server for microsatellite prediction. Bioinformatics 2017, 33, 2583–2585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Marinedrugs 19 00392 g001 550
Figure 1. Structure of α-chitin and β-chitin. Anti-parallel arrangement of chitin chains in α-chitin whereas parallel arrangement for β-chitin.
Figure 1. Structure of α-chitin and β-chitin. Anti-parallel arrangement of chitin chains in α-chitin whereas parallel arrangement for β-chitin.
Marinedrugs 19 00392 g001
Marinedrugs 19 00392 g002 550
Figure 2. Fluorescence observations indicating the presence of chitin and chitosan in T. weissflogii; (A) Cells stained with CBP-eGFP; (B) Chitin microfibrils observation after CBP-eGFP staining; (C) Cells stained with CAP-eGFP. TL: transmission light; GFP: green fluorescent protein. Scale bar = 5 μm.
Figure 2. Fluorescence observations indicating the presence of chitin and chitosan in T. weissflogii; (A) Cells stained with CBP-eGFP; (B) Chitin microfibrils observation after CBP-eGFP staining; (C) Cells stained with CAP-eGFP. TL: transmission light; GFP: green fluorescent protein. Scale bar = 5 μm.
Marinedrugs 19 00392 g002
Marinedrugs 19 00392 g003 550
Figure 3. Results of BUSCO analysis. “C”: Complete; “S”: Single-copy; “D”: Duplicated; “F”: Fragmented; “M”: Missing.
Figure 3. Results of BUSCO analysis. “C”: Complete; “S”: Single-copy; “D”: Duplicated; “F”: Fragmented; “M”: Missing.
Marinedrugs 19 00392 g003
Marinedrugs 19 00392 g004 550
Figure 4. Composition of the T. weissflogii FL transcriptome; (A) Length distribution of unigenes; (B) Length distribution of CDSs; (C) LncRNA numbers predicted by CNCI and CPC; (D) Numbers of top ten TF families.
Figure 4. Composition of the T. weissflogii FL transcriptome; (A) Length distribution of unigenes; (B) Length distribution of CDSs; (C) LncRNA numbers predicted by CNCI and CPC; (D) Numbers of top ten TF families.
Marinedrugs 19 00392 g004
Marinedrugs 19 00392 g005 550
Figure 5. Annotation information of T. weissflogii FL transcriptome unigenes; (A) Annotation summary from the NR, Swiss-Prot, KEGG and KOG databases; (B) Species distribution annotation from the NR database; (C) Functional annotation from the KOG database; (D) Functional annotation from the KEGG database.
Figure 5. Annotation information of T. weissflogii FL transcriptome unigenes; (A) Annotation summary from the NR, Swiss-Prot, KEGG and KOG databases; (B) Species distribution annotation from the NR database; (C) Functional annotation from the KOG database; (D) Functional annotation from the KEGG database.
Marinedrugs 19 00392 g005
Marinedrugs 19 00392 g006 550
Figure 6. AS events and SSRs of T. weissflogii FL transcriptome; (A) Alternatively spliced isoform distribution; (B) AS types identified in the FL transcriptome, including A3 (alternative 3′ splice sites), A5 (alternative 5′ splice sites), AL (alternative last exons) and RI (retained intron); (C) Percentages of SSR motifs; (D) 3D histogram of SSR components. Di- stands for dinucleotide repeats, Tri- for trinucleotide repeats, Tetra- for tetranucleotide repeats, Penta- for pentanucleotide repeats, Hexa- for hexanucleotide repeats.
Figure 6. AS events and SSRs of T. weissflogii FL transcriptome; (A) Alternatively spliced isoform distribution; (B) AS types identified in the FL transcriptome, including A3 (alternative 3′ splice sites), A5 (alternative 5′ splice sites), AL (alternative last exons) and RI (retained intron); (C) Percentages of SSR motifs; (D) 3D histogram of SSR components. Di- stands for dinucleotide repeats, Tri- for trinucleotide repeats, Tetra- for tetranucleotide repeats, Penta- for pentanucleotide repeats, Hexa- for hexanucleotide repeats.
Marinedrugs 19 00392 g006
Marinedrugs 19 00392 g007 550
Figure 7. Hypothetical chitin metabolism pathways in four diatom species. Enzymes are represented by EC numbers: Glutamine--fructose-6-phosphate transaminase (isomerising) (EC:2.6.1.16), N-acetylglucosamine-6-phosphate deacetylase (EC:3.5.1.25), Phosphoacetylglucosamine mutase (EC:5.4.2.3), UDP-N-acetylglucosamine diphosphorylase (EC:2.7.7.23), Chitin synthase (EC:2.4.1.16), Chitinase (EC:3.2.1.14), Chitin deacetylase (EC:3.5.1.41), Beta-N-acetylhexosaminidase (EC:3.2.1.52), N-acetylglucosamine kinase (EC:2.7.1.59). Numbers of identified enzyme-encoding genes are represented in differently coloured boxes. Enzyme-encoding genes unidentified in the T. weissflogi FL transcriptome are indicated by dashed lines. The schematic was adapted from Traller et al. (2016) [14].
Figure 7. Hypothetical chitin metabolism pathways in four diatom species. Enzymes are represented by EC numbers: Glutamine--fructose-6-phosphate transaminase (isomerising) (EC:2.6.1.16), N-acetylglucosamine-6-phosphate deacetylase (EC:3.5.1.25), Phosphoacetylglucosamine mutase (EC:5.4.2.3), UDP-N-acetylglucosamine diphosphorylase (EC:2.7.7.23), Chitin synthase (EC:2.4.1.16), Chitinase (EC:3.2.1.14), Chitin deacetylase (EC:3.5.1.41), Beta-N-acetylhexosaminidase (EC:3.2.1.52), N-acetylglucosamine kinase (EC:2.7.1.59). Numbers of identified enzyme-encoding genes are represented in differently coloured boxes. Enzyme-encoding genes unidentified in the T. weissflogi FL transcriptome are indicated by dashed lines. The schematic was adapted from Traller et al. (2016) [14].
Marinedrugs 19 00392 g007
Table
Table 1. Glycosidic linkage composition of the three Thalassiosira microalgae samples.
Table 1. Glycosidic linkage composition of the three Thalassiosira microalgae samples.
No.LinkageComposition (Mol%)
TwTpTr
14-Glcp15.525.44.3
23-Glcp12.827.856.4
34-GlcNAcp6.90.61.4
4t-Galp6.72.23.1
52,3-Glcp5.51.52.0
6t-Manp5.41.71.7
72-Manp4.92.43.1
84-Xylp3.71.91.4
9t-Xylp3.32.90.8
106-Manp2.9NDND
Tw: Thalassiosira weissflogii; Tp: Thalassiosira pseudonana; Tr: Thalassiosira rotula; ND means “not detected”.
Table
Table 2. Summary of the T. weissflogii FL transcriptome statistics.
Table 2. Summary of the T. weissflogii FL transcriptome statistics.
Statistical DataT. Weissflogii
Raw readsSubread number44,233,932
Average length (bp)1583
N50 (bp)2003
CCSsNumber of reads
Number of CCS bases
CCS read length (mean) (bp)
Number of passes (mean)		1,021,310
2,172,901,990
2127
8
Clustered readsNumber of polished high-quality isoforms110,527
Number of polished low-quality isoforms338
UnigenesTotal number25,412
Total length (bp)51,968,546
Maximum length (bp)11,939
Minimum length (bp)64
Average length (bp)2045.04
N50 length (bp)2417
GC content46.95%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, H.; Bowler, C.; Xing, X.; Bulone, V.; Shao, Z.; Duan, D. Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes. Mar. Drugs 2021, 19, 392. https://doi.org/10.3390/md19070392

AMA Style

Cheng H, Bowler C, Xing X, Bulone V, Shao Z, Duan D. Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes. Marine Drugs. 2021; 19(7):392. https://doi.org/10.3390/md19070392

Chicago/Turabian Style

Cheng, Haomiao, Chris Bowler, Xiaohui Xing, Vincent Bulone, Zhanru Shao, and Delin Duan. 2021. "Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes" Marine Drugs 19, no. 7: 392. https://doi.org/10.3390/md19070392

APA Style

Cheng, H., Bowler, C., Xing, X., Bulone, V., Shao, Z., & Duan, D. (2021). Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes. Marine Drugs, 19(7), 392. https://doi.org/10.3390/md19070392

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop