VcMid Controls Spermatogenesis

We tested the role of VcMid protein in sexual differentiation by generating female transgenic lines with an autosomally integratedVcMID transgene (pVcMID-BH) expressed under its own promoter and fused to a blue fluorescent protein (BFP) and a hemagglutinin (HA) epitope tag at its C-terminus to detect expression (Eve::VcMID-BH) (Figures 2A andS1B). Female transformants carrying an untagged version ofVcMID (Figure S1C) had identical phenotypes as those carrying the tagged version, and all subsequent work was done with tagged strains and untagged transformants as a negative control for detection of VcMid protein.Eve::VcMID-BH lines showed a normal vegetative phenotype and constitutively expressed the mRNA for theVcMID-BH transgene, an expression pattern identical to the endogenousVcMID mRNA in males (Figure S2)[20]. As described above, when wild-type vegetative female gonidia are exposed to sex inducer they undergo modified embryogenesis and develop into sexual spheroids with 32–48 eggs and ∼2,000 sexual somatic cells (Figure 2B). When wild-type vegetative male gonidia are exposed to sex inducer they undergo modified embryogenesis and develop with 128 sperm packets and 128 somatic cells (Figure 2C). When vegetative gonidia fromEve::VcMID-BH lines were exposed to sex inducer they developed into progeny spheroids with a novel pseudo-male sexual phenotype: They produced ∼2,000 somatic cells and 32–48 sexual germ-cell precursors in a pattern similar to that of female eggs; but each of the 32–48 germ-cell precursors in theEve::VcMID-BH lines underwent additional cleavage divisions like male androgonidia to produce sperm packets (Figure 2D). Moreover, the sperm produced inEve::VcMID-BH lines were capable of fertilizing wild-type female eggs to produce characteristically orange-pigmented, thick-walled zygotes (Figure 2E). However, male function was incomplete as fertility defects were noted (see next section). Similar to wild-type sperm (and unlike wild-type female eggs) theEve::VcMID-BH sperm cells were terminally differentiated and could not revert back to vegetative growth if left unfertilized. Another male-specific phenotype exhibited byEve::VcMID-BH lines was frequent spontaneous occurrence of sexual differentiation in vegetative cultures[36], a trait whose underlying basis is not clear, but which appears to be under the control ofVcMID.

Fertility Defects in Females Expressing VcMid

Although some of theEve::VcMID-BH sperm were functional and could fertilize wild-type eggs, the sperm packets and sperm cells from these strains had multiple defects including heterochronic delays in maturation and hatching defects (Figure S3A). The sperm packets were four times larger than those from wild-type males and contained about four times as many sperm cells (256 sperm/packet) (Figure 2F and 2G), some of which were aberrantly formed in contrast with wild-type sperm cells that had uniform morphology (Figures 2H, 2I, andS3B–S3F). Nonetheless the crosses betweenEve::VcMID-BH pseudo-males and wild-type females producedMTF/MTF diploid zygotes (Figure 2E) that could germinate and produce haploid progeny. Forty-one progeny from one such cross were genotyped, half of which (19/41) inherited theVcMID-BH transgene and developed as pseudo-males, and half of which lacked the transgene and developed as normal females (22/41).

Expression and Localization of VcMid Are under Post-transcriptional Control

The absence of a vegetative phenotype inEve::VcMID-BH transgenic lines despite constitutive expression of theVcMID-BH mRNA (Figure S2) suggested that VcMid protein expression or localization might be under posttranscriptional control. Although we could not detect BFP fluorescence inEve::VcMID-BH strains, we could detect the HA epitope tag by Western blotting at all stages of vegetative and sexual development (Figures 3A andS4). Using immunofluorescence (IF) a nuclear-localized signal was detected for VcMid-BH protein in cleaving androgonidia and in mature sperm packets (Figures 3B–3J,S5A–S5R, andS6C–S6H), a result similar to earlier findings of Mid protein localization in sperm nuclei ofPleodorina[21]. However, nuclear VcMid-BH was not detected during early stages of sexual male embryogenesis prior to androgonidia cleavage (Figure S7).

We also examined VcMid-BH expression and localization in gonidia and somatic cells from vegetative spheroids. RNA was prepared from purified gonidial or somatic cells and reverse transcription and PCR (RT-PCR) detectedVcMID-BH mRNA at similar levels in both cell types from transgenic females (Eve::VcMID-BH) and males (AichM::VcMID-BH) (Figure S8B and S8D). The endogenousVcMID transcript from wild-type males was also expressed in both vegetative cell types (Figure S8C and S8D). Whole cell extracts were prepared from purifiedEve::VcMID-BH orAichM::VcMID-BH somatic cells and gonidia and subjected to SDS-PAGE and Western blotting to detect VcMid-BH protein (Figures 3K,S9A, and S9B). In contrast toVcMID-BH mRNA that was present in both vegetative cell types, VcMid-BH protein was only detected in vegetative somatic cells indicating that there is cell-type–specific regulation of VcMid protein synthesis or stability that restricts its accumulation to somatic cells during vegetative growth (Figures 3K andS6B). However, unlike the case for androgonidia and sperm cells, the VcMid-BH protein signal in vegetative somatic cells was excluded from the nucleus and was instead observed only in the cytosol and peri-nuclear region (Figures 3L–3O,S6I–S6P, andS9C–S9J). Together these data suggest that cell-type–limited expression and regulated nuclear localization of VcMid restrict its function to developing and mature sexual male germ cells inV. carteri.

VcMid Is Required for Spermatogenesis and Represses Oogenesis in Males

