Chlamydomonas reinhardtii: biological rationale for genomics

Mutation of CFAP57, a protein required for the asymmetric targeting of a subset of inner dynein arms in Chlamydomonas, causes primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is characterized by chronic airway disease, reduced fertility, and randomization of the left/right body axis. It is caused by defects of motile cilia and sperm flagella. We screened a cohort of affected individuals that lack an obvious axonemal defect for pathogenic variants using whole exome capture, next generation sequencing, and bioinformatic analysis assuming an autosomal recessive trait. We identified one subject with an apparently homozygous nonsense variant [(c.1762C>T), p.(Arg588*)] in the uncharacterized CFAP57 gene. Interestingly, the variant results in the skipping of exon 11 (58 amino acids), which may be due to disruption of an exonic splicing enhancer. In normal human nasal epithelial cells, CFAP57 localizes throughout the ciliary axoneme. Nasal cells from the PCD patient express a shorter, mutant version of CFAP57 and the protein is not incorporated into the axoneme. The missing 58 amino acids include portions of WD repeats that may be important for loading onto the intraflagellar transport (IFT) complexes for transport or docking onto the axoneme. A reduced beat frequency and an alteration in ciliary waveform was observed. Knockdown of CFAP57 in human tracheobronchial epithelial cells (hTECs) recapitulates these findings. Phylogenetic analysis showed that CFAP57 is highly conserved in organisms that assemble motile cilia. CFAP57 is allelic with the BOP2/IDA8/FAP57 gene identified previously in Chlamydomonas reinhardtii. Two independent, insertional fap57 Chlamydomonas mutant strains show reduced swimming velocity and altered waveforms. Tandem mass tag (TMT) mass spectroscopy shows that FAP57 is missing, and the "g" inner dyneins (DHC7 and DHC3) and the "d" inner dynein (DHC2) are reduced, but the FAP57 paralog FBB7 is increased. Together, our data identify a homozygous variant in CFAP57 that causes PCD that is likely due to a defect in the inner dynein arm assembly process.

C2 domain protein MIN1 promotes eyespot organization in Chlamydomonas reinhardtii

Assembly and asymmetric localization of the photosensory eyespot in the biflagellate, unicellular green alga Chlamydomonas reinhardtii requires coordinated organization of photoreceptors in the plasma membrane and pigment granule/thylakoid membrane layers in the chloroplast. min1 (mini-eyed) mutant cells contain abnormally small, disorganized eyespots in which the chloroplast envelope and plasma membrane are no longer apposed. The MIN1 gene, identified here by phenotypic rescue, encodes a protein with an N-terminal C2 domain and a C-terminal LysM domain separated by a transmembrane sequence. This novel domain architecture led to the hypothesis that MIN1 is in the plasma membrane or the chloroplast envelope, where membrane association of the C2 domain promotes proper eyespot organization. Mutation of conserved C2 domain loop residues disrupted association of the MIN1 C2 domain with the chloroplast envelope in moss cells but did not abolish eyespot assembly in Chlamydomonas. In min1 null cells, channelrhodopsin-1 (ChR1) photoreceptor levels were reduced, indicating a role for MIN1 in ChR1 expression and/or stability. However, ChR1 localization was only minimally disturbed during photoautotrophic growth of min1 cells, conditions under which the pigment granule layers are disorganized. The data are consistent with the hypothesis that neither MIN1 nor proper organization of the plastidic components of the eyespot is essential for localization of ChR1.

Mutation of murine adenylate kinase 7 underlies a primary ciliary dyskinesia phenotype

Primary ciliary dyskinesia (PCD) is a genetically and phenotypically heterogeneous disorder, characterized by progressive development of bronchiectasis, inflammation, and features characteristic of chronic obstructive pulmonary disease. We report here that a murine mutation of the evolutionarily conserved adenylate kinase 7 (Ak7) gene results in animals presenting with pathological signs characteristic of PCD, including ultrastructural ciliary defects and decreased ciliary beat frequency in respiratory epithelium. The mutation is associated with hydrocephalus, abnormal spermatogenesis, mucus accumulation in paranasal passages, and a dramatic respiratory pathology upon allergen challenge. Ak7 appears to be a marker for cilia with (9 + 2) microtubular organization. This is suggested by its tissue specificity of expression and also the stringent conservation of Ak7 ortholog structure only in protozoans and metazoans possessing motile (9 + 2) cilia. Collectively, our results indicate an ancestral and crucial role of Ak7 in maintaining ciliary structure and function, and suggest that mutations of the human ortholog may underlie a subset of genetically uncharacterized PCD cases.

Generation of an oligonucleotide array for analysis of gene expression in Chlamydomonas reinhardtii

