Experimental Genome-Wide Determination of RNA Polyadenylation in Chlamydomonas reinhardtii

Nucleotide-resolution Mapping of RNA N6-Methyladenosine (m6A) modifications and comprehensive analysis of global polyadenylation events in mRNA 3' end processing in malaria pathogen Plasmodium falciparum

Plasmodium falciparum is an obligate human parasite of the phylum Apicomplexa and is the causative agent of the most lethal form of human malaria. Although N6-methyladenosine modification is thought to be one of the major post-transcriptional regulatory mechanisms for stage-specific gene expression in apicomplexan parasites, the precise base position of m6A in mRNAs or noncoding RNAs in these parasites remains unknown. Here, we report global nucleotide-resolution mapping of m6A residues in P. falciparum using DART-seq technology, which quantitatively displayed a stage-specific, dynamic distribution pattern with enrichment near mRNA 3' ends. In this process we identified 894, 788, and 1,762 m6A-modified genes in Ring, Trophozoite and Schizont stages respectively, with an average of 5-7 m6A sites per-transcript at the individual gene level. Notably, several genes involved in malaria pathophysiology, such as KAHRP, ETRAMPs, SERA and stress response genes, such as members of Heat Shock Protein (HSP) family are highly enriched in m6A and therefore could be regulated by this RNA modification. Since we observed preferential methylation at the 3' ends of P. falciparum transcripts and because malaria polyadenylation specificity factor PfCPSF30 harbors an m6A reader 'YTH' domain, we reasoned that m6A might play an important role in 3'-end processing of malaria mRNAs. To investigate this, we used two complementary high-throughput RNA 3'-end mapping approaches, which provided an initial framework to explore potential roles of m6A in the regulation of alternative polyadenylation (APA) during malaria development in human hosts.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Phylogenetic analysis and structural prediction reveal the potential functional diversity between green algae SWEET transporters

Sugar-Will-Eventually-be-Exported-Transporters (SWEETs) are an important family of sugar transporters that appear to be ubiquitous in all organisms. Recent research has determined the structure of SWEETs in higher plants, identified specific residues required for monosaccharide or disaccharide transport, and begun to understand the specific functions of individual plant SWEET proteins. However, in green algae (Chlorophyta) these transporters are poorly characterised. This study identified SWEET proteins from across representative Chlorophyta with the aim to characterise their phylogenetic relationships and perform protein structure modelling in order to inform functional prediction. The algal genomes analysed encoded between one and six SWEET proteins, which is much less than a typical higher plant. Phylogenetic analysis identified distinct clusters of over 70 SWEET protein sequences, taken from almost 30 algal genomes. These clusters remain separate from representative higher or non-vascular plant SWEETs, but are close to fungi SWEETs. Subcellular localisation predictions and analysis of conserved amino acid residues revealed variation between SWEET proteins of different clusters, suggesting different functionality. These findings also showed conservation of key residues at the substrate-binding site, indicating a similar mechanism of substrate selectivity and transport to previously characterised higher plant monosaccharide-transporting SWEET proteins. Future work is now required to confirm the predicted sugar transport specificity and determine the functional role of these algal SWEET proteins.

Identification of Alternative Polyadenylation in Cyanidioschyzon merolae Through Long-Read Sequencing of mRNA

Alternative polyadenylation (APA) is widespread among metazoans and has been shown to have important impacts on mRNA stability and protein expression. Beyond a handful of well-studied organisms, however, its existence and consequences have not been well investigated. We therefore turned to the deep-branching red alga, Cyanidioschyzon merolae, to study the biology of polyadenylation in an organism highly diverged from humans and yeast. C. merolae is an acidothermophilic alga that lives in volcanic hot springs. It has a highly reduced genome (16.5 Mbp) and has lost all but 27 of its introns and much of its splicing machinery, suggesting that it has been under substantial pressure to simplify its RNA processing pathways. We used long-read sequencing to assess the key features of C. merolae mRNAs, including splicing status and polyadenylation cleavage site (PAS) usage. Splicing appears to be less efficient in C. merolae compared with yeast, flies, and mammalian cells. A high proportion of transcripts (63%) have at least two distinct PAS's, and 34% appear to utilize three or more sites. The apparent polyadenylation signal UAAA is used in more than 90% of cases, in cells grown in both rich media or limiting nitrogen. Our documentation of APA for the first time in this non-model organism highlights its conservation and likely biological importance of this regulatory step in gene expression.

Transcriptional dynamics in the protozoan parasite Sarcocystis neurona and mammalian host cells after treatment with a specific inhibitor of apicomplexan mRNA polyadenylation