In order to test whetherVcMID is necessary for sexual differentiation of males we developed a new strategy for gene knockdown on the basis of RNAi-inducing hairpin constructs. The hairpin-forming portion of the construct corresponding toVcMID sequences was inserted directly into the 3′ UTR of thenitA selectable marker gene to allow direct and simultaneous selection of both NitA+ and hairpin expression (Figures 4A andS10;Text S1). This strategy has been successful for other loci besidesVcMID (unpublished data), but only the results forVcMID are presented here. We also note that theVcMID knockdown phenotype was gene specific and did not occur when hairpins targeting other loci were introduced intoV. carteri (unpublished data). Two hairpins targetingVcMID (VcMID-hp1 andVcMID-hp2) were introduced into aV. carteri wild-type male strainAichiM to generateAichiM::VcMID-hp1 andAichiM::VcMID-hp2 transgenic lines (seeMaterials and Methods for details). All transgenic male strains had normal vegetative phenotypes, and both hairpin constructs reducedVcMID expression; butAichiM::VcMID-hp1 lines had lowerVcMID transcript levels thanAichiM::VcMID-hp2 lines (Figure S2). We note that unlike wild-type males or pseudomales (see above), vegetative cultures ofAichiM::VcMID-hp1 lines did not undergo spontaneous sexual induction. Sexually inducedAichiM::VcMID-hp1 lines with strong knockdowns showed a novel phenotype: Their early sexual development proceeded as it would for a wild-type male strain and resulted in spheroids that contained 128 small somatic cells and 128 large cells that resembled uncleaved androgonidia (Figure 4B), but the large cells never underwent further cleavage into sperm packets. Instead, many of them could be successfully fertilized with wild-type male sperm to makeMTM/MTM diploid zygospores (Figure 4C). The ability to differentiate as zygospores when fertilized indicates that the presumptive androgonidia inAichiM::VcMID-hp1 strains were converted to functional eggs and that these strains were behaving as pseudo-females. However, unlike normal zygotes from a wild-type cross, 30%–50% of the zygotes fromAichiM×AichiM::VcMID-hp1 pseudo-female crosses died and bleached shortly after fertilization (Figure 4C), a phenotype that depended on addition of exogenous sperm. The surviving zygotes from these crosses produced some viable meiotic progeny, but germination and survival of the progeny were reduced compared with normal wild-type zygotes (Table S1). 29/43 viable progeny inherited theVcMID-hp1 transgene and developed as pseudo-females, while 14/43 lacked the transgene and developed as normal males. The apparent deviation from a 1∶1 inheritance pattern of the transgene was noted but was not pursued further in this study. The high mortality of pseudo-female eggs—whose mating loci are genetically male—suggest thatMTF contains genes that promote female gamete and/or zygote fitness that are absent fromMTM. If left unfertilized, the eggs fromAichiM::VcMID-hp1 lines could de-differentiate and reenter the vegetative reproductive cycle as do unfertilized female eggs. A similar phenotype as our pseudo-male strain was reported previously for a male mutant[41], but the mutant strain is no longer available for characterization.

PartialVcMID Knockdown Generates Self-Fertile Hermaphrodites

AichiM::VcMID-hp2 lines were not as severely knocked down forVcMID expression asAichiM::VcMID-hp1 lines (Figure S2), and had a distinct hermaphrodite phenotype in which sexual spheroids developed with a mixture of normal-looking male sperm packets and pseudo-female eggs (Figure 4D). The hermaphrodite lines exhibited self-fertility as evidenced by zygospores that formed in sexually induced monocultures (Figure 4E). These results indicate thatV. carteri sex determination is highly sensitive to VcMid dosage where either a male or female fate is established depending on the level of VcMid.

Volvox andChlamydomonas Mid Proteins Are Functionally Distinct

In several instances genes fromC. reinhardtii have been shown to function interchangeably with theirV. carteri orthologs in developmental processes such as inversion and asymmetric cell division (reviewed in[42],[43]). Mid proteins have at least two domains: The C-terminal region has a predicted RWP-RK motif DNA binding domain, while the N-terminal region does not show similarity to characterized protein domains from other organisms (Figure S1)[28]. We tested whether either of the domains from CrMid could substitute for those of VcMid to control spermatogenesis inV. carteri. To do so we generated three constructs in which all or part of theCrMID genomic coding region was substituted forVcMID sequences in pVcMID-BH. One construct contained the entireCrMID gene (pCrMID-BH) (Figure 5A) while the other two contained theCrMID N-terminal domain fused to theVcMID DNA binding domain (pMID-V_NC_C-BH) (Figure 5B), or theVcMID N-terminal domain fused to theCrMID DNA binding domain (pMID-C_NV_C-BH) (Figure 5C). All three constructs as well as pVcMID-BH were introduced into wild-type females, and transformants that expressed the different predicted proteins were identified (Figures 5D–5F). UnlikeEve::VcMID-BH transformants that showed a pseudo-male phenotype (Figure 1D), none of the transformants that expressed CrMid or Mid chimeras had this phenotype, but instead always developed as wild-type females (Figure 5G;Table S2). The CrMid-BH protein was expressed as a full-length form and two shorter isoforms (Figure 5D), possibly due to proteolytic cleavage or incorrect pre-mRNA processing. However, lines with either of the two chimeric constructs expressed only full-length predicted proteins at levels comparable to VcMid-BH (Figure 5E and 5F). We conclude from these experiments that CrMid and VcMid are not functionally interchangeable, and that both the N-terminal and RWP-RK domains of VcMid are required to activate spermatogenesis in sexual germ cells.

Deep Homology between Mating Types and Sexes

Although it was reasonable to predict that the route for evolving sexual dimorphism would be through addition of new genetic functions and pathways to a core mating locus as proposed originally by Charlesworth for the evolution of anisogamy[11], we found instead that a single conserved volvocine algal mating locus gene,MID, is largely responsible for controlling male versus female sexual differentiation inV. carteri. It is notable thatVcMID controls multiple sexual traits inV. carteri that have no analogs inC. reinhardtii. The male traits controlled byVcMID include specialized cleavage divisions of androgonidia, sperm packet inversion, sperm packet hatching, specialized sperm cell morphology, and gamete recognition without flagellar adhesion.

The term deep homology is used to describe ancient and conserved genetic mechanisms that control traits that, on the surface, appear disparate or have no obvious homology relationship[44]. Metazoan eye development is a classic example; across phyla with very different eye architecture, it is controlled by the conserved transcription factor, eyeless/Pax6[45]. In the case of volvocine algae Mid proteins control two very different manifestations of sexual reproduction inC. reinhardtii andV. carteri whoseMID orthologs have been diverging for as long as 200 million years[46]. Sexually dimorphic gametes have evolved from isogamous ancestors several times in independent multicellular taxa, and the mechanism could be similar to, or different from, that in volvocine algae where a master regulatory gene acquired the ability to direct sperm-egg dimorphism. Oogamy was likely established very early in the Streptophyte lineage before the split between Charophyte algae and land plants, but the origins and bases of oogamy in Charophytes remain unclear[47]. The recent identification of an anisogamous sexual cycle in a Choanoflagellate—a unicellular relative of metazoans—introduces the potential for investigating the early evolution of gamete size dimorphism in animals[48].

