---
title: "Modeling human monogenic diseases using the choanoflagellate Salpingoeca rosetta"
authors:
- "Seemay Chou"
- "Dennis A. Sun"
doi: "10.57844/arcadia-bp0f-v1xx"
license: "https://creativecommons.org/licenses/by/4.0/"
date: "2025-06-22"
version: 2
canonical_url: "https://thestacks.org/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta"
---
# Modeling human monogenic diseases using the choanoflagellate *Salpingoeca rosetta*
_We applied a decision-making framework for identifying tractable genes from our organism selection dataset for pilot experiments in the choanoflagellate *Salpingoeca rosetta*._
## Abstract
We developed a framework for identifying actionable disease–gene pairs from organism selection data [](https://doi.org/10.57844/arcadia-ugg5-emyd). Here, we apply that framework to genes in the choanoflagellate *Salpingoeca rosetta* and identify seven genes to pursue for experimental testing. We share information about *S. rosetta*, including its life history and currently available genetic tools. We also provide examples of genes we rejected or were particularly excited about in our decision-making process.
Scientists interested in using choanoflagellates or other microbial eukaryotes to model human disease may find our analyses and experimental design process useful in their work. We include scripts to help visualize gene expression in *S. rosetta* and a summary of the state of technology, which can aid in scientific decision-making.
- All **code** related to this pub is available in [this GitHub repository](https://github.com/Arcadia-Science/2025-zoogle-experiments). Code specific to *S. rosetta* is in [this analysis notebook](https://github.com/Arcadia-Science/2025-zoogle-experiments/blob/main/notebooks/Salpingoeca-rosetta/2_salpingoeca-rosetta-expression.ipynb) and [these visualization scripts](https://github.com/Arcadia-Science/2025-zoogle-collabs/tree/main/zoogletools/salpingoeca).
- **Data**, including Foldseek and ProteinCartography results for certain *S. rosetta* proteins, is available on [Zenodo](https://zenodo.org/records/15724261).
# About the organism {#h.ywnpx1mo4cnu}
*Salpingoeca rosetta* is a choanoflagellate and an emerging genetic model organism [](https://doi.org/10.1016/bs.ctdb.2022.01.001). Choanoflagellates are a group of microbial eukaryotes that share a common ancestor with the last common ancestor of animals [](https://doi.org/10.1016/j.cub.2005.02.004). Research groups often study this diverse group of organisms to understand the origins of multicellularity. Choanoflagellates often have transient multicellular life stages generated through cell division [](https://doi.org/10.1016/j.cub.2010.09.014). Cooperativity in these colonies provides insight into the transition from unicellular to multicellular life.
Free-swimming single *S. rosetta* cells have a stereotypical cell shape characterized by a flagellum circumscribed by an array of microvilli referred to as a “collar” [](https://doi.org/10.1016/j.ydbio.2011.06.003) ([Figure 1](#Z03), A). *S. rosetta* cells can divide to form fragile “chains” of related cells attached by cell–cell adhesion [](https://doi.org/10.1016/j.ydbio.2011.06.003). *S. rosetta* also forms spherical colonies referred to as “rosettes” in response to a sulfonolipid produced by the bacteria *Algorphiagus machipongonensis* [](https://doi.org/10.1073/pnas.1605015113). In addition to its multicellular life stages, *S. rosetta* also has diverse unicellular life stages, including slow swimmers, fast swimmers, and thecate cells ([Figure 1](#Z03), A) [](https://doi.org/10.1016/j.ydbio.2011.06.003). *S. rosetta* can even form amoeboid cells in response to physical compression [](https://doi.org/10.7554/elife.61037).
::::::figure{#Z03 type=img src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/asset_e6fc98fc.png' align=center}
:::::image{src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/asset_e6fc98fc.png' width=90% alt='Drawing of cell types in S. rosetta next to its genetic tools. Bar charts show that we trimmed our long list of 41 genes to seven, mostly eliminating candidates where the model wasn’t feasible or there wasn’t an unmet need.'}
:::::
:::::figcaption
**Figure 1:** **Overview of diligence on *Salpingoeca rosetta***.
(A) Illustration of different cell types of *S. rosetta* and schematic of available genetic tools and types of human alleles that we expect to be able to model, given available tools. Cell types include slow swimmers, rosettes, fast swimmers, and thecate cells.
(B) Illustration of manual diligencing outcomes listing the number of genes eliminated by each failure mode. The initial short list was seven genes, all of which we selected for experimental testing.
:::::
::::::
## State of technology {#h.hxekyi8mgixj}
*S. rosetta* has robust genetic tools and genomic resources.
- **Genetic tools**. *S. rosetta* has robust protocols for generating targeted genetic variation. CRISPR-Cas9 mutagenesis is feasible in these organisms through a marker-based approach, and it's possible to generate and maintain stable mutant lines [](https://doi.org/10.7554/elife.56193). Transgenic expression of genes through episomal markers is also possible, with some size limitations [](https://doi.org/10.1091/mbc.e18-08-0514).
- **Phenotypic assays**. *S. rosetta* cells are around 2.5–4 µm in diameter [](https://doi.org/10.1016/j.ydbio.2011.06.003), and can be grown in diverse formats, from flasks to 96-well plates, allowing for robust phenotypic assays across scales. Recent work has shown that it's even possible to screen drug molecules for activity in *S. rosetta* [](https://doi.org/10.1101/2022.08.26.505350). *S. rosetta* cultures are amenable to bulk phenotyping through growth curves and time course imaging. Established protocols for immunofluorescence [](https://doi.org/10.7554/elife.04070), scanning and transmission electron microscopy [](https://doi.org/10.1016/j.ydbio.2011.06.003), and other approaches make it possible to comprehensively identify diverse phenotypes in *S. rosetta* cells.
- **Gene expression**. Bulk RNA-seq data for *S. rosetta* are available for various cell types and stages. We visualized gene expression results from Leon et al. 2025 [](https://doi.org/10.1128/msphere.00917-24), who profiled bulk gene expression across slow swimmers, fast swimmers, rosettes, and thecate cells.
::::::div{.info-box}
You can **generate gene expression visualizations** from *S. rosetta* data by following the examples in [this Jupyter notebook](https://github.com/Arcadia-Science/2025-zoogle-collabs/blob/main/notebooks/Salpingoeca-rosetta/2_salpingoeca-rosetta-expression.ipynb).
::::::
- **Feasible genetic experiments**. Given the wealth of gene-editing tools available in *S. rosetta*, it should be possible to model a variety of different kinds of human disease alleles, including missense, nonsense, and frameshift mutations. The ability to transgenically express genes also makes it possible to test dominant-negative mutations, or to complement mutations with wild-type protein. Conventionally, researchers perform rescue experiments to validate the genetic basis of observed phenotypes by generating CRISPR revertants — mutations that change a previously mutated allele back to the wild-type allele — due to a lack of episomal replication.
# Diligence results {#h.58gy7n44p4do}
We manually diligenced 41 different genes in *Salpingoeca rosetta* to arrive at seven potentially actionable genes for pilot experiments ([Figure 1](#Z03), B).
::::::div{.info-box}
You can see a copy of the **working document** we used to catalog our thoughts for each gene we diligenced in *S. rosetta* [here](https://docs.google.com/spreadsheets/d/1DeJ9X1YGajwvR6zB6bz0RCPz4ClbSYUQDk4jjKj2P2I/edit?usp=sharing).
::::::
::::::div{.info-box}
Download **Foldseek results** for the *S. rosetta* plastin 1 ([F2TWP3](https://www.uniprot.org/uniprotkb/F2TWP3/entry)) and coronin 1a ([F2TZV2](https://www.uniprot.org/uniprotkb/F2TZV2/entry)) homologs, plus **ProteinCartography data** for the *S. rosetta* homolog of the protein encoded by *DDC* ([F2UCR6](https://www.uniprot.org/uniprotkb/F2UCR6/entry)) from [Zenodo](https://zenodo.org/records/15724261) (DOI\: [10.5281/zenodo.15724261](https://doi.org/10.5281/zenodo.15724261)).
::::::
## Example genes we rejected {#h.75u3b76ey9qk}
Here, we describe our specific reasoning for rejecting a handful of select genes from our short list for *S. rosetta* to give you a sense of the variety of reasons we might have passed on a given gene.
### *DDC* (dopamine decarboxylase) {#h.n3xm97d86h8u}
**Lack of homolog confidence**
The *DDC* gene is implicated in deficiency of aromatic-L-amino-acid decarboxylase, which leads to severe postembryonic developmental phenotypes in patients [](https://doi.org/10.1023/a:1005650325003). *S. rosetta* is a top hit for this gene in our organism selection analysis. However, the *S. rosetta* protein (UniProt ID [F2UCR6](https://www.uniprot.org/uniprotkb/F2UCR6/entry)) is putatively annotated as “glutamate decarboxylase,” and FoldSeek searches of this protein don’t return the protein encoded by *DDC* as a top hit.
To understand the broader scope of related gene families, we ran sequence- and structure-based searching and clustering using ProteinCartography [](https://doi.org/10.57844/ARCADIA-A5A6-1068) and observed that the *S. rosetta* protein is found in a different Leiden cluster (LC00) than dopamine decarboxylases from vertebrates (mostly in LC12) and invertebrates (LC05) ([Figure 2](#ZPC); download inputs and results from [Zenodo](https://zenodo.org/records/15724261)). This suggests that the *S. rosetta* protein may not be a dopamine decarboxylase. Due to this uncertainty, we rejected this candidate gene.
::::::figure{#ZPC type=iframe align=center}
:::::iframe{src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/1dAyzLEBEgo_pgsMRLxoNJ8Om_lJZ-zwe' width='100%' height='730px' aria-label='Interactive scatter plot showing that the identified choanoflagellate homolog of the protein encoded by DDC doesn't appear in the same cluster as the protein product of human DDC.'}
:::::
:::::figcaption
**Figure 2:** **Interactive scatter plot that ProteinCartography generated for the choanoflagellate homolog of the protein encoded by *DDC* identified in the organism selection dataset**.
:::::
::::::
### *TREH* (α,α-trehalase) {#h.8awmnwb0t5gd}
**Lack of unmet need**
Trehalose is a disaccharide commonly found in mushrooms and yeast. Patients with a deficiency in α,α-trehalase are unable to process this sugar, resulting in gastrointestinal symptoms such as vomiting and diarrhea [](https://doi.org/10.1017/s0007114500000313). While unpleasant for patients, this disease is treatable by omitting food products containing mushrooms and yeast in patient diets. Due to the low severity of the disease’s symptoms and ease of treatment, we didn’t consider this disease to have a high unmet need, and therefore, we rejected it.
### *HYLS1* (centriolar and ciliogenesis-associated protein HYLS1) {#h.9ur3inrywl89}
**Treatment not possible**
Homozygous loss of function of *HYLS1* causes a lethal malformation syndrome with defects across organ systems, resulting in stillbirth or death shortly after birth [](https://doi.org/10.1093/hmg/ddi157). Due to the lethal and developmental nature of this disease, we didn't believe it would be possible to develop a therapeutic, and therefore rejected this gene.
## Our initial short list {#h.glfmfxf9f8bb}
Here are some of the best-supported case studies from our short list in *S. rosetta*.
::::::figure{#ZS1 type=img src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/1TXHtTsIFyVaD6oMGQQtHpUqOC-hjTn3Z' align=center}
:::::image{src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/1TXHtTsIFyVaD6oMGQQtHpUqOC-hjTn3Z' width=85% alt='Structural alignments of proteins from humans and S. rosetta show high structural similarity due to their overlapping shapes.'}
:::::
:::::figcaption
**Figure 3:** **Structural comparisons of top-selected short-list proteins to human proteins**.
(A) Plastin 1
(B) Coronin 1a
(C) Beta-1,3-N-acetylgalactosaminyltransferase 2
Tan proteins are human structures; blue proteins are predicted *S. rosetta* homolog structures. “Trait distance” is a multivariate Mahalanobis distance between pairs of proteins calculated based on 10 physicochemical properties (see the “Methods” section of our previous pub for more details [](https://doi.org/10.57844/ARCADIA-33B4-4DC5)). “Portfolio rank” is the relative rank of this protein compared to all other proteins found in the organism selection/the Zoogle portal. TM-score is a standard measure of solid-body protein structural similarity [](https://doi.org/10.1093/nar/gki524). In general, TM-scores of 0.5 or above indicate the same fold.
:::::
::::::
### *PLS1* (plastin 1) {#h.ishff7lggu98}
Plastin 1, also known as fimbrin, is an actin-bundling protein implicated in autosomal-dominant hearing loss [](https://doi.org/10.1002/humu.23891). Missense mutations in this protein lead to destabilization of plastin binding to actin and progressive or non-progressive hearing loss with variable age of onset [](https://doi.org/10.1002/humu.23891). While hearing loss in humans is reported as a missense allele, mouse knockouts have shown comparable modeling of human disease. *Pls1* knockout mice show defects in the morphology of stereocilia on ear hair cells [](https://doi.org/10.1093/hmg/ddu417), which are the mechanosensitive cells in the cochlea responsible for transducing sound into an electrical signal. Knockout mice also show an increasing loss of hearing as they age [](https://doi.org/10.1093/hmg/ddu417).
Choanoflagellate cells normally exhibit a stereociliar structure, the collar, which resembles ear hair cell stereocilia [](https://doi.org/10.1007/s10071-023-01776-z). This structure is involved in feeding and possibly in mechanosensation [](https://doi.org/10.1371/journal.pbio.3000226). The choanoflagellate plastin homolog shows high structural similarity (TM-score = 0.9) to plastin 1 and slightly lower structural similarity to two other major human homologs of plastin, LCP1 (also known as plastin 2, TM-score = 0.88) and plastin 3 (TM-score = 0.85) ([Figure 3](#ZS1)) (download Foldseek results from [Zenodo](https://zenodo.org/records/15724261)). The choanoflagellate gene also appears to be expressed in all four cell states in RNA-seq data, with particularly high expression in fast swimmers ([Figure 4](#Z07)).
Given the high degree of structural similarity between the human and choanoflagellate plastins, we were interested in generating knockouts of choanoflagellate plastin and assaying for possible phenotypes. Modeling plastin function in choanoflagellates also poses certain advantages over modeling the same function in mice or in differentiated human stem cells. Protocols for differentiating iPSC-derived hair-like cells take 50–60 days to complete [](https://doi.org/10.1073/pnas.2334503100), while assaying for strong phenotypes in mice can take months, due to an increase in hearing loss signal as the animals age [](https://doi.org/10.1093/hmg/ddu417). In contrast, choanoflagellate cells exhibit collar morphologies at all life stages. Generating knockouts, identifying possible phenotypes, and screening for molecules that rescue plastin function can all be performed with greater speed in choanoflagellates than in existing models. For these reasons, we were particularly excited to model *PLS1* function using *S. rosetta*.
::::::figure{#Z07 type=img src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/asset_055aa353.png' align=center}
:::::image{src='https://thestacks-01.s3.amazonaws.com/publications/resource-modeling-human-monogenic-diseases-salpingoeca-rosetta/asset_055aa353.png' width=100% alt='Illustration of different cell states of S. rosetta, including slow swimmers, rosettes, fast swimmers, and thecate cells. Bulk RNA-seq expression of pls1, coro1a, and b4galt7 across four cell states shown as scatter plots, illustrating that different cell types express different levels of all three genes.'}
:::::
:::::figcaption
**Figure 4:** **Bulk gene expression data for *S. rosetta***.
(A) Illustration of different cell states of *S. rosetta*, including slow swimmers, rosettes, fast swimmers, and thecate cells. Colors used for different cell states are reflected across the remaining charts.
(B–D) Bulk RNA-seq expression of *pls1* (B), *coro1a* (C), and *b4galt7* (D) across four cell states.
(B′–D′) Matrix showing degrees of significance of differential expression between pairs of cell states for *pls1* (B′), *coro1a* (C′), and *b4galt7* (D′).
:::::
::::::
### *CORO1A* (coronin 1a) {#h.pguap761c4u}
Coronin 1a is associated with severe combined immunodeficiency, also known as immunodeficiency 8 (IMD8). Patients with loss of coronin have early-childhood onset of recurrent infections, often associated with Epstein-Barr virus [](https://doi.org/10.1016/j.jaci.2013.01.042). Mice with missense mutations in *Coro1a* show defects in actin localization in the leading edge of T cells [](https://doi.org/10.1038/ni.1662).
Choanoflagellate cells produce a variety of actin-mediated cell structures, such as ruffles, filopodia, exocytic cups, and others [](https://doi.org/10.1371/journal.pbio.3000226)[](https://doi.org/10.1016/j.ydbio.2011.06.003)[](https://doi.org/10.1371/journal.pone.0095577)[](https://doi.org/10.7554/elife.61037). The choanoflagellate homolog of coronin 1a also has a high structural similarity to multiple coronins (the top seven hits in our Foldseek search included CORO6, CORO1C, CORO2B, CORO1A, CORO1B, CORO2A, and CORO7) ([Figure 3](#ZS1)) (download Foldseek results from [Zenodo](https://zenodo.org/records/15724261)). This gene is also expressed across cell states with three different expression classes\: slow swimmers and rosettes, fast swimmers, and thecate cells ([Figure 4](#Z07)). Due to the high number of different phenotypes that could be captured in *S. rosetta*, we felt confident that we could identify some defect in cells with a mutant version of this gene.
### *B3GALNT2* (beta-1,3-N-acetylgalactosaminyltransferase 2) {#h.5o29v8ifqkab}
Missense, nonsense, and frameshift mutations in this glycosylation enzyme are associated with a form of muscular dystrophy caused by congenital brain and eye abnormalities [](https://doi.org/10.1055/s-0038-1651519)[](https://doi.org/10.1016/j.ajhg.2013.01.016). Morpholino knockdown of the zebrafish homolog of beta-1,3-N-acetylgalactosaminyltransferase 2 results in severe morphological defects, including those of the eye, brain, and spinal cord [](https://doi.org/10.1016/j.ajhg.2013.01.016). Mutations in B3GALNT2 also appear to cause hydrocephalus in horses [](https://doi.org/10.1186/s12864-015-1936-z).
The *S. rosetta* homolog of beta-1,3-N-acetylgalactosaminyltransferase 2 is a strong structural match (TM-score = 0.79) to the human protein ([Figure 3](#ZS1)) and is expressed at low levels in all cells and slightly higher levels in thecate cells. ([Figure 4](#Z07)).
Glycosylation plays an important role in the proper development of choanoflagellate rosettes. Recent work has shown that loss of function of glycosylation enzymes can lead to “clumpy” choanoflagellate cells that stick together [](https://doi.org/10.7554/elife.41482). This can be assayed by mixing non-fluorescent and fluorescent choanoflagellate cells and looking for aggregates across strains. We expect that loss of function in this particular glycosylation enzyme may produce similar phenotypes, which can be easily measured using this assay. For these reasons, we decided to pursue this gene.
### Other hits {#h.3uqhex8bxoee}
Several other genes — *CLTC*, *UNC13D*, *GALC*, and *B3GALNT2* — also made our short list ([Table 1](#short-list-table)).
::::::figure{#short-list-table type=table}
| **HGNC gene symbol**
**(UniProt ID)** | **Human protein name** | **Associated human disease (OMIM)** | ***S. rosetta* protein UniProt ID** | ***S. rosetta* identifier** | **Possible *S. rosetta* phenotypes** |
| --- | --- | --- | --- | --- | --- |
| *PLS1*
([Q14651](https://www.uniprot.org/uniprotkb/Q14651/entry)) | Plastin 1 | [Deafness, autosomal dominant 76](https://omim.org/entry/618787) | [F2TWP3](https://www.uniprot.org/uniprotkb/F2TWP3/entry) | PTSG_00515 | Effects on collar morphology |
| *CORO1A*
([P31146](https://www.uniprot.org/uniprotkb/P31146/entry)) | Coronin 1a | [Immunodeficiency 8 with lymphoproliferation](https://omim.org/entry/615401) | [F2TZV2](https://www.uniprot.org/uniprotkb/F2TZV2/entry) | PTSG_01270 | Effects on membrane dynamics |
| *B4GALT7* ([Q9UBV7](https://www.uniprot.org/uniprotkb/Q9UBV7/entry)) | Beta-1,4-galactosyltransferase 7 | [Ehlers-Danlos syndrome, spondylodysplastic type, 1](https://omim.org/entry/130070) | [F2UGX2](https://www.uniprot.org/uniprotkb/F2UGX2/entry) | PTSG_07984 | Effects on glycosylation, cell clumping |
| *CLTC* ([Q00610](https://www.uniprot.org/uniprotkb/Q00610/entry)) | Clathrin heavy chain | [Intellectual developmental disorder, autosomal dominant 56](https://omim.org/entry/617854) | [F2U3P4](https://www.uniprot.org/uniprotkb/F2U3P4/entry) | PTSG_02909 | Effects on membrane dynamics |
| *UNC13D*
([Q70J99](https://www.uniprot.org/uniprotkb/Q70J99/entry)) | Unc-13 homolog D | [Hemophagocytic lymphohistiocytosis, familial, 3](https://omim.org/entry/608898) | [F2U085](https://www.uniprot.org/uniprotkb/F2U085/entry) | PTSG_11723 | Effects on exocytosis |
| *GALC* ([P54803](https://www.uniprot.org/uniprotkb/P54803/entry)) | Galactosylceramidase | [Krabbe disease](https://omim.org/entry/245200) | [F2TVZ9](https://www.uniprot.org/uniprotkb/F2TVZ9/entry) | PTSG_00266 | Effects on glycosylation, lysosomes |
| *B3GALNT2* ([Q8NCR0](https://www.uniprot.org/uniprotkb/Q8NCR0/entry)) | Beta-1,3-N-acetylgalactosaminyltransferase 2 | [Muscular dystrophy-dystroglycanopathy (congenital with brain and eye anomalies), type A, 11](https://omim.org/entry/615181) | [F2UC26](https://www.uniprot.org/uniprotkb/F2UC26/entry) | PTSG_06143 | Effects on glycosylation, cell clumping |
:::::figcaption
**Table 1:** **Summary of the relevant diseases and our initial guess for phenotypes for each short-list gene**.
[HGNC](https://www.genenames.org/)\: HUGO gene nomenclature committee.
:::::
::::::
## Expert feedback and our final short list {#h.4krhg51b0it3}
We sought expert feedback on our short list from David Booth, a choanoflagellate researcher at UCSF, one of the external scientists we consulted as part of this work. While we didn’t have strong phenotypic hypotheses for all seven genes on our short list, given the ease of performing simple phenotypic assays such as growth, morphology, cell state transitions, and cell clumping, we decided to pursue all of our hits. Below is a rough overview of our projected experimental plan in *S. rosetta*.
### Experimental plan {#h.x190amj87qn}
1. **Generate mutants and characterize growth phenotypes**.
2. **Phenotype cell state transitions**. For example, the ability to transition from slow swimmers to rosettes.
3. **Investigate other cellular phenotypes**. For example, cell motility or morphology. Such phenotypes are interesting because they are high-dimensional, quantitative, and easy to assay in high throughput. We can then follow up to perform downstream, finer-grained analyses using lower-throughput phenotyping methods.
4. **Perform genetic rescue experiments using complementation with wild-type protein or synonymous revertants**.
# Methods {#h.90mqlaow6uc2}
For details about the technical analyses described in this pub, see the [companion pub describing our evaluation framework](https://research.arcadiascience.com/pub/resource-modeling-human-monogenic-diseases).
::::::div{.info-box}
All **code** related to this pub is available in [this GitHub repo](https://github.com/Arcadia-Science/2025-zoogle-experiments/tree/v1.0) (DOI\: [10.5281/zenodo.15707938](https://doi.org/10.5281/zenodo.15707938)). Code specific to *S. rosetta* is in [this analysis notebook](https://github.com/Arcadia-Science/2025-zoogle-experiments/blob/main/notebooks/Salpingoeca-rosetta/2_salpingoeca-rosetta-expression.ipynb) and [these visualization scripts](https://github.com/Arcadia-Science/2025-zoogle-experiments/tree/main/zoogletools/salpingoeca).
::::::
## *Salpingoeca* bulk cell type gene expression {#h.8fdnelwk91z4}
We downloaded the *S. rosetta* cell type-specific RNA-seq data from [this FigShare deposition](https://figshare.com/articles/dataset/Table_of_differential_gene_expression_between_cell_types_of_i_Salpingoeca_rosetta_i_/28225106?file=51726065) and used a [custom script](https://github.com/Arcadia-Science/2025-zoogle-experiments/blob/480f329c1424c86c8a2c10ec0bb43d7329d92ff0/notebooks/Salpingoeca-rosetta/2_salpingoeca-rosetta-expression.ipynb) found in our GitHub repository to visualize gene expression counts from bulk data for four cell types\: slow swimmers, fast swimmers, thecate cells, and rosettes. We also visualized the computed pairwise differential expression significance *p*-values, such as in [Figure 4](#Z07).
## AI tool usage {#h.3zvz9jt468vp}
We used Claude to help write code, clean up code, comment our code, draft text that we edited, suggest wording ideas for small phrases and sentence structures, and clarify and streamline text that we wrote. We also used Cursor to access Claude.
## Visualization {#h.ca8vu9ce951z}
We used plotly (v5.17.0) [](https://plot.ly) arcadia-pycolor (v0.6.3) [](https://github.com/Arcadia-Science/arcadia-pycolor) to generate figures before manual adjustment.
# Key takeaways {#h.tj9oadwf3l0c}
## Advantages of *S. rosetta* as a disease model {#h.92gdi0k9ves7}
During our diligence process, we came to realize a variety of advantages that make *S. rosetta* a valuable potential model for human diseases, including\:
- **Multiple cell types**. In contrast to many microbial eukaryotes, which largely have a single cell type, *S. rosetta* has the ability to “differentiate” into different states, each of which may have different phenotypes. This allows for a broader variety of phenotypes that may be identified in mutants.
- **Persistent and dynamic cytoskeletal structures**. *S. rosetta* displays stereocilia and flagella in a variety of cell states. These structures aren’t readily found in all human cell types, requiring potentially expensive differentiation protocols to obtain cells with relevant structures. *S. rosetta* also displays some transient structures, such as filopodia and exocytic cups, which could expand the search space of possible phenotypes.
- **Ease of genetic experimentation**. The state of genetic manipulations in *S. rosetta* is more robust and precise than in many other microbial eukaryotes. The lack of complex multinuclear structures, such as those found in organisms like *Stentor coerulus* or *Tetrahymena thermophila*, means that genetic manipulation is more likely to result in an observable phenotype.
# Next steps {#h.r9ndqas11tdd}
Our diligence process uncovered just seven of many possible disease-causing genes that researchers can feasibly model using *S. rosetta*. Moreover, we sought to identify “low-hanging fruit” among our candidates. Scientists with greater experience with particular types of assays or a willingness to generate targeted mutations might be able to pursue disease–gene pairs that we rejected in our analysis. We encourage others to review our [working diligence document](https://docs.google.com/spreadsheets/d/1DeJ9X1YGajwvR6zB6bz0RCPz4ClbSYUQDk4jjKj2P2I/edit?usp=sharing) and see whether their evaluation of feasibility differs from ours.
We’re funding downstream work by researchers in the [David Booth lab](https://www.biobooth.org/) at UCSF. They’ll share their results openly during their experiments. Stay tuned for more!