In recent years, a class of chemical compounds (benzoxaboroles) that are active against a range of parasites has been shown to target mRNA polyadenylation by inhibiting the activity of CPSF73, the endonucleolytic core of the eukaryotic polyadenylation complex. One particular compound, termed AN3661, is active against several apicomplexan parasites that cause disease in humans. In this study, we report that AN3661 is active against an apicomplexan that causes disease in horses and marine mammals (Sarcocystis neurona), with an approximate IC50 value of 14.99 nM. Consistent with the reported mode of action of AN3661 against other apicomplexans, S. neurona mutants resistant to AN3661 had an alteration in CPSF73 that was identical to a mutation previously documented in AN3661-resistant Toxoplasma gondii and Plasmodium falciparum. AN3661 had a wide-ranging effect on poly(A) site choice in S. neurona, with more than half of all expressed genes showing some alteration in mRNA 3' ends. This was accompanied by changes in the relative expression of more than 25% of S. neurona genes and an overall 5-fold reduction of S. neurona transcripts in infected cells. In contrast, AN3661 had no discernible effect on poly(A) site choice or gene expression in the host cells. These transcriptomic studies indicate that AN3661 is exceedingly specific for the parasite CPSF73 protein, and has the potential to augment other therapies for the control of apicomplexan parasites in domestic animals.

Exploring the Impact of Terminators on Transgene Expression in Chlamydomonas reinhardtii with a Synthetic Biology Approach

Chlamydomonas reinhardtii has many attractive features for use as a model organism for both fundamental studies and as a biotechnological platform. Nonetheless, despite the many molecular tools and resources that have been developed, there are challenges for its successful engineering, in particular to obtain reproducible and high levels of transgene expression. Here we describe a synthetic biology approach to screen several hundred independent transformants using standardised parts to explore different parameters that might affect transgene expression. We focused on terminators and, using a standardised workflow and quantitative outputs, tested 9 different elements representing three different size classes of native terminators to determine their ability to support high level expression of a GFP reporter gene. We found that the optimal size reflected the median size of element found in the C. reinhardtii genome. The behaviour of the terminator parts was similar with different promoters, in different host strains and with different transgenes. This approach is applicable to the systematic testing of other genetic elements, facilitating comparison to determine optimal transgene design.

A marine Chlamydomonas sp. emerging as an algal model

The freshwater microalga Chlamydomonas reinhardtii, which lives in wet soil, has served for decades as a model for numerous biological processes, and many tools have been introduced for this organism. Here, we have established a stable nuclear transformation for its marine counterpart, Chlamydomonas sp. SAG25.89, by fusing specific cis-acting elements from its Actin gene with the gene providing hygromycin resistance and using an elaborated electroporation protocol. Like C. reinhardtii, Chlamydomonas sp. has a high GC content, allowing reporter genes and selection markers to be applicable in both organisms. Chlamydomonas sp. grows purely photoautotrophically and requires ammonia as a nitrogen source because its nuclear genome lacks some of the genes required for nitrogen metabolism. Interestingly, it can grow well under both low and very high salinities (up to 50 g · L^-1 ) rendering it as a model for osmotolerance. We further show that Chlamydomonas sp. grows well from 15 to 28°C, but halts its growth at 32°C. The genome of Chlamydomonas sp. contains some gene homologs the expression of which is regulated according to the ambient temperatures and/or confer thermal acclimation in C. reinhardtii. Thus, knowledge of temperature acclimation can now be compared to the marine species. Furthermore, Chlamydomonas sp. can serve as a model for studying marine microbial interactions and for comparing mechanisms in freshwater and marine environments. Chlamydomonas sp. was previously shown to be immobilized rapidly by a cyclic lipopeptide secreted from the antagonistic bacterium Pseudomonas protegens PF-5, which deflagellates C. reinhardtii.

Profiling Alternative 3' Untranslated Regions in Sorghum using RNA-seq Data