MID Is a Master Regulator of Sex Determination in Volvocine Algae

Our results demonstrate for the first time a function for Mid protein in a volvocine algal sexual cycle outside of the genusChlamydomonas and suggest that Mid could be the master regulator of mating type, gamete size, and gender throughout the volvocine lineage whereMID genes have been identified in most genera[21],[28],[30]. AlthoughMID sequences were previously shown to evolve rapidly, the finding that Mid protein fromC. incerta (now reclassified asC. globosa[49]) can substitute forC. reinhardtii Mid indicates that functional conservation can be retained after speciation[27]. However,C. reinhardtii Mid protein could not substitute for theV. carteri ortholog, which appears to require both its native DNA binding and N-terminal domains to function in sex determination (Figure 5). Future work using cross-species complementation will help clarify whether Mid protein function co-evolved with sexual dimorphism in volvocine algae as our data suggest might be the case.

Transcriptional regulatory network evolution has been studied in various developmental contexts[50]–[53], but very little is known about how regulatory networks are modified or coopted during unicellular to multicellular transitions[54]. The Mid system in volvocine algae represents a new opportunity for understanding how a cell-type specification pathway in a unicellular ancestor evolved to control a complex developmental program in a multicellular descendant. While mating-type differentiation inChlamydomonas appears to involve differential expression of a small number of genes between theplus andminus gametes[14],[22], a larger set of differentially expressed genes might be expected to specify sperm and eggs, which are developmentally very different from each other[37],[55]. An important future goal will be to identify and compare the direct and indirect targets of Mid proteins in bothC. reinhardtii andV. carteri that are predicted to be more numerous and diverse inV. carteri.

Another interesting question is whetherMID-like genes function in sex determination in green algae outside of the volvocine lineage. The molecular bases for sex determination in most green algae and protists are poorly understood, but RWP-RK family proteins are found throughout the green eukaryotic lineage[56], including small Mid-like proteins in Prasinophyte algae[57], and even in distantly related Cryptophyte algae[58]. It remains to be seen whether these Mid-like proteins in non-volvocine species function in sex determination.

Evolution of Mid and Sexual Cycles in the Volvocine Lineage

As described in theIntroduction, volvocine algae exhibit wide diversity in their sexual cycles: There are isogamous, anisogamous, and oogamous species with sexual cycles that can be heterothallic or homothallic. Among homothallic species some are dioecious (producing a mixture of all male and all female sexual offspring) or monoecious (producing sexual offspring containing gametes of both sexes in one individual)[16],[55].

V. carteri is a heterothallic species with genetically determined sexes; but, a remarkable phenotype was produced by a partial RNAi knockdown ofVcMID in males (Figure 4D). Rather than developing as male or female, the partial-knockdown male spheroids became self-fertile monoecious hermaphrodites that produced sperm and eggs in a single individual. On the basis of this observation, it can be inferred that sexual differentiation inV. carteri is bi-stable and highly sensitive to initialMID dosage. Once the Mid-sensitive step of development is initiated (which may coincide with nuclear translocation of VcMid protein), positive and negative feedback loops may be used to lock the sex determination program into a male or female state. This state could be achieved either independently of VcMid concentration, or through positive/negative reinforcement of initial expression states. The phenotype of hermaphroditic development by partialVcMID knockdown may be relevant to the evolution of homothallism, which appears to have arisen in all three major clades ofVolvox[16],[55],[59]. We speculate that naturally evolved homothallic volvocine algae possess aMID gene whose expression is insufficient to specify 100% male gamete production—much like ourVcMID-hp2 strains. Moreover, the timing ofMID expression in homothallic species ofVolvox could play a role in determining monoecious versus dioecious reproductive development: If the Mid-sensitive developmental switch is triggered relatively late in development after germ-cell precursors are formed then a mixture of male and female gametes could develop within a single spheroid as we found withVcMID-hp2 strains (i.e., homothallism, monoecy). In contrast, if the Mid-sensitive step occurred very early in development before individual germ cells were established, then the fate of all germ cells within a mature spheroid might be locked into a male or female program and the resulting population would produce a mixture of all-female or all-male spheroids (i.e., homothallism, dioecy).

The bi-stability of sex determination at intermediate levels ofMID expression inV. carteri is also reminiscent of theiso1 mutant phenotype inC. reinhardtii whereMT− iso1 cells differentiate into a mixture ofplus andminus gametes that iso-agglutinate but cannot self-fertilize[60]. The self-infertility ofiso1 MT− strains is due to the absence ofFUS1, a gene from theMT+ haplotype that is required for gamete fusion[61]. In contrast, there are no essential genes inMTF ofV. carteri that are absolutely required for fertilization and subsequent germination of progeny from matings between pseudo-females (VcMID knockdown strains) and males. This lack of essential femaleMT genes for completing the sexual cycle may have facilitated transitions from heterothallism to homothallism in volvocine algae.

V. carteri Has Evolved New Regulatory Inputs for Mid

In volvocine genera other thanVolvox—including the anisogamous genusPleodorina— sexual differentiation andMID expression are both triggered by the absence of nitrogen (−N)[21],[25],[28],[62]. In contrast,V. carteri f. nagariensis and otherVolvox species use species-specific pheromones called sex inducers to trigger sexual differentiation[16],[55],[63]. A seemingly parsimonious evolutionary route for rewiring input into the Mid pathway inV. carteri would simply placeVcMID transcription under the control of sex inducer instead of nitrogen availability. Unexpectedly, however,VcMID mRNA and VcMid protein are both expressed constitutively at all life cycle stages (Figures 3A andS4)[20], and unlike the case inC. reinhardtii, VcMid appears to be under at least three types of posttranscriptional control: (i) AlthoughVcMID mRNA is present in both vegetative cell types (somatic and gonidial cells), VcMid protein is only translated or stably produced in somatic cells and is absent from vegetative gonidia and vegetative embryos (Figures 3K andS6B). (ii) The VcMid protein produced in somatic cells is excluded from the nucleus (Figures 3L–3O,S6I–S6P, andS9C–S9J). (iii) In response to the presence of sex inducer, VcMid protein accumulates in the nuclei of cleaving androgonidial cells and in the nuclei of sperm cells where it is presumed to function in specifying sperm development (Figure 3B–3J,S5, andS6C–S6H).

