Transcriptional analysis of primary ciliary dyskinesia airway cells reveals a dedicated cilia glutathione pathway.

Primary ciliary dyskinesia (PCD) is a genetic condition that results in dysmotile cilia. The repercussions of cilia dysmotility and gene variants on the multiciliated cell remain poorly understood. We used single-cell RNA sequencing, proteomics, and advanced microscopy to compare primary culture epithelial cells from patients with PCD, their heterozygous mothers, healthy individuals, and induced pluripotent stem (iPS) cells generated from a PCD patient. Transcriptomic analysis revealed unique signatures in PCD airway cells compared to their mothers and healthy individuals. Gene expression in heterozygous mothers' cells diverged from both control and PCD cells, marked by increased inflammatory and cellular stress signatures. Primary and iPS-derived PCD multiciliated cells had increased expression of glutathione-S-transferases, GSTA2 and GSTA1, as well as NRF2 target genes, accompanied by elevated levels of reactive oxygen species (ROS). Immunogold labeling in human cilia and proteomic analysis of the ciliated organism, Chlamydomonas reinhardtii, demonstrated that GSTA2 localizes to motile cilia. Loss of human GSTA2 and C. reinhardtii GSTA resulted in slowed cilia motility pointing to local cilia regulatory roles. Our findings identify cellular responses unique to PCD variants and independent of environmental stress and uncover a dedicated ciliary GSTA2 pathway essential for normal motility that may be a therapeutic target.

Ultrastructure expansion microscopy (U-ExM) of mouse and human kidneys for analysis of subcellular structures.

Ultrastructure expansion microscopy (U-ExM) involves the physical magnification of specimens embedded in hydrogels, which allows for super-resolution imaging of subcellular structures using a conventional diffraction-limited microscope. Methods for expansion microscopy exist for several organisms, organs, and cell types, and used to analyze cellular organelles and substructures in nanoscale resolution. Here, we describe a simple step-by-step U-ExM protocol for the expansion, immunostaining, imaging, and analysis of cytoskeletal and organellar structures in kidney tissue. We detail the critical modified steps to optimize isotropic kidney tissue expansion, and preservation of the renal cell structures of interest. We demonstrate the utility of the approach using several markers of renal cell types, centrioles, cilia, the extracellular matrix, and other cytoskeletal elements. Finally, we show that the approach works well on mouse and human kidney samples that were preserved using different fixation and embedding conditions. Overall, this protocol provides a simple and cost-effective approach to analyze both preclinical and clinical renal samples in high detail, using conventional lab supplies and standard widefield or confocal microscopy.

Loss of an extensive ciliary connectome induces proteostasis and cell fate switching in a severe motile ciliopathy.

Motile cilia have essential cellular functions in development, reproduction, and homeostasis. Genetic causes for motile ciliopathies have been identified, but the consequences on cellular functions beyond impaired motility remain unknown. Variants in CCDC39 and CCDC40 cause severe disease not explained by loss of motility. Using human cells with pathological variants in these genes, Chlamydomonas genetics, cryo-electron microscopy, single cell RNA transcriptomics, and proteomics, we identified perturbations in multiple cilia-independent pathways. Absence of the axonemal CCDC39/CCDC40 heterodimer results in loss of a connectome of over 90 proteins. The undocked connectome activates cell quality control pathways, switches multiciliated cell fate, impairs microtubule architecture, and creates a defective periciliary barrier. Both cilia-dependent and independent defects are likely responsible for the disease severity. Our findings provide a foundation for reconsidering the broad cellular impact of pathologic variants in ciliopathies and suggest new directions for therapies.

Mastigoneme structure reveals insights into the O-linked glycosylation code of native hydroxyproline-rich helices.

Hydroxyproline-rich glycoproteins (HRGPs) are a ubiquitous class of protein in the extracellular matrices and cell walls of plants and algae, yet little is known of their native structures or interactions. Here, we used electron cryomicroscopy (cryo-EM) to determine the structure of the hydroxyproline-rich mastigoneme, an extracellular filament isolated from the cilia of the alga Chlamydomonas reinhardtii. The structure demonstrates that mastigonemes are formed from two HRGPs (a filament of MST1 wrapped around a single copy of MST3) that both have hyperglycosylated poly(hydroxyproline) helices. Within the helices, O-linked glycosylation of the hydroxyproline residues and O-galactosylation of interspersed serine residues create a carbohydrate casing. Analysis of the associated glycans reveals how the pattern of hydroxyproline repetition determines the type and extent of glycosylation. MST3 possesses a PKD2-like transmembrane domain that forms a heteromeric polycystin-like cation channel with PKD2 and SIP, explaining how mastigonemes are tethered to ciliary membranes.

