In silico analysis of the sequence features responsible for alternatively spliced introns in the model green alga Chlamydomonas reinhardtii

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

Transcriptome Analysis Reveals the Involvement of Alternative Splicing in the Nitrogen Starvation Response of Chlamydomonas reinhardtii

Alternative splicing (AS) is a regulatory mechanism of post-transcriptional regulation that plays an important role in plant response to abiotic stresses. However, corresponding research involving the mechanism of AS in the nitrogen starvation response of C. reinhardtii is rare. This study performed a comprehensive and systematic analysis of AS events in C. reinhardtii at nine time points (0 h, 10 m, 30 m, 1 h, 6 h, 8 h, 24 h, and 48 h) under nitrogen starvation. It used STAR and rMATS tools to identify and quantify the probability of the AS event happening through the transcriptome high-throughput sequencing data. A total of 5806 AS events in 3500 genes were identified, and the retained intron and skipped exon were considered the main AS types. The genes related to the AS event in nitrogen starvation were mainly involved in spliceosome and transporter and enriched in the citrate cycle and fatty acid degradation pathways. These results suggested that AS may play an important role in the nitrogen starvation response in C. reinhardtii, and provided insights into post-transcriptional regulation under nitrogen starvation.

Comparative genomics of Chlamydomonas

Despite its role as a reference organism in the plant sciences, the green alga Chlamydomonas reinhardtii entirely lacks genomic resources from closely related species. We present highly contiguous and well-annotated genome assemblies for three unicellular C. reinhardtii relatives: Chlamydomonas incerta, Chlamydomonas schloesseri, and the more distantly related Edaphochlamys debaryana. The three Chlamydomonas genomes are highly syntenous with similar gene contents, although the 129.2 Mb C. incerta and 130.2 Mb C. schloesseri assemblies are more repeat-rich than the 111.1 Mb C. reinhardtii genome. We identify the major centromeric repeat in C. reinhardtii as a LINE transposable element homologous to Zepp (the centromeric repeat in Coccomyxa subellipsoidea) and infer that centromere locations and structure are likely conserved in C. incerta and C. schloesseri. We report extensive rearrangements, but limited gene turnover, between the minus mating type loci of these Chlamydomonas species. We produce an eight-species core-Reinhardtinia whole-genome alignment, which we use to identify several hundred false positive and missing genes in the C. reinhardtii annotation and >260,000 evolutionarily conserved elements in the C. reinhardtii genome. In summary, these resources will enable comparative genomics analyses for C. reinhardtii, significantly extending the analytical toolkit for this emerging model system.

Alternative Splicing During the Chlamydomonas reinhardtii Cell Cycle

Genome-wide analysis of transcriptome data in Chlamydomonas reinhardtii shows periodic patterns in gene expression levels when cultures are grown under alternating light and dark cycles so that G1 of the cell cycle occurs in the light phase and S/M/G0 occurs during the dark phase. However, alternative splicing, a process that enables a greater protein diversity from a limited set of genes, remains largely unexplored by previous transcriptome based studies in C. reinhardtii. In this study, we used existing longitudinal RNA-seq data obtained during the light-dark cycle to investigate the changes in the alternative splicing pattern and found that 3277 genes (19.75% of 17,746 genes) undergo alternative splicing. These splicing events include Alternative 5′ (Alt 5′), Alternative 3′ (Alt 3′) and Exon skipping (ES) events that are referred as alternative site selection (ASS) events and Intron retention (IR) events. By clustering analysis, we identified a subset of events (26 ASS events and 10 IR events) that show periodic changes in the splicing pattern during the cell cycle. About two-thirds of these 36 genes either introduce a pre-termination codon (PTC) or introduce insertions or deletions into functional domains of the proteins, which implicate splicing in altering gene function. These findings suggest that alternative splicing is also regulated during the Chlamydomonas cell cycle, although not as extensively as changes in gene expression. The longitudinal changes in the alternative splicing pattern during the cell cycle captured by this study provides an important resource to investigate alternative splicing in genes of interest during the cell cycle in Chlamydomonas reinhardtii and other eukaryotes.

Introns mediate post-transcriptional enhancement of nuclear gene expression in the green microalga Chlamydomonas reinhardtii

Efficient nuclear transgene expression in the green microalga Chlamydomonas reinhardtii is generally hindered by low transcription rates. Introns can increase transcript abundance by a process called Intron-Mediated Enhancement (IME) in this alga and has been broadly observed in other eukaryotes. However, the mechanisms of IME in microalgae are poorly understood. Here, we identified 33 native introns from highly expressed genes in C. reinhardtii selected from transcriptome studies as well as 13 non-native introns. We investigated their IME capacities and probed the mechanism of action by modification of splice sites, internal sequence motifs, and position within transgenes. Several introns were found to elicit strong IME and found to be broadly applicable in different expression constructs. We determined that IME in C. reinhardtii exclusively occurs from introns within transcribed ORFs regardless of the promoter and is not induced by traditional enhancers of transcription. Our results elucidate some mechanistic details of IME in C. reinhardtii, which are similar to those observed in higher plants yet underly distinctly different induction processes. Our findings narrow the focus of targets responsible for algal IME and provides evidence that introns are underestimated regulators of C. reinhardtii nuclear gene expression.