In depth study will be required to determine how cell-type–regulated production and localization of VcMid are achieved, but it seems reasonable to infer that one or more factors are produced in sexually induced spheroids that promote the translation and/or stability of VcMid and its nuclear localization in androgonidia and sperm. The factor may interact directly with VcMid as a partner and may help specify its localization at promoters of target genes, but this idea remains to be tested. The absence of VcMid protein in vegetative gonidia despite its message accumulating to the same extent as in somatic cells seems puzzling at first glance. However, the block in stability or translation of VcMid in gonidia may have evolved as a failsafe mechanism to prevent sexual differentiation during vegetative embryogenesis. Such a mechanism would be unnecessary in vegetative somatic cells that are already terminally differentiated, but could potentially be important for vegetative gonidial cells because male sexual differentiation is irreversible and would be fatal if it occurred at the wrong time. However, why vegetative somatic cells express VcMid remains a mystery. As noted inResults, we observed no obvious vegetative phase phenotypes inAichiM::VcMID-hp strains whose somatic cells were missing VcMid, or inEve::VcMID-BH strains that contained VcMid in somatic cells, but a definitive conclusion about whether cytoplasmic VcMid has a role in vegetative somatic cells awaits more in depth examination.

Uncoupling of Gender from Sex Chromosome Identity Uncovers Possible Sexually Antagonistic Interactions in theMT Locus

Our results show that the presence or absence ofVcMID is the key determinant of differentiation inV. carteri sexual spheroids; yet theMT locus of this alga has around 70 additional genes, many of which show sex-biased gene expression[20]. Some of theMT genes are present only in the male or only in female haplotype, while the majority are male and female gametolog pairs that are highly diverged in sequence and expression pattern[20]. Our ability to uncouple gender from sex chromosome identity by manipulation ofVcMID expression allowed us to uncover potential contributions of male and femaleMT genes to sexual dimorphism and reproductive fitness. Haploid sex chromosome systems have received less attention than diploid systems, but there are several predictions about them that our results begin to address. Under haploid dioecy, recessive mutations in sex-linked genes are not sheltered from selection as they are in the heterogametic sex of diploid systems, and are therefore expected to degenerate equally and lose only genes required for the opposite sex and not their own[64]. We note, however, that inV. carteri MT there is no evidence for an allele of a gametolog pair having been eliminated from one mating haplotype and retained in the other[20]. Haploid sex chromosomes are predicted to be similar to diploid systems in that both should accumulate sexually antagonistic alleles that benefit one sex, but harm the other[39],[65], and should accumulate repeat sequences, as appears to have occurred in bothVolvox and in the bryophyteMarchantia[20],[66]. Accumulation of sexually antagonistic alleles has not been tested for haploid sex chromosomes, but the extensive divergence betweenV. carteri MT gametologs suggests that this phenomenon may contribute to the developmental and fitness defects we observed in pseudo-male and pseudo-female strains.

Although none of the sex-limitedMT genes besidesVcMID appear to be essential forV. carteri sex determination and completion of the sexual cycle, they clearly impact sexual development and reproductive fitness. A striking phenotype for both the pseudo-male and pseudo-female strains was the patterning of their germ-cell precursors formed during sexual development (Figures 2B–2D and6A–6C). It is clear from these phenotypes that sexual germ cell patterning (i.e., the number and distribution of germ-cell precursor cells and ratio of germ-cell precursors to sexual somatic cells) is separable from germ cell differentiation and is not controlled by the VcMid pathway. Instead, this sexual patterning trait must be controlled by otherMT genes (Figure 6). One candidate for this male-female patterning difference is theMAT3 gene that encodes the retinoblastoma-related homolog inV. carteri[20],[67]. InChlamydomonas, Mat3 protein controls the multiple fission cell cycle by establishing the threshold size at which division can occur and by coupling the extent of cell division to mother cell-size[68],[69]. Although not directly involved in determining gamete cell-size as we had originally predicted[40], it is possible that the male and female alleles ofVcMAT3 dictate the timing of asymmetric cell divisions by coupling embryonic blastomere cell-size to the asymmetric division machinery. Future work will be aimed towards determining the role of male and female VcMat3 gametologs in sexual cell division or other parts of theV. carteri sexual cycle.

In addition to defects in germ cell patterning inVolvox pseudo-males and pseudo-females, these strains show other reproductive defects. These include abnormal sperm cell shape and morphology (Figures 2H, 2I, andS3B–S3F), low efficiency of sperm packet hatching (Figure S3D), and delayed timing of androgonidial cleavage into sperm packets (Figure S3A). In depth study may reveal other defects in pseudo-male sperm related to cytoskeletal organization, gamete recognition, motility, and fertilization dynamics. Although egg cells lack distinct morphological features like sperm, we noted a very high mortality rate in pseudo-female eggs that occurred when we attempted fertilization. The mortality we observed under these conditions could be due to pre-zygotic defects in the eggs caused by their smaller size than wild-type female eggs or by a mismatch between the male mating locus and the female sexual differentiation pathway. The mortality might also be due to zygotic defects that occur when a copy of the female mating locus is absent from the zygote immediately after fertilization. Future work aimed at developing quantitative assays for distinct steps of fertilization will allow us to document in more detail the fitness contributions of male and femaleMT genes to the sexual cycle.

Volvox Strains and Culture Conditions

Eve (Volvox carteri. f. nagariensis UTEX 1885) andAichiM (Volvox carteri. f. nagariensis NIES 398) were obtained from stock centershttp://web.biosci.utexas.edu/utex/ andhttp://mcc.nies.go.jp/, respectively. The strains that were used for transformations are described below and inText S1. Other strains are described inTable S3. Growth of strains was in standard Volvox medium (SVM) or urea-free standard Volvox medium (UF-SVM) with growth conditions described in more detail inText S1.

