Distinct genome‐wide alternative polyadenylation during the response to silicon availability in the marine diatom Thalassiosira pseudonana

The draft genome of Nitzschia closterium f. minutissima and transcriptome analysis reveals novel insights into diatom biosilicification

BACKGROUND

Nitzschia closterium f. minutissima is a commonly available diatom that plays important roles in marine aquaculture. It was originally classified as Nitzschia (Bacillariaceae, Bacillariophyta) but is currently regarded as a heterotypic synonym of Phaeodactylum tricornutum. The aim of this study was to obtain the draft genome of the marine microalga N. closterium f. minutissima to understand its phylogenetic placement and evolutionary specialization. Given that the ornate hierarchical silicified cell walls (frustules) of diatoms have immense applications in nanotechnology for biomedical fields, biosensors and optoelectric devices, transcriptomic data were generated by using reference genome-based read mapping to identify significantly differentially expressed genes and elucidate the molecular processes involved in diatom biosilicification.

RESULTS

In this study, we generated 13.81 Gb of pass reads from the PromethION sequencer. The draft genome of N. closterium f. minutissima has a total length of 29.28 Mb, and contains 28 contigs with an N50 value of 1.23 Mb. The GC content was 48.55%, and approximately 18.36% of the genome assembly contained repeat sequences. Gene annotation revealed 9,132 protein-coding genes. The results of comparative genomic analysis showed that N. closterium f. minutissima was clustered as a sister lineage of Phaeodactylum tricornutum and the divergence time between them was estimated to be approximately 17.2 million years ago (Mya). CAFF analysis demonstrated that 220 gene families that significantly changed were unique to N. closterium f. minutissima and that 154 were specific to P. tricornutum, moreover, only 26 gene families overlapped between these two species. A total of 818 DEGs in response to silicon were identified in N. closterium f. minutissima through RNA sequencing, these genes are involved in various molecular processes such as transcription regulator activity. Several genes encoding proteins, including silicon transporters, heat shock factors, methyltransferases, ankyrin repeat domains, cGMP-mediated signaling pathways-related proteins, cytoskeleton-associated proteins, polyamines, glycoproteins and saturated fatty acids may contribute to the formation of frustules in N. closterium f. minutissima.

CONCLUSIONS

Here, we described a draft genome of N. closterium f. minutissima and compared it with those of eight other diatoms, which provided new insight into its evolutionary features. Transcriptome analysis to identify DEGs in response to silicon will help to elucidate the underlying molecular mechanism of diatom biosilicification in N. closterium f. minutissima.

Alternative polyadenylation profiles of susceptible and resistant rice (Oryza sativa L.) in response to bacterial leaf blight using RNA-seq

BACKGROUND

Alternative polyadenylation (APA) is an important pattern of post-transcriptional regulation of genes widely existing in eukaryotes, involving plant physiological and pathological processes. However, there is a dearth of studies investigating the role of APA profile in rice leaf blight.

RESULTS

In this study, we compared the APA profile of leaf blight-susceptible varieties (CT 9737-613P-M) and resistant varieties (NSIC RC154) following bacterial blight infection. Through gene enrichment analysis, we found that the genes of two varieties typically exhibited distal poly(A) (PA) sites that play different roles in two kinds of rice, indicating differential APA regulatory mechanisms. In this process, many disease-resistance genes displayed multiple transcripts via APA. Moreover, we also found five polyadenylation factors of similar expression patterns of rice, highlighting the critical roles of these five factors in rice response to leaf blight about PA locus diversity.

CONCLUSION

Notably, the present study provides the first dynamic changes of APA in rice in early response to biotic stresses and proposes a possible functional conjecture of APA in plant immune response, which lays the theoretical foundation for in-depth determination of the role of APA events in plant stress response and other life processes.

Differential Expression of Stress Adaptation Genes in a Diatom Ulnaria acus under Different Culture Conditions

Diatoms are a group of unicellular eukaryotes that are essential primary producers in aquatic ecosystems. The dynamic nature of their habitat necessitates a quick and specific response to various stresses. However, the molecular mechanisms of their physiological adaptations are still underexplored. In this work, we study the response of the cosmopolitan freshwater diatom Ulnaria acus (Bacillariophyceae, Fragilariophycidae, Licmophorales, Ulnariaceae, Ulnaria) in relation to a range of stress factors, namely silica deficiency, prolonged cultivation, and interaction with an algicidal bacterium. Fluorescent staining and light microscopy were used to determine the physiological state of cells under these stresses. To explore molecular reactions, we studied the genes involved in the stress response-type III metacaspase (MC), metacaspase-like proteases (MCP), death-specific protein (DSP), delta-1-pyrroline-5-carboxylate dehydrogenase (ALDH12), and glutathione synthetase (GSHS). We have described the structure of these genes, analyzed the predicted amino acid sequences, and measured their expression dynamics in vitro using qRT-PCR. We demonstrated that the expression of UaMC1, UaMC3, and UaDSP increased during the first five days of silicon starvation. On the seventh day, it was replaced with the expression of UaMC2, UaGSHS, and UaALDH. After 45 days of culture, cells stopped growing, and the expression of UaMC1, UaMC2, UaGSHS, and UaDSP increased. Exposure to an algicidal bacterial filtrate induced a higher expression of UaMC1 and UaGSHS. Thus, we can conclude that these proteins are involved in diatoms' adaptions to environmental changes. Further, these data show that the molecular adaptation mechanisms in diatoms depend on the nature and exposure duration of a stress factor.

