Genetic approach to regulated exocytosis using functional complementation in Paramecium: identification of the ND7 gene required for membrane fusion

HCN channels are essential for the escape response of Paramecium

When mechanical stimulation was applied to free swimming Paramecium, forward swimming velocity transiently increased due to activation of the posterior mechanosensory channels. The behavior response, known as "escape response," requires membrane hyperpolarization and the activation of K-channel type adenylate cyclases. Our hypothesis is that this escape response also involves activation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. HCN channels are activated by hyperpolarization and are modulated by cyclic nucleotides such as cAMP and cGMP. They play a critical role in many excitable cells in higher animals. If HCN channels act in Paramecium, this should help to enhance and prolong hyperpolarization, thereby increasing the swimming speed of Paramecium. This study used RNAi to examine the role of the HCN channel 1 in the escape responses by generating hcn1-gene knockdown cells (hcn1-KD). These cells showed reduced mechanically-stimulated escape responses and a lack of cGMP-dependent increases in swimming speed. Electrophysiological experiments demonstrated reduced hyperpolarization upon injection of large negative currents in hcn1-KD cells. This is consistent with a decrease in HCN1 channel activity and changes in the escape response. These findings suggest that HCN1 channels are K+ channels that regulate the escape response of Paramecium by amplifying the hyperpolarizations elicited by posterior mechanical stimulation.

Anoctamin-like protein 1 regulates repolarization in Paramecium behavioral responses

Paramecium exhibits responsive behavior to environmental changes, moving either closer to or further away from stimuli. Electrophysiological experiments have revealed that these behavioral responses are controlled by membrane potentials. Anoctamin, a Ca2+-activated Cl- channel, is involved in the regulation of membrane potential in mammals. However, it remains uncertain whether Cl- channels like anoctamin regulate Paramecium behavior. Herein, replacement of external Cl- ions with acetate ion and application of Cl- channel blocker niflumic acid (NFA, 0.1 μM) increased spontaneous avoiding reactions (sARs). Hence, we hypothesized that anoctamin is involved in the stabilization of membrane potential fluctuation. Paramecium cells in which the anoctamin-like protein 1 gene was knocked down displayed frequent sARs in the culture medium without external stimulation. Treatment of anoctamin-like protein 1-knockdown cells with the Ca2+ chelator BAPTA or Ca-channel blocker nicardipine reversed the increase in sARs. Electrophysiological experiments revealed extension of membrane depolarization when positive currents were applied to anoctamin-like protein 1-knockdown cells. We concluded that anoctamin-like protein 1 works as a Cl-channel and stabilizes the membrane potential oscillation, reducing sARs.

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

Paramecium Genetics, Genomics, and Evolution

The ciliate genus Paramecium served as one of the first model systems in microbial eukaryotic genetics, contributing much to the early understanding of phenomena as diverse as genome rearrangement, cryptic speciation, cytoplasmic inheritance, and endosymbiosis, as well as more recently to the evolution of mating types, introns, and roles of small RNAs in DNA processing. Substantial progress has recently been made in the area of comparative and population genomics. Paramecium species combine some of the lowest known mutation rates with some of the largest known effective populations, along with likely very high recombination rates, thereby harboring a population-genetic environment that promotes an exceptionally efficient capacity for selection. As a consequence, the genomes are extraordinarily streamlined, with very small intergenic regions combined with small numbers of tiny introns. The subject of the bulk of Paramecium research, the ancient Paramecium aurelia species complex, is descended from two whole-genome duplication events that retain high degrees of synteny, thereby providing an exceptional platform for studying the fates of duplicate genes. Despite having a common ancestor dating to several hundred million years ago, the known descendant species are morphologically indistinguishable, raising significant questions about the common view that gene duplications lead to the origins of evolutionary novelties.

Chromatin remodeling is required for sRNA-guided DNA elimination in Paramecium

Small RNAs mediate the silencing of transposable elements and other genomic loci, increasing nucleosome density and preventing undesirable gene expression. The unicellular ciliate Paramecium is a model to study dynamic genome organization in eukaryotic cells, given its unique feature of nuclear dimorphism. Here, the formation of the somatic macronucleus during sexual reproduction requires eliminating thousands of transposon remnants (IESs) and transposable elements scattered throughout the germline micronuclear genome. The elimination process is guided by Piwi-associated small RNAs and leads to precise cleavage at IES boundaries. Here we show that IES recognition and precise excision are facilitated by recruiting ISWI1, a Paramecium homolog of the chromatin remodeler ISWI. ISWI1 knockdown substantially inhibits DNA elimination, quantitatively similar to development-specific sRNA gene knockdowns but with much greater aberrant IES excision at alternative boundaries. We also identify key development-specific sRNA biogenesis and transport proteins, Ptiwi01 and Ptiwi09, as ISWI1 cofactors in our co-immunoprecipitation studies. Nucleosome profiling indicates that increased nucleosome density correlates with the requirement for ISWI1 and other proteins necessary for IES excision. We propose that chromatin remodeling together with small RNAs is essential for efficient and precise DNA elimination in Paramecium.

The evolutionary conserved proteins CEP90, FOPNL, and OFD1 recruit centriolar distal appendage proteins to initiate their assembly

In metazoa, cilia assembly is a cellular process that starts with centriole to basal body maturation, migration to the cell surface, and docking to the plasma membrane. Basal body docking involves the interaction of both the distal end of the basal body and the transition fibers/distal appendages, with the plasma membrane. Mutations in numerous genes involved in basal body docking and transition zone assembly are associated with the most severe ciliopathies, highlighting the importance of these events in cilium biogenesis. In this context, the ciliate Paramecium has been widely used as a model system to study basal body and cilia assembly. However, despite the evolutionary conservation of cilia assembly events across phyla, whether the same molecular players are functionally conserved, is not fully known. Here, we demonstrated that CEP90, FOPNL, and OFD1 are evolutionary conserved proteins crucial for ciliogenesis. Using ultrastructure expansion microscopy, we unveiled that these proteins localize at the distal end of both centrioles/basal bodies in Paramecium and mammalian cells. Moreover, we found that these proteins are recruited early during centriole duplication on the external surface of the procentriole. Functional analysis performed both in Paramecium and mammalian cells demonstrate the requirement of these proteins for distal appendage assembly and basal body docking. Finally, we show that mammalian centrioles require another component, Moonraker (MNR), to recruit OFD1, FOPNL, and CEP90, which will then recruit the distal appendage proteins CEP83, CEP89, and CEP164. Altogether, we propose that this OFD1, FOPNL, and CEP90 functional module is required to determine in mammalian cells the future position of distal appendage proteins.

CEP90, FOPNL and OFD1 form an evolutionary conserved module which promotes the assembly of centriolar distal appendages. This study uses ultrastructure expansion microscopy to reveal the recruitment of this module on early-born procentrioles to in turn recruit centriolar distal appendage proteins, proposing that this dictates the future location of distal appendages.

A novel membrane complex is required for docking and regulated exocytosis of lysosome-related organelles in Tetrahymena thermophila

