Mobile self-splicing group I introns from the psbA gene of Chlamydomonas reinhardtii: highly efficient homing of an exogenous intron containing its own promoter

Promising prospects of nanopore sequencing for algal hologenomics and structural variation discovery

BACKGROUND

The MinION Access Program (MAP, 2014-2016) allowed selected users to test the prospects of long nanopore reads for diverse organisms and applications through the rapid development of improving chemistries. In 2014, faced with a fragmented Illumina assembly for the chloroplast genome of the green algal holobiont Caulerpa ashmeadii, we applied to the MAP to test the prospects of nanopore reads to investigate such intricacies, as well as further explore the hologenome of this species with native and hybrid approaches.

RESULTS

The chloroplast genome could only be resolved as a circular molecule in nanopore assemblies, which also revealed structural variants (i.e. chloroplast polymorphism or heteroplasmy). Signal and Illumina polishing of nanopore-assembled organelle genomes (chloroplast and mitochondrion) reflected the importance of coverage on final quality and current limitations. In hybrid assembly, our modest nanopore data sets showed encouraging results to improve assembly length, contiguity, repeat content, and binning of the larger nuclear and bacterial genomes. Profiling of the holobiont with nanopore or Illumina data unveiled a dominant Rhodospirillaceae (Alphaproteobacteria) species among six putative endosymbionts. While very fragmented, the cumulative hybrid assembly length of C. ashmeadii's nuclear genome reached 24.4 Mbp, including 2.1 Mbp in repeat, ranging closely with GenomeScope's estimate (> 26.3 Mbp, including 4.8 Mbp in repeat).

CONCLUSION

Our findings relying on a very modest number of nanopore R9 reads as compared to current output with newer chemistries demonstrate the promising prospects of the technology for the assembly and profiling of an algal hologenome and resolution of structural variation. The discovery of polymorphic 'chlorotypes' in C. ashmeadii, most likely mediated by homing endonucleases and/or retrohoming by reverse transcriptases, represents the first report of chloroplast heteroplasmy in the siphonous green algae. Improving contiguity of C. ashmeadii's nuclear and bacterial genomes will require deeper nanopore sequencing to greatly increase the coverage of these larger genomic compartments.

Expression of a Synthetic Gene for the Major Cytotoxin (Cyt1Aa) of Bacillus thuringiensis subsp. israelensis in the Chloroplast of Wild-Type Chlamydomonas

Chlamydomonas reinhardtii (Chlamydomonas) strains that are toxic to mosquito larvae because they express chloroplast transgenes that are based on the mosquitocidal proteins of Bacillus thuringiensis subsp. israelensis (Bti) could be very useful in mosquito control. Chlamydomonas has several advantages for this approach, including genetic controls not generally available with industrial algae. The Bti toxin is produced by sporulating bacteria and has been used for mosquito control for >30 years without creating highly resistant mosquito populations. The suite of toxins is four main proteins: three Cry proteins and the cytotoxic Cyt1Aa (27 kDa). Cyt1Aa is not very toxic to mosquitoes by itself, but it prevents the development of resistance. The production of Cyt1Aa in other microbes, however, has been challenging due to its affinity for certain membrane phospholipids. Here we report on the production of recombinant Cyt1Aa (rCyt1A) in the chloroplast of photosynthetic Chlamydomonas at levels of at least 0.3% total protein. Live cell bioassays demonstrated toxicity of the rCyt1Aa Chlamydomonas to larvae of Aedes aegypti. We also expressed the chloroplast cyt1Aa gene in a wild-type Chlamydomonas strain (21 gr) that can grow on nitrate. These results have implications for developing a Chlamydomonas strain that will be toxic to mosquito larvae but will not induce strongly resistant populations.

Small RNA profiling in Chlamydomonas: insights into chloroplast RNA metabolism

