Analysis of Chlamydomonas reinhardtii genome structure using large-scale sequencing of regions on linkage groups I and III

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

Alternatively spliced introns are the ones that are usually spliced but can be occasionally retained in a transcript isoform. They are the most frequently used alternative splice form in plants (~50% of alternative splicing events). Chlamydomonas reinhardtii, a unicellular alga, is a good model to understand alternative splicing (AS) in plants from an evolutionary perspective as it diverged from land plants a billion years ago. Using over 7 million cDNA sequences from both pyrosequencing and Sanger sequencing, we found that a much higher percentage of genes (~20% of multi-exon genes) undergo AS than previously reported (3-5%). We found a full component of SR and SR-like proteins possibly involved in AS. The most prevalent type of AS event (40%) was retention of introns, most of which were supported by multiple cDNA evidence (72%) while only 20% of them have coding capacity. By comparing retained and constitutive introns, we identified sequence features potentially responsible for the retention of introns, in the framework of an "intron definition" model for splicing. We find that retained introns tend to have a weaker 5' splice site, more Gs in their poly-pyrimidine tract and a lesser conservation of nucleotide 'C' at position -3 of the 3' splice site. In addition, the sequence motifs found in the potential branch-point region differed between retained and constitutive introns. Furthermore, the enrichment of G-triplets and C-triplets among the first and last 50 nt of the introns significantly differ between constitutive and retained introns. These could serve as intronic splicing enhancers. All the alternative splice forms can be accessed at http://bioinfolab.miamioh.edu/cgi-bin/PASA_r20140417/cgi-bin/status_report.cgi?db=Chre_AS .

Mapping of the css (chloroplast splicing suppressor) gene(s) to a recombinationally suppressed region of chromosome III in Chlamydomonas reinhardtii

In previous work, three suppressors of defective group I introns (7151, 71N1, 7120) were isolated from a mutant of Chlamydomonas reinhardtii that had a splicing-deficient chloroplast large subunit (LSU) rRNA intron. Genetic analysis indicated that the 7151 and 71N1 suppressor mutations each involved single nuclear loci, and that the 7151 mutation was dominant. Here we present genetic evidence that the 7120 suppressor also involves a single nuclear locus and that the mutation is dominant in vegetative diploids. Moreover, we have employed crosses with the S1D2 strain and molecular markers to map the 7120 and 71N1 suppressors. Based on an analysis of 800 progeny from 7120 × S1D2, the 7120 suppressor is located in a region of ~400 kb on chromosome III that is devoid of recombination. The ~400-kb region contains at least 72 genes, about one-third of which (i.e., 22) are predicted to be organelle targeted. Similar analysis of 71N1 × S1D2 using 400 progeny also pointed to the recombination-deficient region of chromosome III, raising the possibility that these mutations could affect the same gene. These efforts lay the foundation for identifying the css (chloroplast splicing suppressor) gene(s), which promotes splicing of multiple chloroplast group I introns.

Regulation of ciliary motility: conserved protein kinases and phosphatases are targeted and anchored in the ciliary axoneme

Recent evidence has revealed that the dynein motors and highly conserved signaling proteins are localized within the ciliary 9+2 axoneme. One key mechanism for regulation of motility is phosphorylation. Here, we review diverse evidence, from multiple experimental organisms, that ciliary motility is regulated by phosphorylation/dephosphorylation of the dynein arms through kinases and phosphatases that are anchored immediately adjacent to their axonemal substrates.

Improving gene-finding in Chlamydomonas reinhardtii:GreenGenie2

BACKGROUND

The availability of whole-genome sequences allows for the identification of the entire set of protein coding genes as well as their regulatory regions. This can be accomplished using multiple complementary methods that include ESTs, homology searches and ab initio gene predictions. Previously, the Genie gene-finding algorithm was trained on a small set of Chlamydomonas genes and shown to improve the accuracy of gene prediction in this species compared to other available programs. To improve ab initio gene finding in Chlamydomonas, we assemble a new training set consisting of over 2,300 cDNAs by assembling over 167,000 Chlamydomonas EST entries in GenBank using the EST assembly tool PASA.