Plasmid Construction

All plasmid constructs were made using standard molecular cloning methods and manipulations[70]. PCR amplifications used for plasmid construction were done with Phusion Polymerase (Thermo Scientific) according to the manufacturer's guidelines (see also PCR amplification conditions below). Primers used in this study are described inTable S4.

Nuclear Transformation ofV. carteri and Generation of Transgenic Strains

Constructs were introduced intoNitA− female strainE15 or male strainA18. Transformation or co-transformation ofE15 orA18 with nitrate-reductase (NitA) encoding plasmid pVcNR15[71], pVcMID hp1, or pVcMID hp2 was done as previously reported[72] with minor modifications described inText S1. Transformed gonidia cells were selected in UF-SVM.

Induction of Sexual Development and Phenotypic Scoring

For each assay three vegetative juvenile-stage spheroids were placed in a well of a six-well microtiter plate with ∼9 ml SVM per well and 10 µl of sex inducer with a titer of 10⁶, and maintained at 32°C in a 16 h∶8 h light∶dark cycle[34]. The phenotypes of the 40–50 sexual progeny spheroids that resulted from sexual development of gonidia in the three starting spheroids were scored visually under a dissecting microscope and documented using a compound light microscope (Leica DMI6000B, 40× objective, differential interference contrast [DIC] optics) after 3–7 days.

Mating

Mating and zygote germination were performed as previously described with minor modifications[73]. Parental strains were grown to a density of ∼300 unhatched juveniles in 350 ml of SVM at which point sex-inducer was added. In the subsequent cleavage cycle sexual spheroids were produced and allowed to mature. Egg-bearing females were released from their parental spheroid by gentle pipetting 3–5 hours prior to hatching and mixed with males that had their sperm packets released from their vesicles within the parental spheroid by a similar procedure. Matings took place in a glass 150 mm×25 mm petri dish at a density of 10–15 spheroids/ml. The petri dish was placed on a light box within a 30°C growth chamber for 8–16 hours, and then the dishes were wrapped with aluminum foil and left at 32°C for at least one week and typically for three weeks prior to germination.

Zygote Germination

Drawn-out Pasteur pipets were used to manipulate zygotes. In order to remove residual sex-inducer and any other potential inhibitors of germination, zygotes were put into 1.5 ml tubes with ∼1 ml SVM medium and washed as follows: Tubes were vortexed for 10 min, then spun briefly at ∼1,000 rpm in a microcentrifuge to pellet zygotes. The supernatant was removed and the washing step was repeated three times. The washed zygotes were transferred to a sterile glass depression slide, washed briefly with SVM, and allowed to settle to the bottom. Ten to 50 zygotes were transferred to a sterile glass depression well containing SVM with 60 µg/ml carbenicillin, and incubated in a 16 h∶8 h light∶dark 30°C growth chamber. Zygotes germinated after two to six days. The germling colonies were individually transferred to six-well microtiter plates for growth and clonal expansion.

RNA Preparation, cDNA Preparation, and Semi-quantitative RT-PCR

TotalVolvox RNA was prepared as previously described[20]. cDNA was prepared from 5 µg total RNA following the manufacturer's protocol for Thermoscript (Invitrogen) using a 10∶1 mixture of oligo dT and random hexamer for priming and the following cDNA synthesis reaction temperatures: 25°C 10′, 42°C 10′, 50°C 20′, 55°C 20′, 60°C 20′, 85°C 5′, after which the reactions were treated with RNaseH. Reactions were diluted 1∶10 with 10 mM Tris pH 8.0, 1 mM EDTA (TE), and stored at −20°C.S18 andVcMID expression was measured using amplification with primer setsVcMid.f1 andVcMID.r1 and S18.1 and S18.2 as described previously[20] and inTable S4.

Cell Separation

Vegetative gonidia and somatic cell separation was done by mechanical disruption and differential centrifugation. Cultures grown in four 350 ml SVM flasks with ∼5,000 spheroids/flask were collected with a magnetic filter funnel with 25 µm nylon mesh filter, and transferred to a 40 ml Kimble Kontes Dounce homogenizer. Spheroids were broken with a tight-fitting pestle (B type) with six strokes. Broken spheroids were transferred to a 50 ml Falcon tube, and the volume adjusted to ∼40 ml with SVM, after which 2.8 ml Percoll was added and the tubes spun at 200g for 5 min at room temperature. The supernatant containing somatic cells was transferred to a beaker and diluted to 200 ml with SVM. The pellet containing gonidia or embryos was washed two times with 50 ml SVM and gonidia collected after each wash by centrifugation at 200g for 5 minutes. Pure gonidia were then collected in a filter funnel using 10 µm nylon mesh, which allows any remaining somatic cells to pass through. To obtain pure somatic cells the diluted supernatant from the Percoll step above was spun at 460g for ∼3 min to pellet any contaminating gonidia. The low speed spin supernatant was then spun at 3,220g for 5 min to obtain a pure somatic cell pellet, which was washed twice with 50 ml SVM prior to extraction of RNA or protein.

Western Blotting