Gene dosage of independent dynein arm motor preassembly factors influences cilia assembly in Chlamydomonas reinhardtii.

Motile cilia assembly utilizes over 800 structural and cytoplasmic proteins. Variants in approximately 58 genes cause primary ciliary dyskinesia (PCD) in humans, including the dynein arm (pre)assembly factor (DNAAF) gene DNAAF4. In humans, outer dynein arms (ODAs) and inner dynein arms (IDAs) fail to assemble motile cilia when DNAAF4 function is disrupted. In Chlamydomonas reinhardtii, a ciliated unicellular alga, the DNAAF4 ortholog is called PF23. The pf23-1 mutant assembles short cilia and lacks IDAs, but partially retains ODAs. The cilia of a new null allele (pf23-4) completely lack ODAs and IDAs and are even shorter than cilia from pf23-1. In addition, PF23 plays a role in the cytoplasmic modification of IC138, a protein of the two-headed IDA (I1/f). As most PCD variants in humans are recessive, we sought to test if heterozygosity at two genes affects ciliary function using a second-site non-complementation (SSNC) screening approach. We asked if phenotypes were observed in diploids with pairwise heterozygous combinations of 21 well-characterized ciliary mutant Chlamydomonas strains. Vegetative cultures of single and double heterozygous diploid cells did not show SSNC for motility phenotypes. When protein synthesis is inhibited, wild-type Chlamydomonas cells utilize the pool of cytoplasmic proteins to assemble half-length cilia. In this sensitized assay, 8 double heterozygous diploids with pf23 and other DNAAF mutations show SSNC; they assemble shorter cilia than wild-type. In contrast, double heterozygosity of the other 203 strains showed no effect on ciliary assembly. Immunoblots of diploids heterozygous for pf23 and wdr92 or oda8 show that PF23 is reduced by half in these strains, and that PF23 dosage affects phenotype severity. Reductions in PF23 and another DNAAF in diploids affect the ability to assemble ODAs and IDAs and impedes ciliary assembly. Thus, dosage of multiple DNAAFs is an important factor in cilia assembly and regeneration.

Ultrastructure expansion microscopy (U-ExM) of mouse and human kidneys for analysis of subcellular structures.

Ultrastructure expansion microscopy (U-ExM) involves the physical magnification of specimens embedded in hydrogels, which allows for super-resolution imaging of subcellular structures using a conventional diffraction-limited microscope. Methods for expansion microscopy exist for several organisms, organs, and cell types, and used to analyze cellular organelles and substructures in nanoscale resolution. Here, we describe a simple step-by-step U-ExM protocol for the expansion, immunostaining, imaging, and analysis of cytoskeletal and organellar structures in kidney tissue. We detail the critical modified steps to optimize isotropic kidney tissue expansion, and preservation of the renal cell structures of interest. We demonstrate the utility of the approach using several markers of renal cell types, centrioles, cilia, the extracellular matrix, and other cytoskeletal elements. Finally, we show that the approach works well on mouse and human kidney samples that were preserved using different fixation and storage conditions. Overall, this protocol provides a simple and cost-effective approach to analyze both pre-clinical and clinical renal samples in high detail, using conventional lab supplies and standard widefield or confocal microscopy.

Assembly of FAP93 at the proximal axoneme in Chlamydomonas cilia.

To identify proteins specific to the proximal ciliary axoneme, we used iTRAQ to compare short (~2 µm) and full-length (~11 µm) axonemes of Chlamydomonas. Known compoents of the proximal axoneme such as minor dynein heavy chains and LF5 kinase as well as the ciliary tip proteins FAP256 (CEP104) and EB1 were enriched in short axonemes whereas proteins present along the length of the axoneme were of similar abundance in both samples. The iTRAQ analysis revealed that FAP93, a protein of unknown function, and protein phosphatase 2A (PP2A) are enriched in the short axonemes. Consistently, immunoblots show enrichment of FAP93 and PP2A in short axonemes and immunofluorescence confirms the localization of FAP93 and enrichment of PP2A at the proximal axoneme. Ciliary regeneration reveals that FAP93 assembles continuously but more slowly than other axonemal structures and terminates at 1.03 µm in steady-state axonemes. The length of FAP93 assembly correlates with ciliary length, demonstrating ciliary length-dependent assembly of FAP93. Dikaryon rescue experiments show that FAP93 can assemble independently of IFT transport. In addition, FRAP analysis of GFP-tagged FAP93 demonstrates that FAP93 is stably anchored in axoneme. FAP93 may function as a scaffold for assembly of other specific proteins at the proximal axoneme.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Robust acoustic trapping and perturbation of single-cell microswimmers illuminate three-dimensional swimming and ciliary coordination.