RESULTS

The prediction accuracy of our cDNA-trained gene-finder, GreenGenie2, attains 83% sensitivity and 83% specificity for exons on short-sequence predictions. We predict about 12,000 genes in the version v3 Chlamydomonas genome assembly, most of which (78%) are either identical to or significantly overlap the published catalog of Chlamydomonas genes 1. 22% of the published catalog is absent from the GreenGenie2 predictions; there is also a fraction (23%) of GreenGenie2 predictions that are absent from the published gene catalog. Randomly chosen gene models were tested by RT-PCR and most support the GreenGenie2 predictions.

CONCLUSION

These data suggest that training with EST assemblies is highly effective and that GreenGenie2 is a valuable, complementary tool for predicting genes in Chlamydomonas reinhardtii.

Insertional mutagenesis as a tool to study genes/functions in Chlamydomonas

The unicellular alga Chlamydomonas reinhardtii has emerged during the last decades as a model system to understand gene functions, many of them shared by bacteria, fungi, plants, animals and humans. A powerful resource for the research community is the availability of complete collections of stable mutants for studying whole genome function. In the meantime other strategies might be developed; insertional mutagenesis has become currently the best strategy to disrupt and tag nuclear genes in Chlamydomonas allowing forward and reverse genetic approaches. Here, we outline the mutagenesis technique stressing the idea of generating databases for ordered mutant libraries, and also of improving efficient methods for reverse genetics to identify mutants defective in a particular gene.

EST assembly supported by a draft genome sequence: an analysis of the Chlamydomonas reinhardtii transcriptome

Clustering and assembly of expressed sequence tags (ESTs) constitute the basis for most genomewide descriptions of a transcriptome. This approach is limited by the decline in sequence quality toward the end of each EST, impacting both sequence clustering and assembly. Here, we exploit the available draft genome sequence of the unicellular green alga Chlamydomonas reinhardtii to guide clustering and to correct errors in the ESTs. We have grouped all available EST and cDNA sequences into 12 063 ACEGs (assembly of contiguous ESTs based on genome) and generated 15 857 contigs of average length 934 nt. We predict that roughly 3000 of our contigs represent full-length transcripts. Compared to previous assemblies, ACEGs show extended contig length, increased accuracy and a reduction in redundancy. Because our assembly protocol also uses ESTs with no corresponding genomic sequences, it provides sequence information for genes interrupted by sequence gaps. Detailed analysis of randomly sampled ACEGs reveals several hundred putative cases of alternative splicing, many overlapping transcription units and new genes not identified by gene prediction algorithms. Our protocol, although developed for and tailored to the C. reinhardtii dataset, can be exploited by any eukaryotic genome project for which both a draft genome sequence and ESTs are available.

Comparative genomic analysis of fungal genomes reveals intron-rich ancestors

BACKGROUND

Eukaryotic protein-coding genes are interrupted by spliceosomal introns, which are removed from transcripts before protein translation. Many facets of spliceosomal intron evolution, including age, mechanisms of origins, the role of natural selection, and the causes of the vast differences in intron number between eukaryotic species, remain debated. Genome sequencing and comparative analysis has made possible whole genome analysis of intron evolution to address these questions.

RESULTS

We analyzed intron positions in 1,161 sets of orthologous genes across 25 eukaryotic species. We find strong support for an intron-rich fungus-animal ancestor, with more than four introns per kilobase, comparable to the highest known modern intron densities. Indeed, the fungus-animal ancestor is estimated to have had more introns than any of the extant fungi in this study. Thus, subsequent fungal evolution has been characterized by widespread and recurrent intron loss occurring in all fungal clades. These results reconcile three previously proposed methods for estimation of ancestral intron number, which previously gave very different estimates of ancestral intron number for eight eukaryotic species, as well as a fourth more recent method. We do not find a clear inverse correspondence between rates of intron loss and gain, contrary to the predictions of selection-based proposals for interspecific differences in intron number.