In the ciliate Tetrahymena thermophila, lysosome-related organelles called mucocysts accumulate at the cell periphery where they secrete their contents in response to extracellular events, a phenomenon called regulated exocytosis. The molecular bases underlying regulated exocytosis have been extensively described in animals but it is not clear whether similar mechanisms exist in ciliates or their sister lineage, the Apicomplexan parasites, which together belong to the ecologically and medically important superphylum Alveolata. Beginning with a T. thermophila mutant in mucocyst exocytosis, we used a forward genetic approach to uncover MDL1 (Mucocyst Discharge with a LamG domain), a novel gene that is essential for regulated exocytosis of mucocysts. Mdl1p is a 40 kDa membrane glycoprotein that localizes to mucocysts, and specifically to a tip domain that contacts the plasma membrane when the mucocyst is docked. This sub-localization of Mdl1p, which occurs prior to docking, underscores a functional asymmetry in mucocysts that is strikingly similar to that of highly polarized secretory organelles in other Alveolates. A mis-sense mutation in the LamG domain results in mucocysts that dock but only undergo inefficient exocytosis. In contrast, complete knockout of MDL1 largely prevents mucocyst docking itself. Mdl1p is physically associated with 9 other proteins, all of them novel and largely restricted to Alveolates, and sedimentation analysis supports the idea that they form a large complex. Analysis of three other members of this putative complex, called MDD (for Mucocyst Docking and Discharge), shows that they also localize to mucocysts. Negative staining of purified MDD complexes revealed distinct particles with a central channel. Our results uncover a novel macromolecular complex whose subunits are conserved within alveolates but not in other lineages, that is essential for regulated exocytosis in T. thermophila.

An Alveolata secretory machinery adapted to parasite host cell invasion

Apicomplexa are unicellular eukaryotes and obligate intracellular parasites, including Plasmodium (the causative agent of malaria) and Toxoplasma (one of the most widespread zoonotic pathogens). Rhoptries, one of their specialized secretory organelles, undergo regulated exocytosis during invasion¹. Rhoptry proteins are injected directly into the host cell to support invasion and subversion of host immune function². The mechanism by which they are discharged is unclear and appears distinct from those in bacteria, yeast, animals and plants. Here, we show that rhoptry secretion in Apicomplexa shares structural and genetic elements with the exocytic machinery of ciliates, their free-living relatives. Rhoptry exocytosis depends on intramembranous particles in the shape of a rosette embedded into the plasma membrane of the parasite apex. Formation of this rosette requires multiple non-discharge (Nd) proteins conserved and restricted to Ciliata, Dinoflagellata and Apicomplexa that together constitute the superphylum Alveolata. We identified Nd6 at the site of exocytosis in association with an apical vesicle. Sandwiched between the rosette and the tip of the rhoptry, this vesicle appears as a central element of the rhoptry secretion machine. Our results describe a conserved secretion system that was adapted to provide defence for free-living unicellular eukaryotes and host cell injection in intracellular parasites.

A Functional Aqp1 Gene Product Localizes on The Contractile Vacuole Complex in Paramecium multimicronucleatum

In a ciliate Paramecium, the presence of water channels on the membrane of contractile vacuole has long been predicted by both morphological and physiological data, however, to date either the biochemical or the molecular biological data have not been provided. In the present study, to examine the presence of aquaporin in Paramecium, we carried out RT-PCR with degenerated primers designed based on the ParameciumDB, and an aquaporin cDNA (aquaporin 1, aqp1) with a full-length ORF encoding 251 amino acids was obtained from Paramecium multimicronucleatum by using RACE. The deduced amino acid sequence of AQP1 had NPA-NPG motifs, and the prediction of protein secondary structure by CNR5000 and hydropathy plot showed the presence of six putative transmembrane domains and five connecting loops. Phylogenetic analysis results showed that the amino acid sequence of AQP1 was close to that of the Super-aquaporin group. The AQP1-GFP fusion protein clearly demonstrated the subcellular localization of AQP1 on the contractile vacuole complex, except for the decorated spongiome membrane. The functional analyses of aqp1 were done by RNA interference-based gene silencing, using an established feeding method. The aqp1 was found to be crucial for the total fluid output of the cell, the function of contractile vacuole membranes.

Intraflagellar Transport 80 Is Required for Cilia Construction and Maintenance in Paramecium tetraurelia

Intraflagellar transport (IFT) represents a bidirectional dynamic process that carries cargo essential for cilia building and the maintenance of ciliary function, which is important for the locomotion of single cells, intracellular and intercellular signalling transduction. Accumulated evidence has revealed that defects in IFT cause several clinical disorders. Here, we determined the role of IFT80, an IFT-B protein that is mutated in Jeune asphyxiating thoracic dystrophy. Using the RNAi method in the ciliate Paramecium as model, we found that loss of IFT80 prevents cilia biogenesis and causes strong cell lethality. A specific antibody against IFT80 was also prepared in our study, which labelled IFT80 in cilia of Paramecium. GFP fusion experiments were performed to illustrate the dynamic movement of IFT-A and IFT-B proteins in cilia of Paramecium; then, we found that the depletion of IFT80 in cells prevents IFT-A and IFT-B proteins from entering the cilia. Our results showed the distribution change of other IFT proteins in cells that were depleted of IFT80, and we discuss the possible roles of IFT80 in Paramecium.

Roles of Adenylate Cyclases in Ciliary Responses of Paramecium to Mechanical Stimulation

Paramecium shows rapid forward swimming due to increased beat frequency of cilia in normal (forward swimming) direction in response to various kinds of stimuli applied to the cell surface that cause K⁺ -outflow accompanied by a membrane hyperpolarization. Some adenylate cyclases are known to be functional K⁺ channels in the membrane. Using gene-specific knockdown methods, we examined nine paralogues of adenylate cyclases in P. tetraurelia to ascertain whether and how they are involved in the mechanical stimulus-induced hyperpolarization-coupled acceleration of forward swimming. Results demonstrated that knockdown of the adenylate cyclase 1 (ac1)-gene and 2 (ac2)-gene inhibited the acceleration of forward swimming in response to mechanical stimulation of the cell, whereas that spared the acceleration response to external application of 8-Br-cAMP and dilution of extracellular [K⁺ ] induced hyperpolarization. Electrophysiological examination of the knockdown cells revealed that the hyperpolarization-activated inward K⁺ current is smaller than that of a normal cell. Our results suggest that AC1 and AC2 are involved in the mechanical stimulus-induced acceleration of ciliary beat in Paramecium.

Functional diversification of Paramecium Ku80 paralogs safeguards genome integrity during precise programmed DNA elimination

Gene duplication and diversification drive the emergence of novel functions during evolution. Because of whole genome duplications, ciliates from the Paramecium aurelia group constitute a remarkable system to study the evolutionary fate of duplicated genes. Paramecium species harbor two types of nuclei: a germline micronucleus (MIC) and a somatic macronucleus (MAC) that forms from the MIC at each sexual cycle. During MAC development, ~45,000 germline Internal Eliminated Sequences (IES) are excised precisely from the genome through a 'cut-and-close' mechanism. Here, we have studied the P. tetraurelia paralogs of KU80, which encode a key DNA double-strand break repair factor involved in non-homologous end joining. The three KU80 genes have different transcription patterns, KU80a and KU80b being constitutively expressed, while KU80c is specifically induced during MAC development. Immunofluorescence microscopy and high-throughput DNA sequencing revealed that Ku80c stably anchors the PiggyMac (Pgm) endonuclease in the developing MAC and is essential for IES excision genome-wide, providing a molecular explanation for the previously reported Ku-dependent licensing of DNA cleavage at IES ends. Expressing Ku80a under KU80c transcription signals failed to complement a depletion of endogenous Ku80c, indicating that the two paralogous proteins have distinct properties. Domain-swap experiments identified the α/β domain of Ku80c as the major determinant for its specialized function, while its C-terminal part is required for excision of only a small subset of IESs located in IES-dense regions. We conclude that Ku80c has acquired the ability to license Pgm-dependent DNA cleavage, securing precise DNA elimination during programmed rearrangements. The present study thus provides novel evidence for functional diversification of genes issued from a whole-genome duplication.

MKS-NPHP module proteins control ciliary shedding at the transition zone