We report a label-free acoustic microfluidic method to confine single, cilia-driven swimming cells in space without limiting their rotational degrees of freedom. Our platform integrates a surface acoustic wave (SAW) actuator and bulk acoustic wave (BAW) trapping array to enable multiplexed analysis with high spatial resolution and trapping forces that are strong enough to hold individual microswimmers. The hybrid BAW/SAW acoustic tweezers employ high-efficiency mode conversion to achieve submicron image resolution while compensating for parasitic system losses to immersion oil in contact with the microfluidic chip. We use the platform to quantify cilia and cell body motion for wildtype biciliate cells, investigating effects of environmental variables like temperature and viscosity on ciliary beating, synchronization, and three-dimensional helical swimming. We confirm and expand upon the existing understanding of these phenomena, for example determining that increasing viscosity promotes asynchronous beating. Motile cilia are subcellular organelles that propel microorganisms or direct fluid and particulate flow. Thus, cilia are critical to cell survival and human health. The unicellular alga Chlamydomonas reinhardtii is widely used to investigate the mechanisms underlying ciliary beating and coordination. However, freely swimming cells are difficult to image with sufficient resolution to capture cilia motion, necessitating that the cell body be held during experiments. Acoustic confinement is a compelling alternative to use of a micropipette, or to magnetic, electrical, and optical trapping that may modify the cells and affect their behavior. Beyond establishing our approach to studying microswimmers, we demonstrate a unique ability to mechanically perturb cells via rapid acoustic positioning.

Structural specializations of the sperm tail.

Sperm motility is crucial to reproductive success in sexually reproducing organisms. Impaired sperm movement causes male infertility, which is increasing globally. Sperm are powered by a microtubule-based molecular machine-the axoneme-but it is unclear how axonemal microtubules are ornamented to support motility in diverse fertilization environments. Here, we present high-resolution structures of native axonemal doublet microtubules (DMTs) from sea urchin and bovine sperm, representing external and internal fertilizers. We identify >60 proteins decorating sperm DMTs; at least 15 are sperm associated and 16 are linked to infertility. By comparing DMTs across species and cell types, we define core microtubule inner proteins (MIPs) and analyze evolution of the tektin bundle. We identify conserved axonemal microtubule-associated proteins (MAPs) with unique tubulin-binding modes. Additionally, we identify a testis-specific serine/threonine kinase that links DMTs to outer dense fibers in mammalian sperm. Our study provides structural foundations for understanding sperm evolution, motility, and dysfunction at a molecular level.

The effect of Dnaaf5 gene dosage on primary ciliary dyskinesia phenotypes.

DNAAF5 is a dynein motor assembly factor associated with the autosomal heterogenic recessive condition of motile cilia, primary ciliary dyskinesia (PCD). The effects of allele heterozygosity on motile cilia function are unknown. We used CRISPR-Cas9 genome editing in mice to recreate a human missense variant identified in patients with mild PCD and a second, frameshift-null deletion in Dnaaf5. Litters with Dnaaf5 heteroallelic variants showed distinct missense and null gene dosage effects. Homozygosity for the null Dnaaf5 alleles was embryonic lethal. Compound heterozygous animals with the missense and null alleles showed severe disease manifesting as hydrocephalus and early lethality. However, animals homozygous for the missense mutation had improved survival, with partially preserved cilia function and motor assembly observed by ultrastructure analysis. Notably, the same variant alleles exhibited divergent cilia function across different multiciliated tissues. Proteomic analysis of isolated airway cilia from mutant mice revealed reduction in some axonemal regulatory and structural proteins not previously reported in DNAAF5 variants. Transcriptional analysis of mouse and human mutant cells showed increased expression of genes coding for axonemal proteins. These findings suggest allele-specific and tissue-specific molecular requirements for cilia motor assembly that may affect disease phenotypes and clinical trajectory in motile ciliopathies.

The effect of Dnaaf5 gene dosage on primary ciliary dyskinesia phenotypes.

DNAAF5 is a dynein motor assembly factor associated with the autosomal heterogenic recessive condition of motile cilia, primary ciliary dyskinesia (PCD). The effects of allele heterozygosity on motile cilia function are unknown. We used CRISPR-Cas9 genome editing in mice to recreate a human missense variant identified in patients with mild PCD and a second, frameshift null deletion in Dnaaf5 . Litters with Dnaaf5 heteroallelic variants showed distinct missense and null gene dosage effects. Homozygosity for the null Dnaaf5 alleles was embryonic lethal. Compound heterozygous animals with the missense and null alleles showed severe disease manifesting as hydrocephalus and early lethality. However, animals homozygous for the missense mutation had improved survival, with partial preserved cilia function and motor assembly observed by ultrastructure analysis. Notably, the same variant alleles exhibited divergent cilia function across different multiciliated tissues. Proteomic analysis of isolated airway cilia from mutant mice revealed reduction in some axonemal regulatory and structural proteins not previously reported in DNAAF5 variants. While transcriptional analysis of mouse and human mutant cells showed increased expression of genes coding for axonemal proteins. Together, these findings suggest allele-specific and tissue-specific molecular requirements for cilia motor assembly that may affect disease phenotypes and clinical trajectory in motile ciliopathies. Brief Summary: A mouse model of human DNAAF5 primary ciliary dyskinesia variants reveals gene dosage effects of mutant alleles and tissue-specific molecular requirements for cilia motor assembly.