Approximately 1,000 synchronized spheroids at designated stages were hand picked for protein sample preparation. Pelleted spheroids were mixed 1∶1 with Volvox Lysis Buffer (1× PBS supplemented with 1% NP40 [IPEGAL], 1× Sigma Plant Protease Inhibitor Cocktail [catalog number P9599], 5 mM PMSF, 10 mM benzamidine, 5 mM EDTA, 5 mM EGTA). Spheroids and cells were disrupted using a Covaris S220 ultrasonicator according to the manufacturers instructions with the following program settings: PP = 200, DF = 20, CpB = 300, T = 6°C, and t = 300 s in TC 12×12 tubes at 4°C. After lysis the samples were centrifuged at full speed in a microfuge to pellet debris and the supernatant mixed with sample buffer and boiled prior to gel fractionation. SDS-PAGE and Western blotting were performed using standard procedures[70]. SDS-PAGE gels were blotted to Immobilon-P PVDF membranes (Millipore) prior to immunodetection. The rat monoclonal antibody 3F10 (Roche) was used for detection of the HA epitope on immunoblots and for IF. Tubulin was detected with an anti-α-tubulin antibody purchased from Sigma-Aldrich (clone B-5-1-2, catalog number T6074). Antibodies were used at the following dilutions: 250 mg/ml anti-HA, 1∶2,000; and 2 mg/ml anti-α-tubulin, 1∶20,000, all diluted in PBS with 0.05% Tween 20. Blocking was performed with 40 ml of 5% nonfat dry milk in PBS with 0.05% Tween 20 for one hour. Blots were incubated overnight at 4°C min (anti-HA antibody). Blots were washed three times with PBS and 0.05% Tween 20 at room temperature for 10 min, then horse radish peroxidase (HRP)-conjugated goat anti-rat secondary antibody (Thermo Scientific) was used at 1∶2,500 dilution and incubated with blots at room temperature in PBS, 5% nonfat dry milk, 0.05% Tween 20 for 1 h. Blots were then washed three times with PBS and 0.05% Tween 20 at room temperature for 10 min and then briefly rinsed with deionized water. Antigen was detected by chemiluminescence (Luminata Forte Western HRP Substrate, Millipore) using autoradiographic film (HyBlot CL autoradiography film, Denville Scientific Inc.)

Immunofluorescent Staining and Microscopy

Eve::VcMID-BH orEve::VcMID (as a negative control) samples were processed in parallel and imaged under identical conditions. Spheroids were collected using a magnetic funnel with a 25 µm nylon mesh filter (Pall Scientific). Sexual spheroids were fixed with 2% or 4% paraformaldehyde with 1× plant protease inhibitor cocktail (Sigma) (10 µl/ml), 1 µm ALLN (catalog number 208719, VWR International), 1 µm MG132 (catalog number 133407-82-6, Cayman Chemical), 1 mM dithiothreitol (DTT) for 1 h on ice. Spheroids or cells were then washed with PBS and resuspended in cold methanol (−20°C) for five minutes. The methanol wash was repeated two more times and then samples were washed again with PBS and rehydrated in PBS at room temperature for 15 min after which they were adhered to poly-L-lysine–coated cover slips. Cover slips were blocked for 30 min in blocking buffer (5% BSA and 1% cold-water fish gelatin) and incubated for 30 min in the same buffer with 10% (v/v) normal goat sera (Antibodies Incorporated). Cover slips were incubated overnight in anti-HA 1∶500 (1∶300 for somatic cell staining) in 20% blocking buffer at room temperature or at 4°C and washed 6×10 min with 20% blocking buffer in PBS 0.05% Tween 20. After washing out unbound primary antibody cover slips were incubated with AlexaFluor 488 conjugated goat anti-mouse secondary antibodies (Invitrogen) 1∶500 (1∶300 for somatic cells staining) in 20% blocking buffer and incubated for 1 h at room temperature (for 4 h at 4°C for somatic cells staining) in the dark. Following six washes with 20% blocking buffer diluted in PBS and 0.05% Tween20, the cover slips were incubated for 15 min in 2 µg/ml DAPI and then washed for 5 min in PBS. Excess liquid was removed, and the cover slips were mounted with VectaShield (Vector Labs) or Mowiol∶PPD (PPD = p-phenylenediamine 1,4-Benzenediamine hydrochloride [Sigma P1519]) with a 9∶1 ratio. Microscopy was performed with a Leica DMI6000 B using DIC optics or using the following filter cube sets and illumination with a Prior Lumen 200 light source: A4: excitation BP 360/40; dichroic 400; emission BP470/40, L5: excitation BP 480/40; dichroic 505; emission BP 527/30. Where indicated Z stacks were subject to deconvolution using the Leica Advanced Fluorescence Application Suite. Images are representative of results from at least three independent experiments.