Ciliary shedding occurs from unicellular organisms to metazoans. Although required during the cell cycle and during neurogenesis, the process remains poorly understood. In all cellular models, this phenomenon occurs distal to the transition zone (TZ), suggesting conserved molecular mechanisms. The TZ module proteins (Meckel Gruber syndrome [MKS]/Nephronophtysis [NPHP]/Centrosomal protein of 290 kDa [CEP290]/Retinitis pigmentosa GTPase regulator-Interacting Protein 1-Like Protein [RPGRIP1L]) are known to cooperate to establish TZ formation and function. To determine whether they control deciliation, we studied the function of 5 of them (Transmembrane protein 107 [TMEM107], Transmembrane protein 216 [TMEM216], CEP290, RPGRIP1L, and NPHP4) in Paramecium. All proteins are recruited to the TZ of growing cilia and localize with 9-fold symmetry at the level of the most distal part of the TZ. We demonstrate that depletion of the MKS2/TMEM216 and TMEM107 proteins induces constant deciliation of some cilia, while depletion of either NPHP4, CEP290, or RPGRIP1L prevents Ca2+/EtOH deciliation. Our results constitute the first evidence for a role of conserved TZ proteins in deciliation and open new directions for understanding motile cilia physiology.

Functional analysis and subcellular localisation of the conserved transition zone proteins in the ciliate Paramecium tetraurelia demonstrates their involvement in the ciliary shedding process, opening new avenues fir understanding the molecular mechanism of deciliation.

TTC12 Loss-of-Function Mutations Cause Primary Ciliary Dyskinesia and Unveil Distinct Dynein Assembly Mechanisms in Motile Cilia Versus Flagella

Cilia and flagella are evolutionarily conserved organelles whose motility relies on the outer and inner dynein arm complexes (ODAs and IDAs). Defects in ODAs and IDAs result in primary ciliary dyskinesia (PCD), a disease characterized by recurrent airway infections and male infertility. PCD mutations in assembly factors have been shown to cause a combined ODA-IDA defect, affecting both cilia and flagella. We identified four loss-of-function mutations in TTC12, which encodes a cytoplasmic protein, in four independent families in which affected individuals displayed a peculiar PCD phenotype characterized by the absence of ODAs and IDAs in sperm flagella, contrasting with the absence of only IDAs in respiratory cilia. Analyses of both primary cells from individuals carrying TTC12 mutations and human differentiated airway cells invalidated for TTC12 by a CRISPR-Cas9 approach revealed an IDA defect restricted to a subset of single-headed IDAs that are different in flagella and cilia, whereas TTC12 depletion in the ciliate Paramecium tetraurelia recapitulated the sperm phenotype. Overall, our study, which identifies TTC12 as a gene involved in PCD, unveils distinct dynein assembly mechanisms in human motile cilia versus flagella.

A mating-type mutagenesis screen identifies a zinc-finger protein required for specific DNA excision events in Paramecium

In the ciliate Paramecium tetraurelia, functional genes are reconstituted during development of the somatic macronucleus through the precise excision of ∼45 000 single-copy Internal Eliminated Sequences (IESs), thought to be the degenerate remnants of ancient transposon insertions. Like introns, IESs are marked only by a weak consensus at their ends. How such a diverse set of sequences is faithfully recognized and precisely excised remains unclear: specialized small RNAs have been implicated, but in their absence up to ∼60% of IESs are still correctly excised. To get further insight, we designed a mutagenesis screen based on the hypersensitivity of a specific excision event in the mtA gene, which determines mating types. Unlike most IES-containing genes, the active form of mtA is the unexcised one, allowing the recovery of hypomorphic alleles of essential IES recognition/excision factors. Such is the case of one mutation recovered in the Piwi gene PTIWI09, a key player in small RNA-mediated IES recognition. Another mutation identified a novel protein with a C2H2 zinc finger, mtGa, which is required for excision of a small subset of IESs characterized by enrichment in a 5-bp motif. The unexpected implication of a sequence-specific factor establishes a new paradigm for IES recognition and/or excision.

Dicer-like Enzymes with Sequence Cleavage Preferences

Dicer proteins are known to produce small RNAs (sRNAs) from long double-stranded RNA (dsRNA) templates. These sRNAs are bound by Argonaute proteins, which select the guide strand, often with a 5' end sequence bias. However, Dicer proteins have never been shown to have sequence cleavage preferences. In Paramecium development, two classes of sRNAs that are required for DNA elimination are produced by three Dicer-like enzymes: Dcl2, Dcl3, and Dcl5. Through in vitro cleavage assays, we demonstrate that Dcl2 has a strict size preference for 25 nt and a sequence preference for 5' U and 5' AGA, while Dcl3 has a sequence preference for 5' UNG. Dcl5, however, has cleavage preferences for 5' UAG and 3' CUAC/UN, which leads to the production of RNAs precisely matching short excised DNA elements with corresponding end base preferences. Thus, we characterize three Dicer-like enzymes that are involved in Paramecium development and propose a biological role for their sequence-biased cleavage products.

The Rab7 subfamily across Paramecium aurelia species; evidence of high conservation in sequence and function

We examined sequence conservation and signatures of selection in Rab7 proteins across 11 Paramecium aurelia species, and determined the localization patterns of two P. tetraurelia Rab7 paralogs when expressed as GFP fusions in live cells. We found that, while there is a variable number of Rab7 paralogs per genome, Rab7 genes are highly conserved in sequence and appear to be under strong purifying selection across aurelias. Additionally, and surprisingly based on earlier studies, we found that two P. tetraurelia Rab7 proteins have virtually identical localization patterns. Consistent with this, when we examined the gene family of a highly conserved Rab binding partner across aurelias (Rab-Interacting Lysosomal Protein, or RILP), we found that residues in key binding sites in RILPs were absolutely conserved in 13 of 21 proteins, representing genes from 9 of the 11 species examined. Of note, RILP gene number appears to be even more constrained than Rab7 gene number per genome. Abbreviation: WGD: Whole genome duplication.

C11orf70 Mutations Disrupting the Intraflagellar Transport-Dependent Assembly of Multiple Axonemal Dyneins Cause Primary Ciliary Dyskinesia

Primary ciliary dyskinesia (PCD) is a genetically and phenotypically heterogeneous disorder characterized by destructive respiratory disease and laterality abnormalities due to randomized left-right body asymmetry. PCD is mostly caused by mutations affecting the core axoneme structure of motile cilia that is essential for movement. Genes that cause PCD when mutated include a group that encode proteins essential for the assembly of the ciliary dynein motors and the active transport process that delivers them from their cytoplasmic assembly site into the axoneme. We screened a cohort of affected individuals for disease-causing mutations using a targeted next generation sequencing panel and identified two unrelated families (three affected children) with mutations in the uncharacterized C11orf70 gene (official gene name CFAP300). The affected children share a consistent PCD phenotype from early life with laterality defects and immotile respiratory cilia displaying combined loss of inner and outer dynein arms (IDA+ODA). Phylogenetic analysis shows C11orf70 is highly conserved, distributed across species similarly to proteins involved in the intraflagellar transport (IFT)-dependant assembly of axonemal dyneins. Paramecium C11orf70 RNAi knockdown led to combined loss of ciliary IDA+ODA with reduced cilia beating and swim velocity. Tagged C11orf70 in Paramecium and Chlamydomonas localizes mainly in the cytoplasm with a small amount in the ciliary component. IFT139/TTC21B (IFT-A protein) and FLA10 (IFT kinesin) depletion experiments show that its transport within cilia is IFT dependent. During ciliogenesis, C11orf70 accumulates at the ciliary tips in a similar distribution to the IFT-B protein IFT46. In summary, C11orf70 is essential for assembly of dynein arms and C11orf70 mutations cause defective cilia motility and PCD.

Multimerization properties of PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements

