Evolution of C4 photosynthesis in the genus Flaveria: how many and which genes does it take to make C4?

Selective pressure exerted by a massive decline in atmospheric CO(2) levels 55 to 40 million years ago promoted the evolution of a novel, highly efficient mode of photosynthetic carbon assimilation known as C(4) photosynthesis. C(4) species have concurrently evolved multiple times in a broad range of plant families, and this multiple and parallel evolution of the complex C(4) trait indicates a common underlying evolutionary mechanism that might be elucidated by comparative analyses of related C(3) and C(4) species. Here, we use mRNA-Seq analysis of five species within the genus Flaveria, ranging from C(3) to C(3)-C(4) intermediate to C(4) species, to quantify the differences in the transcriptomes of closely related plant species with varying degrees of C(4)-associated characteristics. Single gene analysis defines the C(4) cycle enzymes and transporters more precisely and provides new candidates for yet unknown functions as well as identifies C(4) associated pathways. Molecular evidence for a photorespiratory CO(2) pump prior to the establishment of the C(4) cycle-based CO(2) pump is provided. Cluster analysis defines the upper limit of C(4)-related gene expression changes in mature leaves of Flaveria as 3582 alterations.

An mRNA blueprint for C4 photosynthesis derived from comparative transcriptomics of closely related C3 and C4 species.

C(4) photosynthesis involves alterations to the biochemistry, cell biology, and development of leaves. Together, these modifications increase the efficiency of photosynthesis, and despite the apparent complexity of the pathway, it has evolved at least 45 times independently within the angiosperms. To provide insight into the extent to which gene expression is altered between C(3) and C(4) leaves, and to identify candidates associated with the C(4) pathway, we used massively parallel mRNA sequencing of closely related C(3) (Cleome spinosa) and C(4) (Cleome gynandra) species. Gene annotation was facilitated by the phylogenetic proximity of Cleome and Arabidopsis (Arabidopsis thaliana). Up to 603 transcripts differ in abundance between these C(3) and C(4) leaves. These include 17 transcription factors, putative transport proteins, as well as genes that in Arabidopsis are implicated in chloroplast movement and expansion, plasmodesmatal connectivity, and cell wall modification. These are all characteristics known to alter in a C(4) leaf but that previously had remained undefined at the molecular level. We also document large shifts in overall transcription profiles for selected functional classes. Our approach defines the extent to which transcript abundance in these C(3) and C(4) leaves differs, provides a blueprint for the NAD-malic enzyme C(4) pathway operating in a dicotyledon, and furthermore identifies potential regulators. We anticipate that comparative transcriptomics of closely related species will provide deep insight into the evolution of other complex traits.

Lipid accumulation during the establishment of kleptoplasty in Elysia chlorotica.

The establishment of kleptoplasty (retention of "stolen plastids") in the digestive tissue of the sacoglossan Elysia chlorotica Gould was investigated using transmission electron microscopy. Cellular processes occurring during the initial exposure to plastids were observed in laboratory raised animals ranging from 1-14 days post metamorphosis (dpm). These observations revealed an abundance of lipid droplets (LDs) correlating to plastid abundance. Starvation of animals resulted in LD and plastid decay in animals <5 dpm that had not yet achieved permanent kleptoplasty. Animals allowed to feed on algal prey (Vaucheria litorea C. Agardh) for 7 d or greater retained stable plastids resistant to cellular breakdown. Lipid analysis of algal and animal samples supports that these accumulating LDs may be of plastid origin, as the often algal-derived 20∶5 eicosapentaenoic acid was found in high abundance in the animal tissue. Subsequent culturing of animals in dark conditions revealed a reduced ability to establish permanent kleptoplasty in the absence of photosynthetic processes, coupled with increased mortality. Together, these data support an important role of photosynthetic lipid production in establishing and stabilizing this unique animal kleptoplasty.

Arabidopsis tic62 trol mutant lacking thylakoid-bound ferredoxin-NADP+ oxidoreductase shows distinct metabolic phenotype.