In Chlamydomonas reinhardtii, regulation of chloroplast gene expression is mainly post-transcriptional. It requires nucleus-encoded trans-acting protein factors for maturation/stabilization (M factors) or translation (T factors) of specific target mRNAs. We used long- and small-RNA sequencing to generate a detailed map of the transcriptome. Clusters of sRNAs marked the 5' end of all mature mRNAs. Their absence in M-factor mutants reflects the protection of transcript 5' end by the cognate factor. Enzymatic removal of 5'-triphosphates allowed identifying those cosRNA that mark a transcription start site. We detected another class of sRNAs derived from low abundance transcripts, antisense to mRNAs. The formation of antisense sRNAs required the presence of the complementary mRNA and was stimulated when translation was inhibited by chloramphenicol or lincomycin. We propose that they derive from degradation of double-stranded RNAs generated by pairing of antisense and sense transcripts, a process normally hindered by the traveling of the ribosomes. In addition, chloramphenicol treatment, by freezing ribosomes on the mRNA, caused the accumulation of 32-34 nt ribosome-protected fragments. Using this 'in vivo ribosome footprinting', we identified the function and molecular target of two candidate trans-acting factors.

Toward mosquito control with a green alga: Expression of Cry toxins of Bacillus thuringiensis subsp. israelensis (Bti) in the chloroplast of Chlamydomonas

We are developing Chlamydomonas strains that can be used for safe and sustainable control of mosquitoes, because they produce proteins from Bacillus thuringiensis subsp. israelensis (Bti) in the chloroplast. Chlamydomonas has a number of advantages for this approach, including genetic controls that are not generally available with industrial algae. The Bti toxin has been used for mosquito control for > 30 years and does not engender resistance; it contains three Cry proteins, Cry4Aa (135 kDa), Cry4Ba (128 kDa) and Cry11Aa (72 kDa), and Cyt1Aa (25 kDa). To express the Cry proteins in the chloroplast, the three genes were resynthesized and cry4Aa was truncated to the first 700 amino acids (cry4Aa₇₀₀ ); also, since they can be toxic to host cells, the inducible Cyc6:Nac2-psbD expression system was used. Western blots of total protein from the chloroplast transformants showed accumulation of the intact polypeptides, and the relative expression level was Cry11Aa > Cry4Aa₇₀₀ > Cry4Ba. Quantitative western blots with purified Cry11Aa as a standard showed that Cry11Aa accumulated to 0.35% of total cell protein. Live cell bioassays in dH₂0 demonstrated toxicity of the cry4Aa₇₀₀ and cry11Aa transformants to larvae of Aedes aegypti and Culex quinquefasciatus. These results demonstrate that the Cry proteins that are most toxic to Aedes and Culex mosquitoes, Cry4Aa and Cry11Aa, can be successfully expressed in the chloroplast of Chlamydomonas.

PCR analysis of chloroplast double-strand break (DSB) repair products induced by I-CreII in Chlamydomonas and Arabidopsis

Homing endonuclease I-CreII has been used to study the consequences and repair of a double-strand break (DSB) in the chloroplast genome of Chlamydomonas and Arabidopsis. Since I-CreII is from a mobile psbA intron of Chlamydomonas, it cleaves the psbA gene of an intronless-psbA strain of Chlamydomonas. And it cleaves specifically in the psbA gene of Arabidopsis, which is naturally intronless. We have shown further that most of the repair products of an I-CreII-induced break in chloroplast DNA can be defined by PCR analysis with total nucleic acids and the appropriate primers. Here, we provide protocols for small-scale preparation of nucleic acids from Chlamydomonas and Arabidopsis, as well as guidelines for the subsequent PCR analysis.

Reverse transcription of spliced psbA mRNA in Chlamydomonas spp. and its possible role in evolutionary intron loss

Reverse transcription of mRNA is thought to be an important first step in a model that explains certain evolutionary changes within genes, such as the loss of introns or RNA editing sites. In this model, reverse transcription of mRNA produces cDNA molecules that replace part of the parental gene by homologous recombination. In vivo evidence of reverse transcription of physiologically relevant mRNAs is generally lacking, however, except in genetically engineered cells. Here, we provide in vivo evidence for reverse transcription of the chloroplast psbA mRNA in two naturally occurring species of Chlamydomonas (raudensis and subcaudata) that is based on the presence of spliced cDNAs in both organisms. The psbA cDNAs, which lack the group II intron of the genomic gene, are nearly full length, and the majority of them--though not all--are in the form of RNA-cDNA hybrids. Moreover, the presence in these species of psbA cDNAs is correlated with the loss of an early group I intron from the same psbA gene. The group II intron that interrupts psbA in C. raudensis and C. subcaudata potentially encodes a protein with a reverse transcriptase domain, and the C. raudensis protein was shown to have reverse transcriptase activity in vitro. These results provide strong evidence for reverse transcription of a physiologically important mRNA (psbA) in two species of Chlamydomonas that have also lost an intron from the same gene, possibly through recombination with the cDNA.

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3' overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding-thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.