During sexual processes, the ciliate Paramecium eliminates 25–30% of germline DNA from its somatic genome. DNA elimination includes excision of ∼45 000 short, single-copy internal eliminated sequences (IESs) and depends upon PiggyMac (Pgm), a domesticated piggyBac transposase that is essential for DNA cleavage at IES ends. Pgm carries a core transposase region with a putative catalytic domain containing three conserved aspartic acids, and a downstream cysteine-rich (CR) domain. A C-terminal extension of unknown function is predicted to adopt a coiled-coil (CC) structure. To address the role of the three domains, we designed an in vivo complementation assay by expressing wild-type or mutant Pgm-GFP fusions in cells depleted for their endogenous Pgm. The DDD triad and the CR domain are essential for Pgm activity and mutations in either domain have a dominant-negative effect in wild-type cells. A mutant lacking the CC domain is partially active in the presence of limiting Pgm amounts, but inactive when Pgm is completely absent, suggesting that presence of the mutant protein increases the overall number of active complexes. We conclude that IES excision involves multiple Pgm subunits, of which at least a fraction must contain the CC domain.

Centrin diversity and basal body patterning across evolution: new insights from Paramecium

First discovered in unicellular eukaryotes, centrins play crucial roles in basal body duplication and anchoring mechanisms. While the evolutionary status of the founding members of the family, Centrin2/Vfl2 and Centrin3/cdc31 has long been investigated, the evolutionary origin of other members of the family has received less attention. Using a phylogeny of ciliate centrins, we identify two other centrin families, the ciliary centrins and the centrins present in the contractile filaments (ICL centrins). In this paper, we carry on the functional analysis of still not well-known centrins, the ICL1e subfamily identified in Paramecium, and show their requirement for correct basal body anchoring through interactions with Centrin2 and Centrin3. Using Paramecium as well as a eukaryote-wide sampling of centrins from completely sequenced genomes, we revisited the evolutionary story of centrins. Their phylogeny shows that the centrins associated with the ciliate contractile filaments are widespread in eukaryotic lineages and could be as ancient as Centrin2 and Centrin3.

Summary: Functional and phylogenetic analyses reveal the existence of five centrin families and show that basal body patterning in Paramecium requires a third centrin present in many eukaryote lineages.

Basal body positioning and anchoring in the multiciliated cell Paramecium tetraurelia: roles of OFD1 and VFL3

Background

The development of a ciliary axoneme requires the correct docking of the basal body at cytoplasmic vesicles or plasma membrane. In the multiciliated cell Paramecium, three conserved proteins, FOR20, Centrin 2, and Centrin 3 participate in this process, FOR20 and Centrin 2 being involved in the assembly of the transition zone. We investigated the function of two other evolutionary conserved proteins, OFD1 and VFL3, likely involved in this process.

Results

In Paramecium tetraurelia, a single gene encodes OFD1, while four genes encode four isoforms of VFL3, grouped into two families, VFL3-A and VFL3-B. Depletion of OFD1 and the sole VFL3-A family impairs basal body docking. Loss of OFD1 yields a defective assembly of the basal body distal part. Like FOR20, OFD1 is recruited early during basal body assembly and localizes at the transition zone between axoneme and membrane at the level of the microtubule doublets. While the recruitment of OFD1 and Centrin 2 proceed independently, the localizations of OFD1 and FOR20 at the basal body are interdependent. In contrast, in VFL3-A depleted cells, the unanchored basal bodies harbor a fully organized distal part but display an abnormal distribution of their associated rootlets which mark their rotational asymmetry. VFL3-A, which is required for the recruitment of Centrin 3, is transiently present near the basal bodies at an early step of their duplication. VFL3-A localizes at the junction between the striated rootlet and the basal body.

Conclusion

Our results demonstrate the conserved role of OFD1 in the anchoring mechanisms of motile cilia and establish its relations with FOR20 and Centrin 2. They support the hypothesis of its association with microtubule doublets. They suggest that the primary defect of VFL3 depletion is a loss of the rotational asymmetry of the basal body which specifies the sites of assembly of the appendages which guide the movement of basal bodies toward the cell surface. The localization of VFL3 outside of the basal body suggests that extrinsic factors could control this asymmetry.

Electronic supplementary material

The online version of this article (doi:10.1186/s13630-017-0050-z) contains supplementary material, which is available to authorized users.

Early stages of functional diversification in the Rab GTPase gene family revealed by genomic and localization studies in Paramecium species

Rab GTPase family members in Paramecium have higher retention rates and more-divergent expression levels than other genes after whole-genome duplications, consistent with early steps in functional diversification. Localization analysis also uncovers functionally ­diversifying Rab11 genes.

New gene functions arise within existing gene families as a result of gene duplication and subsequent diversification. To gain insight into the steps that led to the functional diversification of paralogues, we tracked duplicate retention patterns, expression-level divergence, and subcellular markers of functional diversification in the Rab GTPase gene family in three Paramecium aurelia species. After whole-genome duplication, Rab GTPase duplicates are more highly retained than other genes in the genome but appear to be diverging more rapidly in expression levels, consistent with early steps in functional diversification. However, by localizing specific Rab proteins in Paramecium cells, we found that paralogues from the two most recent whole-genome duplications had virtually identical localization patterns, and that less closely related paralogues showed evidence of both conservation and diversification. The functionally conserved paralogues appear to target to compartments associated with both endocytic and phagocytic recycling functions, confirming evolutionary and functional links between the two pathways in a divergent eukaryotic lineage. Because the functionally diversifying paralogues are still closely related to and derived from a clade of functionally conserved Rab11 genes, we were able to pinpoint three specific amino acid residues that may be driving the change in the localization and thus the function in these proteins.

Whole Genome Sequencing Identifies a Novel Factor Required for Secretory Granule Maturation in Tetrahymena thermophila

Unbiased genetic approaches have a unique ability to identify novel genes associated with specific biological pathways. Thanks to next generation sequencing, forward genetic strategies can be expanded to a wider range of model organisms. The formation of secretory granules, called mucocysts, in the ciliate Tetrahymena thermophila relies, in part, on ancestral lysosomal sorting machinery, but is also likely to involve novel factors. In prior work, multiple strains with defects in mucocyst biogenesis were generated by nitrosoguanidine mutagenesis, and characterized using genetic and cell biological approaches, but the genetic lesions themselves were unknown. Here, we show that analyzing one such mutant by whole genome sequencing reveals a novel factor in mucocyst formation. Strain UC620 has both morphological and biochemical defects in mucocyst maturation-a process analogous to dense core granule maturation in animals. Illumina sequencing of a pool of UC620 F2 clones identified a missense mutation in a novel gene called MMA1 (Mucocyst maturation). The defects in UC620 were rescued by expression of a wild-type copy of MMA1, and disrupting MMA1 in an otherwise wild-type strain phenocopies UC620. The product of MMA1, characterized as a CFP-tagged copy, encodes a large soluble cytosolic protein. A small fraction of Mma1p-CFP is pelletable, which may reflect association with endosomes. The gene has no identifiable homologs except in other Tetrahymena species, and therefore represents an evolutionarily recent innovation that is required for granule maturation.

Ku-mediated coupling of DNA cleavage and repair during programmed genome rearrangements in the ciliate Paramecium tetraurelia