Ferredoxin-NADP+ oxidoreductase (FNR), functioning in the last step of the photosynthetic electron transfer chain, exists both as a soluble protein in the chloroplast stroma and tightly attached to chloroplast membranes. Surface plasmon resonance assays showed that the two FNR isoforms, LFNR1 and LFNR2, are bound to the thylakoid membrane via the C-terminal domains of Tic62 and TROL proteins in a pH-dependent manner. The tic62 trol double mutants contained a reduced level of FNR, exclusively found in the soluble stroma. Although the mutant plants showed no visual phenotype or defects in the function of photosystems under any conditions studied, a low ratio of NADPH/NADP+ was detected. Since the CO2 fixation capacity did not differ between the tic62 trol plants and wild-type, it seems that the plants are able to funnel reducing power to most crucial reactions to ensure survival and fitness of the plants. However, the activity of malate dehydrogenase was down-regulated in the mutant plants. Apparently, the plastid metabolism is able to cope with substantial changes in directing the electrons from the light reactions to stromal metabolism and thus only few differences are visible in steady-state metabolite pool sizes of the tic62 trol plants.

In vivo function of Tic22, a protein import component of the intermembrane space of chloroplasts.

Preprotein import into chloroplasts depends on macromolecular machineries in the outer and inner chloroplast envelope membrane (TOC and TIC). It was suggested that both machineries are interconnected by components of the intermembrane space (IMS). That is, amongst others, Tic22, of which two closely related isoforms exist in Arabidopsis thaliana, namely atTic22-III and atTic22-IV. We investigated the function of Tic22 in vivo by analyzing T-DNA insertion lines of the corresponding genes. While the T-DNA insertion in the individual genes caused only slight defects, a double mutant of both isoforms showed retarded growth, a pale phenotype under high-light conditions, a reduced import rate, and a reduction in the photosynthetic performance of the plants. The latter is supported by changes in the metabolite content of mutant plants when compared to wild-type. Thus, our results support the notion that Tic22 is directly involved in chloroplast preprotein import and might point to a particular importance of Tic22 in chloroplast biogenesis at times of high import rates.

Sampling the Arabidopsis transcriptome with massively parallel pyrosequencing.

Massively parallel sequencing of DNA by pyrosequencing technology offers much higher throughput and lower cost than conventional Sanger sequencing. Although extensively used already for sequencing of genomes, relatively few applications of massively parallel pyrosequencing to transcriptome analysis have been reported. To test the ability of this technology to provide unbiased representation of transcripts, we analyzed mRNA from Arabidopsis (Arabidopsis thaliana) seedlings. Two sequencing runs yielded 541,852 expressed sequence tags (ESTs) after quality control. Mapping of the ESTs to the Arabidopsis genome and to The Arabidopsis Information Resource 7.0 cDNA models indicated: (1) massively parallel pyrosequencing detected transcription of 17,449 gene loci providing very deep coverage of the transcriptome. Performing a second sequencing run only increased the number of genes identified by 10%, but increased the overall sequence coverage by 50%. (2) Mapping of the ESTs to their predicted full-length transcripts indicated that all regions of the transcript were well represented regardless of transcript length or expression level. Furthermore, short, medium, and long transcripts were equally represented. (3) Over 16,000 of the ESTs that mapped to the genome were not represented in the existing dbEST database. In some cases, the ESTs provide the first experimental evidence for transcripts derived from predicted genes, and, for at least 60 locations in the genome, pyrosequencing identified likely protein-coding sequences that are not now annotated as genes. Together, the results indicate massively parallel pyrosequencing provides novel information helpful to improve the annotation of the Arabidopsis genome. Furthermore, the unbiased representation of transcripts will be particularly useful for gene discovery and gene expression analysis of nonmodel plants with less complete genomic information. EST sequence accession numbers in GenBank are EH 795234 through EH 995233 and EL 000001 through EL 341852.