CONCLUSION

Our results underscore the high intron density of eukaryotic ancestors and the widespread importance of intron loss through eukaryotic evolution.

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.

Natural history of transposition in the green alga Chlamydomonas reinhardtii: use of the AMT4 locus as an experimental system

The AMT4 locus of the green alga Chlamydomonas reinhardtii, which we mapped to the long arm of chromosome 8, provides a good experimental system for the study of transposition. Most mutations that confer resistance to the toxic ammonium analog methylammonium are in AMT4 and a high proportion of spontaneous mutations are caused by transposon-related events. Among the 15 such events that we have characterized at the molecular level, 9 were associated with insertions of the retrotransposon TOC1, 2 with a small Gulliver-related transposon, and 1 with the Tcr1 transposon. We found that Tcr1 is apparently a foldback transposon with terminal inverted repeats that are much longer and more complex than previously realized. A duplication of Tcr1 yielded a configuration thought to be important for chromosomal evolution. Other mutations in AMT4 were caused by two mobile elements that have not been described before. The sequence of one, which we propose to call the Bill element, indicates that it probably transposes by way of a DNA intermediate and requires functions that it does not encode. The sequence of the other and bioinformatic analysis indicates that it derives from a miniature retrotransposon or TRIM, which we propose to call MRC1 (miniature retrotransposon of Chlamydomonas).

Very little intron loss/gain in Plasmodium: intron loss/gain mutation rates and intron number

We compared intron positions in conserved regions of 3479 orthologous gene pairs from Plasmodium falciparum and Plasmodium yoelii, which likely diverged >or=100 million years ago (Mya). Only 27 out of 2212 positions were specific to one of the two species. Intron presence in related species shows that at least 19 and possibly 26 of the changes are due to intron loss, depending on phylogeny. The implied intron loss and gain rates are much lower than previously estimated for nematodes, arthropods, fungi, and plants, and are comparable only with the rates in vertebrates. That all observed changes were exact, occurring without loss or gain of flanking coding sequence, suggests intron loss via an mRNA intermediate, as does a nonsignificant trend toward loss of introns at adjacent positions. Many of the intron changes occurred in genes encoding proteins involved in nucleic acid-related processes, as previously found for intron gains in nematodes. Two changes occurred in the chloroquine resistance transporter, suggesting a role for positive selection in intron loss in Plasmodium. The dearth of intron loss and gain could be explained by the lack of known transposable elements in Plasmodium, since transposable elements and/or reverse transcriptase are thought to be necessary for both processes. The observed pattern suggests that the availability of stochastic intron loss and gain mutations can be a major determinant of changes in intron number.

Beyond complementation. Map-based cloning in Chlamydomonas reinhardtii

Chlamydomonas reinhardtii is an excellent model system for plant biologists because of its ease of manipulation, facile genetics, and the ability to transform the nuclear, chloroplast, and mitochondrial genomes. Numerous forward genetics studies have been performed in Chlamydomonas, in many cases to elucidate the regulation of photosynthesis. One of the resultant challenges is moving from mutant phenotype to the gene mutation causing that phenotype. To date, complementation has been the primary method for gene cloning, but this is impractical in several situations, for example, when the complemented strain cannot be readily selected or in the case of recessive suppressors that restore photosynthesis. New tools, including a molecular map consisting of 506 markers and an 8X-draft nuclear genome sequence, are now available, making map-based cloning increasingly feasible. Here we discuss advances in map-based cloning developed using the strains mcd4 and mcd5, which carry recessive nuclear suppressors restoring photosynthesis to chloroplast mutants. Tools that have not been previously applied to Chlamydomonas, such as bulked segregant analysis and marker duplexing, are being implemented to increase the speed at which one can go from mutant phenotype to gene. In addition to assessing and applying current resources, we outline anticipated future developments in map-based cloning in the context of the newly extended Chlamydomonas genome initiative.