A gap-free genome assembly of Chlamydomonas reinhardtii and detection of translocations induced by CRISPR-mediated mutagenesis.

Genomic assemblies of the unicellular green alga Chlamydomonas reinhardtii have provided important resources for researchers. However, assembly errors, large gaps, and unplaced scaffolds as well as strain-specific variants currently impede many types of analysis. By combining PacBio HiFi and Oxford Nanopore long-read technologies, we generated a de novo genome assembly for strain CC-5816, derived from crosses of strains CC-125 and CC-124. Multiple methods of evaluating genome completeness and base-pair error rate suggest that the final telomere-to-telomere assembly is highly accurate. The CC-5816 assembly enabled previously difficult analyses that include characterization of the 17 centromeres, rDNA arrays on three chromosomes, and 56 insertions of organellar DNA into the nuclear genome. Using Nanopore sequencing, we identified sites of cytosine (CpG) methylation, which are enriched at centromeres. We analyzed CRISPR-Cas9 insertional mutants in the PF23 gene. Two of the three alleles produced progeny that displayed patterns of meiotic inviability that suggested the presence of a chromosomal aberration. Mapping Nanopore reads from pf23-2 and pf23-3 onto the CC-5816 genome showed that these two strains each carry a translocation that was initiated at the PF23 gene locus on chromosome 11 and joined with chromosomes 5 or 3, respectively. The translocations were verified by demonstrating linkage between loci on the two translocated chromosomes in meiotic progeny. The three pf23 alleles display the expected short-cilia phenotype, and immunoblotting showed that pf23-2 lacks the PF23 protein. Our CC-5816 genome assembly will undoubtedly provide an important tool for the Chlamydomonas research community.

Ciliary central apparatus structure reveals mechanisms of microtubule patterning.

A pair of extensively modified microtubules form the central apparatus (CA) of the axoneme of most motile cilia, where they regulate ciliary motility. The external surfaces of both CA microtubules are patterned asymmetrically with large protein complexes that repeat every 16 or 32 nm. The composition of these projections and the mechanisms that establish asymmetry and longitudinal periodicity are unknown. Here, by determining cryo-EM structures of the CA microtubules, we identify 48 different CA-associated proteins, which in turn reveal mechanisms for asymmetric and periodic protein binding to microtubules. We identify arc-MIPs, a novel class of microtubule inner protein, that bind laterally across protofilaments and remodel tubulin structure and lattice contacts. The binding mechanisms utilized by CA proteins may be generalizable to other microtubule-associated proteins. These structures establish a foundation to elucidate the contributions of individual CA proteins to ciliary motility and ciliopathies.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Structures of radial spokes and associated complexes important for ciliary motility.

In motile cilia, a mechanoregulatory network is responsible for converting the action of thousands of dynein motors bound to doublet microtubules into a single propulsive waveform. Here, we use two complementary cryo-EM strategies to determine structures of the major mechanoregulators that bind ciliary doublet microtubules in Chlamydomonas reinhardtii. We determine structures of isolated radial spoke RS1 and the microtubule-bound RS1, RS2 and the nexin-dynein regulatory complex (N-DRC). From these structures, we identify and build atomic models for 30 proteins, including 23 radial-spoke subunits. We reveal how mechanoregulatory complexes dock to doublet microtubules with regular 96-nm periodicity and communicate with one another. Additionally, we observe a direct and dynamically coupled association between RS2 and the dynein motor inner dynein arm subform c (IDAc), providing a molecular basis for the control of motor activity by mechanical signals. These structures advance our understanding of the role of mechanoregulation in defining the ciliary waveform.

Basal Feet: Walking to the Discovery of a Novel Hybrid Cilium.

Cilia are important cell structures found on nearly all cells. In this issue of Developmental Cell, Mennella and colleagues investigate the molecular architecture of basal foot proteins in cells with primary or motile cilia and discover a hybrid cilium with a unique assembly that regulates polarity in multiciliated cells.

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.

Alternative Splicing During the Chlamydomonasreinhardtii 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.