The availability of genome sequences makes it possible to develop microarrays that can be used for profiling gene expression over developmental time, as organisms respond to environmental challenges, and for comparison between wild-type and mutant strains under various conditions. The desired characteristics of microarrays (intense signals, hybridization specificity and extensive coverage of the transcriptome) were not fully met by the previous Chlamydomonas reinhardtii microarray: probes derived from cDNA sequences (approximately 300 bp) were prone to some nonspecific cross-hybridization and coverage of the transcriptome was only approximately 20%. The near completion of the C. reinhardtii nuclear genome sequence and the availability of extensive cDNA information have made it feasible to improve upon these aspects. After developing a protocol for selecting a high-quality unigene set representing all known expressed sequences, oligonucleotides were designed and a microarray with approximately 10,000 unique array elements (approximately 70 bp) covering 87% of the known transcriptome was developed. This microarray will enable researchers to generate a global view of gene expression in C. reinhardtii. Furthermore, the detailed description of the protocol for selecting a unigene set and the design of oligonucleotides may be of interest for laboratories interested in developing microarrays for organisms whose genome sequences are not yet completed (but are nearing completion).

3D structure of eukaryotic flagella in a quiescent state revealed by cryo-electron tomography

We have used cryo-electron tomography to investigate the 3D structure and macromolecular organization of intact, frozen-hydrated sea urchin sperm flagella in a quiescent state. The tomographic reconstructions provide information at a resolution better than 6 nm about the in situ arrangements of macromolecules that are key for flagellar motility. We have visualized the heptameric rings of the motor domains in the outer dynein arm complex and determined that they lie parallel to the plane that contains the axes of neighboring flagellar microtubules. Both the material associated with the central pair of microtubules and the radial spokes display a plane of symmetry that helps to explain the planar beat pattern of these flagella. Cryo-electron tomography has proven to be a powerful technique for helping us understand the relationships between flagellar structure and function and the design of macromolecular machines in situ.

Comparative genomics of two closely related unicellular thermo-acidophilic red algae, Galdieria sulphuraria and Cyanidioschyzon merolae, reveals the molecular basis of the metabolic flexibility of Galdieria sulphuraria and significant differences in carbohydrate metabolism of both algae

Unicellular algae serve as models for the study and discovery of metabolic pathways, for the functional dissection of cell biological processes such as organellar division and cell motility, and for the identification of novel genes and gene functions. The recent completion of several algal genome sequences and expressed sequence tag collections and the establishment of nuclear and organellar transformation methods has opened the way for functional genomics approaches using algal model systems. The thermo-acidophilic unicellular red alga Galdieria sulphuraria represents a particularly interesting species for a genomics approach owing to its extraordinary metabolic versatility such as heterotrophic and mixotrophic growth on more than 50 different carbon sources and its adaptation to hot acidic environments. However, the ab initio prediction of genes required for unknown metabolic pathways from genome sequences is not trivial. A compelling strategy for gene identification is the comparison of similarly sized genomes of related organisms with different physiologies. Using this approach, candidate genes were identified that are critical to the metabolic versatility of Galdieria. Expressed sequence tags and high-throughput genomic sequence reads covering >70% of the G. sulphuraria genome were compared to the genome of the unicellular, obligate photoautotrophic red alga Cyanidioschyzon merolae. More than 30% of the Galdieria sequences did not relate to any of the Cyanidioschyzon genes. A closer inspection of these sequences revealed a large number of membrane transporters and enzymes of carbohydrate metabolism that are unique to Galdieria. Based on these data, it is proposed that genes involved in the uptake of reduced carbon compounds and enzymes involved in their metabolism are crucial to the metabolic flexibility of G. sulphuraria.

Cilia, primary ciliary dyskinesia and molecular genetics

Primary ciliary dyskinesia (PCD) is a phenotypically and genetically heterogeneous condition in which three genetic mutations have already been identified. The primary defect is in the ultrastructure or function of cilia, highly complex organelles that are structurally related to the flagella of sperm and protozoa. The clinical features of PCD include recurrent sinopulmonary infections, subfertility and laterality defects; the latter due to ciliary dysfunction at the embryological node. Completion of the human genome sequence has accelerated the identification and characterisation of disease genes, and the current molecular strategy in PCD includes candidate gene analysis, positional cloning, model organism analysis and proteomic analysis. The identification of these genes will provide new insights into the molecular mechanisms involved in the assembly and function of cilia and the pathway that determines left-right axis in man. This may also allow the development of new methods for diagnosis, prevention and treatment of PCD.

Mutations in α-tubulin promote basal body maturation and flagellar assembly in the absence of δ-tubulin