Sorghum is an important crop widely used for food, feed, and fuel. Transcriptome-wide studies of 3′ untranslated regions (3′UTR) using regular RNA-seq remain scarce in sorghum, while transcriptomes have been characterized extensively using Illumina short-read sequencing platforms for many sorghum varieties under various conditions or developmental contexts. 3′UTR is a critical regulatory component of genes, controlling the translation, transport, and stability of messenger RNAs. In the present study, we profiled the alternative 3′UTRs at the transcriptome level in three genetically related but phenotypically contrasting lines of sorghum: Rio, BTx406, and R9188. A total of 1,197 transcripts with alternative 3′UTRs were detected using RNA-seq data. Their categorization identified 612 high-confidence alternative 3′UTRs. Importantly, the high-confidence alternative 3′UTR genes significantly overlapped with the genesets that are associated with RNA N⁶-methyladenosine (m⁶A) modification, suggesting a clear indication between alternative 3′UTR and m⁶A methylation in sorghum. Moreover, taking advantage of sorghum genetics, we provided evidence of genotype specificity of alternative 3′UTR usage. In summary, our work exemplifies a transcriptome-wide profiling of alternative 3′UTRs using regular RNA-seq data in non-model crops and gains insights into alternative 3′UTRs and their genotype specificity.

Wide-ranging transcriptome remodelling mediated by alternative polyadenylation in response to abiotic stresses in Sorghum

Alternative polyadenylation (APA) regulates diverse developmental and physiological processes through its effects on gene expression, mRNA stability, translatability, and transport. Sorghum is a major cereal crop in the world and, despite its importance, not much is known about the role of post-transcriptional regulation in mediating responses to abiotic stresses in Sorghum. A genome-wide APA analysis unveiled widespread occurrence of APA in Sorghum in response to drought, heat, and salt stress. Abiotic stress treatments incited changes in poly(A) site choice in a large number of genes. Interestingly, abiotic stresses led to the re-directing of transcriptional output into non-productive pathways defined by the class of poly(A) site utilized. This result revealed APA to be part of a larger global response of Sorghum to abiotic stresses that involves the re-direction of transcriptional output into non-productive transcriptional and translational pathways. Large numbers of stress-inducible poly(A) sites could not be linked with known, annotated genes, suggestive of the existence of numerous unidentified genes whose expression is strongly regulated by abiotic stresses. Furthermore, we uncovered a novel stress-specific cis-element in intronic poly(A) sites used in drought- and heat-stressed plants that might play an important role in non-canonical poly(A) site choice in response to abiotic stresses.

Characterization of mRNA polyadenylation in the apicomplexa

Messenger RNA polyadenylation is a universal aspect of gene expression in eukaryotes. In well-established model organisms, this process is mediated by a conserved complex of 15–20 subunits. To better understand this process in apicomplexans, a group of unicellular parasites that causes serious disease in humans and livestock, a computational and high throughput sequencing study of the polyadenylation complex and poly(A) sites in several species was conducted. BLAST-based searches for orthologs of the human polyadenylation complex yielded clear matches to only two—poly(A) polymerase and CPSF73—of the 19 proteins used as queries in this analysis. As the human subunits that recognize the AAUAAA polyadenylation signal (PAS) were not immediately obvious, a computational analysis of sequences adjacent to experimentally-determined apicomplexan poly(A) sites was conducted. The results of this study showed that there exists in apicomplexans an A-rich region that corresponds in position to the AAUAAA PAS. The set of experimentally-determined sites in one species, Sarcocystis neurona, was further analyzed to evaluate the extent and significance of alternative poly(A) site choice in this organism. The results showed that almost 80% of S. neurona genes possess more than one poly(A) site, and that more than 780 sites showed differential usage in the two developmental stages–extracellular merozoites and intracellular schizonts–studied. These sites affected more than 450 genes, and included a disproportionate number of genes that encode membrane transporters and ribosomal proteins. Taken together, these results reveal that apicomplexan species seem to possess a poly(A) signal analogous to AAUAAA even though genes that may encode obvious counterparts of the AAUAAA-recognizing proteins are absent in these organisms. They also indicate that, as is the case in other eukaryotes, alternative polyadenylation is a widespread phenomenon in S. neurona that has the potential to impact growth and development.

Pipecolic acid confers systemic immunity by regulating free radicals

Pipecolic acid (Pip), a non-proteinaceous product of lysine catabolism, is an important regulator of immunity in plants and humans alike. In plants, Pip accumulates upon pathogen infection and has been associated with systemic acquired resistance (SAR). However, the molecular mechanisms underlying Pip-mediated signaling and its relationship to other known SAR inducers remain unknown. We show that in plants, Pip confers SAR by increasing levels of the free radicals, nitric oxide (NO), and reactive oxygen species (ROS), which act upstream of glycerol-3-phosphate (G3P). Plants defective in NO, ROS, G3P, or salicylic acid (SA) biosynthesis accumulate reduced Pip in their distal uninfected tissues although they contain wild-type-like levels of Pip in their infected leaves. These data indicate that de novo synthesis of Pip in distal tissues is dependent on both SA and G3P and that distal levels of SA and G3P play an important role in SAR. These results also suggest a unique scenario whereby metabolites in a signaling cascade can stimulate each other's biosynthesis depending on their relative levels and their site of action.