Chlamydomonas chloroplasts can use short dispersed repeats and multiple pathways to repair a double-strand break in the genome

Certain group I introns insert into intronless DNA via an endonuclease that creates a double-strand break (DSB). There are two models for intron homing in phage: synthesis-dependent strand annealing (SDSA) and double-strand break repair (DSBR). The Cr.psbA4 intron homes efficiently from a plasmid into the chloroplast psbA gene in Chlamydomonas, but little is known about the mechanism. Analysis of co-transformants selected using a spectinomycin-resistant 16S gene (16S(spec)) provided evidence for both pathways. We also examined the consequences of the donor DNA having only one-sided or no homology with the psbA gene. When there was no homology with the donor DNA, deletions of up to 5 kb involving direct repeats that flank the psbA gene were obtained. Remarkably, repeats as short as 15 bp were used for this repair, which is consistent with the single-strand annealing (SSA) pathway. When the donor had one-sided homology, the DSB in most co-transformants was repaired using two DNAs, the donor and the 16S(spec) plasmid, which, coincidentally, contained a region that is repeated upstream of psbA. DSB repair using two separate DNAs provides further evidence for the SDSA pathway. These data show that the chloroplast can repair a DSB using short dispersed repeats located proximally, distally, or even on separate molecules relative to the DSB. They also provide a rationale for the extensive repertoire of repeated sequences in this genome.

A molecular approach to the characterization of the eukaryotic communities of an extreme acidic environment: methods for DNA extraction and denaturing gradient gel electrophoresis analysis

The diversity of the phytobenthonic community present in six acidophilic microbial mats from Río Tinto (Iberian Pyritic Belt, SW Spain) was analysed by optical microscopy and two molecular techniques, denaturing gradient gel electrophoresis (DGGE) and sequence analysis of 18S rDNA cloned gene fragments. Sixteen DNA isolation protocols as well as two commercial DNA extraction kits were tested and their efficiency compared. Purified DNA extracts were amplified by PCR using universal eukaryotic primers and the PCR products analysed by DGGE. Bead-mill homogenization was found to be superior to the other cell lysis methodologies assayed (sonication or freeze-thawing cycles) as it allowed efficiencies of cell disruption of over 95%. The methods combining bead-mill homogenization in the presence of SDS, treatment with chemical extractants (hexadecylmethylammonium bromide or guanidine isothiocyanate) and phenol extraction resulted in DNA preparations that amplified the same number of bands when analysed by DGGE as the two commercial kits assayed. The phylogenetic affiliations of the DGGE bands were determined by a BLAST search, and nine different species related to the Chlorophyta, Ciliophora, Kinetoplastida, Ascomycota, Streptophyta and Colcochaetales taxonomical groups were identified. Similar levels of diversity were found using cloning procedures. Although not all the species observed under the microscope were detected using molecular techniques, e.g. euglenas, heliozoan, or amoebae, DGGE fingerprints showed rather well the level of diversity present in the samples analysed, with limitations similar to cloning techniques.

Expression, purification, and biochemical characterization of the intron-encoded endonuclease, I-CreII

The ORF of the Cr.psbA4 intron of Chlamydomonas reinhardtii mediates efficient intron homing, and contains an H-N-H and possibly a GIY-YIG motif. The ORF was over-expressed in Escherichia coli without non-native amino acids, but was mostly insoluble. However, co-over-expression of E. coli chaperonins GroEL/GroES solubilized approximately 50% of the protein, which was purified by ion-exchange and heparin-affinity chromatography. Biochemical characterization showed that the protein is a double-strand-specific endonuclease that cleaves fused psbA exon 4-exon 5 DNA, and was named I-CreII. I-CreII has a relatively relaxed divalent metal ion requirement (Mg(2+), Mn(2+), Ca(2+), and Fe(2+) supported cleavage), is insensitive to salt <350 mM, and is stabilized by DNA. Cleavage of target DNA occurs close (4 nt on the top strand) to the intron-insertion site, and leaves 2-nt 3'-OH overhangs, similar to GIY-YIG endonucleases. The boundaries of the recognition sequence span approximately 30 bp, and encompass the cleavage and intron-insertion sites. Cleavage of heterologous psbA DNAs indicates the enzyme can tolerate multiple, but not all, substitutions in the recognition site. This work will facilitate further study of this novel endonuclease, which may also find use in site-specific manipulation of chloroplast DNA.