Comparative profiling of lipid-soluble antioxidants and transcripts reveals two phases of photo-oxidative stress in a xanthophyll-deficient mutant of Chlamydomonas reinhardtii

Excess light can impose severe oxidative stress on photosynthetic organisms. We have characterized high-light responses in wild-type Chlamydomonas reinhardtii and in the npq1 lor1 double mutant. The npq1 lor1 strain lacks two photoprotective carotenoids, lutein and zeaxanthin, and experiences acute photo-oxidative stress upon exposure to excess light. To examine the ability of npq1 lor1 cells to respond to photo-oxidative stress, we measured changes in lipid-soluble antioxidants following a shift from low light to high light in the wild type and the double mutant. The size of the xanthophyll cycle pool increased in both the wild type and mutant during the first 6 h of exposure to high light levels, but then decreased in the mutant during photo-oxidative bleaching. The level of alpha-tocopherol (vitamin E) was constant in the wild type and mutant during the first 6 h; then it increased by three-fold in the wild type but declined in npq1 lor1 cells. We also used cDNA microarrays and RNA gel-blot analysis to monitor differences in gene expression. Both strains showed an initial light-stress response in the form of a transient increase in expression of (1) GPXH, a glutathione peroxidase gene that has been shown to respond specifically to singlet oxygen and lipid peroxidation; (2) SMT1, a gene for a putative sterol C-methyltransferase; and (3) LI818r, a stress-responsive member of the light-harvesting complex superfamily. These transient changes in gene expression in high light were followed by a second series of changes in npq1 lor1, coincident with declines in lipid-soluble antioxidants but preceding detectable photo-oxidative damage to proteins and lipids. Thus, the response of npq1 lor1 to high light is unexpectedly complex, with initial changes in lipid-soluble antioxidants and RNA levels that are associated with acclimation in the wild type and a second wave of changes that accompanies photo-oxidative bleaching.

Comparative genomics identifies a flagellar and basal body proteome that includes the BBS5 human disease gene

Cilia and flagella are microtubule-based structures nucleated by modified centrioles termed basal bodies. These biochemically complex organelles have more than 250 and 150 polypeptides, respectively. To identify the proteins involved in ciliary and basal body biogenesis and function, we undertook a comparative genomics approach that subtracted the nonflagellated proteome of Arabidopsis from the shared proteome of the ciliated/flagellated organisms Chlamydomonas and human. We identified 688 genes that are present exclusively in organisms with flagella and basal bodies and validated these data through a series of in silico, in vitro, and in vivo studies. We then applied this resource to the study of human ciliation disorders and have identified BBS5, a novel gene for Bardet-Biedl syndrome. We show that this novel protein localizes to basal bodies in mouse and C. elegans, is under the regulatory control of daf-19, and is necessary for the generation of both cilia and flagella.

Identification and sequence analysis of six new members of the NIMA-related kinase family in Chlamydomonas

The NIMA kinases are an evolutionarily conserved protein family with enigmatic roles in the regulation of mitosis. We report six new members of this family in Chlamydomonas, in addition to the previously identified NIMA-related kinase, Fa2p. Chlamydomonas NIMA-related kinases (CNKs) 1-6 were sequenced from subclones generated by RT-PCR using information from EST libraries and the recently sequenced Chlamydomonas genome. Phylogenetic and bioinformatic approaches were used to determine the relationships of the six new members with known members of the NIMA-related kinase family. Although humans express at least eleven NIMA-related kinases, the eukaryotic microbes that have been studied to date express only one or two members of the family. Thus, the discovery that Chlamydomonas expresses a total of at least seven NIMA-related kinases is intriguing. Our analyses suggest that members of this family may play roles in the assembly and function of cilia.

Oda5p, a novel axonemal protein required for assembly of the outer dynein arm and an associated adenylate kinase