During somatic differentiation, physiological DNA double-strand breaks (DSB) can drive programmed genome rearrangements (PGR), during which DSB repair pathways are mobilized to safeguard genome integrity. Because of their unique nuclear dimorphism, ciliates are powerful unicellular eukaryotic models to study the mechanisms involved in PGR. At each sexual cycle, the germline nucleus is transmitted to the progeny, but the somatic nucleus, essential for gene expression, is destroyed and a new somatic nucleus differentiates from a copy of the germline nucleus. In Paramecium tetraurelia, the development of the somatic nucleus involves massive PGR, including the precise elimination of at least 45,000 germline sequences (Internal Eliminated Sequences, IES). IES excision proceeds through a cut-and-close mechanism: a domesticated transposase, PiggyMac, is essential for DNA cleavage, and DSB repair at excision sites involves the Ligase IV, a specific component of the non-homologous end-joining (NHEJ) pathway. At the genome-wide level, a huge number of programmed DSBs must be repaired during this process to allow the assembly of functional somatic chromosomes. To understand how DNA cleavage and DSB repair are coordinated during PGR, we have focused on Ku, the earliest actor of NHEJ-mediated repair. Two Ku70 and three Ku80 paralogs are encoded in the genome of P. tetraurelia: Ku70a and Ku80c are produced during sexual processes and localize specifically in the developing new somatic nucleus. Using RNA interference, we show that the development-specific Ku70/Ku80c heterodimer is essential for the recovery of a functional somatic nucleus. Strikingly, at the molecular level, PiggyMac-dependent DNA cleavage is abolished at IES boundaries in cells depleted for Ku80c, resulting in IES retention in the somatic genome. PiggyMac and Ku70a/Ku80c co-purify as a complex when overproduced in a heterologous system. We conclude that Ku has been integrated in the Paramecium DNA cleavage factory, enabling tight coupling between DSB introduction and repair during PGR.

DNA double-strand breaks (DSBs) are potential threats for chromosome stability, but they are usually repaired by two major pathways, homologous recombination or non-homologous end joining (NHEJ). DSBs can also be essential during physiological processes, such as the programmed removal of germline sequences that takes place in various eukaryotes, including ciliates, during somatic differentiation. We use the ciliate Paramecium tetraurelia as a unicellular model to study how DNA breakage and DSB repair are coordinated during programmed genome rearrangements. In this organism, assembly of the somatic genome involves the elimination of ∼25% of germline DNA, including the precise excision of thousands of short Internal Eliminated Sequences (IES) scattered along germline chromosomes. A domesticated piggyBac transposase, PiggyMac, is required for double-strand DNA cleavage at IES ends and IES excision sites are very precisely repaired by the NHEJ pathway. Here, we report that a specialized Ku heterodimer, specifically expressed during programmed genome rearrangements, is an essential partner of PiggyMac and activates DNA cleavage. We propose that incorporation of DSB repair proteins in a pre-cleavage complex constitutes a safe and efficient way for Paramecium to direct thousands of programmed DSBs to the NHEJ pathway and make sure that somatic chromosomes are assembled correctly.

A forward genetic screen reveals essential and non-essential RNAi factors in Paramecium tetraurelia

In most eukaryotes, small RNA-mediated gene silencing pathways form complex interacting networks. In the ciliate Paramecium tetraurelia, at least two RNA interference (RNAi) mechanisms coexist, involving distinct but overlapping sets of protein factors and producing different types of short interfering RNAs (siRNAs). One is specifically triggered by high-copy transgenes, and the other by feeding cells with double-stranded RNA (dsRNA)-producing bacteria. In this study, we designed a forward genetic screen for mutants deficient in dsRNA-induced silencing, and a powerful method to identify the relevant mutations by whole-genome sequencing. We present a set of 47 mutant alleles for five genes, revealing two previously unknown RNAi factors: a novel Paramecium-specific protein (Pds1) and a Cid1-like nucleotidyl transferase. Analyses of allelic diversity distinguish non-essential and essential genes and suggest that the screen is saturated for non-essential, single-copy genes. We show that non-essential genes are specifically involved in dsRNA-induced RNAi while essential ones are also involved in transgene-induced RNAi. One of the latter, the RNA-dependent RNA polymerase RDR2, is further shown to be required for all known types of siRNAs, as well as for sexual reproduction. These results open the way for the dissection of the genetic complexity, interconnection, mechanisms and natural functions of RNAi pathways in P. tetraurelia.

Calcium signaling in closely related protozoan groups (Alveolata): non-parasitic ciliates (Paramecium, Tetrahymena) vs. parasitic Apicomplexa (Plasmodium, Toxoplasma)

The importance of Ca2+-signaling for many subcellular processes is well established in higher eukaryotes, whereas information about protozoa is restricted. Recent genome analyses have stimulated such work also with Alveolates, such as ciliates (Paramecium, Tetrahymena) and their pathogenic close relatives, the Apicomplexa (Plasmodium, Toxoplasma). Here we compare Ca2+ signaling in the two closely related groups. Acidic Ca2+ stores have been characterized in detail in Apicomplexa, but hardly in ciliates. Two-pore channels engaged in Ca2+-release from acidic stores in higher eukaryotes have not been stingently characterized in either group. Both groups are endowed with plasma membrane- and endoplasmic reticulum-type Ca2+-ATPases (PMCA, SERCA), respectively. Only recently was it possible to identify in Paramecium a number of homologs of ryanodine and inositol 1,3,4-trisphosphate receptors (RyR, IP3R) and to localize them to widely different organelles participating in vesicle trafficking. For Apicomplexa, physiological experiments suggest the presence of related channels although their identity remains elusive. In Paramecium, IP3Rs are constitutively active in the contractile vacuole complex; RyR-related channels in alveolar sacs are activated during exocytosis stimulation, whereas in the parasites the homologous structure (inner membrane complex) may no longer function as a Ca2+ store. Scrutinized comparison of the two closely related protozoan phyla may stimulate further work and elucidate adaptation to parasitic life. See also "Conclusions" section.

FOR20, a conserved centrosomal protein, is required for assembly of the transition zone and basal body docking at the cell surface

Within the FOP family of centrosomal proteins, the conserved FOR20 protein has been implicated in the control of primary cilium assembly in human cells. To ascertain its role in ciliogenesis, we have investigated the function of its ortholog, PtFOR20p, in a multiciliated unicellular organism, Paramecium. By a combined functional and cytological analysis, we found that PtFOR20p specifically localizes at basal bodies and is required to build the transition zone, a prerequisite to their maturation and docking at the cell surface, hence to ciliogenesis. We also found that PtCen2p (one of the two basal body specific centrins, ortholog of HsCen2) is required to recruit PtFOR20p at the developing basal body and to control its length. In contrast, the other basal body specific centrin, PtCen3p, is not needed for assembly of the transition zone, but is required downstream, for basal body docking. Comparison of the structural defects induced by depletion of PtFOR20p, PtCen2p or PtCen3p respectively illustrates the dual role of the transition zone in the biogenesis of the basal body and in cilium assembly. The multiple potential roles of the transition zone during basal body biogenesis and the evolutionary conserved function of the FOP proteins in microtubule membrane interactions are discussed.

Sas-4 proteins are required during basal body duplication in Paramecium

This study investigated the role of Sas-4 in basal body duplication in Paramecium and found that Sas-4 proteins are required to assemble and stabilize the germinative disk and cartwheel, which suggests that Sas-4 plays an essential role in basal body duplication.

Centrioles and basal bodies are structurally related organelles composed of nine microtubule (MT) triplets. Studies performed in Caenorhabditis elegans embryos have shown that centriole duplication takes place in sequential way, in which different proteins are recruited in a specific order to assemble a procentriole. ZYG-1 initiates centriole duplication by triggering the recruitment of a complex of SAS-5 and SAS-6, which then recruits the final player, SAS-4, to allow the incorporation of MT singlets. It is thought that a similar mechanism (that also involves additional proteins) is present in other animal cells, but it remains to be investigated whether the same players and their ascribed functions are conserved during basal body duplication in cells that exclusively contain basal bodies. To investigate this question, we have used the multiciliated protist Paramecium tetraurelia. Here we show that in the absence of PtSas4, two types of defects in basal body duplication can be identified. In the majority of cases, the germinative disk and cartwheel, the first structures assembled during duplication, are not detected. In addition, if daughter basal bodies were formed, they invariably had defects in MT recruitment. Our results suggest that PtSas4 has a broader function than its animal orthologues.