• 1.Hoekstra RF (1987) The evolution of sexes. ExperientiaSuppl 55: 59–91. [DOI] [PubMed] [Google Scholar]
• 2.Parker GA (2011) The origin and maintenance of two sexes (anisogamy), and their gamete sizes by gamete competition. Togashi T, Cox PA, editors. The Evolution of Anisogamy. Cambridge University Press. [Google Scholar]
• 3.Ni M, Feretzaki M, Sun S, Wang X, Heitman J (2011) Sex in fungi. Annu Rev Genet45: 405–430. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 4.Umen JG (2011) Evolution of sex and mating loci: an expanded view from Volvocine algae. Curr Opin Microbiol14: 634–641. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 5.Cervantes MD, Hamilton EP, Xiong J, Lawson MJ, Yuan D, et al. (2013) Selecting one of several mating types through gene segment joining and deletion inTetrahymena thermophila. PLoS Biol11: e1001518. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 6.Bloomfield G, Skelton J, Ivens A, Tanaka Y, Kay RR (2010) Sex determination in the social amoebaDictyostelium discoideum. Science330: 1533–1536. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 7.Parker GA, Baker RR, Smith VG (1972) The origin and evolution of gamete dimorphism and the male-female phenomenon. Journal of Theoretical Biology36: 529–553. [DOI] [PubMed] [Google Scholar]
• 8.Randerson JP, Hurst LD (2001) The uncertain evolution of the sexes. Trends Ecol Evol16: 571–579. [Google Scholar]
• 9.Billiard S, López-Villavicencio M, Devier B, Hood ME, Fairhead C, et al. (2011) Having sex, yes, but with whom? Inferences from fungi on the evolution of anisogamy and mating types. Biological reviews of the Cambridge Philosophical Society86: 421–442. [DOI] [PubMed] [Google Scholar]
• 10.Blute M (2012) The evolution of anisogamy: more questions than answers. Biol Theory7: 3–9. [Google Scholar]
• 11.Charlesworth B (1978) The population genetics of anisogamy. J Theor Biol73: 347–357. [DOI] [PubMed] [Google Scholar]
• 12.Kirk DL (2005) A twelve-step program for evolving multicellularity and a division of labor. BioEssays27: 299–310. [DOI] [PubMed] [Google Scholar]
• 13.Bell G (1985) The origin and early evolution of germ cells as illustrated by the Volvocales. Halvorson HO, Monroy A, editors. The origin and evolution of sex. New York: Alan R. Liss Inc. pp. 221–256. [Google Scholar]
• 14.Goodenough U, Lin H, Lee J-H (2007) Sex determination inChlamydomonas. Semin Cell Dev Biol18: 350–361. [DOI] [PubMed] [Google Scholar]
• 15.Nozaki H (1996) Morphology and evolution of sexual reproduction in the Volvocaceae (Chlorophyta). J Plant Res109: 353–361. [Google Scholar]
• 16.Coleman AW (2012) A comparative analysis of the Volvocaceae (Chlorophyta). J Phycol48: 491–513. [DOI] [PubMed] [Google Scholar]
• 17.Knowlton N (1974) A note on the evolution of gamete dimorphism. J Theor Biol46: 283–285. [DOI] [PubMed] [Google Scholar]
• 18.Madsen JD, Waller DM (1983) A note on the evolution of gamete dimorphism in algae. American Naturalist121: 443–447. [Google Scholar]
• 19.Wiese L, Wiese W, Edwards DA (1979) Inducible anisogamy and the evolution of oogamy from isogamy. Ann Bot44: 131–139. [Google Scholar]
• 20.Ferris P, Olson BJSC, De Hoff PL, Douglass S, Casero D, et al. (2010) Evolution of an expanded sex-determining locus inVolvox. Science328: 351–354. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 21.Nozaki H, Mori T, Misumi O, Matsunaga S, Kuroiwa T (2006) Males evolved from the dominant isogametic mating type. Curr Biol16: R1018–R1020. [DOI] [PubMed] [Google Scholar]
• 22.Ning J, Otto TD, Pfander C, Schwach F, Brochet M, et al. (2013) Comparative genomics inChlamydomonas andPlasmodium identifies an ancient nuclear envelope protein family essential for sexual reproduction in protists, fungi, plants, and vertebrates. Genes Dev27: 1198–1215. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 23.Ferris PJ, Goodenough UW (1994) The mating-type locus ofChlamydomonas reinhardtii contains highly rearranged DNA sequences. Cell76: 1135–1145. [DOI] [PubMed] [Google Scholar]
• 24.De Hoff PL, Ferris P, Olson BJSC, Miyagi A, Geng S, et al. (2013) Species and population level molecular profiling reveals cryptic recombination and emergent asymmetry in the dimorphic mating locus ofC. reinhardtii. PLoS Genet9: e1003724. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 25.Ferris PJ, Goodenough UW (1997) Mating type inChlamydomonas is specified by mid, the minus-dominance gene. Genetics146: 859–869. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 26.Lin H, Goodenough UW (2007) Gametogenesis in theChlamydomonas reinhardtii minus mating type is controlled by two genes, MID and MTD1. Genetics176: 913–925. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 27.Ferris PJ, Pavlovic C, Fabry S, Goodenough UW (1997) Rapid evolution of sex-related genes inChlamydomonas. Proc Natl Acad Sci U S A94: 8634–8639. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 28.Hamaji T, Ferris PJ, Coleman AW, Waffenschmidt S, Takahashi F, et al. (2008) Identification of the minus-dominance gene ortholog in the mating-type locus ofGonium pectorale. Genetics178: 283–294. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 29.Setohigashi Y, Hamaji T, Hayama M, Matsuzaki R, Nozaki H (2011) Uniparental inheritance of chloroplast DNA is strict in the isogamous volvocaleanGonium. PLoS ONE6: e19545. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 30.Hamaji T, Ferris PJ, Nishii I, Nishimura Y, Nozaki H (2013) Distribution of the sex-determining gene MID and molecular correspondence of mating types within the isogamous genusGonium (Volvocales, Chlorophyta). PLoS ONE8: e64385. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 31.Kirk DL, Nishii I (2001)Volvox carteri as a model for studying the genetic and cytological control of morphogenesis. Dev Growth Differ43: 621–631. [DOI] [PubMed] [Google Scholar]
• 32.Kirk DL (2001) Germ–soma differentiation inVolvox. Dev Biol238: 213–223. [DOI] [PubMed] [Google Scholar]
• 33.Kochert G, Yates I (1974) Purification and partial characterization of a glycoprotein sexual inducer fromVolvox carteri. Proc Natl Acad Sci U S A [DOI] [PMC free article] [PubMed] [Google Scholar]
• 34.Starr RC, Jaenicke L (1974) Purification and characterization of the hormone initiating sexual morphogenesis inVolvox carteri f.nagariensis Iyengar. Proc Natl Acad Sci U S A71: 1050–1054. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 35.Mages HW, Tschochner H, Sumper M (1988) The sexual inducer ofVolvox carteri. Primary structure deduced from cDNA sequence. FEBS Letters234: 407–410. [DOI] [PubMed] [Google Scholar]
• 36.Starr RC (1969) Structure, reproduction and differentiation in Volvox carteri f. nagariensis Iyengar, strains HK9 & 10. Arch Protistenkd111: 204–222. [Google Scholar]
• 37.Kochert G (1968) Differentiation of reproductive cells inVolvox carteri. J Protozool15: 438–452. [DOI] [PubMed] [Google Scholar]
• 38.Starr RC (1975) Meiosis inVolvox carteri f.nagariensis. Archiv fur Protistenkunde117: 187–191. [Google Scholar]
• 39.Bachtrog D, Kirkpatrick M, Mank JE, McDaniel SF, Pires JC, et al. (2011) Are all sex chromosomes created equal?Trends Genet27: 350–357. [DOI] [PubMed] [Google Scholar]
• 40.Hiraide R, Kawai-Toyooka H, Hamaji T, Matsuzaki R, Kawafune K, et al. (2013) The evolution of male-female sexual dimorphism predates the gender-based divergence of the mating locus gene MAT3/RB. Mol Biol Evol30: 1038–1040. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 41.Starr RC, Jaenicke L (1989) Cell differentiation inVolvox carteri (Chlorophyceae): the use of mutants in understanding patterns and control. Algae as experimental systems. New York: Alan R. Liss Inc. pp. 135–147. [Google Scholar]
• 42.Nishii I, Miller SM (2010)Volvox: simple steps to developmental complexity?Curr Opin Plant Biol13: 646–653. [DOI] [PubMed] [Google Scholar]
• 43.Umen JG, Olson BJSC (2012) Genomics of volvocine algae. Piganeau G, editor. Genomic insights into the biology of algae. New York: Elsevier, Vol. 64. pp. 185–243. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 44.Shubin N, Tabin C, Carroll S (2009) Deep homology and the origins of evolutionary novelty. Nature457: 818–823. [DOI] [PubMed] [Google Scholar]
• 45.Gehring WJ (2005) New perspectives on eye development and the evolution of eyes and photoreceptors. J Hered96: 171–184. [DOI] [PubMed] [Google Scholar]
• 46.Herron MD, Hackett JD, Aylward FO, Michod RE (2009) Triassic origin and early radiation of multicellular volvocine algae. Proc Natl Acad Sci U S A106: 3254–3258. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 47.McCourt RM, Delwiche CF, Karol KG (2004) Charophyte algae and land plant origins. Trends Ecol Evol (Amst)19: 661–666. [DOI] [PubMed] [Google Scholar]
• 48.Levin TC, King N (2013) Evidence for sex and recombination in the choanoflagellateSalpingoeca rosetta. Curr Biol23: 2176–2180. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 49.Nakada T, Shinkawa H, Ito T, Tomita M (2009) Recharacterization ofChlamydomonas reinhardtii and its relatives with new isolates from Japan. J Plant Res123: 67–78. [DOI] [PubMed] [Google Scholar]
• 50.Pires ND, Dolan L (2012) Morphological evolution in land plants: new designs with old genes. Philos Trans R Soc Lond B Biol Sci367: 508–518. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 51.Davidson EH (2006) Gene regulatory networks and the evolution of animal body plans. Science311: 796–800. [DOI] [PubMed] [Google Scholar]
• 52.Wray GA (2003) Transcriptional regulation and the evolution of development. Int J Dev Biol47: 675–684. [PubMed] [Google Scholar]
• 53.Monteiro A (2012) Gene regulatory networks reused to build novel traits. BioEssays34: 181–186. [DOI] [PubMed] [Google Scholar]
• 54.Aravind L, Anantharaman V, Venancio TM (2009) Apprehending multicellularity: regulatory networks, genomics, and evolution. Birth Defects Res C Embryo Today87: 143–164. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 55.Starr R (1968) Cellular differentiation inVolvox. Proc Natl Acad Sci U S A59: 1082–1088. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 56.Schauser L, Wieloch W, Stougaard J (2005) Evolution of NIN-like proteins inArabidopsis, rice, andLotus japonicus. J Mol Evol60: 229–237. [DOI] [PubMed] [Google Scholar]
• 57.Worden AZ, Lee J-H, Mock T, Rouzé P, Simmons MP, et al. (2009) Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotesMicromonas. Science324: 268–272. [DOI] [PubMed] [Google Scholar]
• 58.Curtis BA, Tanifuji G, Burki F, Gruber A, Irimia M, et al. (2012) Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs. Nature492: 59–65. [DOI] [PubMed] [Google Scholar]
• 59.Smith GM (1944) A comparative study of the species ofVolvox. Transactions of the American Microscopical Society63: 265–310. [Google Scholar]
• 60.Campbell AM, Rayala HJ, Goodenough UW (1995) The iso1 gene ofChlamydomonas is involved in sex determination. Mol Biol Cell6: 87–95. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 61.Ferris PJ, Woessner JP, Goodenough UW (1996) A sex recognition glycoprotein is encoded by the plus mating-type gene fus1 ofChlamydomonas reinhardtii. Mol Biol Cell7: 1235–1248. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 62.Nozaki H, Ott FD, Coleman AW (2006) Morphology, molecular phylogeny and taxonomy of two new species of Pleodorina (Volvoceae, Chlorophyceae). J Phycol42: 1072–1080. [Google Scholar]
• 63.Al-Hasani H, Jaenicke L (1992) Characterization of the sex-inducer glycoprotein ofVolvox carteri f.weismannia. Sex Plant Reprod5: 8–12. [Google Scholar]
• 64.Bull J (1978) Sex chromosomes in haploid dioecy: a unique contrast to Muller's theory for diploid dioecy. Am Nat112: 245–250. [Google Scholar]
• 65.Rice W (1984) Sex chromosomes and the evolution of sexual dimorphism. Evolution38: 735–742. [DOI] [PubMed] [Google Scholar]
• 66.Yamato KT, Ishizaki K, Fujisawa M, Okada S, Nakayama S, et al. (2007) Gene organization of the liverwort Y chromosome reveals distinct sex chromosome evolution in a haploid system. Proc Natl Acad Sci U S A104: 6472–6477. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 67.Prochnik SE, Umen J, Nedelcu AM, Hallmann A, Miller SM, et al. (2010) Genomic analysis of organismal complexity in the multicellular green algaVolvox carteri. Science329: 223–226. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 68.Umen JG, Goodenough UW (2001) Control of cell division by a retinoblastoma protein homolog inChlamydomonas. Genes Dev15: 1652–1661. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 69.Fang S-C, de los Reyes C, Umen JG (2006) Cell size checkpoint control by the retinoblastoma tumor suppressor pathway. PLoS Genet2: e167. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 70.Sambrook JJ, Russell DW (2001) Molecular cloning: a laboratory manual. 3rd edition. Cold Spring Harbor (New York): CSHL Press. 1 pp. [Google Scholar]
• 71.Gruber H, Kirzinger S, Schmitt R (1996) Expression of theVolvox gene encoding nitrate reductase: Mutation-dependent activation of cryptic splice sites and intron-enhanced gene expression from a cDNA. Plant Mol Biol31: 1–12. [DOI] [PubMed] [Google Scholar]
• 72.Schiedlmeier B, Schmitt R, Müller W, Kirk MM, Gruber H, et al. (1994) Nuclear transformation ofVolvox carteri. Proc Natl Acad Sci U S A91: 5080–5084. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 73.Huskey RJ, Griffin BE, Cecil PO, Callahan AM (1979) A preliminary genetic investigation ofVolvox carteri. Genetics91: 229–244. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 74.Edgar RC (2004) MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics5: 113. [DOI] [PMC free article] [PubMed] [Google Scholar]

Supplementary Materials

Data Availability Statement

The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.