Although many genetic tools and basic transformation strategies exist for the model microalga Chlamydomonas reinhardtii, high-level genetic engineering with this organism is hindered by its inherent recalcitrance to foreign gene expression and limited knowledge of responsible expression regulators. In this work, we characterized the dynamics of 33 endogenous and 13 non-native introns and their effect on gene expression as artificial insertions into codon optimized transgenes. We found that introns from different origins have the capacity to increase gene expression rates. Intron-mediated enhancement was observed exclusively when these elements were placed in transcripts but not outside of transcribed mRNA regions. Insertion of different endogenous introns into coding sequences was found to positively affect expression rates through a synergy of additive transcription enhancement and exon length reduction, similar to those natively found in the C. reinhardtii genome. Our results indicate that intensive mRNA processing plays an underestimated role in the regulation of native gene expression in C. reinhardtii. In addition to internal sequence motifs, the location of artificially introduced introns greatly affected transgene expression levels. This work is highly valuable to the greater microalgal and synthetic biology research communities and contributes to broadening our understanding of eukaryotic intron-mediated enhancement.

A Musashi Splice Variant and Its Interaction Partners Influence Temperature Acclimation in Chlamydomonas

Microalgae contribute significantly to carbon fixation on Earth. Global warming influences their physiology and growth rates. To understand algal short-term acclimation and adaptation to changes in ambient temperature, it is essential to identify and characterize the molecular components that sense small temperature changes as well as the downstream signaling networks and physiological responses. Here, we used the green biflagellate alga Chlamydomonas reinhardtii as a model system in which to study responses to temperature. We report that an RNA recognition motif (RRM)-containing RNA-binding protein, Musashi, occurs in 25 putative splice variants. These variants bear one, two, and three RRM domains or even lack RRM domains. The most abundant Musashi variant, 12, with a molecular mass of 60 kD, interacts with two clock-relevant members of RNA metabolism, the subunit C3 of the RNA-binding protein CHLAMY1 and the 5'-3' exoribonuclease XRN1. These proteins are able to integrate temperature information by up- or down-regulation of their protein levels in cells grown at low (18°C) or high (28°C) temperature. We further show that the 60-kD Musashi variants with three RRM domains can bind to (UG)₇ repeat-containing RNAs and are up-regulated in cells grown at a higher temperature during early night. Intriguingly, the 60-kD Musashi variant 12, as well as C3 and XRN1, confer thermal acclimation to C. reinhardtii, as shown with mutant lines. Our data suggest that these three proteins of the RNA metabolism machinery are key members of the thermal signaling network in C. reinhardtii.

Genome-wide analysis of the invertase gene family from maize

The recent release of the maize genome (AGPv4) contains annotation errors of invertase genes and therefore the enzymes are bestly curated manually at the protein level in a comprehensible fashion The synthesis, transport and degradation of sucrose are determining factors for biomass allocation and yield of crop plants. Invertase (INV) is a key enzyme of carbon metabolism in both source and sink tissues. Current releases of the maize genome correctly annotates only two vacuolar invertases (ivr1 and ivr2) and four cell wall invertases (incw1, incw2 (mn1), incw3, and incw4). Our comprehensive survey identified 21 INV isogenes for which we propose a standard nomenclature grouped phylogenetically by amino acid similarity: three vacuolar (INVVR), eight cell wall (INVCW), and ten alkaline/neutral (INVAN) isogenes which form separate dendogram branches due to distinct molecular features. The acidic enzymes were curated for the presence of the DPN tripeptide which is coded by one of the smallest exons reported in plants. Particular attention was placed on the molecular role of INV in vascular tissues such as the nodes, internodes, leaf sheath, husk leaves and roots. We report the expression profile of most members of the maize INV family in nine tissues in two developmental stages, R1 and R3. INVCW7, INVVR2, INVAN8, INVAN9, INVAN10, and INVAN3 displayed the highest absolute expressions in most tissues. INVVR3, INVCW5, INVCW8, and INVAN1 showed low mRNA levels. Expressions of most INVs were repressed from stage R1 to R3, except for INVCW7 which increased significantly in all tissues after flowering. The mRNA levels of INVCW7 in the vegetative stem correlated with a higher transport rate of assimilates from leaves to the cob which led to starch accumulation and growth of the female reproductive organs.

The evolution of posttranscriptional regulation

"DNA makes RNA makes protein." After transcription, mRNAs undergo a series of intertwining processes to be finally translated into functional proteins. The "posttranscriptional" regulation (PTR) provides cells an extended option to fine-tune their proteomes. To meet the demands of complex organism development and the appropriate response to environmental stimuli, every step in these processes needs to be finely regulated. Moreover, changes in these regulatory processes are important driving forces underlying the evolution of phenotypic differences across different species. The major PTR mechanisms discussed in this review include the regulation of splicing, polyadenylation, decay, and translation. For alternative splicing and polyadenylation, we mainly discuss their evolutionary dynamics and the genetic changes underlying the regulatory differences in cis-elements versus trans-factors. For mRNA decay and translation, which, together with transcription, determine the cellular RNA or protein abundance, we focus our discussion on how their divergence coordinates with transcriptional changes to shape the evolution of gene expression. Then to highlight the importance of PTR in the evolution of higher complexity, we focus on their roles in two major phenomena during eukaryotic evolution: the evolution of multicellularity and the division of labor between different cell types and tissues; and the emergence of diverse, often highly specialized individual phenotypes, especially those concerning behavior in eusocial insects. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution Translation > Translation Regulation RNA Processing > Splicing Regulation/Alternative Splicing.