Functional specialization of Piwi proteins in Paramecium tetraurelia from post-transcriptional gene silencing to genome remodelling

Proteins of the Argonaute family are small RNA carriers that guide regulatory complexes to their targets. The family comprises two major subclades. Members of the Ago subclade, which are present in most eukaryotic phyla, bind different classes of small RNAs and regulate gene expression at both transcriptional and post-transcriptional levels. Piwi subclade members appear to have been lost in plants and fungi and were mostly studied in metazoa, where they bind piRNAs and have essential roles in sexual reproduction. Their presence in ciliates, unicellular organisms harbouring both germline micronuclei and somatic macronuclei, offers an interesting perspective on the evolution of their functions. Here, we report phylogenetic and functional analyses of the 15 Piwi genes from Paramecium tetraurelia. We show that four constitutively expressed proteins are involved in siRNA pathways that mediate gene silencing throughout the life cycle. Two other proteins, specifically expressed during meiosis, are required for accumulation of scnRNAs during sexual reproduction and for programmed genome rearrangements during development of the somatic macronucleus. Our results indicate that Paramecium Piwi proteins have evolved to perform both vegetative and sexual functions through mechanisms ranging from post-transcriptional mRNA cleavage to epigenetic regulation of genome rearrangements.

Gene expression in a paleopolyploid: a transcriptome resource for the ciliate Paramecium tetraurelia

Background

The genome of Paramecium tetraurelia, a unicellular model that belongs to the ciliate phylum, has been shaped by at least 3 successive whole genome duplications (WGD). These dramatic events, which have also been documented in plants, animals and fungi, are resolved over evolutionary time by the loss of one duplicate for the majority of genes. Thanks to a low rate of large scale genome rearrangement in Paramecium, an unprecedented large number of gene duplicates of different ages have been identified, making this organism an outstanding model to investigate the evolutionary consequences of polyploidization. The most recent WGD, with 51% of pre-duplication genes still in 2 copies, provides a snapshot of a phase of rapid gene loss that is not accessible in more ancient polyploids such as yeast.

Results

We designed a custom oligonucleotide microarray platform for P. tetraurelia genome-wide expression profiling and used the platform to measure gene expression during 1) the sexual cycle of autogamy, 2) growth of new cilia in response to deciliation and 3) biogenesis of secretory granules after massive exocytosis. Genes that are differentially expressed during these time course experiments have expression patterns consistent with a very low rate of subfunctionalization (partition of ancestral functions between duplicated genes) in particular since the most recent polyploidization event.

Conclusions

A public transcriptome resource is now available for Paramecium tetraurelia. The resource has been integrated into the ParameciumDB model organism database, providing searchable access to the data. The microarray platform, freely available through NimbleGen Systems, provides a robust, cost-effective approach for genome-wide expression profiling in P. tetraurelia. The expression data support previous studies showing that at short evolutionary times after a whole genome duplication, gene dosage balance constraints and not functional change are the major determinants of gene retention.

Basal body duplication in Paramecium: the key role of Bld10 in assembly and stability of the cartwheel

Basal bodies which nucleate cilia and flagella, and centrioles which organize centrosomes share the same architecture characterized by the ninefold symmetry of their microtubular shaft. Among the conserved proteins involved in the biogenesis of the canonical 9-triplet centriolar structures, Sas-6 and Bld10 proteins have been shown to play central roles in the early steps of assembly and in establishment/stabilization of the ninefold symmetry. Using fluorescent tagged proteins and RNAi to study the localization and function of these two proteins in Paramecium, we focused on the early effects of their depletion, the consequences of their overexpression and their functional interdependence. We find that both genes are essential and their depletion affects cartwheel assembly and hence basal body duplication. We also show that, contrary to Sas6p, Bld10p is not directly responsible for the establishment of the ninefold symmetry, but is required not only for new basal body assembly and stability but also for Sas6p maintenance at mature basal bodies. Finally, ultrastructural analysis of cells overexpressing either protein revealed two types of early assembly intermediates, hub-like structures and generative discs, suggesting a conserved scaffolding process.

Protein phosphatase 2B (PP2B, calcineurin) in Paramecium: partial characterization reveals that two members of the unusually large catalytic subunit family have distinct roles in calcium-dependent processes

We characterized the calcineurin (CaN) gene family, including the subunits CaNA and CaNB, based upon sequence information obtained from the Paramecium genome project. Paramecium tetraurelia has seven subfamilies of the catalytic CaNA subunit and one subfamily of the regulatory CaNB subunit, with each subfamily having two members of considerable identity on the amino acid level (>or=55% between subfamilies, >or=94% within CaNA subfamilies, and full identity in the CaNB subfamily). Within CaNA subfamily members, the catalytic domain and the CaNB binding region are highly conserved and molecular modeling revealed a three-dimensional structure almost identical to a human ortholog. At 14 members, the size of the CaNA family is unprecedented, and we hypothesized that the different CaNA subfamily members were not strictly redundant and that at least some fulfill different roles in the cell. This was tested by selecting two phylogenetically distinct members of this large family for posttranscriptional silencing by RNA interference. The two targets resulted in differing effects in exocytosis, calcium dynamics, and backward swimming behavior that supported our hypothesis that the large, highly conserved CaNA family members are not strictly redundant and that at least two members have evolved diverse but overlapping functions. In sum, the occurrence of CaN in Paramecium spp., although disputed in the past, has been established on a molecular level. Its role in exocytosis and ciliary beat regulation in a protozoan, as well as in more complex organisms, suggests that these roles for CaN were acquired early in the evolution of this protein family.

Distinct subcellular localization of a group of synaptobrevin-like SNAREs in Paramecium tetraurelia and effects of silencing SNARE-specific chaperone NSF

We have identified new synaptobrevin-like SNAREs and localized the corresponding gene products with green fluorescent protein (GFP)-fusion constructs and specific antibodies at the light and electron microscope (EM) levels. These SNAREs, named Paramecium tetraurelia synaptobrevins 8 to 12 (PtSyb8 to PtSyb12), showed mostly very restricted, specific localization, as they were found predominantly on structures involved in endo- or phagocytosis. In summary, we found PtSyb8 and PtSyb9 associated with the nascent food vacuole, PtSyb10 near the cell surface, at the cytostome, and in close association with ciliary basal bodies, and PtSyb11 on early endosomes and on one side of the cytostome, while PtSyb12 was found in the cytosol. PtSyb4 and PtSyb5 (identified previously) were localized on small vesicles, PtSyb5 probably being engaged in trichocyst (dense core secretory vesicle) processing. PtSyb4 and PtSyb5 are related to each other and are the furthest deviating of all SNAREs identified so far. Because they show no similarity with any other R-SNAREs outside ciliates, they may represent a ciliate-specific adaptation. PtSyb10 forms small domains near ciliary bases, and silencing slows down cell rotation during depolarization-induced ciliary reversal. NSF silencing supports a function of cell surface SNAREs by revealing vesicles along the cell membrane at sites normally devoid of vesicles. The distinct distributions of these SNAREs emphasize the considerable differentiation of membrane trafficking, particularly along the endo-/phagocytic pathway, in this protozoan.

Bug22p, a conserved centrosomal/ciliary protein also present in higher plants, is required for an effective ciliary stroke in Paramecium