Chloroplast RNA processing and stability

Primary chloroplast transcripts are processed in a number of ways, including intron splicing, internal cleavage of polycistronic RNAs, and endonucleolytic or exonucleolytic cleavages at the transcript termini. All chloroplast RNAs are also subject to degradation, although a curious feature of many chloroplast mRNAs is their relative longevity. Some of these processes, e.g., psbA splicing and stability of a number of chloroplast mRNAs, are regulated in response to light-dark cycles or nutrient availability. This review highlights recent advances in our understanding of these processes in the model organism Chlamydomonas reinhardtii, focusing on results since the extensive reviews published in 1998 [Herrin DL et al. 1998 (pp. 183-195), Nickelsen Y 1998 (pp. 151-163), Stern DB and Drager RG 1998 (pp. 164-182), in Rochaix JD et al. (eds) The Molecular Biology of Chloroplasts and Mitochondria in Chlamydomonas. Kluwer Academic Publishers, Dordrecht, The Netherlands]. We also allude to studies with other organisms, and to the potential impact of the Chlamydomonas genome project where appropriate.

Mutagenesis of a light-regulated psbA intron reveals the importance of efficient splicing for photosynthetic growth

The chloroplast-encoded psbA gene encodes the D1 polypeptide of the photosystem II reaction center, which is synthesized at high rates in the light. In Chlamydomonas reinhardtii, the psbA gene contains four self-splicing group I introns whose rates of splicing in vivo are increased at least 6-10-fold by light. However, because psbA is an abundant mRNA, and some chloroplast mRNAs appear to be in great excess of what is needed to sustain translation rates, the developmental significance of light-promoted splicing has not been clear. To address this and other questions, potentially destabilizing substitutions were made in several predicted helices of the fourth psbA intron, Cr.psbA4, and their effects on in vitro and in vivo splicing assessed. Two-nucleotide substitutions in P4 and P7 were necessary to substantially reduce splicing of this intron in vivo, although most mutations reduced self-splicing in vitro. The P7-4,5 mutant, whose splicing was completely blocked, showed no photoautotrophic growth and synthesis of a truncated D1 (exons 1-4) polypeptide from the unspliced mRNA. Most informative was the P4'-3,4 mutant, which exhibited a 45% reduction in spliced psbA mRNA, a 28% reduction in synthesis of full-length D1, and an 18% reduction in photoautotrophic growth. These results indicate that psbA mRNA is not in great excess, and that highly efficient splicing of psbA introns, which is afforded by light conditions, is necessary for optimal photosynthetic growth.

Nuclear genes that promote splicing of group I introns in the chloroplast 23S rRNA and psbA genes in Chlamydomonas reinhardtii

Single nucleotide substitutions were made in the core helices P4, P6, and P7, and in the metal-binding GAAA motif in the J4/5 region of the chloroplast group I rRNA intron of Chlamydomonas reinhardtii, Cr.LSU. In vitro assays showed that these substitutions had surprisingly strong effects on Cr.LSU self-splicing; however, splicing of all but the P6 mutations could be at least partially recovered by increasing the Mg2+ concentration. The mutant constructs were transformed into chloroplasts to replace the wild-type intron; however, only the P4 mutants became homoplasmic, indicating that the other mutations were lethal. The splicing-deficient P4125A mutant, which exhibited slow growth and light sensitivity, was used to isolate suppressor strains that showed a substantial restoration of Cr.LSU splicing. Genetic analysis of the 7151, 7120 and 71N1 suppressors indicated that these mutations are in at least two nuclear genes. The 7151 suppressor mutation, which defines the chloroplast-splicing suppressor (css1) gene, had no obviously altered growth phenotype with the wild-type intron, and was dominant in vegetative diploids containing the mutant intron. All three of the suppressor strains also suppressed a mutation in the P4 region of the fourth psbA intron, Cr.psbA4, indicating that these genes play a role in splicing of multiple group I introns in the chloroplast.