The Art of Exploring Diatom Biosilica Biomaterials: From Biofabrication Perspective

Diatom is a common single-cell microalgae with large species and huge biomass. Diatom biosilica (DB), the shell of diatom, is a natural inorganic material with a micro-nanoporous structure. Its unique hierarchical porous structure gives it great application potential in drug delivery, hemostat materials, and biosensors, etc. However, the structural diversity of DB determines its different biological functions. Screening hundreds of thousands of diatom species for structural features of DB that meet application requirements is like looking for a needle in a seaway. And the chemical modification methods lack effective means to control the micro-nanoporous structure of DB. The formation of DB is a typical biomineralization process, and its structural characteristics are affected by external environmental conditions, genes, and other factors. This allows to manipulate the micro-nanostructure of DB through biological regulation method, thereby transforming the screening mode of the structure function of DB from a needle in a seaway to biofabrication mode. This review focuses on the formation, biological modification, functional activity of DB structure, and its application in biomaterials field, providing regulatory strategies and research idea of DB from the perspective of biofabrication. It will also maximize the possibility of using DB as biological materials.

Thalassiosira pseudonana (Cyclotella nana) (Hustedt) Hasle et Heimdal (Bacillariophyceae): A genetically tractable model organism for studying diatom biology, including biological silica formation

In 2004, Thalassiosira pseudonana was the first eukaryotic marine alga to have its genome sequenced. Since then, this species has quickly emerged as a valuable model species for investigating the molecular underpinnings of essentially all aspects of diatom life, particularly bio-morphogenesis of the cell wall. An important prerequisite for the model status of T. pseudonana is the ongoing development of increasingly precise tools to study the function of gene networks and their encoded proteins in vivo. Here, we briefly review the current toolbox for genetic manipulation, highlight specific examples of its application in studying diatom metabolism, and provide a peek into the role of diatoms in the emerging field of silica biotechnology.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

Divergence in the Regulation of the Salt Tolerant Response Between Arabidopsis thaliana and Its Halophytic Relative Eutrema salsugineum by mRNA Alternative Polyadenylation

Salt tolerance is an important mechanism by which plants can adapt to a saline environment. To understand the process of salt tolerance, we performed global analyses of mRNA alternative polyadenylation (APA), an important regulatory mechanism during eukaryotic gene expression, in Arabidopsis thaliana and its halophytic relative Eutrema salsugineum with regard to their responses to salt stress. Analyses showed that while APA occurs commonly in both Arabidopsis and Eutrema, Eutrema possesses fewer APA genes than Arabidopsis (47% vs. 54%). However, the proportion of APA genes was significantly increased in Arabidopsis under salt stress but not in Eutrema. This indicated that Arabidopsis is more sensitive to salt stress and that Eutrema exhibits an innate response to such conditions. Both species utilized distal poly(A) sites under salt stress; however, only eight genes were found to overlap when their 3' untranslated region (UTR) lengthen genes were compared, thus revealing their distinct responses to salt stress. In Arabidopsis, genes that use distal poly(A) sites were enriched in response to salt stress. However, in Eutrema, the use of poly(A) sites was less affected and fewer genes were enriched. The transcripts with upregulated poly(A) sites in Arabidopsis showed enriched pathways in plant hormone signal transduction, starch and sucrose metabolism, and fatty acid elongation; in Eutrema, biosynthetic pathways (stilbenoid, diarylheptanoid, and gingerol) and metabolic pathways (arginine and proline) showed enrichment. APA was associated with 42% and 29% of the differentially expressed genes (DE genes) in Arabidopsis and Eutrema experiencing salt stress, respectively. Salt specific poly(A) sites and salt-inducible APA events were identified in both species; notably, some salt tolerance-related genes and transcription factor genes exhibited differential APA patterns, such as CIPK21 and LEA4-5. Our results suggest that adapted species exhibit more orderly response at the RNA maturation step under salt stress, while more salt-specific poly(A) sites were activated in Arabidopsis to cope with salinity conditions. Collectively, our findings not only highlight the importance of APA in the regulation of gene expression in response to salt stress, but also provide a new perspective on how salt-sensitive and salt-tolerant species perform differently under stress conditions through transcriptome diversity.

QPAT-seq, a rapid and deduplicatable method for quantification of poly(A) site usages

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples.

HDA6-dependent histone deacetylation regulates mRNA polyadenylation in Arabidopsis

Eukaryotic histone deacetylation, critical for maintaining nucleosome structure and regulating gene expression, is mediated by histone deacetylases (HDACs). Although nucleosomes have been reported to regulate mRNA polyadenylation in humans, the role of HDACs in regulating polyadenylation has not been uncovered. Taking advantage of phenotypic studies on Arabidopsis, HDA6 (one of HDACs) was found to be a critical part of many biological processes. Here, we report that HDA6 affects mRNA polyadenylation in Arabidopsis. Poly(A) sites of up-regulated transcripts are closer to the histone acetylation peaks in hda6 compared to the wild-type Col-0. HDA6 is required for the deacetylation of histones around DNA on nucleosomes, which solely coincides with up-regulated or uniquely presented poly(A) sites in hda6. Furthermore, defective HDA6 results in an overrepresentation of the canonical poly(A) signal (AAUAAA) usage. Chromatin loci for generating AAUAAA-type transcripts have a comparatively low H3K9K14ac around poly(A) sites when compared to other noncanonical poly(A) signal–containing transcripts. These results indicate that HDA6 regulates polyadenylation in a histone deacetylation–dependent manner in Arabidopsis.