Centrioles, cilia, and flagella are ancestral conserved organelles of eukaryotic cells. Among the proteins identified in the proteomics of ciliary proteins in Paramecium, we focus here on a protein, Bug22p, previously detected by cilia and basal-body high-throughput studies but never analyzed per se. Remarkably, this protein is also present in plants, which lack centrioles and cilia. Bug22p sequence alignments revealed consensus positions that distinguish species with centrioles/cilia from plants. In Paramecium, antibody and green fluorescent protein (GFP) fusion labeling localized Bug22p in basal bodies and cilia, and electron microscopy immunolabeling refined the localization to the terminal plate of the basal bodies, the transition zone, and spots along the axoneme, preferentially between the membrane and the microtubules. RNA interference (RNAi) depletion of Bug22p provoked a strong decrease in swimming speed, followed by cell death after a few days. High-speed video microscopy and morphological analysis of Bug22p-depleted cells showed that the protein plays an important role in the efficiency of ciliary movement by participating in the stroke shape and rigidity of cilia. The defects in cell swimming and growth provoked by RNAi can be complemented by expression of human Bug22p. This is the first reported case of complementation by a human gene in a ciliate.

Novel types of Ca2+ release channels participate in the secretory cycle of Paramecium cells

A database search of the Paramecium genome reveals 34 genes related to Ca(2+)-release channels of the inositol-1,4,5-trisphosphate (IP(3)) or ryanodine receptor type (IP(3)R, RyR). Phylogenetic analyses show that these Ca(2+) release channels (CRCs) can be subdivided into six groups (Paramecium tetraurelia CRC-I to CRC-VI), each one with features in part reminiscent of IP(3)Rs and RyRs. We characterize here the P. tetraurelia CRC-IV-1 gene family, whose relationship to IP(3)Rs and RyRs is restricted to their C-terminal channel domain. CRC-IV-1 channels localize to cortical Ca(2+) stores (alveolar sacs) and also to the endoplasmic reticulum. This is in contrast to a recently described true IP(3) channel, a group II member (P. tetraurelia IP(3)R(N)-1), found associated with the contractile vacuole system. Silencing of either one of these CRCs results in reduced exocytosis of dense core vesicles (trichocysts), although for different reasons. Knockdown of P. tetraurelia IP(3)R(N) affects trichocyst biogenesis, while CRC-IV-1 channels are involved in signal transduction since silenced cells show an impaired release of Ca(2+) from cortical stores in response to exocytotic stimuli. Our discovery of a range of CRCs in Paramecium indicates that protozoans already have evolved multiple ways for the use of Ca(2+) as signaling molecule.

Biogenesis of Dense-Core Secretory Granules

Dense core granules (DCGs) are vesicular organelles derived from outbound traffic through the eukaryotic secretory pathway. As DCGs are formed, the secretory pathway can also give rise to other types of vesicles, such as those bound for endosomes, lysosomes, and the cell surface. DCGs differ from these other vesicular carriers in both content and function, storing highly concentrated cores’ of condensed cargo in vesicles that are stably maintained within the cell until a specific extracellular stimulus causes their fusion with the plasma membrane. These unique features are imparted by the activities of membrane and lumenal proteins that are specifically delivered to the vesicles during synthesis. This chapter will describe the DCG biogenesis pathway, beginning with the sorting of DCG proteins from proteins that are destined for other types of vesicle carriers. In the trans-Golgi network (TGN), sorting occurs as DCG proteins aggregate, causing physical separation from non-DCG proteins. Recent work addresses the nature of interactions that produce these aggregates, as well as potentially important interactions with membranes and membrane proteins. DCG proteins are released from the TGN in vesicles called immature secretory granules (ISGs). The mechanism of ISG formation is largely unclear but is not believed to rely on the assembly of vesicle coats like those observed in other secretory pathways. The required cytosolic factors are now beginning to be identified using in vitro systems with purified cellular components. ISG transformation into a mature fusion-competent, stimulus-dependent DCG occurs as endoproteolytic processing of many DCG proteins causes continued condensation of the lumenal contents. At the same time, proteins that fail to be incorporated into the condensing core are removed by a coat-mediated budding mechanism, which also serves to remove excess membrane and membrane proteins from the maturing vesicle. This chapter will summarize the work leading to our current view of granule synthesis, and will discuss questions that need to be addressed in order to gain a more complete understanding of the pathway.

Rapid downregulation of the Ca2+-signal after exocytosis stimulation in Paramecium cells: essential role of a centrin-rich filamentous cortical network, the infraciliary lattice

We analysed in Paramecium tetraurelia cells the role of the infraciliary lattice, a cytoskeletal network containing numerous centrin isoforms tightly bound to large binding proteins, in the re-establishment of Ca2+ homeostasis following exocytosis stimulation. The wild type strain d4-2 has been compared with the mutant cell line Delta-PtCenBP1 which is devoid of the infraciliary lattice ("Delta-PtCenBP1" cells). Exocytosis is known to involve the mobilization of cortical Ca2+-stores and a superimposed Ca2+-influx and was analysed using Fura Red ratio imaging. No difference in the initial signal generation was found between wild type and Delta-PtCenBP1 cells. In contrast, decay time was greatly increased in Delta-PtCenBP1 cells particularly when stimulated, e.g., in presence of 1mM extracellular Ca2+, [Ca2+]o. Apparent halftimes of f/f0 decrease were 8.5 s in wild type and approximately 125 s in Delta-PtCenBP1 cells, requiring approximately 30 s and approximately 180 s, respectively, to re-establish intracellular [Ca2+] homeostasis. Lowering [Ca2+]o to 0.1 and 0.01 mM caused an acceleration of intracellular [Ca2+] decay to t(1/2)=33 s and 28 s, respectively, in Delta-PtCenBP1 cells as compared to 8.1 and 5.6, respectively, for wild type cells. We conclude that, in Paramecium cells, the infraciliary lattice is the most efficient endogenous Ca2+ buffering system allowing the rapid downregulation of Ca2+ signals after exocytosis stimulation.

Maternal noncoding transcripts antagonize the targeting of DNA elimination by scanRNAs in Paramecium tetraurelia

The germline genome of ciliates is extensively rearranged during the development of a new somatic macronucleus from the germline micronucleus, after sexual events. In Paramecium tetraurelia, single-copy internal eliminated sequences (IESs) are precisely excised from coding sequences and intergenic regions. For a subset of IESs, introduction of the IES sequence into the maternal macronucleus specifically inhibits excision of the homologous IES in the developing zygotic macronucleus, suggesting that epigenetic regulation of excision involves a global comparison of germline and somatic genomes. ScanRNAs (scnRNAs) produced during micronuclear meiosis by a developmentally regulated RNAi pathway have been proposed to mediate this transnuclear cross-talk. In this study, microinjection experiments provide direct evidence that 25-nucleotide (nt) scnRNAs promote IES excision. We further show that noncoding RNAs are produced from the somatic maternal genome, both during vegetative growth and during sexual events. Maternal inhibition of IES excision is abolished when maternal somatic transcripts containing an IES are targeted for degradation by a distinct RNAi pathway involving 23-nt siRNAs. The results strongly support a scnRNA/macronuclear RNA scanning model in which a natural genomic subtraction, occurring during meiosis between deletion-inducing scnRNAs and antagonistic transcripts from the maternal macronucleus, regulates rearrangements of the zygotic genome.

Translational control of intron splicing in eukaryotes

Most eukaryotic genes are interrupted by non-coding introns that must be accurately removed from pre-messenger RNAs to produce translatable mRNAs. Splicing is guided locally by short conserved sequences, but genes typically contain many potential splice sites, and the mechanisms specifying the correct sites remain poorly understood. In most organisms, short introns recognized by the intron definition mechanism cannot be efficiently predicted solely on the basis of sequence motifs. In multicellular eukaryotes, long introns are recognized through exon definition and most genes produce multiple mRNA variants through alternative splicing. The nonsense-mediated mRNA decay (NMD) pathway may further shape the observed sets of variants by selectively degrading those containing premature termination codons, which are frequently produced in mammals. Here we show that the tiny introns of the ciliate Paramecium tetraurelia are under strong selective pressure to cause premature termination of mRNA translation in the event of intron retention, and that the same bias is observed among the short introns of plants, fungi and animals. By knocking down the two P. tetraurelia genes encoding UPF1, a protein that is crucial in NMD, we show that the intrinsic efficiency of splicing varies widely among introns and that NMD activity can significantly reduce the fraction of unspliced mRNAs. The results suggest that, independently of alternative splicing, species with large intron numbers universally rely on NMD to compensate for suboptimal splicing efficiency and accuracy.