In vivo chlorophyll fluorescence screening allows the isolation of a Chlamydomonas mutant defective for NDUFAF3, an assembly factor involved in mitochondrial complex I assembly

The qualitative screening method used to select complex I mutants in the microalga Chlamydomonas, based on reduced growth under heterotrophic conditions, is not suitable for high-throughput screening. In order to develop a fast screening method based on measurements of chlorophyll fluorescence, we first demonstrated that complex I mutants displayed decreased photosystem II efficiency in the genetic background of a photosynthetic mutation leading to reduced formation of the electrochemical proton gradient in the chloroplast (pgrl1 mutation). In contrast, single mutants (complex I and pgrl1 mutants) could not be distinguished from the wild type by their photosystem II efficiency under the conditions tested. We next performed insertional mutagenesis on the pgrl1 mutant. Out of about 3000 hygromycin-resistant insertional transformants, 46 had decreased photosystem II efficiency and three were complex I mutants. One of the mutants was tagged and whole genome sequencing identified the resistance cassette in NDUFAF3, a homolog of the human NDUFAF3 gene, encoding for an assembly factor involved in complex I assembly. Complemented strains showed restored complex I activity and assembly. Overall, we describe here a screening method which is fast and particularly suited for the identification of Chlamydomonas complex I mutants.

Noncanonical Alternative Polyadenylation Contributes to Gene Regulation in Response to Hypoxia

Stresses from various environmental challenges continually confront plants, and their responses are important for growth and survival. One molecular response to such challenges involves the alternative polyadenylation of mRNA. In plants, it is unclear how stress affects the production and fate of alternative mRNA isoforms. Using a genome-scale approach, we show that in Arabidopsis thaliana, hypoxia leads to increases in the number of mRNA isoforms with polyadenylated 3' ends that map to 5'-untranslated regions (UTRs), introns, and protein-coding regions. RNAs with 3' ends within protein-coding regions and introns were less stable than mRNAs that end at 3'-UTR poly(A) sites. Additionally, these RNA isoforms were underrepresented in polysomes isolated from control and hypoxic plants. By contrast, mRNA isoforms with 3' ends that lie within annotated 5'-UTRs were overrepresented in polysomes and were as stable as canonical mRNA isoforms. These results indicate that the generation of noncanonical mRNA isoforms is an important feature of the abiotic stress response. The finding that several noncanonical mRNA isoforms are relatively unstable suggests that the production of non-stop and intronic mRNA isoforms may represent a form of negative regulation in plants, providing a conceptual link with mechanisms that generate these isoforms (such as alternative polyadenylation) and RNA surveillance.

PlantAPA: A Portal for Visualization and Analysis of Alternative Polyadenylation in Plants

Alternative polyadenylation (APA) is an important layer of gene regulation that produces mRNAs that have different 3' ends and/or encode diverse protein isoforms. Up to 70% of annotated genes in plants undergo APA. Increasing numbers of poly(A) sites collected in various plant species demand new methods and tools to access and mine these data. We have created an open-access web service called PlantAPA (http://bmi.xmu.edu.cn/plantapa) to visualize and analyze genome-wide poly(A) sites in plants. PlantAPA provides various interactive and dynamic graphics and seamlessly integrates a genome browser that can profile heterogeneous cleavage sites and quantify expression patterns of poly(A) sites across different conditions. Particularly, through PlantAPA, users can analyze poly(A) sites in extended 3' UTR regions, intergenic regions, and ambiguous regions owing to alternative transcription or RNA processing. In addition, it also provides tools for analyzing poly(A) site selections, 3' UTR lengthening or shortening, non-canonical APA site switching, and differential gene expression between conditions, making it more powerful for the study of APA-mediated gene expression regulation. More importantly, PlantAPA offers a bioinformatics pipeline that allows users to upload their own short reads or ESTs for poly(A) site extraction, enabling users to further explore poly(A) site selection using stored PlantAPA poly(A) sites together with their own poly(A) site datasets. To date, PlantAPA hosts the largest database of APA sites in plants, including Oryza sativa, Arabidopsis thaliana, Medicago truncatula, and Chlamydomonas reinhardtii. As a user-friendly web service, PlantAPA will be a valuable addition to the community of biologists studying APA mechanisms and gene expression regulation in plants.