Of the uncloned ODA genes required for outer dynein arm assembly in Chlamydomonas, ODA5 and ODA10 are of particular interest because they do not encode known subunits of the outer arm or the outer dynein arm-docking complex (ODA-DC), and because genetic studies suggest their products interact. Beginning with a tagged oda5 allele, we isolated genomic and cDNA clones of the wild-type gene. ODA5 predicts a novel, 66-kDa coiled-coil protein. Immunoblotting indicates Oda5p is an axonemal component that assembles onto the axoneme independently of the outer arm and ODA-DC and is uniquely missing in oda5 and oda10 axonemes. Oda5p is released from the axoneme by extraction with 0.6 M KCl, but the soluble Oda5p does not cosediment with the outer dynein arm/ODA-DC in sucrose gradients. Quantitative mass spectrometry by using isotope coded affinity tagging revealed that a previously unidentified adenylate kinase is reduced 35-50% in oda5 flagella. Direct enzymatic assays demonstrated a comparable reduction in adenylate kinase activity in oda5 flagella, and also in oda10 flagella, but not in flagella of other oda mutants. We propose that Oda5p is part of a novel axonemal complex that is required for outer arm assembly and anchors adenylate kinase in proximity to the arm.

Genomic analysis of mutants affecting xanthophyll biosynthesis and regulation of photosynthetic light harvesting in Chlamydomonas reinhardtii

When the absorption of light energy exceeds the capacity for its utilization in photosynthesis, regulation of light harvesting is critical in order for photosynthetic organisms to minimize photo-oxidative damage. Thermal dissipation of excess absorbed light energy, measured as non-photochemical quenching (NPQ) of chlorophyll fluorescence, is induced rapidly in response to excess light conditions, and it is known that xanthophylls such as zeaxanthin and lutein, the transthylakoid pH gradient, and the PsbS protein are involved in this mechanism. Although mutants affecting NPQ and the biosynthesis of zeaxanthin and lutein were originally isolated and characterized at the physiological level in the unicellular green alga Chlamydomonas reinhardtii, the molecular basis of several of these mutants, such as npq1 and lor1, has not been determined previously. The recent sequencing of the C. reinhardtii nuclear genome has facilitated the search for C. reinhardtii homologs of plant genes involved in xanthophyll biosynthesis and regulation of light harvesting. Here we report the identification of C. reinhardtii genes encoding PsbS and lycopene varepsilon-cyclase, and we show that the lor1 mutation, which affects lutein synthesis, is located within the lycopene varepsilon-cyclase gene. In contrast, no homolog of the plant violaxanthin de-epoxidase (VDE) gene was found. Molecular markers were used to map the npq1 mutation, which affects VDE activity, as a first step toward the map-based cloning of the NPQ1 gene.

Molecular markers for rapidly identifying candidate genes in Chlamydomonas reinhardtii. Ery1 and ery2 encode chloroplast ribosomal proteins

To take advantage of available expressed sequence tags and genomic sequence, we have developed 64 PCR-based molecular markers in Chlamydomonas reinhardtii that map to the 17 linkage groups. These markers will allow the rapid association of a candidate gene sequence with previously identified mutations. As proof of principle, we have identified the genes encoded by the ERY1 and ERY2 loci. Mendelian mutations that confer resistance to erythromycin define three unlinked nuclear loci in C. reinhardtii. Candidate genes ribosomal protein L4 (RPL4) and L22 (RPL22) are tightly linked to the ERY1 locus and ERY2 locus, respectively. Genomic DNA for RPL4 from wild type and five mutant ery1 alleles was amplified and sequenced and three different point mutations were found. Two different glycine residues (G(102) and G(112)) are replaced by aspartic acid and both are in the unstructured region of RPL4 that lines the peptide exit tunnel of the chloroplast ribosome. The other two alleles change a splice site acceptor site. Genomic DNA for RPL22 from wild type and three mutant ery2 alleles was amplified and sequenced and revealed three different point mutations. Two alleles have premature stop codons and one allele changes a splice site acceptor site.