Analysis of sequence variability in the macronuclear DNA of Paramecium tetraurelia: a somatic view of the germline

Ciliates are the only unicellular eukaryotes known to separate germinal and somatic functions. Diploid but silent micronuclei transmit the genetic information to the next sexual generation. Polyploid macronuclei express the genetic information from a streamlined version of the genome but are replaced at each sexual generation. The macronuclear genome of Paramecium tetraurelia was recently sequenced by a shotgun approach, providing access to the gene repertoire. The 72-Mb assembly represents a consensus sequence for the somatic DNA, which is produced after sexual events by reproducible rearrangements of the zygotic genome involving elimination of repeated sequences, precise excision of unique-copy internal eliminated sequences (IES), and amplification of the cellular genes to high copy number. We report use of the shotgun sequencing data (>10(6) reads representing 13 x coverage of a completely homozygous clone) to evaluate variability in the somatic DNA produced by these developmental genome rearrangements. Although DNA amplification appears uniform, both of the DNA elimination processes produce sequence heterogeneity. The variability that arises from IES excision allowed identification of hundreds of putative new IESs, compared to 42 that were previously known, and revealed cases of erroneous excision of segments of coding sequences. We demonstrate that IESs in coding regions are under selective pressure to introduce premature termination of translation in case of excision failure.

Functional diversification of centrins and cell morphological complexity

In addition to their key role in the duplication of microtubule organising centres (MTOCs), centrins are major constituents of diverse MTOC-associated contractile arrays. A centrin partner, Sfi1p, has been characterised in yeast as a large protein carrying multiple centrin-binding sites, suggesting a model for centrin-mediated Ca2+-induced contractility and for the duplication of MTOCs. In vivo validation of this model has been obtained in Paramecium, which possesses an extended contractile array - the infraciliary lattice (ICL) - essentially composed of centrins and a huge Sfi1p-like protein, PtCenBP1p, which is essential for ICL assembly and contractility. The high molecular diversity revealed here by the proteomic analysis of the ICL, including ten subfamilies of centrins and two subfamilies of Sf1p-like proteins, led us to address the question of the functional redundancy, either between the centrin-binding proteins or between the centrin subfamilies. We show that all are essential for ICL biogenesis. The two centrin-binding protein subfamilies and nine of the centrin subfamilies are ICL specific and play a role in its molecular and supramolecular architecture. The tenth and most conserved centrin subfamily is present at three cortical locations (ICL, basal bodies and contractile vacuole pores) and might play a role in coordinating duplication and positioning of cortical organelles.

A role for convergent evolution in the secretory life of cells

The role of convergent evolution in biological adaptation is increasingly appreciated. Many clear examples have been described at the level of individual proteins and for organismal morphology, and convergent mechanisms have even been invoked to account for similar community structures that are shared between ecosystems. At the cellular level, an important area that has received scant attention is the potential influence of convergent evolution on complex subcellular features, such as organelles. Here, we show that existing data strongly argue that convergent evolution underlies the similar properties of specialized secretory vesicles, called dense core granules, in the animal and ciliate lineages. We discuss both the criteria for judging convergent evolution and the contribution that such evolutionary analysis can make to improve our understanding of processes in cell biology. The elucidation of these underlying evolutionary relationships is vital because cellular structures that are assumed to be analogous, owing to shared features, might in fact be governed by different molecular mechanisms.

A broad spectrum of actin paralogs inParamecium tetraureliacells display differential localization and function

To localize the different actin paralogs found in Paramecium and to disclose functional implications, we used overexpression of GFP-fusion proteins and antibody labeling, as well as gene silencing. Several isoforms are associated with food vacuoles of different stages. GFP-actin either forms a tail at the lee side of the organelle, or it is vesicle bound in a homogenous or in a speckled arrangement, thus reflecting an actin-based mosaic of the phagosome surface appropriate for association and/or dissociation of other vesicles upon travel through the cell. Several paralogs occur in cilia. A set of actins is found in the cell cortex where actin outlines the regular surface pattern. Labeling of defined structures of the oral cavity is due to other types of actin, whereas yet more types are distributed in a pattern suggesting association with the numerous Golgi fields. A substantial fraction of actins is associated with cytoskeletal elements that are known to be composed of other proteins. Silencing of the respective actin genes or gene subfamilies entails inhibitory effects on organelles compatible with localization studies. Knock down of the actin found in the cleavage furrow abolishes cell division, whereas silencing of other actin genes alters vitality, cell shape and swimming behavior.

Nowa1p and Nowa2p: novel putative RNA binding proteins involved in trans-nuclear crosstalk in Paramecium tetraurelia

BACKGROUND

The germline genome of ciliates is extensively rearranged during development of a new somatic macronucleus from the germline micronucleus, a process that follows sexual events. In Paramecium tetraurelia, single-copy internal eliminated sequences (IESs) and multicopy transposons are eliminated, whereas cellular genes are amplified to approximately 800 n. For a subset of IESs, introduction of the IES sequence into the maternal (prezygotic) macronucleus specifically inhibits excision of the homologous IES in the developing zygotic macronucleus. This and other homology-dependent maternal effects have suggested that rearrangement patterns are epigenetically determined by an RNA-mediated, trans-nuclear comparison, involving the RNA interference pathway, of germline and somatic genomes.

RESULTS

We report the identification of novel developmentally regulated RNA binding proteins, Nowa1p and Nowa2p, which are required for the survival of sexual progeny. Green fluorescent protein (GFP) fusions show that Nowa1p accumulates into the maternal macronucleus shortly before meiosis of germline micronuclei and is later transported to developing macronuclei. Nowa1p/2p depletion impairs the elimination of transposons and of those IESs that are controlled by maternal effects, confirming the existence of distinct IES classes.

CONCLUSIONS

The results indicate that Nowa proteins are essential components of the trans-nuclear-crosstalk mechanism that is responsible for epigenetic programming of genome rearrangements. We discuss implications for the current models of genome scanning in ciliates, a process related to the formation of heterochromatin by RNA interference in other eukaryotes.

A Multigene Family Encoding R-SNAREs in the Ciliate Paramecium tetraurelia

SNARE proteins (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) mediate membrane interactions and are conventionally divided into Q-SNAREs and R-SNAREs according to the possession of a glutamine or arginine residue at the core of their SNARE domain. Here, we describe a set of R-SNAREs from the ciliate Paramecium tetraurelia consisting of seven families encoded by 12 genes that are expressed simultaneously. The complexity of the endomembrane system in Paramecium can explain this high number of genes. All P. tetraurelia synaptobrevins (PtSybs) possess a SNARE domain and show homology to the Longin family of R-SNAREs such as Ykt6, Sec22 and tetanus toxin-insensitive VAMP (TI-VAMP). We localized four exemplary PtSyb subfamilies with GFP constructs and antibodies on the light and electron microscopic level. PtSyb1-1, PtSyb1-2 and PtSyb3-1 were found in the endoplasmic reticulum, whereas PtSyb2 is localized exclusively in the contractile vacuole complex. PtSyb6 was found cytosolic but also resides in regularly arranged structures at the cell cortex (parasomal sacs), the cytoproct and oral apparatus, probably representing endocytotic compartments. With gene silencing, we showed that the R-SNARE of the contractile vacuole complex, PtSyb2, functions to maintain structural integrity as well as functionality of the osmoregulatory system but also affects cell division.