We have isolated suppressors of the deletion allele of δ-tubulin, uni3-1, in the biflagellate green alga Chlamydomonas reinhardtii. The deletion of δ-tubulin produces cells that assemble zero, one or two flagella and have basal bodies composed primarily of doublet rather than triplet microtubules. Flagellar number is completely restored in the suppressed strains. Most of the uni3-1 suppressors map to the TUA2 locus, which encodes α2-tubulin. Twelve independent tua2 mutations were sequenced. Amino acids D205 or A208, which are nearly invariant residues in α-tubulin, were altered. The tua2 mutations on their own have a second phenotype - they make the cells colchicine supersensitive. Colchicine supersensitivity itself is not needed for suppression and colchicine cannot phenocopy the suppression. The suppressors partially restore the assembly of triplet microtubules. These results suggest that the δ-tubulin plays two roles: it is needed for extension or stability of the triplet microtubule and also for early maturation of basal bodies. We suggest that the mutant α-tubulin promotes the early maturation of the basal body in the absence of δ-tubulin, perhaps through interactions with other partners, and this allows assembly of the flagella.

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.

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

Chlamydomonas reinhardtii is a unicellular green alga that has been used as a model organism for the study of flagella and basal bodies as well as photosynthesis. This report analyzes finished genomic DNA sequence for 0.5% of the nuclear genome. We have used three gene prediction programs as well as EST and protein homology data to estimate the total number of genes in Chlamydomonas to be between 12,000 and 16,400. Chlamydomonas appears to have many more genes than any other unicellular organism sequenced to date. Twenty-seven percent of the predicted genes have significant identity to both ESTs and to known proteins in other organisms, 32% of the predicted genes have significant identity to ESTs alone, and 14% have significant similarity to known proteins in other organisms. For gene prediction in Chlamydomonas, GreenGenie appeared to have the highest sensitivity and specificity at the exon level, scoring 71% and 82%. respectively. Two new alternative splicing events were predicted by aligning Chlamydomonas ESTs to the genomic sequence. Finally recombination differs between the two sequenced contigs. The 350-Kb of the Linkage group III contig is devoid of recombination, while the Linkage group I contig is 30 map units long over 33-kb.

Molecular map of the Chlamydomonas reinhardtii nuclear genome

We have prepared a molecular map of the Chlamydomonas reinhardtii genome anchored to the genetic map. The map consists of 264 markers, including sequence-tagged sites (STS), scored by use of PCR and agarose gel electrophoresis, and restriction fragment length polymorphism markers, scored by use of Southern blot hybridization. All molecular markers tested map to one of the 17 known linkage groups of C. reinhardtii. The map covers approximately 1,000 centimorgans (cM). Any position on the C. reinhardtii genetic map is, on average, within 2 cM of a mapped molecular marker. This molecular map, in combination with the ongoing mapping of bacterial artificial chromosome (BAC) clones and the forthcoming sequence of the C. reinhardtii nuclear genome, should greatly facilitate isolation of genes of interest by using positional cloning methods. In addition, the presence of easily assayed STS markers on each arm of each linkage group should be very useful in mapping new mutations in preparation for positional cloning.

In vivo localization of centrin in the green alga Chlamydomonas reinhardtii

The green alga Chlamydomonas reinhardtii has been used as a model system to study flagellar assembly, centriole assembly, and cell cycle events. These processes are dynamic. Therefore, protein targeting and protein-protein interactions should be evaluated in vivo. To be able to study dynamic processes in C. reinhardtii in vivo, we have explored the use of the green fluorescent protein (GFP). A construct containing a fusion of centrin and GFP was incorporated into the genome as a single copy. The selected clone shows expression in 25-50% of the cells. Centrin-GFP was targeted in vivo to the nuclear basal body connectors and the distal connecting fibers. At the electron microscopic level, it was also localized to the flagellar transitional regions. EM data of transformants indicate that there are some abnormalities in the centrin-containing structures. The transitional region consists of only the transverse septum or has lesions in the H-piece. The distal connecting fibers are thinner and their characteristic crossbands seem to be incomplete. Deflagellation is not affected since more than 95% of the cells deflagellate. Also basal body segregation is not affected since cells with an abnormal flagellar number were not detected. Functional studies of the centrin-GFP fusion show the characteristic calcium-induced mobility shift in SDS-PAGE. Immunofluorescence revealed that during cell division, centrin-GFP remains associated with the basal bodies. In vivo localization of the fusion protein during cell division shows that in metaphase centrin-GFP appears as two opposing spots located close to the spindle poles. The distance between the spots increases as the cells progress through anaphase and then decreases during telophase. GFP is a useful tool to study dynamic processes in the cytoskeleton of C. reinhardtii.